sem.r

options <- commandArgs(trailingOnly = TRUE)
nobs <- as.numeric(options[1])

library(lavaan)
model <- readLines("./tmp/model.lav")
elems <- readLines("./tmp/elems.lav")

cov <- getCov(elems[1], names=strsplit(elems[2], " ", fixed = TRUE)[[1]])
fit <- sem(model, sample.cov=cov, sample.nobs=nobs)
summary(fit, fit.measures=TRUE, standardized=TRUE)
fitMeasures(fit, fit.measures = "all", baseline.model = NULL)