Re: [R] variance functions in glmmPQL or glm?

From: Spencer Graves <spencer.graves_at_pdf.com>
Date: Sun 01 Oct 2006 - 23:14:01 GMT

      What are you trying to do with varIdent with logistic regression with variance components? Have you considered the weights argument in glm or glmmPQL? It sounds to me like you are trying to estimate 'size' for a binomial distribution, and I have trouble visualizing a context where that would be reasonable that would NOT be handled by the random effects portion of glmmPQL{MASS} or lmer{lme4}.

      If you'd like more help from this listserve, please provide commented, minimal, self-contained, reproducible code, as suggested in the posting guide "www.R-project.org/posting-guide.html".

      Hope this helps, even if it doesn't answer your question. 
      Spencer Graves

Gretchen wrote:
> Hello R users-
> I am new to R, and tried searching the archives and literature for an answer
> to this - please be patient if I missed something obvious.
>
> I am fitting a logistic regression model, and would like to include variance
> functions (specifically the varIdent function). I cannot figure out how to
> do this either in glmmPQL (or something similar) for the model with random
> effects, or in glm for the model without random effects. Is it possible to
> fit a varIdent function in a generalized linear model? If so, what are the
> appropriate packages/functions to use?
>
> Any help would be appreciated.
> Thank you,
> Gretchen Anderson
>
>
>
> M.Sc. Candidate
>
> Dept. of Fisheries and Wildlife
>
> Michigan State University
>
> 13 Natural Resources Building
>
> East Lansing, MI 48824
>
> Phone: (517) 353-0731
>
> E-Mail: ande1125@msu.edu
>
> ______________________________________________
> R-help@stat.math.ethz.ch 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.
>



R-help@stat.math.ethz.ch 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. Received on Mon Oct 02 09:28:15 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Mon 02 Oct 2006 - 00:30:10 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.