Discussion :

[CupOfAppleTea]

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 :

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

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 :

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

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 ?

[arm3366]

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

Code :

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

?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.

[CupOfAppleTea]

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".

[dev_ggy]

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.

[dev_ggy]

Monsieur,

Code :

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

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.

[CupOfAppleTea]

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

[CupOfAppleTea]

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 :

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

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])'

[arm3366]

Bonjour,

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

Code :

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

is.data.frame(a)

BC

[CupOfAppleTea]

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