Scribd
Upload
33 views12 pages

SAS Ejemplos

1. This document contains examples and code for analyzing experimental data where the assumptions of ANOVA are and are not met. It includes tests for normality and homogeneity of variance, and shows how to perform ANOVA, Tukey's test, and Kruskal-Wallis test for non-normal data. 2. The first example analyzes data that meets ANOVA assumptions and performs ANOVA, Tukey's test, and LSD test. The second example analyzes non-normal data using Kruskal-Wallis test instead of ANOVA. 3. The last section shows the code to perform a Kruskal-Wallis test on diameter data for different treatments of the fungus Lasiodiplodia citric

Uploaded by

Germán Morales
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
33 views12 pages

SAS Ejemplos

1. This document contains examples and code for analyzing experimental data where the assumptions of ANOVA are and are not met. It includes tests for normality and homogeneity of variance, and shows how to perform ANOVA, Tukey's test, and Kruskal-Wallis test for non-normal data. 2. The first example analyzes data that meets ANOVA assumptions and performs ANOVA, Tukey's test, and LSD test. The second example analyzes non-normal data using Kruskal-Wallis test instead of ANOVA. 3. The last section shows the code to perform a Kruskal-Wallis test on diameter data for different treatments of the fungus Lasiodiplodia citric

Uploaded by

Germán Morales
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 12

1.

DCA
Ejp 1. Se cumplen los supuestos del ANOVA /
TUKEY Y LSD
DATA Variable;/*Ejercicio 1*/
INPUT Trat cuartatt cuartaT;
DATALINES;
1 1.41 1.73
1 1.41 1.73
1 1.41 1.73
1 2.83 3.00
2 2.24 2.45
2 1.41 1.73
2 2.45 2.65
2 3.00 3.16
3 2.83 3.00
3 2.83 3.00
3 2.45 2.65
3 3.00 3.16
4 2.45 2.65
4 1.41 1.73
4 2.83 3.00
4 2.65 2.83
5 0.00 1.00
5 0.00 1.00
5 0.00 1.00
5 0.00 1.00
;
ODS SELECT TestsForNormality;
proc univariate data=Variable normaltest;
var cuartatt cuartaT;
by Trat;
run;
proc print;
PROC GLM;
CLASS Trat;
MODEL cuartatt cuartaT=Trat;
MEANS Trat/HOVTEST=BARTLETT;
MEANS Trat/HOVTEST=LEVENE;
OUTPUT OUT=DOS R=ERROR;
PROC UNIVARIATE DATA=DOS NORMAL;
VAR cuartatt cuartaT;
RUN;
1. DCA
COMPARACIÓN DE MEDIAS POR TUKEY
DATA Variable;/*Ejp 1 cumplimiento de supuestos del ANOVA*/
INPUT Trat Rep PreEval Primera;
DATALINES;
1 1 4 4
1 2 5 1
1 3 7 2
1 4 2 1
2 1 4 4
2 2 4 4
2 3 6 4
2 4 2 2
3 1 4 4
3 2 3 4
3 3 5 7
3 4 1 4
4 1 1 1
4 2 3 5
4 3 1 3
4 4 3 2
5 1 4 0
5 2 5 0
5 3 3 0
5 4 3 0
;
PROC GLM;
CLASS TRAT;
MODEL PreEval Primera=TRAT;
MEANS TRAT/tukey;
RUN;

COMPARACIÓN DE MEDIAS POR LSD


