# NOT RUN {
data(industry)
indus.code <- list(
Consum~ecq(Job,0,5),
Pollution~ecq(Job,1,8)+ecq(Consum,1,7)
)
indus.mod <- dlsem(indus.code,group="Region",exogenous=c("Population","GDP"),data=industry,
log=TRUE)
# the lag shape of the causal effect associated to specific paths
lagPlot(indus.mod,path="Job*Pollution")
lagPlot(indus.mod,path="Job*Consum*Pollution")
# the lag shape of an overall causal effect
lagPlot(indus.mod,from="Job",to="Pollution")
# }
Run the code above in your browser using DataLab