df=data.frame(
Y=c(101,101,100,99,201,201,200,199,51,51,50,49),
X2=c("1","1","1",'1','2','2','2','2','3','3','3','3'),
X=c(1,1,1,1,2,2,2,2,3,3,3,3),
R=c(1,1,2,3,1,1,2,3,1,1,2,3)
)
# BLUP
# https://mattzheng.blog.csdn.net/article/details/51636011?utm_medium=distribute.pc_relevant.none-task-blog-2%7Edefault%7EBlogCommendFromBaidu%7Edefault-14.control&depth_1-utm_source=distribute.pc_relevant.none-task-blog-2%7Edefault%7EBlogCommendFromBaidu%7Edefault-14.control
library(lme4)
fm1 <- lmer(Y ~ 1 + (1|X) + (1|R), df)
summary(fm1)
ranef(fm1)
# BLUE
# https://cloud.tencent.com/developer/article/1549810
library(lme4)
fm2 <- lmer(Y ~ X2 + (1|R), df)
summary(fm2)
fixef(fm2)
