experiments/Orthodont.R

In davidbolin/ngme: Linear Mixed Effects Models With Flexible Distributions

library(ngme)
data(Orthodont)
set.seed(123)
NormalMVD <- ngme( fixed       = distance ~ age,
                   random      = ~ 1|id,
                   data        = Orthodont,
                   reffects    = "tdist",
                   use.process = F,
                   silent      = T,
                   nIter       = 2,
                   controls    = list(estimate.fisher = FALSE,
                                      subsample.type  = 0,
                                      nSim  =2,
                                      nBurnin = 2))
print(NormalMVD$fixed_est)