mod <- psem(
lm(rich ~ cover, data = keeley),
lm(cover ~ firesev, data = keeley),
lm(firesev ~ age, data = keeley),
data = keeley
)
# Get log-likelihood based AIC
AIC(mod)
# Get d-sep based AIC
AIC(mod, AIC.type = "dsep")
Run the code above in your browser using DataLab