2017-03-07 4 views
0

質問1は、ベクトル、他のベクトルの各要素に関数を適用

私はpnにかかる機能を持って考えます。私は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

私はnpの組み合わせこれに属するプロット追跡するpaste0()を使用して、プロットにnpの値を入れしようとしていますが、私は使用している場合、動作するようには思えません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

それは0.05のちょうど特異値です。 – ozarka

+1

'mapply'が並行して反復するので、他の組み合わせは作成されません。それらを生成するために 'expand.grid'または' rep'を使い、次に組み合わせを 'mapply'に渡します。 – alistaire

答えて

0
df<-expand.grid(p,n) 
mapply(Fun = ci.auto, df$Var1,df$Var2,alpha) 
関連する問題