Skip to content
Prev 76694 / 398502 Next

Multivariate Skew Normal distribution

On Thu, 1 Sep 2005 13:09:00 +0000 (GMT), Caio Lucidius Naberezny
Azevedo wrote:
CLNA> 
CLNA>  Could anyone tell me if there is any package (or function) that
CLNA>  generates values from a multivariate skew normal distribution?

try the following

library(sn)
location <- c(20,4) # e.g. for a two-dimensional variable
dispers  <- matrix(c(3,-2,-2,2), 2,2)
skew <- c(10,-6)
rmsn(n=10, xi=location, Omega=dispers, alpha=skew) # for skew-normal data
rmst(n=10, xi=location, Omega=dispers, alpha=skew, df=5) # for skew-t data

see also help(rsn) and help(rst) for the univariate case

for more information, see also (as already pointed out in the list):
  http://azzalini.stat.unipd.it/SN

best wishes,
Adelchi Azzalini