Forum de Statistiques
Vous souhaitez réagir à ce message ? Créez un compte en quelques clics ou connectez-vous pour continuer.
Les posteurs les plus actifs de la semaine

Le Deal du moment : -32%
-32% sur le purificateur d’air Dyson Pure Cool ...
Voir le deal
229 €

Courbe de Roc.

Aller en bas

Courbe de Roc. Empty Courbe de Roc.

Message par FMarwen le Mar 9 Juil 2013 - 10:00

Bonjour,

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
a la page 27 !
Est ce que quelqu'un peux me aider ?
Merci beaucoup d'avance.

FMarwen

Nombre de messages : 226
Date d'inscription : 17/04/2013

Revenir en haut Aller en bas

Courbe de Roc. Empty Re: Courbe de Roc.

Message par Ayana le Mar 9 Juil 2013 - 11:31

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
Ayana
Ayana

Nombre de messages : 539
Localisation : Londres
Date d'inscription : 18/08/2009

Revenir en haut Aller en bas

Courbe de Roc. Empty Re: Courbe de Roc.

Message par FMarwen le Mar 9 Juil 2013 - 11:41

autrement les LTpred sont les ajustements de modèles ?

FMarwen

Nombre de messages : 226
Date d'inscription : 17/04/2013

Revenir en haut Aller en bas

Courbe de Roc. Empty Re: Courbe de Roc.

Message par Ayana le Mar 9 Juil 2013 - 11:47

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
Ayana

Nombre de messages : 539
Localisation : Londres
Date d'inscription : 18/08/2009

Revenir en haut Aller en bas

Courbe de Roc. Empty Re: Courbe de Roc.

Message par FMarwen le Mar 9 Juil 2013 - 12:00

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

Revenir en haut Aller en bas

Courbe de Roc. Empty Re: Courbe de Roc.

Message par Ayana le Mar 9 Juil 2013 - 12:49

Oulà! C'est à quel niveau que ça marche pas, et tu cherches à faire quoi?
Ayana
Ayana

Nombre de messages : 539
Localisation : Londres
Date d'inscription : 18/08/2009

Revenir en haut Aller en bas

Courbe de Roc. Empty Re: Courbe de Roc.

Message par FMarwen le Mar 9 Juil 2013 - 12:56

Bon merci beaucoup @ Ayana Smile , j'ai trouvé le problème (il faut installer le package {base}),sur ce lien
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

Revenir en haut Aller en bas

Courbe de Roc. Empty Re: Courbe de Roc.

Message par FMarwen le Mar 9 Juil 2013 - 13:37

questions de plus :
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

Revenir en haut Aller en bas

Courbe de Roc. Empty Re: Courbe de Roc.

Message par Contenu sponsorisé


Contenu sponsorisé


Revenir en haut Aller en bas

Revenir en haut


 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum