prop.test отображать доверительный интервал в процентах вместо десятичного числа в R? - PullRequest
0 голосов
/ 28 мая 2020

Мне интересно, есть ли способ изменить параметры вывода функции prop.test в R, чтобы он отображал доверительный интервал уже в процентах, а не в десятичном? Например, я пытаюсь найти 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%] вместо десятичных знаков? Спасибо!

Ответы [ 2 ]

1 голос
/ 28 мая 2020

Наверное, будет проще:

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 
0 голосов
/ 28 мая 2020

Непросто. Форматом печати управляет функция print.htest(), которая задокументирована в ?print.htest: похоже, она не предлагает никаких параметров, кроме количества цифр и префикса для компонента «метод».

При желании вы можете самостоятельно взломать функцию. (Возможно, последует более подробная информация; выгрузите 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
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...