Analytical estimates of errors, associated with random statistical noise in magnetic field and other data to which the minimum/maximum variance analysis technique is commonly applied, are derived from first principles. A systematic expansion procedure is used in which the expansion parameter is proportional to the noise amplitude and inversely proportional to the square root of the number of vector data samples, K. The two special cases where the signal-to-noise ratio is large and small are considered for arbitrary noise distributions. The ideal case of small errors and isotropic Gaussian noise allows determination of uncertainty cones of elliptic cross section for all three eigenvectors, x1, x2, x3, and uncertainty intervals for all three eigenvalues, &lgr;1, &lgr;2, &lgr;3, of the variance matrix, all in terms of the eigenvalues themselves and the number of data points. Denoting the angular standard deviation (in radians) of vector xi toward or away from vector xj by Δϕij these results are Δϕij=¿[&lgr;3(&lgr;i+&lgr;j-&lgr;3)/(K-1)(&lgr;i-&lgr;j)2>1/2 and Δ&lgr;i=¿[2&lgr;3(2&lgr;i-&lgr;3)/(K-1)>1/2 where &lgr;3 is the smallest eigenvalue. The predictions from these formulas compare favorably with results from simple bootstrap numerical error analysis. Applications to maximum variance analysis of electric field data and to the study of anisotropic plasma pressure tensors are discussed as special cases. ¿ 1998 American Geophysical Union