isIndep(x, var1, var2, given = NULL, conf = 0.95)
dlsem
.NULL
, marginal independence is checked.dlsem
. Default is 0.95.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)
isIndep(myfit,"Job","Pollution",given=c("Consum"))
Run the code above in your browser using DataLab