<div dir="ltr"><div class="gmail_default" style="font-family:'courier new',monospace;color:rgb(0,0,102)"><span style="color:rgb(0,0,0);font-family:arial,sans-serif;font-size:12.727272033691406px">Amigos de R, </span><div style="color:rgb(0,0,0);font-family:arial,sans-serif;font-size:12.727272033691406px">
<br></div><div style="color:rgb(0,0,0);font-family:arial,sans-serif;font-size:12.727272033691406px">Estamos tentando implementar em código R uma rotina para suavizar uma curva ROC que utiliza modelos mistos (Haobo Ren , Xiao-Hua Zhou & Hua Liang (2004): A Flexible Method for Estimating the ROC Curve, Journal of Applied Statistics, 31:7, 773-784).</div>
<div style="color:rgb(0,0,0);font-family:arial,sans-serif;font-size:12.727272033691406px"><br></div><div style="color:rgb(0,0,0);font-family:arial,sans-serif;font-size:12.727272033691406px">O problema é que na hora de rodar a função lme erros diversos aparecem. Aparentemente não estamos sabendo colocar os argumentos de forma apropriada na função, pois o texto original não especifica uma variável de agrupamento.</div>
<div style="color:rgb(0,0,0);font-family:arial,sans-serif;font-size:12.727272033691406px"><br></div><div style="color:rgb(0,0,0);font-family:arial,sans-serif;font-size:12.727272033691406px">Os objetos são:</div><div style="color:rgb(0,0,0);font-family:arial,sans-serif;font-size:12.727272033691406px">
<br></div><div style="color:rgb(0,0,0);font-family:arial,sans-serif;font-size:12.727272033691406px"><div>Rempi <- c(0.75, 0.80, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.90, 0.95, 0.95, 0.95, 1.00, 1.00, 1.00, 1.00)</div>
<div><br></div><div>X <- matrix(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, </div><div>1, 1, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, </div><div>0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 0.0025, 0.01, </div>
<div>0.0225, 0.04, 0.0625, 0.09, 0.1225, 0.16, 0.2025, 0.25, 0.3025, </div><div>0.36, 0.4225, 0.49, 0.5625, 0.64, 0.7225, 0.81, 0.9025, 0.000125, </div><div>0.001, 0.003375, 0.008, 0.015625, 0.027, 0.042875, 0.064, 0.091125, </div>
<div>0.125, 0.166375, 0.216, 0.274625, 0.343, 0.421875, 0.512, 0.614125, </div><div>0.729, 0.857375),nrow=19,byrow=F)</div><div><br></div><div>Z <- matrix(c(0, 0, 2e-04, 0.0012, 0.0037, 0.0086, 0.0166, 0.0284, </div><div>
0.0447, 0.0664, 0.0942, 0.1288, 0.171, 0.2214, 0.281, 0.3504, </div><div>0.4304, 0.5217, 0.625, 0, 0, 0, 0, 3e-04, 0.0015, 0.0045, 0.0099, </div><div>0.0186, 0.0313, 0.0486, 0.0715, 0.1005, 0.1366, 0.1804, 0.2326, </div>
<div>
0.2941, 0.3655, 0.4477, 0, 0, 0, 0, 0, 0, 4e-04, 0.002, 0.0054, </div><div>0.0114, 0.0208, 0.0343, 0.0527, 0.0768, 0.1072, 0.1447, 0.1901, </div><div>0.2441, 0.3075, 0, 0, 0, 0, 0, 0, 0, 0, 6e-04, 0.0025, 0.0063, </div><div>
0.013, 0.0231, 0.0376, 0.0571, 0.0823, 0.1141, 0.1531, 0.2002, </div><div>0, 0, 0, 0, 0, 0, 0, 0, 0, 1e-04, 9e-04, 0.003, 0.0074, 0.0147, </div><div>0.0257, 0.0411, 0.0616, 0.0881, 0.1213, 0, 0, 0, 0, 0, 0, 0, </div><div>
0, 0, 0, 0, 2e-04, 0.0012, 0.0037, 0.0086, 0.0166, 0.0284, 0.0447, </div><div>0.0664, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3e-04, 0.0015, </div><div>0.0045, 0.0099, 0.0186, 0.0313, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, </div><div>
0, 0, 0, 0, 0, 4e-04, 0.002, 0.0054, 0.0114, 0, 0, 0, 0, 0, 0, </div><div>0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6e-04, 0.0025, 0, 0, 0, 0, 0, </div><div>0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1e-04), nrow=19,byrow=F)</div><div>
<br></div><div>grouping <- 1:19 # isto é o que não faz sentido... só escrevi para tentar rodar</div><div><br></div><div>Em que a expressão do modelo é:</div><div><br></div><div>Rempi = X*Beta + Z*b + e</div><div><br></div>
<div>Tentamos:</div><div><br></div><div><div>> lme(fixed = Rempi ~ X, random = ~ Z|grouping)</div><div><br></div><div>Erro em MEEM(object, conLin, control$niterEM) : <br></div><div> Singularity in backsolve at level 0, block 1</div>
<div>Além disso: Mensagens de aviso perdidas:</div><div>In lme.formula(fixed = Rempi ~ X, random = ~Z | grouping) :</div><div> Fewer observations than random effects in all level 1 groups</div></div><div><br></div><div style>
Alguma luz pra poder fazer esse modelo dar certo?</div><div style><br></div><div style>Abraço forte,</div></div></div><div><div dir="ltr"><span style="color:rgb(0,0,102);font-family:'courier new',monospace"><br></span></div>
<div dir="ltr"><span style="color:rgb(0,0,102);font-family:'courier new',monospace">Dr. Pedro Emmanuel A. A. do Brasil</span><br></div><div dir="ltr"><div><font face="'courier new', monospace" style="background-color:rgb(255,255,255)" color="#000066">Curriculum Lattes:
<span style="text-align:left"><a href="http://lattes.cnpq.br/6597654894290806" target="_blank">http://lattes.cnpq.br/6597654894290806</a></span> <br>ResearchGate.net: <a href="https://www.researchgate.net/profile/Pedro_Brasil2/" target="_blank">https://www.researchgate.net/profile/Pedro_Brasil2/</a><br>
</font></div><div><font face="'courier new', monospace" style="background-color:rgb(255,255,255)" color="#000066">Instituto de Pesquisa Clínica Evandro Chagas<br>Fundação Oswaldo Cruz<br>Rio de Janeiro - Brasil<br>
Av. Brasil 4365, <br>CEP 21040-360,<br>Tel 55 21 3865-9648<br>email: <a href="mailto:pedro.brasil@ipec.fiocruz.br" target="_blank">pedro.brasil@ipec.fiocruz.br</a><br>email: <a href="mailto:emmanuel.brasil@gmail.com" target="_blank">emmanuel.brasil@gmail.com</a><br>
<br>---Apoio aos softwares livres<br><a href="http://www.zotero.org" target="_blank">www.zotero.org</a> - gerenciamento de referências bibliográficas. <br><a href="http://www.broffice.org" target="_blank">www.broffice.org</a> ou <a href="http://www.libreoffice.org/" target="_blank">www.libreoffice.org</a> - textos, planilhas ou apresentações.<br>
<a href="http://www.epidata.dk" target="_blank">www.epidata.dk</a> - entrada de dados.<br><a href="http://www.r-project.org" target="_blank">www.r-project.org</a> - análise de dados.<br><a href="http://www.ubuntu.com" target="_blank">www.ubuntu.com</a> - sistema operacional</font></div>
</div></div>
</div>