#-----------------------------------------------------------------------------------
## Nested Model with no intercept
#-----------------------------------------------------------------------------------
data(Table1.2)
Table1.2$Batch <- factor(x = Table1.2$Batch)
library(nlme)
Exam2.B.6fm1 <- lme(
fixed = Y ~ X
, data = Table1.2
, random = list(Batch = pdDiag(~1), X = pdDiag(~1))
, method = c("REML", "ML")[1]
)
Exam2.B.6fm1
library(broom.mixed)
tidy(Exam2.B.6fm1)
Run the code above in your browser using DataLab