Question

我正在使用for循环对我的数据子集执行可重复性分析。我正在努力寻找一个将for循环结果存储在数据帧中的函数，因为可重复性分析输出给出了R，标准误，置信区间和p值。我只想存储R和CI值。

这是我的代码：

p<-10000

for(i in 1:P){
  newdf<-df[df$ID %in% sample(unique(df$ID), 16), ]
  m1<-rpt(Behaviour~Temperature+(1|ID),grname="ID",data=newdf,datatype="Gaussian",nboot=1000,npermut=1000)
}

有人可以帮忙吗？

Answer 1

您可以创建一个数字矢量来存储R和CI值：

p<-10000
R_value <- numeric(length = p)
CI_value <- numeric(length = p)

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

从for循环存储多个结果

1 个答案: