私は推論パッケージでブートストラップしています。 対象となる統計は平均であり、サンプルデータは3列と5行のティブルで提供されます。私の本当のチブルは86行と40列です。すべての列に対して、tibble "test_tibble"の列 "x"について以下に示すように、ブートストラップシミュレーションを実行したいと思います。
<前>ウィズウィズ私は今、自分のティブルの他のコラムに対して同じことをする方法を探しています。私はこのようなforループを試しました:
<前>ウィズウィズ残念ながら、これは次のエラーを返します
<前>ウィズウィズ「応答」の引数として手動で入力せずに、列x、y、zを反復する方法はありますか?これは、40列ではかなり退屈です。
library(infer)
library(tidyverse)
test_tibble <- tibble(x = 1:5, y = 6:10, z = 11:15)
# A tibble: 5 x 3
x y z
<int> <int> <int>
1 1 6 11
2 2 7 12
3 3 8 13
4 4 9 14
5 5 10 15
specify(test_tibble, response = x) %>%
generate(reps = 100, type = "bootstrap") %>%
calculate(stat = "mean") %>%
summarise(
lower_CI = quantile(probs = 0.025, stat),
upper_CI = quantile(probs = 0.975, stat)
)
# A tibble: 1 x 2
lower_CI upper_CI
<dbl> <dbl>
1 2.10 4
- 解決した方法 # 1
これはトリッキーな質問で、トリッキーな答えです。
response
を見てくださいspecify
の引数 ドキュメントの機能:これを念頭に置いて、プロセスを自動化するようにコードを変更し、元のデータフレームに列を1つ追加し、
を使用しました 1の列を説明変数として使用して、同じ結果を取得する引数。 <前>ウィズウィズそれが役に立てば幸い
formula