# MultinormalDistribution not having the right mean [closed]

I am using the function MultinormalDistribution to generate a random vector using the following code:

RandomReal[
MultinormalDistribution[
c03 ConstantArray[
1, {Binomial[n, 2], Binomial[n, 2]}] + (1 -
c03) IdentityMatrix[Binomial[n, 2]]]]


Here c03 = 0.959016 and n = 100 so the covariance matrix is very large. To check the resulting random vector, I plugged it into Histogram to see if it made sense. The histogram is a bell curve like it should be, but its mean is not the theoretically predicted zero. In fact, every time I run my code I get a histogram with the same shape but with a different mean which ranges as far as -1 to 1. However, if I generate a random vector of the same length using

RandomVariate[NormalDistribution[], Binomial[n,2]]


it has the correct mean of about zero without fail. Am I using MultinormalDistribution incorrectly? Are these errors because the covariance matrix is so big? I really appreciate any thoughts.

## closed as off-topic by MarcoB, m_goldberg, Henrik Schumacher, Sektor, rhermansMay 30 '18 at 17:43

This question appears to be off-topic. The users who voted to close gave this specific reason:

• "This question arises due to a simple mistake such as a trivial syntax error, incorrect capitalization, spelling mistake, or other typographical error and is unlikely to help any future visitors, or else it is easily found in the documentation." – MarcoB, m_goldberg, Henrik Schumacher, Sektor, rhermans
If this question can be reworded to fit the rules in the help center, please edit the question.

• Lots of missing stuff. Can you give typical values for c03, c13, and undirectedlist? – J. M. will be back soon Nov 3 '17 at 3:44
• In my case c03 = c13 = 0.959016 so actually the sparse array with undirectedlist doesn't matter. However undirectedlist is a list of pairs of indices of length 6Binomial[n,3] = 970200 and is such that the sparse array is symmetric. – TimF Nov 3 '17 at 4:40
• Then why did you not give the simpler example to begin with? Please include this additional information in your post by editing it. In any event: RandomVariate[MultinormalDistribution[ToeplitzMatrix[Prepend[ConstantArray[c03, Binomial[n, 2] - 1], 1]]]] – J. M. will be back soon Nov 3 '17 at 4:42
• Also, your guess about the ill-conditioning seems to be correct: With[{c03 = 0.959016, n = 100}, LinearAlgebraMatrixConditionNumber[ToeplitzMatrix[Prepend[ConstantArray[c03, Binomial[n, 2] - 1], 1]]]] gives 231612.`. – J. M. will be back soon Nov 3 '17 at 4:47
• I increased the precision in the computation and now the results are making sense. Thanks! – TimF Nov 3 '17 at 16:13