Arne Henningsen
2016-Aug-27 08:39 UTC
[R] Hickman models with two binary dependent variables in R
See also: http://r-forge.r-project.org/forum/forum.php?thread_id=31866&forum_id=844&group_id=256 On 26 August 2016 at 16:11, PIKAL Petr <petr.pikal at precheza.cz> wrote:> Hi > > See in line > >> -----Original Message----- >> From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Faradj >> Koliev >> Sent: Thursday, August 25, 2016 12:32 PM >> To: r-help at r-project.org >> Subject: [R] Hickman models with two binary dependent variables in R >> >> 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) > > From help page > > The endogenous variable of the argument 'selection' must have exactly two levels (e.g. 'FALSE' and 'TRUE', or '0' and '1'). By default the levels are sorted in increasing order ('FALSE' is before 'TRUE', and '0' is before '1'). This also applies for the binary outcome equation. For continuous-oucome cases, the dependent variable(s) should be numeric. > > seems to me that both equatio0ns can have binary values. > > >> >> 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 > > What does it mean. Be more specific. > > And provide some data (preferably by dput). > Or at least result of str(yourdata) to show that they are appropriate to the function. > And do not post in HTML. > > Cheers > Petr > >> use heckprob command for this. Anyone who knows more than me? >> >> Thanks! >> >> >> >> >> >> [[alternative HTML version deleted]] >> >> ______________________________________________ >> R-help at 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. > > ________________________________ > Tento e-mail a jak?koliv k n?mu p?ipojen? dokumenty jsou d?v?rn? a jsou ur?eny pouze jeho adres?t?m. > Jestli?e jste obdr?el(a) tento e-mail omylem, informujte laskav? neprodlen? jeho odes?latele. Obsah tohoto emailu i s p??lohami a jeho kopie vyma?te ze sv?ho syst?mu. > Nejste-li zam??len?m adres?tem tohoto emailu, nejste opr?vn?ni tento email jakkoliv u??vat, roz?i?ovat, kop?rovat ?i zve?ej?ovat. > Odes?latel e-mailu neodpov?d? za eventu?ln? ?kodu zp?sobenou modifikacemi ?i zpo?d?n?m p?enosu e-mailu. > > V p??pad?, ?e je tento e-mail sou??st? obchodn?ho jedn?n?: > - vyhrazuje si odes?latel pr?vo ukon?it kdykoliv jedn?n? o uzav?en? smlouvy, a to z jak?hokoliv d?vodu i bez uveden? d?vodu. > - a obsahuje-li nab?dku, je adres?t opr?vn?n nab?dku bezodkladn? p?ijmout; Odes?latel tohoto e-mailu (nab?dky) vylu?uje p?ijet? nab?dky ze strany p??jemce s dodatkem ?i odchylkou. > - trv? odes?latel na tom, ?e p??slu?n? smlouva je uzav?ena teprve v?slovn?m dosa?en?m shody na v?ech jej?ch n?le?itostech. > - odes?latel tohoto emailu informuje, ?e nen? opr?vn?n uzav?rat za spole?nost ??dn? smlouvy s v?jimkou p??pad?, kdy k tomu byl p?semn? zmocn?n nebo p?semn? pov??en a takov? pov??en? nebo pln? moc byly adres?tovi tohoto emailu p??padn? osob?, kterou adres?t zastupuje, p?edlo?eny nebo jejich existence je adres?tovi ?i osob? j?m zastoupen? zn?m?. > > This e-mail and any documents attached to it may be confidential and are intended only for its intended recipients. > If you received this e-mail by mistake, please immediately inform its sender. Delete the contents of this e-mail with all attachments and its copies from your system. > If you are not the intended recipient of this e-mail, you are not authorized to use, disseminate, copy or disclose this e-mail in any manner. > The sender of this e-mail shall not be liable for any possible damage caused by modifications of the e-mail or by delay with transfer of the email. > > In case that this e-mail forms part of business dealings: > - the sender reserves the right to end negotiations about entering into a contract in any time, for any reason, and without stating any reasoning. > - if the e-mail contains an offer, the recipient is entitled to immediately accept such offer; The sender of this e-mail (offer) excludes any acceptance of the offer on the part of the recipient containing any amendment or variation. > - the sender insists on that the respective contract is concluded only upon an express mutual agreement on all its aspects. > - the sender of this e-mail informs that he/she is not authorized to enter into any contracts on behalf of the company except for cases in which he/she is expressly authorized to do so in writing, and such authorization or power of attorney is submitted to the recipient or the person represented by the recipient, or the existence of such authorization is known to the recipient of the person represented by the recipient. > ______________________________________________ > R-help at 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.-- Arne Henningsen http://www.arne-henningsen.name
Faradj Koliev
2016-Aug-27 14:11 UTC
[R] Hickman models with two binary dependent variables in R
Dear Arne, Many thanks for this, It actually worked with heckit() command as well, do I need to use selection()? Also, I would be really grateful if you can suggest a package that would allow for estimation of heckman models with two ordered variables (0-1-2). Can sampleSelection handle this? Warm regards, Faradj> 27 aug. 2016 kl. 10:39 skrev Arne Henningsen <arne.henningsen at gmail.com>: > > See also: > > http://r-forge.r-project.org/forum/forum.php?thread_id=31866&forum_id=844&group_id=256 > > > > On 26 August 2016 at 16:11, PIKAL Petr <petr.pikal at precheza.cz> wrote: >> Hi >> >> See in line >> >>> -----Original Message----- >>> From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Faradj >>> Koliev >>> Sent: Thursday, August 25, 2016 12:32 PM >>> To: r-help at r-project.org >>> Subject: [R] Hickman models with two binary dependent variables in R >>> >>> 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) >> >> From help page >> >> The endogenous variable of the argument 'selection' must have exactly two levels (e.g. 'FALSE' and 'TRUE', or '0' and '1'). By default the levels are sorted in increasing order ('FALSE' is before 'TRUE', and '0' is before '1'). This also applies for the binary outcome equation. For continuous-oucome cases, the dependent variable(s) should be numeric. >> >> seems to me that both equatio0ns can have binary values. >> >> >>> >>> 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 >> >> What does it mean. Be more specific. >> >> And provide some data (preferably by dput). >> Or at least result of str(yourdata) to show that they are appropriate to the function. >> And do not post in HTML. >> >> Cheers >> Petr >> >>> use heckprob command for this. Anyone who knows more than me? >>> >>> Thanks! >>> >>> >>> >>> >>> >>> [[alternative HTML version deleted]] >>> >>> ______________________________________________ >>> R-help at 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. >> >> ________________________________ >> Tento e-mail a jak?koliv k n?mu p?ipojen? dokumenty jsou d?v?rn? a jsou ur?eny pouze jeho adres?t?m. >> Jestli?e jste obdr?el(a) tento e-mail omylem, informujte laskav? neprodlen? jeho odes?latele. Obsah tohoto emailu i s p??lohami a jeho kopie vyma?te ze sv?ho syst?mu. >> Nejste-li zam??len?m adres?tem tohoto emailu, nejste opr?vn?ni tento email jakkoliv u??vat, roz?i?ovat, kop?rovat ?i zve?ej?ovat. >> Odes?latel e-mailu neodpov?d? za eventu?ln? ?kodu zp?sobenou modifikacemi ?i zpo?d?n?m p?enosu e-mailu. >> >> V p??pad?, ?e je tento e-mail sou??st? obchodn?ho jedn?n?: >> - vyhrazuje si odes?latel pr?vo ukon?it kdykoliv jedn?n? o uzav?en? smlouvy, a to z jak?hokoliv d?vodu i bez uveden? d?vodu. >> - a obsahuje-li nab?dku, je adres?t opr?vn?n nab?dku bezodkladn? p?ijmout; Odes?latel tohoto e-mailu (nab?dky) vylu?uje p?ijet? nab?dky ze strany p??jemce s dodatkem ?i odchylkou. >> - trv? odes?latel na tom, ?e p??slu?n? smlouva je uzav?ena teprve v?slovn?m dosa?en?m shody na v?ech jej?ch n?le?itostech. >> - odes?latel tohoto emailu informuje, ?e nen? opr?vn?n uzav?rat za spole?nost ??dn? smlouvy s v?jimkou p??pad?, kdy k tomu byl p?semn? zmocn?n nebo p?semn? pov??en a takov? pov??en? nebo pln? moc byly adres?tovi tohoto emailu p??padn? osob?, kterou adres?t zastupuje, p?edlo?eny nebo jejich existence je adres?tovi ?i osob? j?m zastoupen? zn?m?. >> >> This e-mail and any documents attached to it may be confidential and are intended only for its intended recipients. >> If you received this e-mail by mistake, please immediately inform its sender. Delete the contents of this e-mail with all attachments and its copies from your system. >> If you are not the intended recipient of this e-mail, you are not authorized to use, disseminate, copy or disclose this e-mail in any manner. >> The sender of this e-mail shall not be liable for any possible damage caused by modifications of the e-mail or by delay with transfer of the email. >> >> In case that this e-mail forms part of business dealings: >> - the sender reserves the right to end negotiations about entering into a contract in any time, for any reason, and without stating any reasoning. >> - if the e-mail contains an offer, the recipient is entitled to immediately accept such offer; The sender of this e-mail (offer) excludes any acceptance of the offer on the part of the recipient containing any amendment or variation. >> - the sender insists on that the respective contract is concluded only upon an express mutual agreement on all its aspects. >> - the sender of this e-mail informs that he/she is not authorized to enter into any contracts on behalf of the company except for cases in which he/she is expressly authorized to do so in writing, and such authorization or power of attorney is submitted to the recipient or the person represented by the recipient, or the existence of such authorization is known to the recipient of the person represented by the recipient. >> ______________________________________________ >> R-help at 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. > > > > -- > Arne Henningsen > http://www.arne-henningsen.name[[alternative HTML version deleted]]
Arne Henningsen
2016-Aug-29 12:43 UTC
[R] Hickman models with two binary dependent variables in R
Dear Faradj On 27 August 2016 at 16:11, Faradj Koliev <faradj.g at gmail.com> wrote:> It actually worked with heckit() command as well, do I need to use > selection()?I suggest that you use the selection() command/function only. heckit() is just a wrapper to selection() with the only difference that heckit() estimates the model with the 2-step method by default, while selection() estimates the model by the maximum likelihood method by default (unless specified otherwise with argument "method"). I see no good reason for estimating the sample-selection model by the 2-step method instead of the maximum likelihood method.> Also, I would be really grateful if you can suggest a package that would > allow for estimation of heckman models with two ordered variables (0-1-2). > Can sampleSelection handle this?Do you mean that the dependent variable of the selection equation is an ordered variable with three levels (0-1-2), that the dependent variable of the outcome equation is an ordered variable with three levels (0-1-2), or that both of the dependent variables are ordered variable with three levels (0-1-2) each? In any case: no, this is not implemented in the sampleSelection package. Sorry! Anyway, I guess that it is not too complicated to derive the likelihood function and implement the estimation yourself, e.g., using the maxLik package. If you do this, I would be happy to help you to implement this feature in the sampleSelection package. Best wishes, Arne -- Arne Henningsen http://www.arne-henningsen.name