data(agres)
mycode <- list(
GVA~quec(NPATENT,0,4),
PPI~quec(NPATENT,0,6)+quec(GVA,0,10),
ENTR_INCOME~quec(NPATENT,0,3)+quec(GVA,1,10)
)
mod0 <- dlsem(mycode,group="COUNTRY",context=c("GDP","FARM_SIZE"),
data=agres,uniroot.check=TRUE,imputation=FALSE,log=TRUE)
pathAnal(mod0,from="NPATENT",to="PPI",lag=c(0,5,10,15),cumul=TRUE,nitt=10000)
Run the code above in your browser using DataLab