我有一个rpy2脚本,给我95%置信区间的覆盖率。
from rpy2.robjects.packages import importr
binom = importr('binom')
from rpy2 import robjects
robjects.r('''library(binom)
p = seq(0,1,.01)
coverage = binom.coverage(p, 10, method="bayes", type = "central")$coverage
''')
我想得到99%的置信区间的覆盖率,请问:
也许是这种方式?
robjects.r('''library(binom)
p = seq(0,1,.01)
coverage = binom.coverage(p, 10, conf.level = 0.99,
method="bayes", type = "central")$coverage
''')
请看 文件.