Trabalho de casa? Não fazemos trabalhos de casa. Ao menos podia ter tentado correr esse código...
Rui Barradas Em 05-10-2012 14:56, Diego Pujoni escreveu: > Olá pessoal, estou realizando uma ANOVA com medidas repetidas e estou > utilizando a função "Anova" do pacote "car". > > 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"). > >> Pa2 > Day Type Replicate logbiovolume > 1 0 c A 19.34 > 2 0 c B 18.27 > 3 0 c C 18.56 > 4 0 t A 18.41 > 5 0 t B 18.68 > 6 0 t C 18.86 > 7 2 c A 18.81 > 8 2 c B 18.84 > 9 2 c C 18.52 > 10 2 t A 18.29 > 11 2 t B 17.91 > 12 2 t C 17.67 > 13 4 c A 19.16 > 14 4 c B 18.85 > 15 4 c C 19.36 > 16 4 t A 19.05 > 17 4 t B 19.09 > 18 4 t C 18.26 > . > . > . > > Pa2.teste = within(Pa2,{group = factor(Type) > time = factor(Day) > id = factor(Replicate)}) > matrix = with(Pa2.teste,cbind(Pa2[,VAR][group=="c"],Pa2[,VAR][group=="t"])) > matrix > [,1] [,2] > [1,] 19.34 18.41 > [2,] 18.27 18.68 > [3,] 18.56 18.86 > [4,] 18.81 18.29 > [5,] 18.84 17.91 > [6,] 18.52 17.67 > [7,] 19.16 19.05 > [8,] 18.85 19.09 > [9,] 19.36 18.26 > [10,] 19.63 18.96 > [11,] 19.94 18.06 > [12,] 19.54 18.37 > [13,] 19.98 17.96 > [14,] 20.99 17.93 > [15,] 20.45 17.74 > [16,] 21.12 17.60 > [17,] 21.66 17.33 > [18,] 21.51 18.12 > model <- lm(matrix ~ 1) > design <- factor(c("c","t")) > > options(contrasts=c("contr.sum", "contr.poly")) > aov <- Anova(model, idata=data.frame(design), idesign=~design, type="III") > summary(aov, multivariate=F) > > Univariate Type III Repeated-Measures ANOVA Assuming Sphericity > > SS num Df Error SS den Df F Pr(>F) > (Intercept) 12951.2 1 6.3312 17 34775.336 < 2.2e-16 *** > design 19.1 1 17.3901 17 18.697 0.0004606 *** > --- > Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 > > > 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 análise. 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)? > > Muito Obrigado > > Diego PJ > > [[alternative HTML version deleted]] > > > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. [[alternative HTML version deleted]]
______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.