On 09/27/2011 07:53 AM, majesty wrote:
Dear subscribers,

I am looking for a function which would allow me to model the dependent
variable as the number of successes in a series of Bernoulli trials. My data
looks like this

ID  TRIALS     SUCCESSESS   INDEP1 INDEP2 INDEP3
1      4444        0                0.273     0.055  0.156
2      98170      74               0.123     0.456  0.789
3      145486    30               0.124     0.235  0.007
4      147149    49               0.888     0.357  0.321
5      60585      11               0.484     0.235  0.235
6      198953    43               0.295     0.123  0.856

I want to find out how independent variables influence the number of
successful trials (dependent variable)
I had tried to use glm formula

regression.glm<- glm( SUCCESSESS ~ INDEP1 + INDEP2 + INDEP3, data = data,
family = binomial, weights= TRIALS)

Try

regression.glm<- glm(cbind(SUCCESSESS,TRIALS) ~ INDEP1 + INDEP2 + INDEP3, data = data,family = binomial)

as is specified in the details of ?glm.

--
Patrick Breheny
Assistant Professor
Department of Biostatistics
Department of Statistics
University of Kentucky

______________________________________________
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.

Reply via email to