[求助]一个循环问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

一组数据(名称：re）：
34    37    68    47    77    73
79    73    91    76    70    99
85    25    104    84    42    56
106    23    80    99    75    107
101    96    81    90    97    89
98    91    51    37    73    101
81    64    82    100    25    86
97    19    83    97    50    95
...
...
每行分别做bootstrap抽出6个值，进行10000次，得到每行均值的检验p值。
for ( i in 1:nrow(re))
{
  x<-as.numeric(re[i,])
  for(j in 1:10000)
    {
    data=sample(x,6,replace=TRUE,prob=NULL)
    meansamp=mean(data)
    {
    meanbootsamp[j]=meansamp
    }
    }
   pboot <- length(which(meanbootsamp>=mean(x)))/1000
}


运行后，结果有问题，请帮忙看看。
多谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Bootstrap Bootstra boots Trap boot 求助

具体是什么问题呢？ 最后一行总的抽样次数1000应该是10000吧？

另外，

meansamp=mean(data)

    {
    meanbootsamp[j]=meansamp
    }

明显有错。

需事先定义 meanbootsamp向量，也用不着最里层的{}啊！！

哦，最后1000是笔误。
在学习抽样理论，网上找了一个R文件，就照葫芦画瓢，没画好。

re<-read.delim(file="re.txt")
meanbootsamp<-rep(NA,nrow(re))
for ( i in 1:nrow(re)){
x<-as.numeric(re[i,])
for(j in 1:10000){
    data=sample(x,6,replace=TRUE,prob=NULL)
    meansamp=mean(data)
    meanbootsamp[j]=meansamp
    }
    pboot <- length(which(meanbootsamp[j]>=mean(x)))/10000
    }
对循环还是有点晕，麻烦帮看一下，对meanbootsamp和pboot下标的定义好像有点问题。

你定义了meanbootsamp向量，pboot向量也要事先定义啊！