Question

我正在运行逻辑回归，并希望确保获得95％的置信区间。代码：

# Dissable scientific notation. 
# From: stackoverflow.com/questions/25946047
options(scipen=999)



###############################################################################


OR_CI_round_number<-5  # How many decimal point to keep after rounding OR and CI.


dfAPI <- haven::read_dta(
  file = "https://stats.idre.ucla.edu/stat/stata/faq/eyestudy.dta") 


dfAPI$carrot<-factor(dfAPI$carrot)
dfAPI$carrot<-relevel(dfAPI$carrot, ref = "1")


glmAPI = glm(lenses ~ carrot, data= dfAPI, family=(binomial(link = "log")))
#glmAPI
#summary(glmAPI)
round(exp(cbind(RR = coef(glmAPI), confint(glmAPI))), OR_CI_round_number)
round(exp(cbind(RR = coef(glmAPI), confint(glmAPI, level = 0.95))), OR_CI_round_number)

结果：

> round(exp(cbind(RR = coef(glmAPI), confint(glmAPI, level = 0.95))), OR_CI_round_number)
Waiting for profiling to be done...
                 RR   2.5 %  97.5 %
(Intercept) 0.41176 0.28349 0.54870
carrot0     1.58601 1.09250 2.40172

> round(exp(cbind(RR = coef(glmAPI), confint(glmAPI))), OR_CI_round_number)
Waiting for profiling to be done...
                 RR   2.5 %  97.5 %
(Intercept) 0.41176 0.28349 0.54870
carrot0     1.58601 1.09250 2.40172

我问的原因是因为我得到RR 2.5 % 97.5 %。据我了解，它们表示置信区间为95％的上下边界。这是正确的吗？

Answer 1

是的，正确的，这些是您的边界范围，在2.5％左边和97.5％右边是显着性水平。

此代码是否给我置信区间？ 95％

1 个答案: