在 R 中填充表答案

【问题标题】：Populating a table in R在 R 中填充表
【发布时间】：2016-12-30 17:56:22
【问题描述】：

我正在尝试制作一个表格，其中包含针对不同初始概率的几种方法确定的样本量。我希望它看起来像下面的样子

Prob.    .80    .81    .82 ...
Wald     xxx    yyy    zzz
AC       qqq    rrr    sss
Jeff.    ddd    uuu    iii

字母是我想根据间隔类型和概率值填充的实际输出。为了获得单独的价值，我可以使用以下代码，但我很难将它们组合在一起。

library(binom)
n_wald <- ciss.wald(p, d, alpha = a)
n_agricoull <- ciss.agresticoull(p, d, alpha = a)
n_jeffreys <- ciss.binom(p, d, alpha = a, prior.shape1 = 0.5,
 prior.shape2 = 0.5, method = "bayes")

我需要帮助将样本量计算转换为初始概率的每个增量的表格。基本上，我正在尝试使用第二段代码来填写表格的“xxx”。

【问题讨论】：

那么...您的问题具体是什么？
@Alex 我需要帮助将样本量计算转换为初始概率的每个增量的表格。基本上，我正在尝试使用第二段代码来填写表格的“xxx”。更好？
@a.powell 最好修改您的问题，而不是在 cmets 中进行澄清。
@Imo 我这样做了。希望它有助于回答这个问题。

标签： r populate confidence-interval

【解决方案1】：

您引用的函数实际上在binomSamSize 包中。它们是矢量化的，这使得这更容易一些，除非您为多个参数提供多个值。

我使用dplyr 只是为了更容易将这些推送到数据框中。

library(binomSamSize)
library(dplyr)
p0_vec <- seq(0.8, 0.9, by = 0.02)

rbind(
  ciss.wald(p0=p0_vec, 
            alpha=0.1, 
            d=0.05),
  ciss.agresticoull(p0=p0_vec, 
                    alpha=0.1, 
                    d=0.05),
  vapply(p0_vec,
         ciss.binom,
         numeric(1),
         alpha=0.1, 
         d=0.05,
         prior.shape1 = 0.5,
         prior.shape2 = 0.5, 
         method = "bayes")
) %>%
  as.data.frame() %>%
  bind_cols(
    data.frame(type = c("Wald", "AC", "Jeff"), 
               stringsAsFactors = FALSE),
    .
  ) %>%
setNames(c("type", p0_vec))

【讨论】：