lam <- c(0.9, 0.9, 0.9, 0.9, 0.9, 0.9)

