topic binomial unequal group sizes in Statistical Procedures

binomial unequal group sizes

DDK — Fri, 16 Oct 2015 16:20:54 GMT

Dear all,

If you have binomial data for which the group sizes differ, how do you correct for it in proc glimmix?

For example:

proc glimmix data=name;

class var1 var2;

model y/n=var1 var2 / dist=binomial link=logit;

random randomvariable;

run;

y is for example the number of affected individuals per group that differs in size (n). I have been told that you can correct for this by adding an overdispersion parameter 'random _residual_;' but I am not too confident about this. Anyone have any ideas? Thanks in advance.

Re: binomial unequal group sizes

lvm — Fri, 16 Oct 2015 19:03:40 GMT

You don't need to correct for this -- it is automatically taken care of. Do not use the 'random _residual_" statement for this purpose (this is for overdispersion: high variability than one would obtain with a binomial distribution).

Re: binomial unequal group sizes

JacobSimonsen — Fri, 16 Oct 2015 19:26:52 GMT

Agree. And further, since the random-statement is not neccessary, it may be better to use a procedure (proc logistics) that is specialized for logistic regression.

Re: binomial unequal group sizes

DDK — Fri, 16 Oct 2015 21:27:29 GMT

Thanks, thought it was a bit weird to correct for unequal sizes using an overdispersion parameter. Good to hear a confirmation that it is not necessary to do so.