[軟體熟悉度]:
入門
[問題敘述]:
我目前想要做出多組數據的sp^2來估計群體的總變異數,目前遇到的問題是如果單跑
sqrt((((D[3+(i*4500-4500),5])^2)+((D[4+(i*4500-4500),5])^2))/(k-1-(k/2-1)))
這行的話出來的sp^2值是正確的
但如果用i=5000,k=4500時的for迴圈做大量運算時sp^2就會跑掉
請問各位先進大大們有解法嗎,小弟被困住好多天了,謝謝各位大大
[程式範例]:
for(i in c(1:5000)){
for (k in seq(4,4500,by = 2)) {
D[4+i*4500-4500,6] <- sqrt((((D[3+(i*4500-4500),5])^2)+((D[4+(i*4500-4500),5])
^2))/(k-1-(k/2-1)))
}
}