Discussion :

[JechercheT]

Bonjour,
Le code suivant :

Code :

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

lm(DDenc~DTCR)

me retourne l'erreur :

les longueurs des variables diffèrent (trouvé pour 'DTCR')

sachant que la longueur de DTCR est 31 et celle de DDenc est 37.

J'ai essayé:

Code :

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

lm(DDenc~DTCR, na.action="na.exclude")

Toujours le même message !
Cela veut dire que la régression linéaire sous R ne se fait que par des variables de mêmes longueurs ou il y a une solution ?

[faubry]

cela veut dire que la regression linéaire sous R ne se fait que par des variables de meme longueurs ou il y a une solution??

En effet et par définition même de l'analyse linéaire. Il faut qu'à chaque valeur de DDenc corresponde une valeur de DTCR, ce qui n'est pas le cas si les longueurs des deux vecteurs différent. Il faut obligatoirement des NA aux positions non définies. Exemple :

Code :

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

1
2
DDenc <- c( 1, 2, 3)
DTCR <- c( 2, 5)

L'analyse ne peut pas se faire. Par contre, si c'est pour DDenc = 2 que DTCR n'est pas défini, alors on peut analyser :

Code :

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

1
2
DDenc <- c( 1, 2, 3)
DTCR <- c( 2, NA, 5)

En effet, les deux vecteurs ont alors la même longueur et l'analyse ignorera la deuxième valeur de DDenc et de DTCR. Sinon, comment lm peut-il deviner quelles valeurs de DDenc n'a pas de DTCR correspondant ?

C'est aussi vrai si c'était DDenc le plus court.

[JechercheT]

Et comment je peux insérer ces NA sachant que les deux variables ainsi mentionnées sont, déjà, des output de R et je sais pas exactement quelles sont les cellules vides?

[faubry]

Les deux variables proviennent soit de la lecture de fichiers (par exemple, par la fonction read.table), soit de calculs, R ne créant pas des output ex-nihilo. Si tu cherches une relation du type DDenc ~ DTCR, cela signifie que tu connais la fonction au sens mathématique du terme qui est l'ensemble des couples qui associent une valeur de DDenc à DTCR (produit cartésien). Le NA (Not Available) signifie simplement que tu ne connais pas la valeur exacte, soit qu'elle n'a pas été rapportée, soit qu'elle n'est pas calculable pour diverses raisons, mais qu'elle existe dans la réalité. En d'autres termes, tu as un modèle, cette fonction. C'est donc à toi de les insérer à la bonne place. Si c'est le code R que tu veux, voici un exemple

Code :

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

1
2
3
4
 
x <- 1:11
# Je veux inserer un NA entre le 6 et le 7
x.NA <- c( head(x, 6), NA, tail( x, -6))

[JechercheT]

Le DTCR=diff(TCR) la différence première de TCR, le DDenc est la différence seconde de Enc.
Je voudrais savoir s'il existe une relation linéaire entre les deux variables en sorties (DTCR et DDenc).
Mais je ne sais pas quel est le rang des valeurs manquantes pour ajouter ces NA.