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