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