[R] please Help me on a repeated measures anova

soileil soileil at msn.com
Thu Dec 30 11:38:58 CET 2010


So far I've done this:


> tableau = read.table("Classeur1.txt", h=T)
> tableau
    Temps.s. Modalite  AchE
1         15       ER   0,1
2         15       ER 0,112
3         15    T+0,5 0,078
4         15    T+0,5 0,053
5         15      TED 0,107
6         15      TED 0,117
7         15      TED 0,124
8         15      T+1 0,109
9         15      T+1 0,118
10        15   SED+ER 0,158
11        15   SER+ED 0,119
12        15   SED+ED 0,101
13        15   SED+ED 0,129
14        30       ER  0,11
15        30       ER 0,134
16        30    T+0,5 0,082
17        30    T+0,5   0,1
18        30      TED 0,126
19        30      TED 0,135
20        30      TED 0,145
21        30      T+1 0,119
22        30      T+1 0,134
23        30   SED+ER 0,175
24        30   SER+ED  0,14
25        30   SED+ED 0,113
26        30   SED+ED 0,135
27        30       ER 0,123
28        30       ER 0,153
29        30    T+0,5 0,088
30        30    T+0,5 0,109
31        30      TED 0,141
32        30      TED 0,153
33        30      TED 0,163
34        30      T+1 0,134
35        30      T+1 0,153
36        30   SED+ER 0,194
37        30   SER+ED 0,157
38        30   SED+ED  0,18
39        30   SED+ED  0,14
40        60       ER 0,135
41        60       ER 0,174
42        60    T+0,5 0,094
43        60    T+0,5  0,12
44        60      TED 0,161
45        60      TED 0,169
46        60      TED 0,187
47        60      T+1 0,148
48        60      T+1  0,17
49        60   SED+ER 0,208
50        60   SER+ED 0,174
51        60   SED+ED  0,14
52        60   SED+ED 0,146
53        75       ER 0,147
54        75       ER 0,192
55        75    T+0,5 0,101
56        75    T+0,5 0,127
57        75      TED 0,172
58        75      TED 0,183
59        75      TED 0,208
60        75      T+1 0,163
61        75      T+1 0,184
62        75   SED+ER 0,226
63        75   SER+ED 0,192
64        75   SED+ED 0,154
65        75   SED+ED 0,153
66        90       ER 0,159
67        90       ER 0,208
68        90    T+0,5 0,108
69        90    T+0,5 0,136
70        90      TED 0,184
71        90      TED 0,201
72        90      TED 0,227
73        90      T+1 0,174
74        90      T+1 0,197
75        90   SED+ER  0,24
76        90   SER+ED 0,208
77        90   SED+ED 0,166
78        90   SED+ED 0,159
79       105       ER 0,171
80       105       ER 0,226
81       105    T+0,5 0,113
82       105    T+0,5 0,145
83       105      TED   0,2
84       105      TED 0,218
85       105      TED 0,244
86       105      T+1 0,187
87       105      T+1 0,214
88       105   SED+ER 0,259
89       105   SER+ED 0,225
90       105   SED+ED 0,179
91       105   SED+ED 0,165
92       120       ER 0,182
93       120       ER 0,251
94       120    T+0,5  0,12
95       120    T+0,5 0,154
96       120      TED 0,213
97       120      TED 0,229
98       120      TED 0,259
99       120      T+1 0,202
100      120      T+1 0,228
101      120   SED+ER 0,268
102      120   SER+ED  0,24
103      120   SED+ED  0,19
104      120   SED+ED 0,172
> attach (tableau)


> tableau <- lm (AchE ~ Temps.s. * Modalite, random = ~ 1 | individu)
Erreur dans storage.mode(y) <- "double" : 
  la modification du mode de stockage d'un objet 'factor' n'est pas
autorisée
De plus : Messages d'avis :
1: In model.response(mf, "numeric") :
  l'utilisation de type="numeric" avec une réponse de type facteur sera
ignorée
2: In lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
  les arguments surnuméraires random sont ignorés.

What is my mistake?
-- 
View this message in context: http://r.789695.n4.nabble.com/please-Help-me-on-a-repeated-measures-anova-tp3159868p3168134.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list