* Create SAS data set from soil data;
* Using shortcuts to refer to variables;
* Read data file 'soil_data.txt';

DATA soilsamp;
*INFILE 'z:\My Documents\teaching\stat_517\soil_data.txt';
FILENAME webpage URL 'http://people.stat.sc.edu/hitchcock/soil_data.txt';
INFILE webpage;
INPUT Region :$9. Temp Iron1 Iron2 Iron3 Iron4;
ARRAY ironct (4) Iron1-Iron4;
ARRAY logironct (4);
DO i = 1 to 4;
 IF ironct(i) = 0 THEN ironct(i) = 0.01;
    logironct(i) = log(ironct(i));
END;

AvgIron = MEAN(OF Iron1-Iron4);
AvgIronAlt = MEAN(OF Iron:);  /* Another shortcut approach */
MaxAllNumbers = MAX(OF _NUMERIC_);

PROC PRINT DATA = soilsamp;
title 'Some Statistics for the Soil Samples';
RUN;