Power-domain non-orthogonal multiple access (NOMA) assigns different power levels for near and far users in order to discriminate their signals by employing successive interference cancellation (SIC) at the near user. In this context, multiple-input-single-output NOMA (MISO-NOMA), where the base station (BS) is equipped with multiple antennas while each mobile user has a single antenna receiver, is shown to have a better overall performance by using the knowledge of instantaneous channel state information (CSI). However, this requires prior estimation of CSI using pilot transmission, which increases the transmission overhead. Moreover, its performance is severely degraded in the presence of CSI estimation errors. In this work, we provide statistical beamforming solutions for downlink power-domain NOMA that utilize only knowledge of statistical CSI, thus reducing the transmission overhead significantly. First, we derive the outage probabilities for both near and far users in the multi-user NOMA system without imposing strong assumptions, such as Gaussian or Chi-square distribution. This is done by employing the exact characterization of the ratio of indefinite quadratic form (IQF). Second, this work proposes two techniques to obtain the optimal solution for beam vectors which rely on the derived outage probabilities. Specifically, these two methods are based on (1) minimization of total beam power while constraining the outage probabilities to the QoS threshold, and (2) minimization of outage probabilities while constraining the total beam power. These proposed methods are non-convex function of beam vectors and, hence, are solved using numerical optimization via sequential quadratic programming (SQP). Since the proposed methods do not require pilot transmission for channel estimation, they inherit better spectral efficiency. Our results validate the theoretical findings and prove the supremacy of the proposed method.