# NOT RUN {
data( rat.blood )
cons <- list(order = "simple", decreasing = FALSE, node = 1 )
clme.out <- clme(mcv ~ time + temp + sex + (1|id), data = rat.blood ,
constraints = cons, seed = 42, nsim = 0)
BIC( clme.out )
BIC( clme.out, k=log( nobs(clme.out)/(2*pi) ) )
# }
Run the code above in your browser using DataLab