search for: x_hijk

Dear all, I have a somewhat unusual linear mixed model that I can't seem to code in lme. It's only unusual in that one random effect is applied only to some of the observations (I have an indicator variable that specifies which observations have this random effect). The model is: X_hijk = alpha_h + h * b_i + r_(ij) + e_hijk , where h = 0 or 1 (indicator) i = 1, ..., N j = 1, ..., n_i k = 1, ..., K alpha is fixed, and the rest are random. I'm willing to assume b, r, and e are mutually independent and normal with var(b) = sigma^2_b, var(r) = sigma^2_r, and var(e) = sigm...

Dear all, I have a question about mixed effect model in R. The data set has 5 variables, X(response),subject, times, repeat, indicator The model is X_hijk=a_h+Z_h*b_i+r(ij)+e_hijk , where h=0,1(indicator), i=1,...,n(subject), j=1,...,n_i(times within subject; nested effect),k=1,2,3(repeat). Z_h=1 if h=1 =0 if h=0 b_i~N(0,c^2) random effect of subject r(ij)~N(0,d^2) random effect of times within subject e_hijk~N(0,e^2) error term which is independe...