Thank you, this is indeed the case in my data. Case solved, then!

 Dr. Iasonas Lamprianou
Department of Social and Political Sciences
University of Cyprus




________________________________
From: Dennis Murphy <djmuser@gmail.com>

Cc: r-sig-mixed-models@r-project.org
Sent: Thu, 2 June, 2011 16:07:46
Subject: Re: [R-sig-ME] glmer problems

Hi:

If the regions are intrinsically nested within country, then the
second specification of your model is preferred. 'Intrinsic nesting'
in this case means that the region names only occur in combination
with exactly one country. If that is the case in your data, the second
model specification should work.

HTH,
Dennis


rote:
> Dear all. I am running this model but I get an error.
>
>> gmm1a <- glmer(workstatus ~ 1+(1|cntry) + (1|cntry:region),family=binomial,
>>data=retirement2)
> Warning in region:cntry :
>  numerical expression has 7011 elements: only the first used
> Warning in region:cntry :
>  numerical expression has 7011 elements: only the first used
> Warning in inherits(x, "factor") : NAs introduced by coercion
> Error in region:cntry : NA/NaN argument
>
>
>
>
> if I run a slightly different model, then it works
>
>> gmm1a <- glmer(workstatus ~ 1+(1|cntry) + (1|region),family=binomial,
>>data=retirement2)
>
>
> If I run the command below, it seems that there are no missing data either in
> the cntry or in the region variables. Can anyone offer any hints please?
>
> sum(table(retirement2$region, retirement2$cntry))
>
>
>
> Dr. Iasonas Lamprianou
> Department of Social and Political Sciences
> University of Cyprus
>
>
>
>
> ________________________________

	[[alternative HTML version deleted]]


