<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
</head>
<body dir="auto">
<div><span></span></div>
<div>
<div>Dear Graeme,</div>
<div><br>
</div>
<div>Thank you very much! I will be wating for your email and for the update<br>
<br>
Best!</div>
<div><br>
</div>
<div>Juan Felipe<br>
<div><br>
</div>
</div>
<div><br>
On 2/04/2015, at 3:35 p.m., Graeme Blair <<a href="mailto:graeme.blair@columbia.edu">graeme.blair@columbia.edu</a>> wrote:<br>
<br>
</div>
<blockquote type="cite">
<div>
<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>
</div>
</blockquote>
</div>
</body>
</html>