Les posteurs les plus actifs de la semaine
Aucun utilisateur |
Sujets les plus vus
Courbe de Roc.
2 participants
Page 1 sur 1
Courbe de Roc.
Bonjour,
Je suis en train de lire document pdf sur Courbe de Roc sur ce lien :
le problème que j'arrive pas à comprendre l'objet
Est ce que quelqu'un peux me aider ?
Merci beaucoup d'avance.
Je suis en train de lire document pdf sur Courbe de Roc sur ce lien :
http://www.linguist.univ-paris-diderot.fr/~bcrabbe/LingExp/cours7.pdf a écrit:
le problème que j'arrive pas à comprendre l'objet
- Code:
ndata$LTpred
Est ce que quelqu'un peux me aider ?
Merci beaucoup d'avance.
FMarwen- Nombre de messages : 226
Date d'inscription : 17/04/2013
Re: Courbe de Roc.
Bonjour,
Dans le code, ndata$LTpred désigne la variable des prédictions, le plus souvent issue du modèle logistique. Dans cet exemple, les prédictions se trouvent dans une table appelée ndata
Ayana
Dans le code, ndata$LTpred désigne la variable des prédictions, le plus souvent issue du modèle logistique. Dans cet exemple, les prédictions se trouvent dans une table appelée ndata
Ayana
Ayana- Nombre de messages : 550
Localisation : Londres
Date d'inscription : 18/08/2009
Re: Courbe de Roc.
autrement les LTpred sont les ajustements de modèles ?
FMarwen- Nombre de messages : 226
Date d'inscription : 17/04/2013
Re: Courbe de Roc.
Oui. Si tu utilises un modèle logistique voici le code :
- Code:
modele<-glm(groupe~variable,data=tab,family=binomial)
PS<-fitted.values(modele)
tab$PS<-PS
pred <- prediction(tab$PS, tab$groupe)
Ayana- Nombre de messages : 550
Localisation : Londres
Date d'inscription : 18/08/2009
Re: Courbe de Roc.
Désolé je te dérange ,mais ça marche pas ,Voici mon code :
- Code:
logit <- glm(L_T ~ HD + D13 + D_Re + Z, family=binomial(logit), data=MergeTestP1)
summary(logit)
ndata <- data.frame(HD=MergeTestP1$HD,D13=MergeTestP1$D13,D_Re=MergeTestP1$D_Re,Z=MergeTestP1$Z)
fix(ndata)
#Predicted Probabilities
ndata$predict<- predict(logit,ndata,type="response")
#Decision rule (theta = 0.7)
ndata$decision<- ifelse(ndata$predict > 0.2,"1","0")
#Add reference column
ndata$Ref <- MergeTestP1$L_T
#Confusion matrix
t <- table(ndata$Ref,ndata$decision)
t
sum(t)
#Accurracy, TPR,FPR
acc <- (t[1]+t[4]) / sum(t)
TPR <- t[1] / (t[1]+t[2])
FPR <- t[3] / (t[3]+t[4])
library(ROCR)
# Build a prediction object
ndata$fit<-fitted.values(logit)
preds <- prediction(ndata$fit,ndata$Ref)
#Compute the conf matrix for every value of theta
perf <- performance(preds,"TPR","FPR")
#Plot the ROC Curve
plot(perf,colorize=T)
#Compute AUC
> auc <- performance(preds,"auc")
FMarwen- Nombre de messages : 226
Date d'inscription : 17/04/2013
Re: Courbe de Roc.
Oulà! C'est à quel niveau que ça marche pas, et tu cherches à faire quoi?
Ayana- Nombre de messages : 550
Localisation : Londres
Date d'inscription : 18/08/2009
Re: Courbe de Roc.
Bon merci beaucoup @ Ayana , j'ai trouvé le problème (il faut installer le package {base}),sur ce lien
Peux être utile pour les visiteurs de la page.
http://freakonometrics.blog.free.fr/index.php?post/2012/01/11/Construire-une-courbe-ROC a écrit:
Peux être utile pour les visiteurs de la page.
FMarwen- Nombre de messages : 226
Date d'inscription : 17/04/2013
Re: Courbe de Roc.
questions de plus :
Est ce que possible de tracer deux courbe de Roc sur le même graphe ?
Merci beaucoup d'avance .
Est ce que possible de tracer deux courbe de Roc sur le même graphe ?
Merci beaucoup d'avance .
FMarwen- Nombre de messages : 226
Date d'inscription : 17/04/2013
Sujets similaires
» Courbe de gauss sous R
» Tracer des "error bar" sur une courbe
» courbe de référence
» Courbe par colonne .
» Taux d'incidence sous R
» Tracer des "error bar" sur une courbe
» courbe de référence
» Courbe par colonne .
» Taux d'incidence sous R
Page 1 sur 1
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum