lam <- c(0.8, 0.8, 0.8, 0.8)
