# NOT RUN {
mSim <- lvm(Y~G+X1+X2)
addvar(mSim) <- ~Z1+Z2+Z3+Z4+Z5+Z6
df.data <- lava::sim(mSim, 1e2)
mBase <- lvm(Y~G)
addvar(mBase) <- ~X1+X2+Z1+Z2+Z3+Z4+Z5+Z6
e.lvm <- estimate(mBase, data = df.data)
res.x <- modelsearch2(e.lvm, statistic = "score", method.p.adjust = "holm", nStep = 2)
res.y <- modelsearch2(getStep(res.x, slot = "sequenceModel"),
statistic = "score", method.p.adjust = "holm")
res.xy <- merge(res.x,res.y)
modelsearch2(e.lvm, statistic = "score", method.p.adjust = "holm")
# }
Run the code above in your browser using DataLab