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)
isIndep(mod0,"NPATENT","ENTR_INCOME",given=c("GVA","PPI"))
Run the code above in your browser using DataLab