Docsity
Docsity

Prepare-se para as provas
Prepare-se para as provas

Estude fácil! Tem muito documento disponível na Docsity


Ganhe pontos para baixar
Ganhe pontos para baixar

Ganhe pontos ajudando outros esrudantes ou compre um plano Premium


Guias e Dicas
Guias e Dicas


ANOVA e teste de comparacao de medias, Trabalhos de Estatística

anova de varios dados com testes parametrcicos

Tipologia: Trabalhos

2020

Compartilhado em 27/06/2020

nisce-silva-2
nisce-silva-2 🇧🇷

2 documentos

1 / 7

Toggle sidebar

Esta página não é visível na pré-visualização

Não perca as partes importantes!

bg1
Orlenisce Silva Passos Aragão
classe=factor(peixe$classe)
> levels(classe)
[1] "Bagre" "Carpa" "Cascudo" "Dourado" "Pacu" "Pintado" "Tilapia"
> tecnica=factor(peixe$tecnica)
> levels(tecnica)
[1] "Bow" "Forma" "Textura"
> desempenho=factor(peixe$desempenho)
> levels(desempenho)
[1] "70" "72" "74" "76" "78" "80" "81" "82" "83" "84" "85" "86" "87" "90" "92" "93"
"99"
> anova=aov(desempenho~classe+tecnica, data=peixe)
> anova1=aov(desempenho~classe*tecnica, data=peixe)
> #tabela ANOVA
> summary(anova)
Df Sum Sq Mean Sq F value Pr(>F)
classe 6 464.5 77.41 4.969 0.008914 **
tecnica 2 578.4 289.19 18.562 0.000212 ***
Residuals 12 187.0 15.58
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> summary(anova1)
Df Sum Sq Mean Sq
classe 6 464.5 77.41
tecnica 2 578.4 289.19
classe:tecnica 12 187.0 15.58
> coef(anova)
(Intercept) classeCarpa classeCascudo classeDourado classePacu classePintado
79.476190 6.000000 8.666667 10.333333 11.333333 12.666667
classeTilapia tecnicaForma tecnicaTextura
15.666667 -12.714286 -4.714286
> res=residuals(anova)
> res
1 2 3 4 5 6 7
3.2380952 -0.7619048 -1.4285714 -1.0952381 -0.0952381 0.5714286 -0.4285714
8 9 10 11 12 13 14
1.2380952 6.2380952 -0.4285714 -1.0952381 -1.0952381 -1.4285714 -3.4285714
15 16 17 18 19 20 21
-9.4761905 -0.4761905 1.8571429 2.1904762 1.1904762 0.8571429 3.8571429
> fitted(anova)
1 2 3 4 5 6 7 8
66.76190 72.76190 75.42857 77.09524 78.09524 79.42857 82.42857 74.76190
9 10 11 12 13 14 15 16
80.76190 83.42857 85.09524 86.09524 87.42857 90.42857 79.47619 85.47619
17 18 19 20 21
88.14286 89.80952 90.80952 92.14286 95.14286
pf3
pf4
pf5

Pré-visualização parcial do texto

Baixe ANOVA e teste de comparacao de medias e outras Trabalhos em PDF para Estatística, somente na Docsity!

Orlenisce Silva Passos Aragão classe=factor(peixe$classe)

