給定一個向量，將一個函數應用於另一個向量的每個元素

Question

問題1

我有一個接受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返回如下圖所示的圖表（ 編輯：信譽太低，無法上傳img，但我顯示僅獲得6個圖表 ）。 根據我的ablines位置，我認為它嘗試的是所有p值而不是單個n值！

問題2

我正在嘗試使用paste0()將n和p值放入繪圖中，以跟蹤哪個繪圖屬於n和p哪個組合，但是當我使用mapply時，它似乎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)

Answer 1

df<-expand.grid(p,n)
mapply(Fun = ci.auto, df$Var1,df$Var2,alpha)