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

INPUT PROGRAM .
compute #n = 16.
LOOP SAMP=1 to 2000.
LOOP V = 1 to #n.
COMPUTE x=rv.beta(1,5).
END CASE.
LEAVE SAMP.
END LOOP.
END LOOP.
END FILE.
END INPUT PROGRAM .

AGGREGATE
  /OUTFILE=*
  /BREAK=samp
  /x_1 = MEAN(x) /x_2 = SD(x) /N1=N.


COMPUTE variance = x_2**2.
EXECUTE .




  
IGRAPH /VIEWNAME='Scatterplot' /X1 = VAR(x_1) TYPE = SCALE /Y = VAR(x_2) TYPE
  = SCALE /COORDINATE = VERTICAL  /FITLINE METHOD = REGRESSION LINEAR LINE =
  TOTAL SPIKE=OFF /X1LENGTH=3.0 /YLENGTH=3.0 /X2LENGTH=3.0 /CHARTLOOK='NONE'
 /SCATTER COINCIDENT = NONE.
EXE.
COMPUTE TRUEM = 1/6.
COMPUTE T = (X_1 - TRUEM)*SQRT(N1)/X_2.
COMPUTE C = 0.
IF (ABS(T) > CDF.T(.975,N1-1))C = 1 .
EXECUTE .
DESCRIPTIVES
  VARIABLES=c X_1 TRUEM
  /STATISTICS=MEAN .