data(agres)
mycode <- list(
GVA~quec(NPATENT,1,15),
PPI~quec(NPATENT,0,13)+quec(GVA,0,14),
ENTR_INCOME~quec(NPATENT,0,14)+quec(GVA,1,14)
)
myfit <- dlsem(mycode,group="COUNTRY",exogenous=c("GDP","FARM_SIZE"),
data=agres,uniroot.check=TRUE,imputation=FALSE,log=TRUE)
# the lag shape associated to specific paths
lagPlot(myfit,path="NPATENT*PPI")
lagPlot(myfit,path="NPATENT*GVA*PPI")
# the lag shape associated to an overall causal effect
lagPlot(myfit,from="NPATENT",to="PPI")
Run the code above in your browser using DataLab