Regression logarithmique

Questions sur les fonctions statistiques de R

Modérateur : Groupe des modérateurs

LA haumont
Messages : 41
Enregistré le : 16 Mai 2006, 09:33

Regression logarithmique

Messagepar LA haumont » 25 Juil 2006, 13:17

bonjour,
on me demande de réaliser une régression logarithmique mais je ne toruve quasi rien dessus. pas de commandes , d'explication!!
vous auriez des pistes?
merci

Tillard
Messages : 87
Enregistré le : 17 Déc 2004, 10:32

regression logarithmique ???

Messagepar Tillard » 25 Juil 2006, 17:19

Bonjour
jamais entendu parler de ça

regression : y = ax + b
regression logarithmique: log(y) = ax + b ou y = alog(x) + b

???????????

il faut regarder la distribution de vos données (y et x)

cordialement
Emmanuel Tillard
UMR ERRC (Elevage des Ruminants en Regions Chaudes)
CIRAD - St PIERRE (La Réunion)
tel: 02 62 49 92 54

LA haumont
Messages : 41
Enregistré le : 16 Mai 2006, 09:33

Messagepar LA haumont » 26 Juil 2006, 08:52

en fait je crois que c'est même: log(y) = alog(x) + b mais le problème c'est que j'ai des Na dans une de mes colonnes et donc ca pose problème!!

LA haumont
Messages : 41
Enregistré le : 16 Mai 2006, 09:33

Messagepar LA haumont » 02 Aoû 2006, 12:59

Bonjour,
après de nombreux essais je ne trouve pas comment faire ma régression logarithmique. En effet, lorsque je passe ma variable en log elle devient en character

Code : Tout sélectionner

> diff2angiolog <- log(regroup2angio2$diff2angio)
> fix(diff2angiolog)
> diff2angiolog2 <- gsub("-Inf", "0.0000",diff2angiolog, fixed=TRUE)
> fix(diff2angiolog2)
> mode(diff2angiolog2))
Erreur : erreur de syntaxe dans "mode(diff2angiolog2))"
> mode(diff2angiolog2)
[1] "character"


ce qui m'empêche par la suite de faire ma régression.
Quelqu'un voit comment je pourrais m'en sortir?
merci beaucoup

Tillard
Messages : 87
Enregistré le : 17 Déc 2004, 10:32

regression logarithmique

Messagepar Tillard » 03 Aoû 2006, 08:07

Bonjour

Les transformations log posent un problème lorsque qu'il y a des valeurs nulles

Voir le livre de Venables and Ripley:
Venables WN, Ripley BD. 1994. Modern applied statistics with S-Plus. Springer.
Un exemple de transformation log est donné page 171 et 172

le modele est: log(y + A) = B0 + B1(Age) + B2(Sex) + B3(Eth) + B4(Lrn)
La question: comment choisir A de facon a maximiser la vraisemblance marginale ?

Code : Tout sélectionner

library(MASS)
quine$Days  #il y a des 0 dans Days

?logtrans
logtrans(Days ~ Age+Sex+Eth+Lrn, quine, alpha = seq(0.75, 6.5, le=20))

#Le graphe produit montre clairement que la vraisemblance marginale est maximale pour A = 2.2 environ

#donc, le modèle peut s'écrire:
#log(y + 2.2) = B0 + B1(Age) + B2(Sex) + B3(Eth) + B4(Lrn)

summary(lm(log(Days + 2.2) ~ Age+Sex+Eth+Lrn, quine))


Cordialement
Emmanuel Tillard

UMR ERRC (Elevage des Ruminants en Regions Chaudes)

CIRAD - St PIERRE (La Réunion)

tel: 02 62 49 92 54

LA haumont
Messages : 41
Enregistré le : 16 Mai 2006, 09:33

Messagepar LA haumont » 03 Aoû 2006, 09:15

merci beaucoup j'obtiens quelques résultats.
Par contre, je ne comprends pas a quoi correspondent exactement les éléments dans seq. En faisant divers modèles, j'obtiens toujours le même graph! Et la valeur du A on la regarde bien sur l'axe des abscisses? Il ne s'agit pas d'une valeur sûre? plutot approximative?
Aussi je ne vois pas ce que représente exactement l'axe des ordonnés , j'ai des chiffres allant de -5680 a -5780!!!

Tillard
Messages : 87
Enregistré le : 17 Déc 2004, 10:32

regression logarithmique

Messagepar Tillard » 03 Aoû 2006, 12:56

seq représente une fourchette de valeur plausible pour A
on arrive a cerner le probleme apres quelques tatonnements
A se lie bien sur l'axe des x, au niveau ou la courbe présente un maximum

pour avoir la valeur exacte:

Code : Tout sélectionner

val.logtrans <- logtrans(Days ~ Age+Sex+Eth+Lrn, quine, alpha = seq(0.75, 6.5, le=20))

max.val.logtrans <- approx(val.logtrans$y, val.logtrans$x, xout=max(val.logtrans$y)
max.val.logtrans


pouvez vous envoyer vos données ou un extrait pour qu'on puisse refaire l'analyse
Emmanuel Tillard

UMR ERRC (Elevage des Ruminants en Regions Chaudes)

CIRAD - St PIERRE (La Réunion)

tel: 02 62 49 92 54

Vincent Andres
Messages : 13
Enregistré le : 17 Juil 2006, 09:27

Re: Regression logarithmique

Messagepar Vincent Andres » 04 Aoû 2006, 08:30

LA haumont a écrit :bonjour,
on me demande de réaliser une régression logarithmique mais je ne toruve quasi rien dessus. pas de commandes , d'explication!!
vous auriez des pistes?
merci

Peut-être "on" demande-t-il tout bêtement une régression non pas sur les valeurs brutes, mais sur leurs variations (ou "rendements") exprimés en log returns ?
C'est quand même assez courant de travailler sur les returns.
C'est juste une piste !


Retourner vers « Archives : Fonctions statistiques »

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité