我正在R中工作,使用MASS
创建多元分布矢量。
我正在尝试为Y创建1000个值,这些值必须由代码的第二行到最后一行给出的公式来计算。当我运行此命令时,R表示存在50个或更多警告,其中第一个警告是:
错误:“()”中出现意外的')'
有人可以告诉我我的代码是否正确吗?
library("MASS")
Y <- c()
for(i in 1:1000){
set.seed(i)
U <- rnorm(n=40, mean=0, sd=0.04)
X <- mvrnorm(n=40, mu=mu, Sigma=S)
Y[i] <- 1.1 - 2.4*X[,1] + 1.4*X[,2] + U
}