<div dir="ltr">Dear Juan,<div><br></div><div>Thanks for your patience. This error was caused by a change in another software package list depends on, and we have now resolved the problem. I will send you a test version of the package off-list, and we will release an update to CRAN including the fix soon.</div><div><br></div><div>Best,</div><div>Graeme<br><br><div class="gmail_quote">On Fri, Feb 27, 2015 at 1:52 PM Juan Felipe Riaño Rodríguez <<a href="mailto:jf.riano48@uniandes.edu.co">jf.riano48@uniandes.edu.co</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi!<br>
<br>
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>
<br>
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>
<br>
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>
<br>
Thank you very much,<br>
<br>
<br>
<br>
> data(race)<br>
><br>
> 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>
Error in logistic(x.treatment.curr %*% fit.control.coef) :<br>
  object 'fit.control.coef' not found<br>
><br>
> summary(ceiling.results, boundary.proportions = T)<br>
Error in summary(ceiling.results, boundary.proportions = T) :<br>
  object 'ceiling.results' not found<br>
><br>
> # Fit standard design ML model with floor effects alone<br>
> # Replicates Columns 5-6<br>
><br>
> data(race)<br>
><br>
> floor.results <- ictreg(y ~ age + college + male + south, treat = "treat",<br>
+                       J = 3, data = affirm, method = "ml", fit.start = "glm",<br>
+                       floor = TRUE, floor.fit = "bayesglm",<br>
+                       floor.formula = ~ age + college + male + south)<br>
Error in ictreg(y ~ age + college + male + south, treat = "treat", J = 3,  :<br>
  log-likelihood is not monotonically increasing.<br>
><br>
> summary(floor.results, boundary.proportions = T)<br>
Error in summary(floor.results, boundary.proportions = T) :<br>
  object 'floor.results' not found<br>
><br>
______________________________<u></u>_________________<br>
listpackage-discuss mailing list<br>
<a href="mailto:listpackage-discuss@lists.r-forge.r-project.org" target="_blank">listpackage-discuss@lists.r-<u></u>forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/listpackage-discuss" target="_blank">https://lists.r-forge.r-<u></u>project.org/cgi-bin/mailman/<u></u>listinfo/listpackage-discuss</a><br>
</blockquote></div></div></div>