x***2
发帖数: 946 | 1 我在function里面生成multivariate norm变量时出错。
test<-function(rep,n1,n2,mu1,mu2,x,v1,v2,vx,rho){
...
vx1y<-rho*sqrt(v1*vx)
vx2y<-rho*sqrt(v2*vx)
xy1<-rmvnorm(n=n1,mean=c(mu1,x),sigma=matrix(c(v1,vx1y,vx1y,vx),2,2))
xy2<-rmvnorm(n=n2,mean=c(mu2,x),sigma=matrix(c(v2,vx2y,vx2y,vx),2,2))
在这里就报错，错误信息
Error in rmvnorm(n = n1, mean = c(mu1, x), sigma = matrix(c(v1, vx1y, :
mean and sigma have non-conforming size
如果我改成数字，function就没问题。
# xy1<-rmvnorm(n=n1,mean=c(0,0),sigma=matrix(c(1,0.1,0.1,1),2,2))
# xy2<-rmvnorm(n=n2,mean=c(0,0),sigma=matrix(c(1,0.1,0.1,1),2,2))
请问这是什么原因？谢谢点拨。 | f*******m
发帖数: 94 | 2 在run这个function之前先清除R里面的记录：
rm(list=ls(all=TRUE))
应该没有问题 |
|
