[ADMB Users] Values for random effects

H. Skaug hskaug at gmail.com
Fri Nov 12 23:20:50 PST 2010


Hi

Since U is a categorical variable (factor) with 7 levels. I suspect
that you want a random intercept, so that you should have used:

....  random = ~ 1, group="U", ....

This will give you one single variance parameter.

Does this help, or do I miss your point?

Hans

On Fri, Nov 12, 2010 at 8:08 PM, Schweizer, Peter E.
<schweizerpe at ornl.gov> wrote:
> Hello admb users,
>
>
>
> I just recently  ‘discovered’ glmm.admb and  would be interested to learn if
> admb in R can provide a  single value for random effects, similar to a
> ‘working correlation’ in SAS when using GEE.
>
>
>
> As of now I have
>
> M1<-glmm.admb( Y ~  V1 + V2 + V3 +V4 +V5, random = ~ U, group="U", data=A1,
> family="nbinom")
>
> In this case, U has 7 groups and in the model summary for above  model M1 I
> am getting (7) values for the diagonal matrix of the random effects. What I
> would like to obtain is a value that expresses the random effect for all
> groups in ONE single value, similar to a ‘working correlation’ in SAS when
> using GEE.
>
>
>
> Does this sound reasonable?
>
>
>
> Kind regards,
>
>
>
>
>
> Pete
>
> _______________________________________________
> Users mailing list
> Users at admb-project.org
> http://lists.admb-project.org/mailman/listinfo/users
>
>



More information about the Users mailing list