* this assumes you  have a data set with one y.

INPUT PROGRAM .
compute #n = 16.
LOOP SAMP=1 to 3000.
LOOP V = 1 to #n.
COMPUTE x=rv.chi(1).
END CASE.
LEAVE SAMP.
END LOOP.
END LOOP.
END FILE.
END INPUT PROGRAM .

AGGREGATE
  /OUTFILE=*
  /BREAK=samp
  /mean = MEAN(x) /std = SD(x) /N1=N.


CORRELATIONS
  /VARIABLES=std mean
  /PRINT=TWOTAIL NOSIG
  /MISSING=PAIRWISE .