at distribution のシミュレーションを実行した後、R で t テーブルから行を作成するにはどうすればよいですか? 本質的に、信頼レベルを入力として使用するのではなく、正規分布からランダムなサンプルから計算された値のみを使用して qt 関数を実行したいと考えています。
at(自由度 "df") 分布で 5000 回シミュレーションを実行しました。これらの値を使用して、df に対応する t テーブルの線を再作成します。それを理解しようとするのに苦労しています。
どんな助けでも大歓迎です！
分位関数は、私が使用した t テーブルの値と同等の値を提供します。
quants <- 1-c(0.2, 0.1, 0.05, 0.025, 0.01, 0.005, 0.001)
n <- 15
sims <- replicate( 5000, {x <- rnorm(n); sqrt(n)*mean(x)/sd(x)} )
(t.table <- quantile(sims, quants) )
## compare
rbind( t.table, qt( quants, n-1 ) )
pt(t.table, n-1)
これらは左側の領域であり、一般的な t テーブルは右側の領域を示しています。これと一致するようにラベルを変更できますnames(t.table) <- 1-quants。特定のテーブルに一致するように変更quantsします (より厳密な一致が必要な場合は、シミュレーションの数を増やします)。
この場合のシミュレーションは、法線のマトリックス全体を 1 ステップでシミュレートしてから、 などのツールを使用することで、速度を上げることができますrowMeans。