data(industry)
mycode <- list(
Consum~quec(Job,0,6),
Pollution~quec(Job,1,11)+quec(Consum,1,6)
)
myfit <- dlsem(mycode,group="Region",exogenous=c("Population","GDP"),
data=industry,uniroot.check=TRUE,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