Hi everyone, How do I run Heckman models in R with two binary dependent variables?
sampleSelection package in R works with standard heckman models ( binary DV for the selection equation and continuous DV for the outcome equation). In my case dependent variables are both binary (actually ordered but I didn’t find anything on that) So using sampleSelection package one could do this by running: SelectionEquation <- binaryDV1 ~ x1+x2+x3+x4 OutcomeEquation <- binaryDV2~o7+x1+x4+x5 HeckmanModel <- heckit(SelectionEquation,OutcomeEquation, data=mydata, method="2step") The problem is that heckit() doesn’t work here. I think in STATA one could use heckprob command for this. Anyone who knows more than me? Thanks! [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.