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

sens <-seq(from=.1, to=.5, by=0.05)
enroll <-seq(from=20, to=200, by=20)
df <-expand.grid(sens=sens,enroll=enroll)
df <-transform(df,succes=sens*enroll)

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

Меня интересует только верхний предел 95% доверительного интервала биномиального теста. Я хочу, чтобы это единственное число было добавлено к фрейму данных в виде столбца с именем «upper.limit»

Я думал о чем-то вроде

binom.test(succes,enroll)$conf.int    

увы, conf.int дает что-то вроде

[1] 0,1266556 0,2918427
атр ( "conf.level")
[1] 0,95

Все, что я хочу, это просто 0.2918427

Кроме того, у меня есть ощущение, что где-то должен быть do.call и, может быть, даже неудачно, но я не знаю, как это пройдет через весь фрейм данных. Или, может быть, я должен использовать Plyr?

Ясно, что моя голова кружится. Пожалуйста, остановите это.