0
私はp
とn
にかかる機能を持って考えます。私はp=c(0.05,0.1,0.25,0.5,0.9,0.95)
とn=c(5,10,30,50,100)
のすべての組み合わせをグラフにする必要があります。合計で6 * 5 = 30のプロットです。私はmapply
を使って試したところ、以下のプロットが返ってきました(EDIT:imgをアップロードするには評判が低すぎますが、私は6プロットだけをと表示しています)。私のablines
がどこにあるかに基づいて、私はそれがn
の単一の値にだけp
のすべての値を試していると思います!
質問2
私はn
とp
の組み合わせこれに属するプロット追跡するpaste0()
を使用して、プロットにn
とp
の値を入れしようとしていますが、私は使用している場合、動作するようには思えませんmapply
(それは私が手動でプロットし、単一のグラフのために正常に動作します)
ci.auto <- function(p,n,alpha){
# Repeat the process nsim times
nsim <- 10000
ci.mat <- replicate(nsim, conf.int(p, n, alpha))
# Graphing the first 100 intervals
matplot(rbind(1:100, 1:100), ci.mat[, 1:100], type = "l", lty = 1,
xlab = paste0("sample number n=",n), ylab = "confidence interval")
abline(h = p)
text(2.5, 0.8, paste0("p=",p), col = "red")
# Proportion of times the interval is correct
mean((p >= ci.mat[1,])*(p <= ci.mat[2,]))
}
par(mfrow = c(5,6))
mapply(ci.auto, p = c(0.05,0.1,0.25,0.5,0.9,0.95), n = c(5,10,30,50,100), alpha)
それは0.05のちょうど特異値です。 – ozarka
'mapply'が並行して反復するので、他の組み合わせは作成されません。それらを生成するために 'expand.grid'または' rep'を使い、次に組み合わせを 'mapply'に渡します。 – alistaire