We propose a cost efficient cross-subcarrier pre-coder design (CSPD) for massive multiple-input multiple-output orthogonal frequency division multiplexing (MIMO-OFDM) downlink with imperfect channel state information (CSI). To reduce the high computational complexity caused by individual precoder design for each subcarrier, we design transform domain precoding vectors (TDPVs), from which the precoders for a set of subcarriers can be obtained through a transform. The number of TDPVs is much less than that of subcarriers, and the number of total parameters to be designed can be reduced significantly. The main objective is to maximize an upper bound of the ergodic sum-rate by exploiting the a posteriori beam-based statistical channel model. We provide a concave minorizing function of the upper bound of the ergodic sum-rate and then derive the stationary points of a concave quadratic optimization problem with this minorizing function. To reduce the dimension of the matrix inversion in the stationary points, we propose an algorithm by using block coordinate descent (BCD) method with power allocation.