modelsearch output

[fishpm]

Hi Klaus,
The function modelsearch in lava (v 1.4.7) appears to return each additional model path 3 times. Adjusted p-values appear sensitive to this redundancy.

m <- lvm(c(v1, v2, v3, v4) ~ u)
latent(m) <- ~u
dtempo <- sim(m, 1e2)
e <- estimate(m, dtempo)
esearch <- modelsearch(e)
esearch

Thanks,
Patrick

[kkholst]

The current implementation looks at all possible linear extensions of the model (both directed yx, xy, and un-directed/correlation x~~y). In some cases adding a direct path will not be equivalent (i.e., not lead to similar likelihood) to the model where an undirected path/correlation is added between two nodes. To only search among the undirected paths (the old behaviour) you can use the syntax:
modelsearch(e, type="cor").