### 多元正态分布N(mu,Sigma) ###
mu<-runif(3)
Sigma=diag(1:3)
n<-4

# 不用包 #
#cholesky,正定！
LL<-t(chol(Sigma))
Z<-matrix(rnorm(n*length(mu)),ncol=n)
mu+LL%*%Z

# 用包 #
library(MASS)
mvrnorm(n,mu,Sigma) #eigen,半正定！

library(mvtnorm)
rmvnorm(n,mu,Sigma,method=c("eigen","svd","chol")[1])#三种方法可选！推荐【1】。

### 多元t分布t(mu,Sigma) with df v ###
v<-4#任意自由度
# 不用包 #
WW<-sqrt(v/rchisq(n,v))
LL<-t(chol(Sigma))
Z<-matrix(rnorm(n*length(mu)),ncol=n)
mu+t(WW*t(LL%*%Z))

# 用包 #
library(mixAK)
rMVT(n,v,mu,1,Sigma)

library(mvtnorm)
rmvt(n,Sigma,v,mu,type=c("shifted","Kshirsagar")[1])#推荐【2】

### 多元F分布 ###
#!!!没听说过，求大神告诉_(:з」∠)_
#有文献的最好。

嗯，谢谢了，我试试

请问您知道多元t分布的性质吗

这个具体的可以查文献看看，多元也就是一元的推广，很多性质都是可以类比的，理解不太复杂，证的话就难搞了。

哦，谢谢了