2016-06-17 3 views
-2

Я пытаюсь создать таблицу, которая проведет определение размера выборки из нескольких методов для разных начальных вероятностей. Я хотел бы, чтобы выглядеть следующийЗаполнение таблицы в R

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» для таблицы.

+0

Итак ... что, в частности, ваш вопрос? –

+0

@Alex Мне нужна помощь в преобразовании вычислений размера выборки в таблицу для каждого приращения начальной вероятности. В основном, я пытаюсь использовать вторую часть кода, чтобы заполнить «xxx» для таблицы. Лучше? –

+0

@ a.powell Лучше изменить свой вопрос, чем дать пояснения в комментариях. – lmo

ответ

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)) 
Смежные вопросы