[R-br] comando diagnosis

[Alexandro (Yahoo)]

#Boa tarde a todos, numa regressão logística múltipla, após a escolha de um ponto de corte,
>
>#criei a tabela de classificação conforme CMR abaixo, e usei o 
>#comando dignosis do pacote DiagnosisMed
>#para calcular a sensibilidade e especificidade !
>
>
>require(DiagnosisMed)
>
>
>tab<-matrix(c(64,112,17,406),nrow=2,ncol=2,byrow=T,
>dimnames=list(Classificado=c(" 1 "," 0 "),
>Observado=c(" 1 "," 0 ")));tab
>
>
>diagnosis(tab)[]
>
>
>
>#Porém os valores estão trocados, o certo é 
>#Sensitividade =64 /(17+64)  =79,01%
>#Especificidade=406/(112+406)=78,378%
>#Tentei usar o comando:
>
>
>
>
>diagnosis(t(tab))[]
>
>
>
>
>#Mas não deu certo. Sendo assim:
>#1- o que ocorreu ?
>#2- Tem alguma função nos pacotes epiR ou ROCR que
>#calcula a sensibilidade e especificidade ?
>#3- para achar o ponto de corte "ideal" tenho que fazer várias tabelas de classificação
>#e plotar a sens. x espec. 
>#No "satanic", tem-se um comando do proc logistic que
>#mostra todas as sensibilidade e especificidade para diferentes pontos de cortes.
>#Tem algum comando similar no R para regressão logística múltipla ?
>#no caso univariado tem-se:
>y=c(0,1,0,0,1,1,1,0,0,0)
>x=c(12,19,14,13,14,15,16,17,17,18)
>require(DiagnosisMed);arearoc<-DiagnosisMed::ROC(y,x,Full=TRUE) 
>
>
>
>
>###########   MUITO OBRIGADO   ###########
>#Alex
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111013/fcd80362/attachment.html>