Oi Alexandro,<br><br>Se vc usar uma tabela como objeto a partir de onde os dados devem ser calculados esta tem que ser no formato:<br><br>VN FN<br>FP VP<br><br>Assim, vc deve usar:<br>tab<-matrix(c(406,17,112,64),nrow=2,ncol=2,byrow=T)<br>
diagnosis(tab,plot=T)<br><br>Abs,<br><br>D<br><br><br><div class="gmail_quote">2011/10/13 Alexandro (Yahoo) <span dir="ltr"><<a href="mailto:vl.alexandro@yahoo.com.br">vl.alexandro@yahoo.com.br</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div style="color:#000;background-color:#fff;font-family:times new roman, new york, times, serif;font-size:12pt"><div><blockquote style="border-left-width:2px;border-left-style:solid;border-left-color:rgb(16, 16, 255);margin-left:5px;padding-left:5px">
<div style="font-family:'times new roman', 'new york', times, serif"><div style="font-family:'times new roman', 'new york', times, serif"><font><span style="font-size:16px">#Boa tarde a todos, numa regressão logística múltipla, após a escolha de um ponto de corte,</span><br>
</font><div style="font-family:'times new roman', 'new york', times, serif;font-size:12pt"><div><div style="color:rgb(0, 0, 0);background-color:rgb(255, 255, 255);font-size:12pt;font-family:'times new roman', 'new york', times, serif">
<div><div><font face="'times new roman', 'new york', times,
serif">#criei a tabela de classificação conforme CMR abaixo, e usei o </font></div><div><font face="'times new roman', 'new york', times, serif">#comando dignosis do pacote DiagnosisMed</font></div><div>
<font face="'times new roman', 'new york', times, serif">#para calcular a sensibilidade e especificidade !</font></div><div><font face="'times new roman', 'new york', times, serif"><br></font></div>
<div><font face="'times new roman', 'new
york', times, serif">require(DiagnosisMed)</font></div><div><font face="'times new roman', 'new york', times, serif"><br></font></div><div><font face="'times new roman', 'new york', times, serif">tab<-matrix(c(64,112,17,406),nrow=2,ncol=2,byrow=T,</font></div>
<div><font face="'times new roman', 'new york', times, serif">dimnames=list(Classificado=c(" 1 "," 0 "),</font></div><div><font face="'times new roman', 'new york', times, serif">Observado=c(" 1 "," 0 ")));tab</font></div>
<div><span style="font-family:times, serif"><br></span></div><div><span style="font-family:times, serif">diagnosis(tab)[]</span><br></div><div><font face="'times new roman', 'new york', times,
serif"><br></font></div><div><font face="'times new roman', 'new york', times, serif">#Porém os valores estão trocados, o certo é </font></div><div><font face="'times new roman', 'new york', times, serif">#Sensitividade =64 /(17+64) =79,01%</font></div>
<div><font face="'times new roman', 'new york', times, serif">#Especificidade=406/(112+406)=78,378%</font></div><div><font face="'times new roman', 'new york', times, serif">#Tentei usar o comando:</font></div>
<div><font face="'times new roman', 'new york', times, serif"><br></font></div><div><font face="'times new roman', 'new york', times, serif"><br></font></div><div><font face="'times new roman', 'new york', times,
serif">diagnosis(t(tab))[]</font></div><div><font face="'times new roman', 'new york', times,
serif"><br></font></div><div><font face="'times new roman', 'new york', times, serif"><br></font></div><div><font face="'times new roman', 'new york', times, serif">#Mas não deu certo. Sendo assim:</font></div>
<div><font face="'times new roman', 'new york', times, serif">#1- o que ocorreu ?</font></div><div><font face="'times new roman', 'new york', times, serif">#2- Tem alguma função nos pacotes epiR ou ROCR que</font></div>
<div><font face="'times new roman', 'new york', times, serif">#calcula a sensibilidade e especificidade ?</font></div><div><font style="background-color:rgb(255, 255, 255)" color="#7F7F00" face="'times new roman', 'new york', times, serif"><b>#3- para achar o ponto de corte "ideal" tenho que fazer várias
tabelas de classificação</b></font></div><div><font style="background-color:rgb(255, 255, 255)" color="#7F7F00" face="'times new roman', 'new york', times, serif"><b>#e plotar a sens. x espec. </b></font></div>
<div><font style="background-color:rgb(255, 255, 255)" color="#7F7F00" face="'times new roman', 'new york', times, serif"><b>#No "satanic", tem-se um comando do proc logistic que</b></font></div>
<div><font style="background-color:rgb(255, 255, 255)" color="#7F7F00" face="'times new roman', 'new york', times, serif"><b>#mostra todas as sensibilidade e especificidade para diferentes pontos de cortes.</b></font></div>
<div><font style="background-color:rgb(255, 255, 255)" color="#7F7F00" face="'times new roman', 'new york', times, serif"><b>#Tem algum comando similar no R para regressão
logística múltipla ?</b></font></div><div><font face="'times new roman', 'new
york', times, serif">#no caso univariado tem-se:</font></div><div><font face="'times new roman', 'new york', times, serif">y=c(0,1,0,0,1,1,1,0,0,0)</font></div><div><font face="'times new roman', 'new york', times, serif">x=c(12,19,14,13,14,15,16,17,17,18)</font></div>
<div><font face="'times new roman', 'new york', times, serif">require(DiagnosisMed);arearoc<-DiagnosisMed::ROC(y,x,Full=TRUE) </font></div><div style="font-size:12pt;font-family:times, serif"><br></div>
<div style="font-size:12pt;font-family:times, serif"><br></div><div style="font-size:12pt;font-family:times, serif">########### <font color="#C00000">MUITO OBRIGADO</font> ###########</div><div style="font-size:12pt;font-family:times, serif">
#Alex</div></div></div></div></div></div></div></blockquote></div></div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br>