[R-sig-ME] modelling proportions, with aggregated data, and the new/old lme4

Ben Bolker bbolker at gmail.com
Mon Mar 19 02:40:53 CET 2012


Joerg Luedicke <joerg.luedicke at ...> writes:

> I would certainly check out a Poisson model with the number of
> successes as outcome and successes+failures as an offset.

  That seems odd to me; the Poisson+offset model
should be appropriate when p<<1 (i.e. for very small p,
the Poisson variance mu=n*p is approximately the same as
the binomial variance n*p*(1-p); in this case p is not small.

  Of course, I could be wrong.

> 
> On Sun, Mar 18, 2012 at 11:51 AM, Rolf Turner <r.turner at ...> wrote:
> > On 19/03/12 07:10, Ben Bolker wrote:
> >
> > <SNIP>
> >
> >
> >>> head(dat)
> >>>   successes failures id sex subdist
> >>> 1       560      726  1   F       4
> >>> 2       844      510  1   M       4

  [snip]

> >>>
> >>> In community #1, which is in subdistrict #4, there are 560 women
> >>
> >>   you mean 510, right?
> >
> > <SNIP>
> >
> > Sure looks like 560 to me.  Time for a trek to the optometrist, Ben?

  Looked at the wrong line (failures in men rather than successes in women).




More information about the R-sig-mixed-models mailing list