不存储值的For循环列表

不存储值的For循环列表,r,for-loop,subset,repeat,R,For Loop,Subset,Repeat,我有个人在不同温度下反复执行的行为数据集,例如: ID Test Behaviour Temperature A12.4.2 ONE 8.64 4 A12.4.2 TWO 7.63 5 A6.3.3 ONE 1.81 3 A6.3.3 TWO 2.47 9 B12.4.1 ONE

我有个人在不同温度下反复执行的行为数据集,例如:

  ID          Test   Behaviour     Temperature
 A12.4.2      ONE     8.64            4
 A12.4.2      TWO     7.63            5
 A6.3.3       ONE     1.81            3
 A6.3.3       TWO     2.47            9
 B12.4.1      ONE     1.17            12
 B12.4.1      TWO     3.96            2
 E9.4.2       ONE     13.04           13
 E9.4.2       TWO     9.51            6
 ...

我使用下面的代码随机地将该数据集子集,然后对子集运行重复性分析,最后通过重复性分析生成R值和CI值

P<-10000

R_value<-numeric(length=P)
CI_value<-numeric(length=P)
for(i in 1:P){
  newdata<-Data[Data$ID %in% sample(unique(Data$ID), 16), ]
  m1<-rptR::rpt(((Behaviour))~Temperature+(1|ID),grname="ID",data=newdata,datatype="Gaussian",nboot=1000,npermut=1000)
  R_value[i] <- m1$R
  CI_value[i] <- m1$CI
}

P

您的代码正在运行。

示例中有一条错误消息,我已将其更改为
示例(唯一(数据$ID),4)
。然后它运行。您可能还可以添加replace,比如so
sample(unique(Data$ID),16,replace=TRUE)
,这样也可以。我还减少了
rboot
npermut
中的数字

库(rptR)
数据标识测试行为温度
#>1 A12.4.2一个8.64 4
#>2 A12.4.2两个7.63 5
#>3 A6.3.3一个1.81 3
#>4 A6.3.3两个2.47 9
#>5 B12.4.1一个1.17 12
#>6 B12.4.1两个3.96 2
#>7 E9.4.2一个13.04 13
#>8 E9.4.2两个9.51 6
P[[3]]
#> [1] 0.8324396
#> 
#> [[4]]
#> [1] 0.8324396
#> 
#> [[5]]
#> [1] 0.8324396
#> 
#> [[6]]
#> [1] 0.8324396
#> 
#> [[7]]
#> [1] 0.8324396
#> 
#> [[8]]
#> [1] 0.8324396
#> 
#> [[9]]
#> [1] 0.8324396
#> 
#> [[10]]
#> [1] 0.8324396
CI_值
#>  [1] 0.95 0.95 0.95 0.95 0.95 0.95 0.95 0.95 0.95 0.95

数据是什么样子的?。而且,这看起来像是几个小时前你的问题的重复?!这回答了你的问题吗?嗨,Wimpel,我的问题与几小时前发布的问题不同,因为我已经尝试了建议给我的代码,虽然它解决了我的一些问题,但我对这段代码有了新的问题。Lizzie,最好采用循环中的代码,而不是运行循环,设置
I Hi Allan,按照你的建议去做,我确实得到了预期的结果。