" Exercise 9.2

 Statistical Methods in Biology: Design and Analysis of Experiments and Regression
 by S.J. Welham, S.A. Gezan, S.J. Clark & A. Mead (2014)
 Chapman & Hall/CRC Press, Boca Raton, Florida. ISBN: 978-1-4398-0878-8
 Data from E. Wright, Rothamsted Research
 Version 1, 11/07/2014 "

" Set working directory - change to location of your data file "
\SET [WORKINGDIRECTORY='d:/stats4biol/data/']

" Read data from working directory "
FILEREAD [name='SIZE.DAT'; IMETHOD=read] FGROUPS=no,4(yes),no

" ANOVA of colony sizes "
BLOCKSTRUCTURE (Occasion*Incubator)/Dish
TREATMENTSTRUCTURE Temperature
ANOVA [FPROBABILITY=yes] Y=Size

" Residual plots "
APLOT [RMETHOD=standardized] METHOD=fittedvalues,normal,histogram,absresidual

" Get log10 transformation of colony sizes "
CALCULATE log10Size = LOG10(Size)

" ANOVA of logged colony sizes "
BLOCKSTRUCTURE (Occasion*Incubator)/Dish
TREATMENTSTRUCTURE Temperature
ANOVA [FPROBABILITY=yes; PSE=differences,lsd] Y=log10Size

" Residual plots "
APLOT [RMETHOD=standardized] METHOD=fittedvalues,normal,histogram,absresidual

" Treatment means with LSD "
AGRAPH [METHOD=data; PSE=LSD] XFACTOR=Temperature; YTITLE='Log~_{10}(Colony size)'

" Extract treatment means and print on back-transformed scale "
AKEEP TERMS=Temperature; MEANS=means
PRINT STRUCTURE=(10**means)

" End of File "