DATA Variable;/*Ejp 1 cumplimiento de supuestos del ANOVA*/
INPUT Trat Rep PreEval Primera;
DATALINES;
1 1 4 4
1 2 5 1
1 3 7 2
1 4 2 1
2 1 4 4
2 2 4 4
2 3 6 4
2 4 2 2
3 1 4 4
3 2 3 4
3 3 5 7
3 4 1 4
4 1 1 1
4 2 3 5
4 3 1 3
4 4 3 2
5 1 4 0
5 2 5 0
5 3 3 0
5 4 3 0
;
PROC GLM;
CLASS TRAT;
MODEL PreEval Primera=TRAT;
MEANS TRAT/LSD;
RUN;
1. DCA
Ejp 2 No cumplimiento de supuestos de ANOVA /
KRUSKAL WALLIS
DATA Variable;/*Diámetro de colonias Lasiodiplodia citricola_DCA_No cumple
supuestos*/
INPUT Trat Rep NTransf Transf;
DATALINES;
1 1 0.00 0.000
1 2 0.00 0.000
1 3 0.00 0.000
1 4 0.00 0.000
1 5 0.00 0.000
2 1 0.00 0.000
2 2 0.00 0.000
2 3 0.00 0.000
2 4 0.00 0.000
2 5 0.00 0.000
3 1 0.00 0.000
3 2 0.00 0.000
3 3 0.00 0.000
3 4 0.00 0.000
3 5 0.00 0.000
4 1 0.00 0.000
4 2 0.00 0.000
4 3 0.00 0.000
4 4 0.00 0.000
4 5 0.00 0.000
5 1 0.00 0.000
5 2 0.00 0.000
5 3 0.00 0.000
5 4 0.00 0.000
5 5 0.00 0.000
6 1 0.00 0.000
6 2 0.00 0.000
6 3 0.00 0.000
6 4 0.00 0.000
6 5 0.00 0.000
7 1 0.00 0.000
7 2 0.00 0.000
7 3 0.00 0.000
7 4 0.00 0.000
7 5 0.00 0.000
8 1 0.00 0.000
8 2 0.00 0.000
8 3 0.00 0.000
8 4 0.00 0.000
8 5 0.00 0.000
9 1 0.00 0.000
9 2 0.00 0.000
9 3 0.00 0.000
9 4 0.00 0.000
9 5 0.00 0.000
10 1 0.00 0.000
10 2 0.00 0.000
10 3 0.00 0.000
10 4 0.00 0.000
10 5 0.00 0.000
11 1 0.00 0.000
11 2 0.00 0.000
11 3 0.00 0.000
11 4 0.00 0.000
1. DCA
11 5 0.00 0.000
12 1 0.00 0.000
12 2 0.00 0.000
12 3 0.00 0.000
12 4 0.00 0.000
12 5 0.00 0.000
13 1 0.00 0.000
13 2 0.00 0.000
13 3 0.00 0.000
13 4 0.00 0.000
13 5 0.00 0.000
14 1 0.00 0.000
14 2 0.00 0.000
14 3 0.00 0.000
14 4 0.00 0.000
14 5 0.00 0.000
15 1 0.00 0.000
15 2 0.00 0.000
15 3 0.00 0.000
15 4 0.00 0.000
15 5 0.00 0.000
16 1 2.00 1.414
16 2 2.00 1.414
16 3 2.10 1.449
16 4 1.60 1.265
16 5 2.00 1.414
17 1 1.90 1.378
17 2 2.00 1.414
17 3 1.70 1.304
17 4 2.20 1.483
17 5 2.00 1.414
18 1 2.00 1.414
18 2 2.10 1.449
18 3 2.10 1.449
18 4 2.00 1.414
18 5 1.80 1.342
19 1 0.00 0.000
19 2 0.00 0.000
19 3 0.00 0.000
19 4 0.00 0.000
19 5 0.00 0.000
20 1 0.00 0.000
20 2 0.00 0.000
20 3 0.00 0.000
20 4 0.00 0.000
20 5 0.00 0.000
21 1 0.00 0.000
21 2 0.00 0.000
21 3 0.00 0.000
21 4 0.00 0.000
21 5 0.00 0.000
22 1 0.00 0.000
22 2 0.00 0.000
22 3 0.00 0.000
22 4 0.00 0.000
22 5 0.00 0.000
23 1 0.00 0.000
23 2 0.00 0.000
23 3 0.00 0.000
23 4 0.00 0.000
23 5 0.00 0.000
24 1 0.00 0.000
1. DCA
24 2 0.00 0.000
24 3 0.00 0.000
24 4 0.00 0.000
24 5 0.00 0.000
25 1 0.00 0.000
25 2 0.00 0.000
25 3 0.00 0.000
25 4 0.00 0.000
25 5 0.00 0.000
26 1 0.00 0.000
26 2 0.00 0.000
26 3 0.00 0.000
26 4 0.00 0.000
26 5 0.00 0.000
27 1 0.00 0.000
27 2 0.00 0.000
27 3 0.00 0.000
27 4 0.00 0.000
27 5 0.00 0.000
28 1 8.20 2.864
28 2 7.80 2.793
28 3 8.00 2.828
28 4 8.00 2.828
28 5 8.00 2.828
;
ODS SELECT TestsForNormality;
proc univariate data=Variable normaltest;
var NTransf Transf;
by Trat;
run;
proc print;
PROC GLM;
CLASS Trat;
MODEL NTransf Transf=Trat;
MEANS Trat/HOVTEST=BARTLETT;
MEANS Trat/HOVTEST=LEVENE;
OUTPUT OUT=DOS R=ERROR;
PROC UNIVARIATE DATA=DOS NORMAL;
VAR NTransf Transf;
RUN;
1. DCA
Prueba de Kruskal WAllis
DATA Variable;/*Diámetro de colonias Lasiodiplodia citricola*/
INPUT TRAT Rep Eval4;
DATALINES;
1 1 0.00
1 2 0.00
1 3 0.00
1 4 0.00
1 5 0.00
2 1 0.00
2 2 0.00
2 3 0.00
2 4 0.00
2 5 0.00
3 1 0.00
3 2 0.00
3 3 0.00
3 4 0.00
3 5 0.00
4 1 0.00
4 2 0.00
4 3 0.00
4 4 0.00
4 5 0.00
5 1 0.00
5 2 0.00
5 3 0.00
5 4 0.00
5 5 0.00
6 1 0.00
6 2 0.00
6 3 0.00
6 4 0.00
6 5 0.00
7 1 0.00
7 2 0.00
7 3 0.00
7 4 0.00
7 5 0.00
8 1 0.00
8 2 0.00
8 3 0.00
8 4 0.00
8 5 0.00
9 1 0.00
9 2 0.00
9 3 0.00
9 4 0.00
9 5 0.00
10 1 0.00
10 2 0.00
10 3 0.00
10 4 0.00
10 5 0.00
11 1 0.00
11 2 0.00
11 3 0.00
11 4 0.00
11 5 0.00
12 1 0.00
12 2 0.00
12 3 0.00
1. DCA
12 4 0.00
12 5 0.00
13 1 0.00
13 2 0.00
13 3 0.00
13 4 0.00
13 5 0.00
14 1 0.00
14 2 0.00
14 3 0.00
14 4 0.00
14 5 0.00
15 1 0.00
15 2 0.00
15 3 0.00
15 4 0.00
15 5 0.00
16 1 2.00
16 2 2.00
16 3 2.10
16 4 1.60
16 5 2.00
17 1 1.90
17 2 2.00
17 3 1.70
17 4 2.20
17 5 2.00
18 1 2.00
18 2 2.10
18 3 2.10
18 4 2.00
18 5 1.80
19 1 0.00
19 2 0.00
19 3 0.00
19 4 0.00
19 5 0.00
20 1 0.00
20 2 0.00
20 3 0.00
20 4 0.00
20 5 0.00
21 1 0.00
21 2 0.00
21 3 0.00
21 4 0.00
21 5 0.00
22 1 0.00
22 2 0.00
22 3 0.00
22 4 0.00
22 5 0.00
23 1 0.00
23 2 0.00
23 3 0.00
23 4 0.00
23 5 0.00
24 1 0.00
24 2 0.00
24 3 0.00
24 4 0.00
24 5 0.00
1. DCA
25 1 0.00
25 2 0.00
25 3 0.00
25 4 0.00
25 5 0.00
26 1 0.00
26 2 0.00
26 3 0.00
26 4 0.00
26 5 0.00
27 1 0.00
27 2 0.00
27 3 0.00
27 4 0.00
27 5 0.00
28 1 8.20
28 2 7.80
28 3 8.00
28 4 8.00
28 5 8.00
;
PROC RANK TIES=MEAN;
VAR Eval4;
RANKs RANK1;
PROC PRINT;
PROC GLM;
CLASS TRAT;
MODEL RANK1=TRAT;
MEANS TRAT/tukey;
RUN;
PROC PRINT;
PROC FREQ DATA=Variable;
TABLE TRAT*Eval4/CMH2 SCORES=RANK NOPRINT;
RUN;
2. DBCA
Ejemplo 3_ DBCA_ Transformación efectiva /
Prueba de TUKEY Y LSD
DATA Variable;/*Ejp 3_Severidad antracnosis_ DBCA_No se cumple
supuestos_Tukey*/
INPUT Trat Bloq$ NTransf Transf;
DATALINES;
1 1 3.33 1.826
1 2 8.33 2.887
1 3 6.67 2.582
1 4 3.33 1.826
2 1 3.33 1.826
2 2 5.00 2.236
2 3 5.00 2.236
2 4 1.67 1.291
3 1 6.67 2.582
3 2 3.33 1.826
3 3 1.67 1.291
3 4 3.33 1.826
4 1 3.33 1.826
4 2 6.67 2.582
4 3 5.00 2.236
4 4 1.67 1.291
5 1 10.00 3.162
5 2 13.33 3.651
5 3 10.00 3.162
5 4 11.67 3.416
;
ODS SELECT TestsForNormality;
proc univariate data=Variable normaltest;
var NTransf Transf;
by Trat;
run;
proc print;
PROC GLM;
CLASS Trat;
MODEL NTransf Transf=Trat;
MEANS Trat/HOVTEST=BARTLETT;
MEANS Trat/HOVTEST=LEVENE;
OUTPUT OUT=DOS R=ERROR;
PROC UNIVARIATE DATA=DOS NORMAL;
VAR NTransf Transf;
RUN;

