sas_anova

DATA A;
INPUT CLONE	$ 1-4 PARCELA	VOL;
DATALINES;
H13	 1	256.91
H13	 2	211.92
H13	 3	244.17
H13	 4	228.54
H13	 5	201.30
H13	 6	228.33
C206 1	228.81
C206 2	231.91
C206 3	159.34
C206 4	174.82
C206 5	210.49
C206 6	220.50
HG11 1	206.36
HG11 2	178.68
HG11 3	197.69
HG11 4	279.63
HG11 5	206.17
HG11 6	209.80
AC15 1	201.72
AC15 2	228.92
AC15 3	211.29
AC15 4	209.54
AC15 5	234.19
AC15 6	166.38
C105 1	219.71
C105 2	240.36
C105 3	235.48
C105 4	221.17
C105 5	184.62
C105 6	261.64
;
/*proc print data=a; run;*/

data A;
set A;
if CLONE= 'H13'  then BM = VOL * 0.48;
if CLONE= "C206" then BM = VOL * 0.52;
if CLONE= "HG11" then BM = VOL * 0.41;
if CLONE= "AC15" then BM = VOL * 0.59;
if CLONE= "C105" then BM = VOL * 0.5;
run;
proc print data=A; run;

/*Testar normalidade */
PROC UNIVARIATE NORMAL DATA=A;
VAR BM;
RUN;

/*Testar homegeneidade de variancias*/
PROC TRANSREG DATA=A PLOTS=NONE;
MODEL BOXCOX(BM)=CLASS(CLONE);
RUN;

/*ANOVA*/
PROC GLM DATA=A ;
CLASS CLONE;
model BM=CLONE;
MEANS CLONE / TUKEY lines;
run;

/*Grafico IC*/
/*proc mixed data=A method=type3;*/
/*class CLONE;*/
/*model BM = CLONE;*/
/*store out1;  */
/*run;         */
/*ods graphics on;*/
/*ods html style=statistical sge=on; */
/*proc plm restore=out1;             */
/*lsmeans CLONE/ adjust=tukey plot=meanplot cl lines; */
/*ods exclude diffs diffplot;       */
/*run; title; run;*/