<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal>Hi,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Please, could somebody to help me with the errors below. I tried to replicate two analyses from list manual v6 but I received some errors.<o:p></o:p></p><p class=MsoNormal>Thank for your time and help.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>All the best!<o:p></o:p></p><p class=MsoNormal>Mircea<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>data(affirm)<o:p></o:p></p><p class=MsoNormal>ceiling.results <- ictreg(y ~ age + college + male + south, treat = "treat", J = 3, data = affirm, method = "ml", fit.start = "nls", ceiling = TRUE, ceiling.fit = "bayesglm", ceiling.formula = ~ age + college + male + south)<o:p></o:p></p><p class=MsoNormal><span style='background:yellow;mso-highlight:yellow'>Error in cbind(y, 1 - y) : object 'y' not found</span><o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>data(race)<o:p></o:p></p><p class=MsoNormal>mle.estimates <- ictreg(y ~ male + college + age + south, data = race)<o:p></o:p></p><p class=MsoNormal>Note: number of control items set to 3 <o:p></o:p></p><p class=MsoNormal>draws <- mvrnorm(n = 3, mu = coef(mle.estimates), Sigma = vcov(mle.estimates) * 9)<o:p></o:p></p><p class=MsoNormal>bayesDraws.1 <- ictregBayes(y ~ male + college + age + south, data = race, delta.start = draws[1, 1:5], psi.start = draws[1, 6:10], burnin = 10000, n.draws = 100000, delta.tune = diag(.002, 5), psi.tune = diag(.00025, 5), constrained.single = "full")<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>*** Starting posterior sampling... *** <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>10 percent done.<o:p></o:p></p><p class=MsoNormal> Metropolis acceptance ratios<o:p></o:p></p><p class=MsoNormal> Sensitive item model: 0.384<o:p></o:p></p><p class=MsoNormal> Control items model: 0.461<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>20 percent done.<o:p></o:p></p><p class=MsoNormal> Metropolis acceptance ratios<o:p></o:p></p><p class=MsoNormal> Sensitive item model: 0.38<o:p></o:p></p><p class=MsoNormal> Control items model: 0.464<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>30 percent done.<o:p></o:p></p><p class=MsoNormal> Metropolis acceptance ratios<o:p></o:p></p><p class=MsoNormal> Sensitive item model: 0.376<o:p></o:p></p><p class=MsoNormal> Control items model: 0.463<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>40 percent done.<o:p></o:p></p><p class=MsoNormal> Metropolis acceptance ratios<o:p></o:p></p><p class=MsoNormal> Sensitive item model: 0.377<o:p></o:p></p><p class=MsoNormal> Control items model: 0.464<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>50 percent done.<o:p></o:p></p><p class=MsoNormal> Metropolis acceptance ratios<o:p></o:p></p><p class=MsoNormal> Sensitive item model: 0.376<o:p></o:p></p><p class=MsoNormal> Control items model: 0.464<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>60 percent done.<o:p></o:p></p><p class=MsoNormal> Metropolis acceptance ratios<o:p></o:p></p><p class=MsoNormal> Sensitive item model: 0.375<o:p></o:p></p><p class=MsoNormal> Control items model: 0.465<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>70 percent done.<o:p></o:p></p><p class=MsoNormal> Metropolis acceptance ratios<o:p></o:p></p><p class=MsoNormal> Sensitive item model: 0.375<o:p></o:p></p><p class=MsoNormal> Control items model: 0.465<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>80 percent done.<o:p></o:p></p><p class=MsoNormal> Metropolis acceptance ratios<o:p></o:p></p><p class=MsoNormal> Sensitive item model: 0.374<o:p></o:p></p><p class=MsoNormal> Control items model: 0.465<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>90 percent done.<o:p></o:p></p><p class=MsoNormal> Metropolis acceptance ratios<o:p></o:p></p><p class=MsoNormal> Sensitive item model: 0.375<o:p></o:p></p><p class=MsoNormal> Control items model: 0.464<o:p></o:p></p><p class=MsoNormal><span style='background:yellow;mso-highlight:yellow'>Error in gzfile(file, "wb") : cannot open the connection<o:p></o:p></span></p><p class=MsoNormal><span style='background:yellow;mso-highlight:yellow'>In addition: Warning messages:<o:p></o:p></span></p><p class=MsoNormal><span style='background:yellow;mso-highlight:yellow'>1: In if (probit) { :<o:p></o:p></span></p><p class=MsoNormal><span style='background:yellow;mso-highlight:yellow'> the condition has length > 1 and only the first element will be used<o:p></o:p></span></p><p class=MsoNormal><span style='background:yellow;mso-highlight:yellow'>2: In if (robit) { :<o:p></o:p></span></p><p class=MsoNormal><span style='background:yellow;mso-highlight:yellow'> the condition has length > 1 and only the first element will be used<o:p></o:p></span></p><p class=MsoNormal><span style='background:yellow;mso-highlight:yellow'>3: In gzfile(file, "wb") :<o:p></o:p></span></p><p class=MsoNormal><span style='background:yellow;mso-highlight:yellow'> cannot open compressed file 'C:\Users\user\Documents/desktop/res.RData', probable reason 'No such file or directory'</span><o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal style='margin-bottom:12.0pt'><span lang=RO style='font-size:12.0pt;color:#002060'>_________________________________________________________</span><span style='font-size:12.0pt;color:#002060'><o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt;color:#002060'>Mircea Comsa, Assoc. Prof.<br>Faculty of Sociology and Social Work<br>Babes-Bolyai University, Romania<br>Bd. 21 Decembrie 1989, nr.128-130 RO-400604 Cluj-Napoca<br>Tel: + 40-0744-133.871<br>Fax: + 40-264-424.674<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>