levels(classe) [1] "Bagre" "Carpa" "Cascudo" "Dourado" "Pacu" "Pintado" "Tilapia" tecnica=factor(peixe$tecnica) levels(tecnica) [1] "Bow" "Forma" "Textura" desempenho=factor(peixe$desempenho) levels(desempenho) [1] "70" "72" "74" "76" "78" "80" "81" "82" "83" "84" "85" "86" "87" "90" "92" "93" "99" anova=aov(desempenho~classe+tecnica, data=peixe) anova1=aov(desempenho~classetecnica, data=peixe) #tabela ANOVA summary(anova) Df Sum Sq Mean Sq F value Pr(>F) classe 6 464.5 77.41 4.969 0.008914 ** tecnica 2 578.4 289.19 18.562 0.000212 *** Residuals 12 187.0 15. Signif. codes: 0 ‘’ 0.001 ‘’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 summary(anova1) Df Sum Sq Mean Sq classe 6 464.5 77. tecnica 2 578.4 289. classe:tecnica 12 187.0 15. coef(anova) (Intercept) classeCarpa classeCascudo classeDourado classePacu classePintado 79.476190 6.000000 8.666667 10.333333 11.333333 12. classeTilapia tecnicaForma tecnicaTextura 15.666667 -12.714286 -4. res=residuals(anova) res 1 2 3 4 5 6 7 3.2380952 -0.7619048 -1.4285714 -1.0952381 -0.0952381 0.5714286 -0. 8 9 10 11 12 13 14 1.2380952 6.2380952 -0.4285714 -1.0952381 -1.0952381 -1.4285714 -3. 15 16 17 18 19 20 21 -9.4761905 -0.4761905 1.8571429 2.1904762 1.1904762 0.8571429 3. fitted(anova) 1 2 3 4 5 6 7 8 66.76190 72.76190 75.42857 77.09524 78.09524 79.42857 82.42857 74. 9 10 11 12 13 14 15 16 80.76190 83.42857 85.09524 86.09524 87.42857 90.42857 79.47619 85. 17 18 19 20 21 88.14286 89.80952 90.80952 92.14286 95.

boxplot(desempenho~classe+tecnica, data=peixe) plot(anova) plot(residuals(anova), type="l")

levels(densidade) [1] "4200" "4300" "4400" "4500" "4600" "4800" "5000" "5800" "6000" "6300" "7000" "7500" "7700" [14] "8000" "8500" "9000" "9500" "10000" "11000" "12000" "12500" "13000" "14000" "14500" "15000" "16000" [27] "17500" "18000" "19000" "25000" anova=aov(densidade~ponto+estacao, data=bifatorial) #tabela ANOVA summary(anova) Df Sum Sq Mean Sq F value Pr(>F) ponto 1 409921714 409921714 30.053 6.88e-07 *** estacao 3 267793333 89264444 6.544 0.000601 *** Residuals 67 913884952 13640074


Signif. codes: 0 ‘’ 0.001 ‘’ 0.01 ‘’ 0.05 ‘.’ 0.1 ‘ ’ 1

coef(anova) (Intercept) ponto estacaooutono estacaoprimavera estacaoverao 13362.222 -1397.143 3177.778 -1333.333 -1733. res=residuals(anova) res 1 2 3 4 5 6 7 8 9 -5731.7460 -5231.7460 2768.2540 -142.8571 -5142.8571 -4142.8571 -2965.0794 - 2965.0794 -2465. 10 11 12 13 14 15 16 17 18 -1631.7460 -631.7460 -5631.7460 -4334.6032 -2834.6032 8665.3968 3754. 254.2857 11254. 19 20 21 22 23 24 25 26 27 -1567.9365 4432.0635 2432.0635 -234.6032 -3234.6032 -2234.6032 -2937.4603 - 2637.4603 11562. 28 29 30 31 32 33 34 35 36 5651.4286 2151.4286 5651.4286 2829.2063 6829.2063 3829.2063 5162. 4662.5397 4162. 37 38 39 40 41 42 43 44 45 -1740.3175 -1240.3175 1459.6825 -2951.4286 2048.5714 -3451.4286 -273.6508 - 1773.6508 -1473. 46 47 48 49 50 51 52 53 54 -640.3175 -140.3175 -140.3175 -443.1746 -343.1746 -143.1746 -4554.2857 - 2054.2857 -1854. 55 56 57 58 59 60 61 62 63 -1876.5079 -1876.5079 -1376.5079 -443.1746 -443.1746 -543.1746 953. 1053.9683 1153. 64 65 66 67 68 69 70 71 72 -3657.1429 -3157.1429 342.8571 -679.3651 -479.3651 -579.3651 653. 753.9683 553. fitted(anova) #não precisa ser feito 1 2 3 4 5 6 7 8 9 10

boxplot(densidade~ponto+estacao, data=bifatorial)

plot(anova)

shapiro.test(res) #do residuo Shapiro-Wilk normality test data: res W = 0.92299, p-value = 0.