Certain aspects of data generation are studied through multivariate autoregressive (AR) models. The main emphasis is on the preservation of certain desired moments and the effect of initial values on these moments. The problem of preservation of moments is approached in a nontraditional way by starting with the initial values. For this purpose, general AR processes with a random start and with time-varying parameters are introduced to lay a foundation for the analysis of all types of AR processes, including the periodic cases. It is shown that an AR process with a random start and with parameters obtained from the moment equations is capable of generating jointly multivariate normal vectors with any specified means and covariance matrices, and with any specified autocovariance matrices up to a given lag. With a random start, there is no transition period involved for achieving these moments. A simple solution is proposed for matrix equations of the form BB(T) = A which appear in the moment equations. The aggregation properties of general AR process are also studied.