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... 阅读全帖 |
|
