n <- c(12,24)
mu<-c(6.573,6.5)
sigma<-sqrt(0.25)
Diseased.Data<-round(rnorm(n[1], mu[1], sigma), 4)
Healthy.Data<-round(rnorm(n[2], mu[2], sigma), 4)
g <- c(2,3,4)
for(i in 1:3){
pool.dis.data <- matrix(NA,n[1]/g[i],n[1]/g[i])
for(j in n[1]/g[i]){
pool.dis.data <- replicate(n[1]/g[i],mean(sample(Diseased.Data,g[i])))
}
}
当我运行此代码时,我只得到g的最后一个元素的答案。我需要的是每个列矩阵应该包含g中每个元素的元素,例如矩阵应该是这样的 m&lt; -cbind(a1,a2,a3,a4,a5,a6),m1&lt; -cbind(b1,b2,b3,b4), m2&lt; - cbind(c1,c2,c3)