Discussion :

[mauriceletendu]

Bonjour,

J'ai crée un modèle de régression linéaire de température en fonction de l'altitude, ma variable contenant le modèle s'appelle "modele"

J'ai crée ce modèle à partir d'une centaine de données.

Maintenant j'ai une matrice contenant une quantité assez importante d'altitude. Je voudrais calculer les températures à partir du modèle de regression.

Comment faire ? En sachant que je peux convertir ma matrice en autre structure, c'est pas très important pour le moment.

En clair je veux ça :
Données d'entrée:Altitude ----> mon modele -----> Données de sortie:température

Quelqu'un sait comment faire ?

Merci

[Emmanuel R]

Bonjour,

Utilisez la fonction predict.lm. Pour obtenir de l'aide à ce sujet, tapez : ?predict.lm
En lisant l'aide et en ajoutant bien le paramètre type="response", vous devriez arriver à faire ce que vous souhaitez.

Bon courage

[mauriceletendu]

Je m'en sort pas, voici ce que je fais:

j'ai 2 matrices: une "temp" avec 204 températures et beaucoup de valeurs manquantes
Et une matrice "alt" de la même dimension de "temp" mais avec toutes les altitudes.

En clair j'ai 204 points où j'ai mes températures, je veux faire un modèle de regression linéaire en fonction de mes altitudes pour pouvoir interpoler mes données sur toute la matrice.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
modele <- lm(temp[temp>0]~alt[temp>0], na.action=na.omit)  # je veux mes température en fonction de mes altitudes
result <- predict.lm(modele,alt,type="response", na.action=na.omit) # j'ai essayé plusieurs commandes avec predict aussi (ici alt est transformé en data.frame)

Mais dans chaque cas, j'ai l'erreur:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
Message d'avis :
'newdata' had 29036739 rows but variables found have 204 rows

Et ça me calcule seulement 204 valeurs.... alors que je rentre des millions d'altitude (je veux en sortie mes millions de températures, calculées en fonction de mes altitudes...)

De plus, R m'oblige à convertir ma matrice alt en data.frame....

[Sengar]

Tu peux nous montrer tes matrices? (au moins l'entête avec head) et si possible le résultat que tu souhaiterais obtenir sur les premières lignes par exemple?

Quand je vois ton modèle sans en savoir plus sur la structure des données, je ne comprend pas ce que tu fais.

D'autre part, s'il s'agit juste de faire une interpolation, tu peux directement utiliser une fonction comme approxfun().

[mauriceletendu]

En fait je cherche juste à utiliser le modèle de régression linéaire.

Le modèle de régression c'est y = ax+b

Dans le modèle, il y a les valeurs des coefficients a et b.

Je veux savoir comment les utiliser.

Imaginons que mon modèle soit fait à partir de 100 données.

maintenant j'ai ceci par exemple:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

x <- c(2,451,256,21,45,87,35)

Je veux calculer mes y.

Comment faire ?

Dans un premier temps j'ai fait ceci:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
b <- modele$coefficients[1]
a <- modele$coefficients[2]
y <- (a*x + b)

Mais il doit y avoir plus simple non ?

[A. D.]

Bonjour,

La réponse vous a déjà été donnée :

Utilisez la fonction predict.lm. Pour obtenir de l'aide à ce sujet, tapez : ?predict.lm
En lisant l'aide et en ajoutant bien le paramètre type="response", vous devriez arriver à faire ce que vous souhaitez.

Petite précision : l'élément à spécifier pour le paramètre "newdata" doit être de type "dataframe" (comme cela est indiqué dans la documentation de la fonction "predict.lm").

Bonne continuation


Cordialement,


A.D.