* SAS code for adjusted CA trend statistic;
 options linesize=65;

 data vcdp;
   input conc x nrisk nriskBP @@;
   nsurv = nrisk - x;
   nsurvBP = nriskBP - x;
   cards;
    0  4 50 41.8370   25  9 50 40.9661
   50 11 50 38.4904  100  7 50 22.8135

 data BPtrend; set vcdp;
   tumor="yes"; count=x; output;
   tumor=" no"; count=nsurvBP; output;
   keep conc tumor count;

 proc freq order=data;
  table conc*tumor / chisq nocol nopercent;
  weight count;

 run;