DATA a;

INFILE 'qfdata.dat';

INPUT sex year pop dac dlungca dbrca  agefrom ageto;

PROC MEANS;

RUN;

/*

* data manipulation / analysis statements;

if agefrom >= 50 and ageto <= 85;

age = (agefrom-45)/5;

py=pop*1; pm=pop*12;

arate= dac / py ; lnarate=log(arate+0.00000001);
brate=dbrca/py ;  lnbrate=log(brate+0.00000001);
lrate=dlungca/py ;lnlrate=log(lrate+0.00000001);

proc glm;
 class age;
 model lnbrate  = age year / solution;
weight dbrca;

proc logistic ; model dac / py  = age year;
proc logistic ; model dac / pm  = age year;

proc glm; model lnarate  = age year; weight dac;

proc glm; model lnlrate  = age year; weight dlungca;


proc plot;
plot lrate*age =   year;

*/