r - 模拟适用于 n= 100 但不适用于 n = 300

Question

是否有原因这适用于一个样本量但在增加时不起作用？已经研究了一段时间，n=100 运行良好。它是一种模拟，它为 X1、X2 和误差生成数据，然后使用这些值来计算 Y 值。

###Change to sample size n2=300
n1<-100
n2<- 300
n3 <- 1000 
ITER <- 500

#set beta coefficients
B0 <- 1
B1 <- 1
B2 <- 1
B3_1 <-1
B3_2 <- 2
B3_3 <- 5
#set gamma values
c1 <-.1
c2 <- .25
c3 <- .75
sim10 <- matrix(data = 0, nrow = ITER, ncol = 4)

colnames(sim10) <- c("Y", "X1", "X2", "E")

# simulate X1s,
for (i in 1:ITER){
sim10[,2] <-runif(n2, -3, 3)
}

# simulate X2s
for(i in 1:ITER){
sim10[,3] <-  (c1*sim10[,2])+ rnorm(n2,0,1)
}

# simulate error
for(i in 1:ITER){
sim10[,4] <- dnorm(300, 0,1)  
}

# y's
for(i in 1:ITER)
sim10[,1] <- B0+ B1*sim10[,2] + B2*sim10[,3] + B3_1*(sim10[,3]^2) + 
sim10[,4]

标签： rsimulation