Programa em SAS para Aplicar ANOVA e Teste de Comparações Múltiplas de Tukey
data imc_dat;
input cat $ imc corr kcal;
cards;
AT 20.2 60.7 3200
AT 21.3 54.8 3100
AT 19.3 49.6 2800
AT 21.1 52.3 3300
SEM 22.4 14.9 2600
SEM 21.9 17.8 2700
SEM 23.8 18.6 3200
SEM 24.1 15.1 3300
SE 27.3 2.5 2700
SE 23.4 4.3 2300
SE 25.2 2.3 2600
SE 26.4 2.6 3200
PR 26.2 4.1 2600
PR 24.2 2.1 2700
PR 25.4 1.9 2650
;
proc print;
run;
proc anova;
class cat;
model imc corr kcal = cat;
means cat / tukey lines;
run;
Resultado do Programa de MANOVA:
Arquivo para Download (mht abrir no Wodr)
Arquivo Word de Resultados
- Análise Multivariada Exemplos:
var imc corr kcal; ==> Kruskal Wallis
Arquivo Word para Download:
run;
Resultado do Programa de MANOVA:
Arquivo para Download (mht abrir no Wodr)
Arquivo Word de Resultados
Obs
|
cat
|
imc
|
corr
|
kcal
|
1
|
AT
|
20.2
|
60.7
|
3200
|
2
|
AT
|
21.3
|
54.8
|
3100
|
3
|
AT
|
19.3
|
49.6
|
2800
|
4
|
AT
|
21.1
|
52.3
|
3300
|
5
|
SEM
|
22.4
|
14.9
|
2600
|
6
|
SEM
|
21.9
|
17.8
|
2700
|
7
|
SEM
|
23.8
|
18.6
|
3200
|
8
|
SEM
|
24.1
|
15.1
|
3300
|
9
|
SE
|
27.3
|
2.5
|
2700
|
10
|
SE
|
23.4
|
4.3
|
2300
|
11
|
SE
|
25.2
|
2.3
|
2600
|
12
|
SE
|
26.4
|
2.6
|
3200
|
13
|
PR
|
26.2
|
4.1
|
2600
|
14
|
PR
|
24.2
|
2.1
|
2700
|
15
|
PR
|
25.4
|
1.9
|
2650
|
The ANOVA Procedure
Class
Level Information
|
||
Class
|
Levels
|
Values
|
cat
|
4
|
AT PR SE SEM
|
Number of Observations Read
|
15
|
Number of Observations Used
|
15
|
The ANOVA Procedure
Dependent Variable: imc
Source
|
DF
|
Sum
of Squares
|
Mean
Square
|
F
Value
|
Pr > F
|
Model
|
3
|
63.99233333
|
21.33077778
|
14.23
|
0.0004
|
Error
|
11
|
16.49166667
|
1.49924242
|
||
Corrected Total
|
14
|
80.48400000
|
R-Square
|
Coeff
Var
|
Root
MSE
|
imc Mean
|
0.795094
|
5.214802
|
1.224436
|
23.48000
|
Source
|
DF
|
Anova
SS
|
Mean
Square
|
F
Value
|
Pr > F
|
cat
|
3
|
63.99233333
|
21.33077778
|
14.23
|
0.0004
|
Gráfico Chamado: Box and Wisker Plot
The ANOVA Procedure
Dependent Variable: corr
Source
|
DF
|
Sum
of Squares
|
Mean
Square
|
F
Value
|
Pr > F
|
Model
|
3
|
6829.158500
|
2276.386167
|
300.25
|
<.0001
|
Error
|
11
|
83.397500
|
7.581591
|
||
Corrected Total
|
14
|
6912.556000
|
R-Square
|
Coeff
Var
|
Root
MSE
|
corr Mean
|
0.987935
|
13.60410
|
2.753469
|
20.24000
|
Source
|
DF
|
Anova
SS
|
Mean
Square
|
F
Value
|
Pr > F
|
cat
|
3
|
6829.158500
|
2276.386167
|
300.25
|
<.0001
|
The ANOVA Procedure
Dependent Variable: kcal
Source
|
DF
|
Sum
of Squares
|
Mean
Square
|
F
Value
|
Pr > F
|
Model
|
3
|
497333.333
|
165777.778
|
1.95
|
0.1801
|
Error
|
11
|
935000.000
|
85000.000
|
||
Corrected Total
|
14
|
1432333.333
|
R-Square
|
Coeff
Var
|
Root
MSE
|
kcal Mean
|
0.347219
|
10.18210
|
291.5476
|
2863.333
|
Source
|
DF
|
Anova
SS
|
Mean
Square
|
F
Value
|
Pr > F
|
cat
|
3
|
497333.3333
|
165777.7778
|
1.95
|
0.1801
|
Sabemos que
existem diferenças entre as categorias para Índice de Massa Corporal e para Quilômetros
caminhados ou corridos por semana. Para as Quilocalorias não achamos
diferenças entre as categorias.
The ANOVA Procedure
The ANOVA Procedure
Tukey's
Studentized Range (HSD) Test for imc
Note:
|
This test controls
the Type I experimentwise error rate, but it generally has a higher Type II
error rate than REGWQ.
|
Alpha
|
0.05 = 5%
|
Error Degrees of Freedom
|
11
|
Error Mean Square
|
1.499242
|
Critical Value of Studentized Range
|
4.25603
|
Minimum Significant
Difference
|
2.712
|
Harmonic Mean of Cell Sizes
|
3.692308
|
Note:
|
Cell sizes are not
equal.
|
Means with the same letter
are not significantly different. |
||||
Tukey Grouping
|
Mean
|
N
|
cat
|
|
A
|
25.5750
|
4
|
SE
|
|
A
|
||||
A
|
25.2667
|
3
|
PR
|
|
A
|
||||
B
|
A
|
23.0500
|
4
|
SEM
|
B
|
||||
B
|
20.4750
|
4
|
AT
|
|
The ANOVA Procedure
The ANOVA Procedure
Tukey's
Studentized Range (HSD) Test for corr
Note:
|
This test controls
the Type I experimentwise error rate, but it generally has a higher Type II
error rate than REGWQ.
|
Alpha
|
0.05
|
Error Degrees of Freedom
|
11
|
Error Mean Square
|
7.581591
|
Critical Value of Studentized Range
|
4.25603
|
Minimum Significant
Difference
|
6.0987
|
Harmonic Mean of Cell Sizes
|
3.692308
|
Note:
|
Cell sizes are not
equal.
|
Means with the same letter
are not significantly different. |
|||
Tukey Grouping
|
Mean
|
N
|
cat
|
A
|
54.350
|
4
|
AT
|
B
|
16.600
|
4
|
SEM
|
C
|
2.925
|
4
|
SE
|
C
|
|||
C
|
2.700
|
3
|
PR
|
The ANOVA Procedure
The ANOVA Procedure
Tukey's
Studentized Range (HSD) Test for kcal
Note:
|
This test controls
the Type I experimentwise error rate, but it generally has a higher Type II
error rate than REGWQ.
|
Alpha
|
0.05
|
Error Degrees of Freedom
|
11
|
Error Mean Square
|
85000
|
Critical Value of Studentized Range
|
4.25603
|
Minimum Significant
Difference
|
645.75
|
Harmonic Mean of Cell Sizes
|
3.692308
|
Note:
|
Cell sizes are not
equal.
|
Means with the same letter
are not significantly different. |
|||
Tukey Grouping
|
Mean
|
N
|
cat
|
A
|
3100.0
|
4
|
AT
|
A
|
|||
A
|
2950.0
|
4
|
SEM
|
A
|
|||
A
|
2700.0
|
4
|
SE
|
A
|
|||
A
|
2650.0
|
3
|
PR
|
- Análise Multivariada Exemplos:
 |
