options nodate nonumber ps=54 ls=76;

data one;
infile 'H:\public_html\data\atltime1004.dat';
input city $ 22-24 flttime 42-45 distance 49-52;
if flttime > 0;
if city ne 'HNL';
if city ne 'SJU';
if city ne 'STT';
run;

proc sort;
by city;
run;
/*
proc gchart;
by city;
vbar flttime / space=0;
run;
quit;
*/

data two;
set one;
by city;
retain sumtime sumtime2 ntime;
if first.city then do;
sumtime=0;
sumtime2=0;
ntime=0;
end;
sumtime=sumtime+flttime;
sumtime2=sumtime2+(flttime**2);
ntime=ntime+1;
if last.city then do;
meantime=sumtime/ntime;
stdtime=sqrt(((sumtime2-(ntime*meantime**2))/(ntime-1)));
vartime=stdtime**2;
mmbeta=(stdtime**2)/meantime;
mmalpha=(meantime/stdtime)**2;
output;
end;
drop flttime;
run;

proc reg;
model vartime=meantime / noint;
model vartime=meantime;
run;
quit;

symbol1 c=black v=: i=rl;

proc gplot;
plot vartime*meantime;
run;
quit;

proc gplot;
plot meantime*distance;
run;
quit;

proc sort;
by distance;

proc print;
run;

quit;


proc reg data=one;
model flttime = distance;
output out=atlout p=yhat r=e;
run;

symbol1 c=black v=:;

proc gplot;
plot e*yhat / vref=0;
run;

proc gchart;
vbar e / space=0 midpoints=-40 to 40 by 2.5;
run;

proc genmod data=one;
model flttime = distance / link=power(-1) dist=gamma;
run;

quit;