[R] zero cells in one variable in logistic regression

Vito Muggeo (UniPa) vito.muggeo at unipa.it
Mon Jul 13 16:03:41 CEST 2009


dear anna,
if you are not interested in point estimate and SE of the parameter of 
the aforementioned categorical variable, I believe the conventional 
glm(..,family=binomial) is correct. In particular, the returned deviance 
is reliable and also it is the relevant likelihood ratio test..

hope this helps,
vito


David Winsemius ha scritto:
> 
> On Jul 13, 2009, at 5:37 AM, anna.bucharova wrote:
> 
>>
>> Dear all.
>> I am sort of beginner with R. I do logistic regression with binomial
>> response variable and several continuous and categorical variables. In 
>> one
>> categorical variable, zero cell occures (2x2 table looks like
>> 7 - 0
>> 23 - 25
>> This leads to overestimating of odds ratio and inflated confidence 
>> interval
>> for odds for given variable. The variable is significant in univariate 
>> test.
>> I do not necessarilly need odd ratio, but I need the explained 
>> deviance by
>> this variable and I really want to keep this variable in the model. It
>> probably matters for explained deviance. How to treat this problem?
>> Thanks for help, Anna Bucharova
>> -- 
> 
> You might consider glmrob in package:robustbase. See 
> http://www.jstatsoft.org/v10/i04
> 
> David Winsemius, MD
> Heritage Laboratories
> West Hartford, CT
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide 
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
> 

-- 
====================================
Vito M.R. Muggeo
Dip.to Sc Statist e Matem `Vianelli'
Università di Palermo
viale delle Scienze, edificio 13
90128 Palermo - ITALY
tel: 091 6626240
fax: 091 485726/485612
http://dssm.unipa.it/vmuggeo




More information about the R-help mailing list