* SAS code for C-A trend statistic;
 options linesize=65;

 data tBHQ;
  input conc nlabel @@;
  total = 2000;
  notlabel = total - nlabel;
  pi = nlabel/total;

  label="Yes"; count=nlabel; output;
  label=" No"; count=notlabel; output;

  cards;
 0 8 120 9 240 9 480 18 600 39 720 44

 proc freq order=data;
  table conc*label / chisq;
  weight count;

 run;