расчет доверительного интервала коэффициента регрессии Пуассона

В Моем R-коде регрессия Пуассона выглядит следующим образом:

poissmod <- glm(aerobics$y ~ factor(aerobics$x1) + factor(aerobics$x2) + aerobics$x3 + aerobics$x4, family = poisson)
poissmod

Теперь мне нужно вычислить доверительный интервал для фактора aerobics$x1 (в модели без aerobics$x1, поскольку это несущественно).

Это может показаться очень простым, но я не знаком с R и нигде не могу найти ответ ...

Кто-нибудь может мне помочь?

Заранее большое спасибо!


r regression confidence-interval poisson
person user3387899    schedule 06.03.2014    source источник

Ответы (1)


arrow_upward
2
arrow_downward

См., Например, функция confint в пакете MASS (http://stat.ethz.ch/R-manual/R-devel/library/MASS/html/confint.html):

ldose <- rep(0:5, 2)
numdead <- c(1, 4, 9, 13, 18, 20, 0, 2, 6, 10, 12, 16)
sex <- factor(rep(c("M", "F"), c(6, 6)))
SF <- cbind(numdead, numalive = 20 - numdead)
budworm.lg0 <- glm(SF ~ sex + ldose - 1, family = binomial)
confint(budworm.lg0)
confint(budworm.lg0, "ldose")

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

Вот еще один пример из документации пакета stats для регрессии Пуассона (https://stat.ethz.ch/R-manual/R-devel/library/stats/html/confint.html):

## from example(glm)
counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3, 1, 9); treatment <- gl(3, 3)
glm.D93 <- glm(counts ~ outcome + treatment, family = poisson())
confint(glm.D93) # needs MASS to be present on the system
confint.default(glm.D93)  # based on asymptotic normality
person majom    schedule 06.03.2014
comment
Спасибо за ответ. Я уже пробовал это сделать, указав следующий R-код: confint (poissmod2, 'aerobics $ x1', level = 0.95). Однако я получаю сообщение об ошибке: аргумент имеет нулевую длину ... - person user3387899; 06.03.2014
comment
Добро пожаловать в SO. Команды в основном работают. Если с вашим кодом есть какие-то другие проблемы, я мог бы помочь, если бы вы предоставили воспроизводимый пример. См. stackoverflow.com / questions / 5963269 / - person majom; 06.03.2014
comment
Ну, код, похоже, у меня не работает ... Я все время получаю ошибку: «аргумент имеет нулевую длину». Значит, с моим кодом все еще что-то не так? - person user3387899; 06.03.2014
comment
Было бы полезно, если бы вы могли опубликовать воспроизводимый пример игрушки, где возникает эта ошибка. - person majom; 06.03.2014
comment
Извините, я плохо знаком с R и не понимаю, что вы имеете в виду, говоря о воспроизводимой игрушке :-) - person user3387899; 06.03.2014
comment
Пожалуйста, прочтите сообщение, на которое я указал в своем первом комментарии. В нем очень подробно описано, что вам следует делать. - person majom; 06.03.2014
comment
@ user3387899: confint(poissmod2) не дает вам значения для x1? Или попробуйте aerobics <- transform(aerobics, x1=factor(x1), x2=factor(x2)); poissmod <- glm(y ~ x1 + x2 + x3 + x4, data = aerobics, family = poisson); confint(poissmod, "x1") - person Ben Bolker; 10.03.2014