<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class="">Hi!</div><div class=""><br class=""></div><div class="">I was wondering if this could be a simple matter of the R versions. Originally, which R version was used to run the paper’s tables?</div><div class=""><br class=""></div><div class="">Thank you very much,</div><div class=""><br class=""></div><div apple-content-edited="true" class="">
Juan Felipe Riaño Rodriguez<br class="">Centro de Estudios sobre Desarrollo Económico - CEDE,<br class="">Facultad de Economía - Universidad de los Andes,<br class="">Calle 19A No. 1-37 Este - Bloque W (W-815),<br class="">Tel. +57-1-339-4949 ext 3409
</div>
<br class=""><div><blockquote type="cite" class=""><div class="">El 27/02/2015, a las 1:54 p.m., Graeme Blair <<a href="mailto:graeme.blair@columbia.edu" class="">graeme.blair@columbia.edu</a>> escribió:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><span id="mailbox-conversation" class=""><div class="">Hi Juan,</div>
<div class=""><br class=""></div>
<div class="">Thanks for reporting this issue. We will look into it and update you and the list as soon as it is resolved.</div>
<div class=""><br class=""></div>
<div class="">Best,</div>
<div class="">Graeme</div></span><div class="mailbox_signature">
<br class=""><div style="margin-top: 0px; margin-bottom: 0px; line-height: normal; font-family: 'Helvetica Neue'; color: rgb(50, 51, 51);" class=""><br class=""></div><p style="margin-top: 0px; margin-bottom: 0px; line-height: normal; font-family: 'Helvetica Neue'; color: rgb(50, 51, 51);" class=""></p>
</div>
<br class=""><br class=""><div class="gmail_quote"><p class="">On Fri, Feb 27, 2015 at 1:52 PM, Juan Felipe Riaño Rodríguez <span dir="ltr" class=""><<a href="mailto:jf.riano48@uniandes.edu.co" target="_blank" class="">jf.riano48@uniandes.edu.co</a>></span> wrote:<br class=""></p><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><p class="">Hi!
<br class=""><br class="">Working with my data to verify ceiling and floor liars, I used your replication code as a guide and got error messages about the Log likelihood not being monotonically increasing. But then I realized that I also get this error with your own code on your data, that is simply running the replication of your paper.
<br class=""><br class="">In particular when trying to replicate Table 7, columns 3-6, as below, I get the error messages shown below (if I change the fit.start parameter, the error that pops out may change from object not found to non non monotonically increasing Log likelihood, but I always get an error anyways).
<br class=""><br class="">In short, before concluding there is problem with my data specifically, given this finding I would like to point it out and make sure there is not a more general problem with the code (perhaps something got messed up in an software update or that of package dependencies?).
<br class=""><br class="">Thank you very much,
<br class=""><br class=""><br class=""><br class="">> data(race)
<br class="">>
<br class="">> ceiling.results <- ictreg(y ~ age + college + male + south, treat = "treat", J = 3, data = affirm, method = "ml", fit.start = "ml", ceiling = TRUE, ceiling.fit = "bayesglm", ceiling.formula = ~ age + college + male + south)
<br class="">Error in logistic(x.treatment.curr %*% fit.control.coef) :
<br class=""> object 'fit.control.coef' not found
<br class="">>
<br class="">> summary(ceiling.results, boundary.proportions = T)
<br class="">Error in summary(ceiling.results, boundary.proportions = T) :
<br class=""> object 'ceiling.results' not found
<br class="">>
<br class="">> # Fit standard design ML model with floor effects alone
<br class="">> # Replicates Columns 5-6
<br class="">>
<br class="">> data(race)
<br class="">>
<br class="">> floor.results <- ictreg(y ~ age + college + male + south, treat = "treat",
<br class="">+ J = 3, data = affirm, method = "ml", fit.start = "glm",
<br class="">+ floor = TRUE, floor.fit = "bayesglm",
<br class="">+ floor.formula = ~ age + college + male + south)
<br class="">Error in ictreg(y ~ age + college + male + south, treat = "treat", J = 3, :
<br class=""> log-likelihood is not monotonically increasing.
<br class="">>
<br class="">> summary(floor.results, boundary.proportions = T)
<br class="">Error in summary(floor.results, boundary.proportions = T) :
<br class=""> object 'floor.results' not found
<br class="">>
<br class="">_______________________________________________
<br class="">listpackage-discuss mailing list
<br class=""><a href="mailto:listpackage-discuss@lists.r-forge.r-project.org" class="">listpackage-discuss@lists.r-forge.r-project.org</a>
<br class=""><a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/listpackage-discuss" class="">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/listpackage-discuss</a>
<br class=""></p></blockquote></div><br class=""></div></blockquote></div><br class=""></body></html>