次のサンプルデータセットがあります
Sites<-c("A","A","B","B","B","C","C","C")
TR<-c(0.1,0.4,0.4,0.3,0.3,0.6,0.5,0.5)
df<-data.frame(Sites,TR)
変数「サイト」で行をサブセット化し、各サイト トラップ レート (TR) をブートストラップしたいと思います。次のコードを使用してみました
require(iterators)
sites<-isplit(df,df$Sites)
nextElem(sites)
require(foreach)
mT<-foreach(Sites=sites)%do%
{
meanTraps <- vector()
for(i in 1:boots){
meanTraps[i] <- mean(sample(sites$value$TR, nrow(sites$value), replace = TRUE))
}
}
しかし、これは単に返されます
list()
理想的には、サイトごとに meanTraps を含むデータ フレームが必要です。誰でも私のエラーを修正するのを助けることができますか?