2016-10-25 7 views
0

rのデータフレームの行に単純なブートストラップを実行しようとしています。ここまでは私がこれまでに取り組んだことはありますが、私は行き止まりに打ち勝っています。rのデータフレームの行ごとにブートストラップを行う

x1 <- c(1:5) 
x2 <- c(6:10) 
y <- runif(5) 
z <- as.data.frame(rbind(x1, x2, y)) 

trial <- 10 
avg <- rep(0, trial) 
for(i in 1:trial){ 
    ind <- sample(ncol(z), size = ncol(z), replace = TRUE) 
    z.boot <- z[ind, ] 
    mean[i] <- mean(z.boot) 
} 
mean 

理想的には、私がやりたいことは3行目の重みで1行目と2行目のためのブートストラップ加重平均を得ることですが、私も私のループを動作させることはできません。これを行うにはよりよい方法が必要です。すべてのヘルプは...私はかなりの加重平均についてのあなたのポイントを得ることはありません...しかし、あなたは多分、ここからそれをうまくすることができ、これを試してください

+0

あなたは 'のNcoI(Zからのサンプル) 'が、行ごとにサブセットを作成します。また、 'data.frame'(定義されていない)の'平均値 'も見つけようとします。あなたは何をしようとしているのですか? – nicola

+0

別のオプションは、重みを確率として使用し、中央値を探すことです。本当に、私はデータフレームの行に適用されたこの手法を学ぼうとしています。 – mike

答えて

0

を高く評価されています

n= seq(100, 500, 50)  
bootdata=list() 
for (i in 1:length(n)) { 
    bootdata[[i]]=data[sample(nrow(data), n[i], replace=TRUE), ] 
} 
bootdata 
str(bootdata[[1]]) 
関連する問題