search for: rmst

Hi R-users,   I have this code below and use some code from sn package for multivariate skew t (rmst) to generate a set of random numbers.   I try to do a looping to produce ‘rand’ for different values of degrees of freedom (degf).  Actually, I’m not sure how to specify the attribute for the rand.  As we know, rand will give us a set of results + all the parameters used, so just wonder whether we...

...e "vector" is the matrix row, on which I evaluate my likelihood and the rest in parametrized in a standard way, just as the help pages give it. In large dimensions, I get a zero value of the density (which is probably due to numerical issues). I tried the following dummy example t1 <- rmst(1,mu,P,alpha, nu) t2 <- dmst(t1, mu, alpha,nu) and t2 remains to be zero. Can anyone help me on this one? thanks in advance, Konrad -- "We are what we pretend to be, so we must be careful about what we pretend to be" Kurt Vonnegut Jr. "Mother Night"...

> -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Caio Lucidius > Naberezny Azevedo > Sent: 01 September 2005 12:09 > To: Help mailing list - R > Subject: [R] Multivariate Skew Normal distribution > > > Hi all, > > Could anyone tell me if there is any package (or function)