Dear all, I want to do a logistic regression. So far I've only found out how to do that in R, in a dataset of complete cases. I'd like to do logistic regression via max likelihood, using all the study cases (complete and incomplete). Can you help?
I'm using glm() with family=binomial(logit). If any covariate in a study case is missing then the study case is dropped, i.e. it is doing a complete cases analysis. As a lot of study cases are being dropped, I'd rather it did maximum likelihood using all the study cases. I tried setting glm()'s na.action to NULL, but then it complained about NA's present in the study cases. I've about 1000 unmatched study cases and less than 10 covariates so could use unconditional ML estimation (as opposed to conditional ML estimation). regards Desmond -- Desmond Campbell UCL Genetics Institute d.campb...@ucl.ac.uk Tel. ext. 020 31084006, int. 54006 ______________________________________________ R-help@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.
