* Create SAS data set from soil data;
* Creating arrays to help work with many 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 Iron2 Iron3 Iron4;
ARRAY logironct (4);
DO i = 1 to 4;
 IF ironct(i) = 0 THEN ironct(i) = 0.01;
    logironct(i) = log(ironct(i));
END;


PROC PRINT DATA = soilsamp;
title 'Iron Readings for Various Soil Samples';
RUN;