<font color="#000066"><font face="courier new,monospace">Amigos de R,</font></font><div><font color="#000066" face="'courier new', monospace"><br></font></div><div><font color="#000066" face="'courier new', monospace">Pacientes foram expostos a experimentos para verificar o tempo de resposta a determinadas exposições. Para cada exposição (que são 3, mas prefiro aborda-las individualmente) verificou-se se o sujeito acertou a resposta desejada e quanto tempo levou para responder. A hipótese do estudo é que quando o sujeito acerta a resposta o tempo de resposta é menor do que quando erra. O problema é que cada paciente é exposto diversas vezes. Então eu tenho um banco com 10 pacientes, mas com mais de 300 linhas, semelhante a simulação abaixo. A minha questão é que eu imagino que haja dependencia dos pacientes, ou seja, que as respostas sejam mais parecidas dentro de cada sujeito do que entre cada sujeito. Provavelmente a variação teria dois componentes, o de dentro do sujeito e entre os sujeitos. POr isso, eu imaginei que utilizar um modelo aleatório com o paciente como o elemento aleatório fosse mais interpretável que um modelo linear simples. </font></div>
<div><font color="#000066" face="'courier new', monospace"><br></font></div><div><font color="#000066" face="'courier new', monospace">Mas infelizmente nao estou conseguindo fazer a função funcionar. Olhei os exemplos , li algumas vezes a documentção... mas ainda não foi.</font></div>
<div><font color="#000066" face="'courier new', monospace"><br></font></div><div><font color="#000066" face="'courier new', monospace">Qualquer sugestão de como fazer funcionar é bem vinda!<br></font><div>
<font color="#000066"><font face="courier new,monospace"><br></font></font></div><div><font color="#000066"><font face="courier new,monospace"><div>z <- rnorm(200, mean = 120 , sd = 10)</div><div>x <- rnorm(100, mean = 90, sd = 20)</div>
<div>data <- as.data.frame(rbind(cbind(x,1,1:10),cbind(z,0,1:10)))</div><div>names(data) <- c('tempo','desfecho','pac')</div><div><br></div><div>summary(lm(tempo ~ desfecho + pac,data))</div>
<div><br></div></font></font></div><div><font color="#000066"><font face="courier new,monospace">library('nlm')</font></font></div><div><font color="#000066"><font face="courier new,monospace"><div>> lme(tempo ~ desfecho + pac,data)</div>
<div>Erro em getGroups.data.frame(dataMix, groups) : </div><div> Invalid formula for groups</div><div>> lme(tempo ~ desfecho,data,random= pac)</div><div>Erro em reStruct(random, REML = REML, data = NULL) : </div><div>
objeto 'pac' não encontrado</div><div>> lme(tempo ~ desfecho,data,random= ~ pac)</div><div>Erro em getGroups.data.frame(dataMix, groups) : </div><div> Invalid formula for groups</div><div>> lme(tempo ~ desfecho,data,random= ~desfecho|pac)</div>
<div>Erro em lme.formula(tempo ~ desfecho, data, random = ~desfecho | pac) : </div><div> nlminb problem, convergence error code = 1</div><div> message = iteration limit reached without convergence (10)</div></font></font></div>
<div><font color="#000066"><font face="courier new,monospace"><br></font></font></div><div><font color="#000066"><font face="courier new,monospace">Abraço forte,</font></font></div><div><font color="#000066"><font face="courier new,monospace"><br clear="all">
</font></font><font face="'courier new', monospace" style="background-color:rgb(255,255,255)" color="#000066">Dr. Pedro Emmanuel A. A. do Brasil</font><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>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>
<br>
</div></div>