<span id="mailbox-conversation"><div>Hi Juan,</div>

<div><br></div>

<div>Thanks for reporting this issue. We will look into it and update you and the list as soon as it is resolved.</div>

<div><br></div>

<div>Best,</div>

<div>Graeme</div></span><div class="mailbox_signature">

<br><p style="margin-top: 0px; margin-bottom: 0px; line-height: normal; font-family: 'Helvetica Neue'; color: rgb(50, 51, 51);"><br></p>

<p style="margin-top: 0px; margin-bottom: 0px; line-height: normal; font-family: 'Helvetica Neue'; color: rgb(50, 51, 51);"></p>

</div>

<div class="gmail_quote">On Fri, Feb 27, 2015 at 1:52 PM, Juan Felipe Riaño Rodríguez <span dir="ltr"><<a href="mailto:jf.riano48@uniandes.edu.co" target="_blank">jf.riano48@uniandes.edu.co</a>> wrote: <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. 

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,  : 

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>_______________________________________________

<br>listpackage-discuss mailing list

<br>listpackage-discuss@lists.r-forge.r-project.org

<br>https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/listpackage-discuss

<br></p></blockquote></div><br>