head(data_Bauer_Curran_2005)
# \donttest{
HSBmod <-nlme::lme(MathAch ~ Sector + CSES + CSES:Sector,
data = data_Bauer_Curran_2005,
random = ~1 + CSES|School, method = "ML")
summary(HSBmod)
REGIONS_OF_SIGNIFICANCE(model=HSBmod,
plot_title='Johnson-Neyman Regions of Significance',
Xaxis_label='Child SES',
Yaxis_label='Slopes of School Sector on Math achievement')
# }
Run the code above in your browser using DataLab