compute p1 =p.
execute.


SORT CASES BY
  p (a) adjrsq (d) .

compute filter1 = 0.
execute.
Matrix.
Get x /file = * /variables = all /names = names.
* print x.
compute nr = nrow(x).
compute nc = ncol(x).
compute x(1,nc) = 1.
compute x(2,nc) = 1.
compute x(nr,nc)=1.
loop i = 3 to nr-1.
do if (x(i,nc-1) <>x(i-1,nc-1)).
compute x(i,nc) = 1.
compute x(i+1,nc) = 1.

end if.
end loop.
save x /outfile = */names = names.
end matrix.


FILTER OFF.
USE ALL.
SELECT IF(filter1=1).
EXECUTE .




GRAPH
  /SCATTERPLOT(OVERLAY)=p p  WITH p1 mallows (PAIR) BY model (IDENTIFY)
  /MISSING=LISTWISE
  /TITLE= "Mallow's Cp Versus P".

GRAPH
  /SCATTERPLOT(BIVAR)=p WITH adjrsq BY model
  /MISSING=LISTWISE
  /TITLE= 'Adjusted Rsq Versus P'.
GRAPH
  /SCATTERPLOT(BIVAR)=p WITH mse BY model
  /MISSING=LISTWISE
  /TITLE= 'MSE Versus P'.