已知生成多维正态分布数据可以由
x<-mvrnorm("定义样本量","给定均值向量"," 给定协方差矩阵")产生。
我只知道产生一维的卡方分布数据可以由
y<-rchisq("定义样本量","定义自由度“，"定义非中心参数")
问题：如何产生满足给定均值和协方差的服从多维卡方分布的样本数据？比如有x1,x2,x3均服从相同自由度的卡方分布，但是两两不独立，协方差矩阵为
    [,1] [,2] [,3]
[1,]  1.0  0.5  0.0
[2,]  0.5  1.0  0.2
[3,]  0.0  0.2  1.0
该如何生成满足此条件下的数据？
（注：我的目标是生成多维的服从非对称分布的数据，所以我拿卡方分布来做数据模拟）
恳请各位坛友给予指导，谢谢！