<div>Olá pessoal, estou realizando uma ANOVA com medidas repetidas e estou utilizando a função "Anova" do pacote "car".</div><div><br></div><div>Medi o biovolume de algas a cada dois dias durante 10 dias (no banco de dados abaixo só coloquei até o 4° dia). Tenho 2 tratamentos ("c","t") e o experimento foi realizado em tréplicas ("A","B","C"). </div>
<div><br></div><div>> Pa2</div><div>   Day Type Replicate logbiovolume</div><div>1    0    c         A        19.34</div><div>2    0    c         B        18.27</div><div>3    0    c         C        18.56</div><div>4    0    t         A        18.41</div>
<div>5    0    t         B        18.68</div><div>6    0    t         C        18.86</div><div>7    2    c         A        18.81</div><div>8    2    c         B        18.84</div><div>9    2    c         C        18.52</div>
<div>10   2    t         A        18.29</div><div>11   2    t         B        17.91</div><div>12   2    t         C        17.67</div><div>13   4    c         A        19.16</div><div>14   4    c         B        18.85</div>
<div>15   4    c         C        19.36</div><div>16   4    t         A        19.05</div><div>17   4    t         B        19.09</div><div>18   4    t         C        18.26</div><div><br></div><div>Pa2.teste = within(Pa2,{group = factor(Type)</div>
<div>                   time = factor(Day)</div><div>                   id = factor(Replicate)})</div><div>matrix = with(Pa2.teste,cbind(Pa2[,VAR][group=="c"],Pa2[,VAR][group=="t"]))</div><div>matrix</div>
<div>       [,1]  [,2]</div><div> [1,] 19.34 18.41</div><div> [2,] 18.27 18.68</div><div> [3,] 18.56 18.86</div><div> [4,] 18.81 18.29</div><div> [5,] 18.84 17.91</div><div> [6,] 18.52 17.67</div><div> [7,] 19.16 19.05</div>
<div> [8,] 18.85 19.09</div><div> [9,] 19.36 18.26</div><div>[10,] 19.63 18.96</div><div>[11,] 19.94 18.06</div><div>[12,] 19.54 18.37</div><div>[13,] 19.98 17.96</div><div>[14,] 20.99 17.93</div><div>[15,] 20.45 17.74</div>
<div>[16,] 21.12 17.60</div><div>[17,] 21.66 17.33</div><div>[18,] 21.51 18.12</div><div> model <- lm(matrix ~ 1)</div><div> design <- factor(c("c","t"))</div><div> </div><div> options(contrasts=c("contr.sum", "contr.poly"))</div>
<div> aov <- Anova(model, idata=data.frame(design), idesign=~design, type="III")</div><div> summary(aov, multivariate=F)</div><div><br></div><div>Univariate Type III Repeated-Measures ANOVA Assuming Sphericity</div>
<div><br></div><div>                 SS num Df Error SS den Df         F    Pr(>F)    </div><div>(Intercept) 12951.2      1   6.3312     17 34775.336 < 2.2e-16 ***</div><div>design         19.1      1  17.3901     17    18.697 0.0004606 ***</div>
<div>---</div><div>Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 </div><div> </div><div><br></div><div>O problema é que eu acho que esta função não está levando em consideração os dias, nem as réplicas. Como faço para introduzir isto na função. Vocês conhecem alguma função correspondente não paramétrica para este teste? Tipo um teste de Friedman com dois grupos (tratamento e réplica) e um bloco (tempo)? Pelo que entendi, a função friedman.test só faz com um grupo apenas.</div>
<div><br></div><div>Muito Obrigado</div><div> </div>                                               Diego PJ<br>