* PropVal4 Example, source code file crossval.sas;

* Compute the crossvalidation R-squared value, which is
  R2 - 1 - SSE/SST;

options linesize=70 nodate pageno=1;

data crossval;
   infile "c:/datasets/actual-predicted-mean.txt" firstobs=2;
   input y yhat ymean;
   title "CrossVal Output";
   ressq = (y - yhat)**2;
   devsq = (y - ymean)**2;

proc print;

proc means sum;
   var ressq devsq;

run;
quit;