mvrnorm problem
As Peter Dalgaard has already pointed out, you both need to read the help page. If you supply a matrix where the help page asks for a vector, you are likely to get troubles. I will coerce it for the next version, but I would like to point out that in about 12 years of providing mvrnorm (often as here without credit), this it the first time anyone has made this error in public. Anticipating what people will try is not easy.
On Mon, 2 Feb 2004 Ted.Harding at nessie.mcc.ac.uk wrote:
On 02-Feb-04 Stuart V Jordan wrote:
I am trying to simulate draws from a multivariate normal using mvrnorm,
and
am getting the following error message:
Error in mu + eS$vectors %*% diag(sqrt(pmax(ev, 0)), p) %*% t(X) :
non-conformable arrays
[...] Hmmm ... using the same B and V as giben by Stuart Jordan, I get:
mvrnorm(1,B,V)
[1] -179.8332342 0.9632282 2.5687489 -2.2337125 47.2717626 [6] -3745.3844310 0.9632282 0.2839965 -0.1500585 0.5804460 [11] 5.1919420 -4.7381677 -0.2382119
mvrnorm(2,B,V)
Error in mu + eS$vectors %*% diag(sqrt(pmax(ev, 0)), p) %*% t(X) :
non-conformable arrays
?????
Ted.
PS By the way, [near] singularity seems to have nothing to do
with it here:
If I brutally make V non-singular:
diag(V) <- diag(V)+10000000
it goes the same way.
--------------------------------------------------------------------
E-Mail: (Ted Harding) <Ted.Harding at nessie.mcc.ac.uk>
Fax-to-email: +44 (0)870 167 1972
Date: 02-Feb-04 Time: 22:17:13
------------------------------ XFMail ------------------------------
______________________________________________ R-help at stat.math.ethz.ch mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595