***Instructions: Change "var1 var2...varx" in line 8 to your desired variables.***





OMS 


 /select tables


 /if commands = ['Correlations'] SUBTYPES = ['Correlations']


 /destination format = sav OUTFILE = 'alpha1.sav' .


 CORRELATIONS 


	/VARIABLES= var1 var2...varx


 /STATISTICS XPROD


 /MISSING=PAIRWISE .


OMSEND .


GET FILE "alpha1.sav".


SELECT IF (Var2='Covariance').


SAVE OUTFILE 'alpha2.sav'


  /DROP = Command_ Subtype_ Label_ Var1 Var2 .


GET FILE 'alpha2.sav'.


MATRIX .


GET X /FILE=* /VARIABLES ALL.


COMPUTE k = nrow(x) .


COMPUTE SUMVAR=TRACE(X).


COMPUTE SUMCOV=(MSUM(X)-SUMVAR)/2.


COMPUTE VBAR = SUMVAR/k .


COMPUTE CBAR = SUMCOV/((k*k-k)/2) .


PRINT VBAR.


PRINT CBAR.


COMPUTE PAIRALPH = (k*CBAR)/(VBAR+(k-1)*CBAR) .


PRINT pairalph .


END MATRIX .