syntaxe_nombreuses variables à expliquer_LM

Mon objectif est de procéder à une analyse canonique discriminante.
Pour faire simple, je dispose du jeu de données suivant:
_Un tableau de 680 colonnes.
_ la première colonne identifie les individus sur lesquels j'ai mesuré chacune des variables.
_ les colonnes 2, 3 et 4 contiennent des variables explicatives de type facteur.
_ les colonnes 5 à 680 contiennent les variables à expliquer.

Je dois d'abord faire un modèle linéaire que je rentrerai ensuite dans la formule de l'analyse canonique discriminante.
Je dois donc, dans mon modèle linéaire, considérer mes 3 variables explicatives ainsi que mes 676 variables à expliquer.
La formule du modèle linéaire étant du type : y~x je dois rentrer 676 variables dans y. Ce que je cherche, c'est une formule me permettant de ne pas avoir à spécifier mes 676 variables une par une. Je voudrais pouvoir les rentrer en une seule fois. Un truc du type:
(avec "a" mon tableau de données)

Code:

---

mod=lm(a[,5:680]~a$Var1*a$Var2*a$Var3)

---

Seulement, "a[,5:680]" ne fonctionne pas. J'obtiens le message d'erreur suivant:

Code:

---

1 2	Error in model.frame.default(formula = a[, 5:680] ~ flatitude * ftype * : invalid type (list) for variable 'a[, 5:680]'

---

Connaissez-vous un moyen de rentrer toutes mes colonnes d'un coup sans avoir à les spécifier une par une ?

Bonjour,
Je pense qu'il faudrait commencer à regarder la syntaxe et le fonctionnement de la fonction lm en tapant :

Code:

---

?lm

---

Ce n'est pas la fonction adaptée à ta problématique à mon avis. Essaye plutôt d'explorer la fonction lda pour ton cas.

Je souhaite en fait procéder à une analyse canonique discriminante, via la fonction "candisc" :
http://127.0.0.1:29294/library/candi...l/candisc.html

Comme vous pourrez le voir sur cette page internet dans la rubrique "exemples", la fonction "candisc" nécessite en entrée un modèle linéaire.

C'est pourquoi je me pose le problème que je vous ai soumis dans mon précédent message et dont je n'ai pas trouvé la réponse sur les pages d'aide concernant la fonction "lm".

Citation:

---

http://127.0.0.1:29294/library/candi...l/candisc.html

---

Ce lien est un lien vers un serveur local. Nous n'y avons pas accès.

Bien cordialement.

Monsieur,

Code:

---

1 2 3 4 5

data(Grass) grass.mod <- lm(cbind(N1,N9,N27,N81,N243) ~ Block + Species, data=Grass) Anova(grass.mod,test="Wilks") grass.can1 <-candisc(grass.mod, term="Species") plot(grass.can1, type="n")

---

Ci-dessus le code pour l'exemple de la fonction candisc.

dans l'exemple de la fonction que vous citez candisc on fait une opération supplémentaire que vous je faite pas avec la fonction cbind.

Si vous la faite qu'obtenez vous ?

bien cordialement.

Dsl, il s'agit de la fonction "candisc" du package "candisc".
Mais j'ai bien peur qu'elle ne soulève pas les foules... :p

Ayant, 676 variables à expliquer, je ne souhaite pas énumérer chacune des colonnes correspondantes avec un cbind().
C'est pourquoi je souhaiterais pouvoir trouver un moyen de "rentrer" ces 676 variables d'un seul coup comme avec un "a[,5:680]" (=colonne 5 à 680 de mon tableau "a") qui ne fonctionne pas:

Code:

---

1 2 3

> abs.mod=lm(cbind(a[,5:680])~flatitude*ftype*forientation,data=a) Error in model.frame.default(formula = cbind(a[, 5:680]) ~ flatitude * : invalid type (list) for variable 'cbind(a[, 5:680])'

---

Bonjour,

assure toi que ton fichier "a" est bien un data.frame

Code:

---

is.data.frame(a)

---

BC

C'est bien un data.frame en effet, et j'obtiens le même message d'erreur que précédemment.

Bonsoir,

ce n'est pas la solution la plus "jolie" mais cela devrait tourner:

Code:

---

1 2 3

test <- a[,5] for(i in 6:680){test <- cbind(test, a[,i])} lm(test ~ flatitude * ftype * forientation, data = a)

---

A améliorer avec un apply éventuellement...

Bonne continuation :)

Merci, ça fonctionne! Par contre je ne comprends pas pourquoi ça fonctionne de cette manière et pas de l'autre? (celle que je tentais d'utiliser)
J'aimerais comprendre...