/* Tables Relating the Price of Wheat and Other Farm Products Since 1890.
Katherine Bennet Davis. Journal of Political Economy. Vol. 6 #3. June 1898 
All prices are "per ct. 100"*/

options nodate nonumber ps=54 ls=76;

data one;
length city $8.;
do city='CHICAGO', 'NEW YORK';
do year=1890 to 1897;
do commodity='Wheat', 'Corn', 'Oats', 'Beef', 'Pork', 'Lard';
input prcper100 @@;
output;
end; end; end;
cards;
84.8 78.7 88.3 56.7 63.6 60.6
91.9 117.0 111.8 64.6 58.0 63.2
74.9 89.1 87.7 53.5 67.1 74.3
63.8 76.7 82.0 69.0 100.4 98.8
54.2 87.1 90.3 60.6 76.5 69.7
59.1 80.5 69.9 67.3 62.8 62.5
60.7 52.0 51.7 60.3 42.7 43.7
76.7 51.3 52.4 65.9 47.8 41.6
74.2 65.4 75.9 63.3 70.0 59.7
82.6 95.7 92.4 75.9 65.7 62.3
68.5 73.3 72.9 62.4 66.4 72.5
55.8 67.7 71.9 74.3 105.3 97.3
46.1 69.1 74.7 74.2 82.5 73.1
51.3 64.8 57.8 73.6 69.5 61.3
58.9 46.1 48.6 68.3 51.6 44.0
72.0 43.3 48.4 70.1 51.0 41.6
;


proc print;
run;

proc glm;
class city year commodity;
model prcper100 = city year commodity city*year city*commodity year*commodity;
run;

data twoc;
set one;
if city='CHICAGO';
prcper100c=prcper100;
drop prcper100 city;
run;

data twon;
set one;
if city='NEW YORK';
prcper100n=prcper100;
drop prcper100 city;
run;

data two; 
merge twoc twon;
run;

proc sort;
by commodity year;
run;

symbol1 c=black l=1 v=: i=join;
symbol2 c=black l=2 v=: i=join;

proc gplot;
by commodity;
plot prcper100c*year=1 prcper100n*year=2 / overlay frame;
run;

quit;




quit;