[R] modeling logit(y/n) using lrm

Kevin E. Thorpe kevin.thorpe at utoronto.ca
Fri Jun 16 19:39:15 CEST 2006

Hamilton, Cody wrote:
> I have a dataset at a hospital level (as opposed to the patient level)
> that contains number of patients experiencing events (call this number
> y), and the number of patients eligible for such events (call this
> number n).  I am trying to model logit(y/n) = XBeta.  In SAS this can be
> done in PROC LOGISTIC or GENMOD with a model statement such as: model
> y/n = <predictors>;.  Can this be done using lrm from the Hmisc library
> without restructuring the dataset so that for each hospital there is one
> row with y = 1 and one row with y = 0 and then using the weight option
> in lrm to weight these two responses by the number of 'successes' and
> 'failures' for that hospital, respectively?  I would like to avoid the
> restructuring, and I understand that the use of the weight function is
> not compatible with a lot of the validation functions available in Hmisc
> (validate, bootcov, etc.).

I don't know about lrm, but for glm you can do
glm(cbind(y,m)~ ...) where y is number of successes and
m is the number of failures.

So, you might try that.

> Cody Hamilton, Ph.D
> Institute for Health Care Research and Improvement
> Baylor Health Care System
> (214) 265-3618

Kevin E. Thorpe
Biostatistician/Trialist, Knowledge Translation Program
Assistant Professor, Department of Public Health Sciences
Faculty of Medicine, University of Toronto
email: kevin.thorpe at utoronto.ca  Tel: 416.946.8081  Fax: 416.946.3297

More information about the R-help mailing list