[R] Help to conduct a random factor analysis with binomial response

Leif Egil Loe l.e.loe at bio.uio.no
Tue Nov 20 09:28:13 CET 2001


Dear users of the R mailing list,

I am a ph.d. student in biology working on red deer in Norway, who would 
like to conduct an analysis with random factor where the response is 
binomially distributed. This cannot be conducted in S-plus, and I was told 
by others that it may be possible in R. However, I soon got into trouble 
which I hope you can help me to solve.

My model is on this form:
cbind(success,1-success) ~ variable 1 (numeric) + variable 2 (numeric) + 
variable 3 (factor) + variable 4 (factor). I have only one random factor, 
variable 5 (factor).

We assume that the analysis to be used is a non-linear mixed effects model, 
and hence use the nlme function in R. However, we failed to specify this 
model in a correct way. My question is therefore:

What is the syntax for this specific model (we got especially confused by 
the "start" argument)? If anyone could write the entire formula for this 
specific model we would be very greatful!


Thank you!


With regrads, Leif Egil Loe


.................................................................................................................................................
Leif Egil Loe                                        Privatadresse:
Avdeling for Zoologi                               c/o Hafting
Biologisk institutt                                 Hauges vei 31
Universitetet i Oslo                               N-0871 Oslo
Postboks 1050 Blindern
N-0316 Oslo

Tlf.    +47 22 85 72 94
Fax.  +47 22 85 46 05
Mob. +47 91 77 33 40
E-post. l.e.loe at bio.uio.no

.................................................................................................................................................

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list