[R-br] comando diagnosis

[Daniel C Bezerra]

Oi Alexandro,

Se vc usar uma tabela como objeto a partir de onde os dados devem ser
calculados esta tem que ser no formato:

VN FN
FP VP

Assim, vc deve usar:
tab<-matrix(c(406,17,112,64),nrow=2,ncol=2,byrow=T)
diagnosis(tab,plot=T)

Abs,

D


2011/10/13 Alexandro (Yahoo) <vl.alexandro em yahoo.com.br>

>  #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
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111013/c20cc3f3/attachment.html>