我以前从未使用过R,并且我正在为类运行方程式,因此什么也没有发生/它花了很长时间才能通过它。
对R缺乏经验的人,我尝试重置它。
> x4 = c(rgamma(5800,3))
> x4ave = mean(x4)
> x4sd = sd(x4)
> hist(x4)
> samp54 = sample(x4,54,replace = F)
> samp54ave = mean(samp54)
> x4distr = c()
> for(i in 1:550){samp=sample(x4,54,replace=F);x4distr=c(x4distr,mean(samp))}
> x4distrterm = x4distr[220]
> x4distrterm
[1] 2.916715
> x4distr2 = c()
> for(i in 1:35000){samp=sample(x4,250,replace=F);x4distr2=c(x4distr2,mean,(samp))}
我希望我的x4distr2可以填满,但是它不会停止运行数小时。
答案 0 :(得分:1)
您有两个错误:
1)您正在将mean
(函数)与输出而不是mean
的结果连接起来。
2)您应该预先分配内存来存储结果,否则每次存储新结果时都要重新分配。
一个简单的解决方法:
x4distr2 = numeric(35000)
for(i in 1:35000){samp=sample(x4,250,replace=F);x4distr2[i]=mean(samp)}