####=========================================####
#### For CRAN time limitations most lines in the
#### examples are silenced with one '#' mark,
#### remove them and run the examples using
#### command + shift + C |OR| control + shift + C
####=========================================####
# you need to install the orthopolynom library to do random regression models
# library(orthopolynom)
# data(DT_legendre)
# DT <- DT_legendre
# mRR2<-mmer(Y~ 1 + Xf
# , random=~ vsr(usr(leg(X,1)),SUBJECT)
# , rcov=~vsr(units)
# , data=DT)
# summary(mRR2)$varcomp
#
# mRR2b<-mmec(Y~ 1 + Xf
# , random=~ vsc(usc(leg(X,1)),isc(SUBJECT))
# , rcov=~vsc(isc(units))
# , data=DT)
# summary(mRR2b)$varcomp
Run the code above in your browser using DataLab