Comparar resultados ANOVA com Kruskal-Wallis.
Criar programa de KruskalWallis para os dados da Manova com base nos Slides a Seguir
data imc_dat;
input cat $ imc corr kcal;
cards;
AT 20.2 60.7 3200
AT 21.3 54.8 3100
AT 19.3 49.6 2800
AT 21.1 52.3 3300
SEM 22.4 14.9 2600
SEM 21.9 17.8 2700
SEM 23.8 18.6 3200
SEM 24.1 15.1 3300
SE 27.3 2.5 2700
SE 23.4 4.3 2300
SE 25.2 2.3 2600
SE 26.4 2.6 3200
PR 26.2 4.1 2600
PR 24.2 2.1 2700
PR 25.4 1.9 2650
;
proc print;
run;
proc npar1way data=imc_dat wilcoxon dscf;
class cat;
var imc corr kcal = cat;
run;
var imc corr kcal = cat; ==> ANOVA, GLM, MANOVA
Tem um erro nesse programa, onde esta? Compare com os slides apresentados antes do programa SAS (Dois slides de Kruskal Wallis com trechos de programa SAS).
var imc corr kcal; ==> Kruskal Wallis
Resultados de Kruskal Wallis
Resultados ANOVA dados brutos, transformados e da Estat. Robusta. Tese doutorado do Gabriel.
Colaborar e Orientar Exercícios
Nenhum comentário:
Postar um comentário