> > Thanks to Spencer Graves for provinding this clarification about pmst. > As the author of pmst, I was really the one expected to answer the > query, but I was on travel in the last two weeks and did not read > this query. > > As a complement to what already explained, the reason of the sharp change > from k=19 to k=20 is as follows: pmst (package sn) calls pmt (package > mnormt) with k increased by 1, and pmt calls a Fortran routine (written by > Alan Genz), which issues an error if k>20 or k<1. Hence, the effective > maximum number of dimensions allowed by pmst is 19.
Ok, thanks a lot to both of You for help - with Your explanations, You really clarified my view of the problem (and saved me a lot of time I would've spent on digging through my own code). best wishes, Konrad Banachewicz best wishes, > > Adelchi Azzalini > > > Also, have you asked about this directly to the maintainers of > the > > "sn", "mnormt" and "mvtnorm" packages? They might have other > suggestions. > > > > Hope this helps. > > Spencer Graves > -- > "We are what we pretend to be, so we must be careful about what we pretend > to be" > > Kurt Vonnegut Jr. "Mother Night" [[alternative HTML version deleted]] ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html