Skip to contents

Model searching

Source: R/modelsearch.R
modelsearch.Rd

Performs Wald or score tests

Usage

modelsearch(x, k = 1, dir = "forward", type = "all", ...)

Arguments

x

lvmfit-object

k

Number of parameters to test simultaneously. For equivalence the number of additional associations to be added instead of rel.

dir

Direction to do model search. "forward" := add associations/arrows to model/graph (score tests), "backward" := remove associations/arrows from model/graph (wald test)

type

If equal to 'correlation' only consider score tests for covariance parameters. If equal to 'regression' go through direct effects only (default 'all' is to do both)

...

Additional arguments to be passed to the low level functions

Value

Matrix of test-statistics and p-values

See also

compare, equivalence

Author

Klaus K. Holst

Examples


m <- lvm();
regression(m) <- c(y1,y2,y3) ~ eta; latent(m) <- ~eta
regression(m) <- eta ~ x
m0 <- m; regression(m0) <- y2 ~ x
dd <- sim(m0,100)[,manifest(m0)]
e <- estimate(m,dd);
modelsearch(e,messages=0)
#>  Score: S P(S>s)    Index  holm      BH       
#>  3.221    0.07272   y1~~x  0.686     0.07272  
#>  3.221    0.07272   y1~x   0.686     0.07272  
#>  3.221    0.07272   x~y1   0.686     0.07272  
#>  3.221    0.07272   y2~~y3 0.686     0.07272  
#>  3.221    0.07272   y2~y3  0.686     0.07272  
#>  3.221    0.07272   y3~y2  0.686     0.07272  
#>  3.618    0.05716   y1~~y2 0.686     0.07272  
#>  3.618    0.05716   y1~y2  0.686     0.07272  
#>  3.618    0.05716   y2~y1  0.686     0.07272  
#>  3.618    0.05716   y3~~x  0.686     0.07272  
#>  3.618    0.05716   y3~x   0.686     0.07272  
#>  3.618    0.05716   x~y3   0.686     0.07272  
#>  22.9     1.702e-06 y2~~x  3.064e-05 5.107e-06
#>  22.9     1.702e-06 y2~x   3.064e-05 5.107e-06
#>  22.9     1.702e-06 x~y2   3.064e-05 5.107e-06
#>  22.9     1.702e-06 y1~~y3 3.064e-05 5.107e-06
#>  22.9     1.702e-06 y1~y3  3.064e-05 5.107e-06
#>  22.9     1.702e-06 y3~y1  3.064e-05 5.107e-06
modelsearch(e,messages=0,type="cor")
#>  Score: S P(S>s)    Index  holm      BH       
#>  3.221    0.07272   y1~~x  0.2287    0.07272  
#>  3.221    0.07272   y2~~y3 0.2287    0.07272  
#>  3.618    0.05716   y1~~y2 0.2287    0.07272  
#>  3.618    0.05716   y3~~x  0.2287    0.07272  
#>  22.9     1.702e-06 y2~~x  1.021e-05 5.107e-06
#>  22.9     1.702e-06 y1~~y3 1.021e-05 5.107e-06