data(keeley)
mod <- psem(
lm(rich ~ cover, data=keeley),
lm(cover ~ firesev, data=keeley),
lm(firesev ~ age, data=keeley),
data = keeley
)
plot(mod)
### More customized plot
plot(mod, node_attrs = list(
shape = "rectangle", color = "black",
fillcolor = "orange", x = 3, y=1:4))
Run the code above in your browser using DataLab