Comparación de medias por Tukey

DATA Variable;/*Ejp 3_Severidad antracnosis_ DBCA_No se cumple


supuestos_Tukey*/
INPUT Bloq Trat NTransf Transf;
DATALINES;
1 1 3.33 1.83
1 2 3.33 1.83
1 3 6.67 2.58
1 4 3.33 1.83
1 5 10.00 3.16
2 1 8.33 2.89
2 2 5.00 2.24
2 3 3.33 1.83
2 4 6.67 2.58
2. DBCA
2 5 13.33 3.65
3 1 6.67 2.58
3 2 5.00 2.24
3 3 1.67 1.29
3 4 5.00 2.24
3 5 10.00 3.16
4 1 3.33 1.83
4 2 1.67 1.29
4 3 3.33 1.83
4 4 1.67 1.29
4 5 11.67 3.42
;
PROC GLM;
CLASS Bloq TRAT;
CLASS TRAT;
MODEL NTransf Transf=Bloq TRAT;
MEANS TRAT/tukey;
RUN;

Ejemplo 4_ DBCA_ Transformación no efectiva/


Prueba de Friedman
DATA Variable;/*Ejp 4_DBCA_No se cumple supuestos_Friedman*/
INPUT Trat Bloq$ NTransf Transf;
DATALINES;
1 I 0.2 1.095
1 II 0.2 1.095
1 III 0.1 1.049
1 IV 0 1.000
2 I 0 1.000
2 II 0 1.000
2 III 0 1.000
2 IV 0 1.000
3 I 0 1.000
3 II 0 1.000
3 III 0 1.000
3 IV 0 1.000
4 I 0 1.000
4 II 0 1.000
4 III 0 1.000
4 IV 0 1.000
5 I 0.8 1.342
5 II 1.2 1.483
5 III 1.2 1.483
5 IV 1.3 1.517
;
ODS SELECT TestsForNormality;
proc univariate data=Variable normaltest;
var NTransf Transf;
by Trat;
run;
proc print;
PROC GLM;
CLASS Trat;
MODEL NTransf Transf=Trat;
MEANS Trat/HOVTEST=BARTLETT;
2. DBCA
MEANS Trat/HOVTEST=LEVENE;
OUTPUT OUT=DOS R=ERROR;
PROC UNIVARIATE DATA=DOS NORMAL;
VAR NTransf Transf;
RUN;
2. DBCA
Prueba de Friedman
DATA Nro_Larvas;/*FITO-35 Insecticida Famectin vs Spodoptera
frugiperda en maiz*/
INPUT Bloq$ Trat Segunda Tercera;
DATALINES;
I 1 0.2 0.2
I 2 0 0.1
I 3 0 0
I 4 0 0
I 5 0.8 1.4
II 1 0.2 0.2
II 2 0 0
II 3 0 0.1
II 4 0 0.1
II 5 1.2 1.4
III 1 0.1 0.2
III 2 0 0.1
III 3 0 0.1
III 4 0 0.1
III 5 1.2 1.3
IV 1 0 0.1
IV 2 0 0.1
IV 3 0 0
IV 4 0 0.1
IV 5 1.3 1.5
;
PROC RANK TIES=MEAN;
VAR Segunda Tercera;
RANKs RANKEv2 RANKEv3;
BY Bloq;
PROC PRINT;
PROC GLM;
CLASS Bloq TRAT;
CLASS TRAT;
MODEL RANKEv2 RANKEv3= TRAT Bloq;
LSMEANS TRAT/TDIFF ADJUST=TUKEY;
RUN;
PROC FREQ DATA=Nro_Larvas;
TABLE Bloq*TRAT*Segunda/CMH2 SCORES=RANK NOPRINT;
TABLE Bloq*TRAT*Tercera/CMH2 SCORES=RANK NOPRINT;
RUN;

You might also like