求助~如何用R生成服从多元正态分布的数据？？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

library(MASS)
options(digits=3)
set.seed(1234)
mean <- c(230.7,146.7,3.6)
sigma <- matrix(c(15360.8,6721.2,-47.1,6721.2,4700.9,-16.5,-47.1,-16.5,0.3),nrow=1,ncol=3)
mydata <- mvrnorm(n=500,mean,sigma)    #总弹出来“ incompatible arguments”不兼容的参数
mydata <- as.data.frame(mydata)
names(mydata) <- c('y','x1','x2')
dim(mydata)


“mydata <- mvrnorm(n=500,mean,sigma)”这句总出问题，总弹出来“ incompatible arguments”不兼容的参数。
哪位大神能帮忙解答一下~~~谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：多元正态分布 正态分布 如何用 Library options 正态分布 library 如何

sigma应该是一个3*3的矩阵，将生成sigma的那句命令里的“nrow = 1"改成“nrow = 3”。

hugebear 发表于 2015-4-29 21:46
sigma应该是一个3*3的矩阵，将生成sigma的那句命令里的“nrow = 1"改成“nrow = 3”。

谢谢你啊~问题解决了~