2

R または Stata ソフトウェアでブートストラップを行う際に助けが必要です。政策の有効性などについて「はい」と「いいえ」と答えた人の割合を計算したい

Stataにはこのコードがあります

bs "summarize y1" "r(mean)", reps(200) size(770)

r(mean)割合を見積もるには、何の値が必要ですか?

また、Rには次のコードがあります。

test <- function (q13){
    test13 <- table(q13)
    rel_freq <- test13/sum(test13)
    return(rel_freq)
      }

results <- boot(data=q13, statistic=test,
                R=200)

コードを修正するにはどうすればよいですか? エラーが発生します

統計 (データ、オリジナル、...) のエラー: 未使用の引数 (オリジナル)

4

1 に答える 1

3

Stataではproportion、変数に 2 つ以上のカテゴリがある場合に使用できます。

//サンプルデータ

sysuse auto, clear
keep if (headroom==2.0 |headroom==2.5)
gen prop=.
replace prop=0 if headroom==2.0
replace prop=1 if headroom==2.5

// 0 は「はい」、1 は「いいえ」とします

set seed 123
bootstrap _b, reps(100):proportion prop

@Nick に従って更新: バイナリ変数の場合、次で十分です

bootstrap r(mean), reps(100): summarize prop, meanonly

................................................................... ................................................................... ................................................................... ...................................

Rでは、bootパッケージと mtcars データを使用して次のように実行できます。

library(boot)
set.seed(123)
x<-mtcars$vs
myprop<-function(x,i){
sum(x[i]==0)/length(x)
}

bootprop <- boot(x,myprop,100)
于 2013-08-14T18:36:23.467 に答える