[软件熟悉度]:
入门
[问题叙述]:
我目前想要做出多组数据的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)))
}
}