droites de régression

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Roland Chiavassa
Messages : 24
Enregistré le : 29 Juin 2008, 09:13

droites de régression

Messagepar Roland Chiavassa » 13 Sep 2008, 13:05

Bonjour,
J'ai voulu, à titre pédagogique, montrer à mes étudiants, différents types de droite de

régression à savoir la droite de régression au sens de moindres carrés, la droite de

régression au sens des moindres rectangles, la droite LAD (least absolute deviation), la

droite LMS (least median of square). Pour les trois dernières, j'ai minimisé le critère

correspondant à l'aide de la fonction optim de R.
Ma question est la suivante : existe-t-il dans R une fonction, un package qui permet de

trouver directement les coefficients de ces droites. Ce qui me permettrait de vérifier mes

résultats !
Avec mes remerciements
Roland

Renaud Lancelot
Messages : 2484
Enregistré le : 16 Déc 2004, 08:01
Contact :

Messagepar Renaud Lancelot » 13 Sep 2008, 15:20

R dispose de fonctions de recherche puissantes. Faire par exemple

Code : Tout sélectionner

RSiteSearch("least absolute deviation")


Renaud

François Bonnot
Messages : 537
Enregistré le : 10 Nov 2004, 15:19
Contact :

Messagepar François Bonnot » 15 Sep 2008, 07:23

Bonjour,
J'ai été amené à écrire récemment une fonction pour calculer la droite des moindres rectangles (inutile d'utiliser optim), la voici pour information.

Code : Tout sélectionner

frect <- function(x,y){ # droite des moindres rectangles
  b <- sqrt(var(y)/var(x))*sign(cor(x,y))
  c(mean(y)-b*mean(x),b)
}

x<-rnorm(1000,5,1)
y<-100+7*x+10*rnorm(1000)

plot(x,y)
abline(lm(y~x))
res <- lm(x~y)$coeff
abline(-res[1]/res[2],1/res[2])
abline(frect(x,y),col="red")


FB

jean lobry
Messages : 733
Enregistré le : 17 Jan 2008, 20:00
Contact :

Messagepar jean lobry » 15 Sep 2008, 15:50

Bonjour,

je ne connais pas de paquet R où toutes ces fonctions seraient disponibles (ce qui ne signifie pas qu'il n'existe pas). Pour des illustrations pédagogiques il est intéressant de signaler que le jeu de "données" d'Anscombe de 1973 est disponible en standard dans R. Essayer :

Code : Tout sélectionner

example(anscombe)


Amicalement,

Jean


Retourner vers « Questions en cours »

Qui est en ligne

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

cron