un question qui va peut être vous paraître simple.
Code : Tout sélectionner
> # modèle complet
> a3m0 <- ASRF3~ty*co*ta*an
> a3aov0 <- aov(a3m0, data=fcb)
> monAIC(a3aov0)
n k nsk AIC AICc
a3aovmodele 226 17 13.29412 300.3598 303.3021
> a3m1 <- ASRF3~ty*co*ta*an-ty:co:ta:an
> a3aov1 <- aov(a3m1, data=fcb)
> a3m2a <- ASRF3~ty*co*ta*an-ty:co:ta:an-ty:co:ta # retrait de - ty:co:ta
> a3aov2a <- aov(a3m2a, data=fcb)
etc.etc...
ici je déclare la réponse d'étude comme étant ASRF3.
hors, je dois réaliser le même type d'analyses un grand nombre de fois, sur plusieurs réponses.
j'aurais donc voulu pouvoir déclarer une fois ma réponse et ne pas modifier manuellement dans un éditeur de texte tous les numéros de réponse (ASRF5, ASRF6, etc...)
j'ai essayé :
Code : Tout sélectionner
> Rep <- ASRF3
Erreur : objet "ASRF3" non trouvé
(en effet, ASRF3 n'est pas un objet mais un nom de colonne)
Code : Tout sélectionner
> Rep <- "ASRF3"
> a3m0 <- Rep~ty*co*ta*an
> a3aov0 <- aov(a3m0, data=fcb)
Erreur dans model.frame(formula, rownames, variables, varnames, extras, extranames, :
les longueurs des variables diffèrent (trouvé pour 'ty')
> Rep <- 'ASRF3'
> a3m0 <- Rep~ty*co*ta*an
> a3aov0 <- aov(a3m0, data=fcb)
Erreur dans model.frame(formula, rownames, variables, varnames, extras, extranames, :
les longueurs des variables diffèrent (trouvé pour 'ty')
Du coup, en raison de cette erreur (dont je ne comprend pas bien la signification), les calculs suivant ne peuvent être effectués...
N'est il pas possible de déclarer spécifiquement ce type d'information ?