[R-sig-ME] Random intercept model- unbalanced cluster
bbolker @ending from gm@il@com
Mon Oct 29 22:23:16 CET 2018
In principle lme4 shouldn't have problems with a subset of groups that
have only one observation (although clearly the model will get more
fragile/unreliable the less information is available about within vs
among group variation ...). I'd expect the random effects for groups
with only one observation to be strongly shrunk toward the population
mean ... if in doubt, it can be very useful to simulate a situation
similar to your real data set to see what happens in cases where you
know the real answer ...
On 2018-10-29 5:13 p.m., Yashree Mehta wrote:
> Or is there an alternative method of modeling this subset of households who
> only own one plot?
> thank you,
> On Thu, Oct 25, 2018 at 6:00 PM Yashree Mehta <yashree19 using gmail.com> wrote:
>> I am working with a random intercept model on a cluster dataset (Repeated
>> measurements of plots per household). I have the usual "X" vector
>> of covariates and one id variable which will make up the random
>> intercept. For example,
>> Response variable: Production of maize
>> Covariates: Size, input quantities, soil fertility dummies etc..
>> ID variable: Household_ID
>> However, about 40% of the households own one plot. The number of plots per
>> household ranges from 1 to 13.
>> When I estimated the random intercept model using lmer, I can extract a
>> random intercept for all households, irrespective of their number of plots.
>> How does lmer treat these households with just 1 plot? Also, is it
>> theoretically correct to include these observations ?
>> Thank you,
> [[alternative HTML version deleted]]
> R-sig-mixed-models using r-project.org mailing list
More information about the R-sig-mixed-models