Two-Level Hierarchical Linear Models Using R

https://stat.utexas.edu/images/SSC/documents/SoftwareTutorials/MultilevelModeling.pdf

R:

• Intercept-only Model (Unconditional Model) 仅截距

library(lme4)
library(languageR)
lmer(popular ~ 1 + (1|class))
pvals.fnc(my_intonly)

• Random Intercept with One Fixed Level-1 Factor (Non-Random Slope) 随机截距

my_1fixed <- lmer(popular ~ 1 + c_extrav + (1|class))
pvals.fnc(my_1fixed)

• Random Intercept and Slope for One Level-1 Factor 随机截距+随机斜率

my_1rnd <- lmer(popular ~ 1 + c_extrav + (1 + c_extrav|class))
p.values.lmer(my_1rnd)

• Random Slope for Two Level-1 Factors 随机斜率（2个第一层变量斜率随机）

my_2fixedrnd <- lmer(popular ~ 1 + c_extrav + c_sex + (1 + c_extrav + c_sex |class))
p.values.lmer(my_2rnd)

• One Level-2 Factor and Two Random Level-1 Factors (No Interactions) 第二层变量+2个随机第一层变量

my_model5 <- lmer(popular ~ c_extrav + c_sex + c_texp + (1 + c_extrav +
c_sex|class))
p.values.lmer(my_model5)

• One Level-2 Factor and Two Random Level-1 Factors with Interaction 第二层变量+2个随机第一层变量+交互项

my_model6 <- lmer(popular ~ c_extrav + c_sex + c_texp + c_extrav*c_texp +
c_sex*c_texp + (1 + c_extrav + c_sex|class))
p.values.lmer(my_model6)