prop.test 在 R 中以百分比而不是小数形式显示置信区间？

Question

我想知道是否有办法更改 R 中 prop.test 函数的输出设置，以便它以百分比而不是小数形式显示置信区间？例如，我试图找到西方移民中患有糖尿病的比例的 95% 置信区间。这是我的代码和输出：

    sum(Immigrant_West$DIABETES)= 8, nrow(Immigrant_West)=144
prop.test(x=sum(Immigrant_West$DIABETES),n=nrow(Immigrant_West),conf.level = .95,correct=TRUE)

    > 1-sample proportions test with continuity correction
    data:  sum(Immigrant_West$DIABETES) out of nrow(Immigrant_West), null probability 0.5 
   X-squared = 112, df = 1, p-value <2e-16
alternative hypothesis: true p is not equal to 0.5
95 percent confidence interval:
 0.02606 0.11017
sample estimates:
      p 
0.05556

那么有没有办法将置信区间输出更改为显示 [2.606%, 11.017%] 而不是小数？谢谢！

Answer 1

不容易。打印格式由 print.htest() 函数控制，该函数记录在 ?print.htest 中：它似乎没有提供除 "method" 组件的位数和前缀之外的任何选项.

如果你愿意，你可以自己破解这个功能。（可能会有更多细节；转储 stats:::print.htest 到文件，然后编辑它，然后 source() 它。

正如@CarlWitthoft 所建议的：

p <- prop.test(x=8,n=144)
str(p) ## see what's there
p$estimate <- p$estimate*100
p$conf.int <- p$conf.int*100

Answer 2

这可能会更简单：

prp.out <- prop.test(x=8, n=144, conf.level=.95, correct=TRUE)
prp.out$conf.int <- prp.out$conf.int * 100
prp.out
# 
#   1-sample proportions test with continuity correction
# 
# data:  8 out of 144, null probability 0.5
# X-squared = 112.01, df = 1, p-value < 2.2e-16
# alternative hypothesis: true p is not equal to 0.5
# 95 percent confidence interval:
#   2.606172 11.016593
# sample estimates:
#          p 
# 0.05555556

prop.test 在 R 中以百分比而不是小数形式显示置信区间？

prop.test display the confidence interval in percentage terms instead of decimal in R?

testing

r

confidence-interval