compute cov matrix
data calc ng=1
matrices
r sy 16 16
s di 16 16
compute s*r*s' /
matrix r file=mtmm.cor
matrix s file=mtmm.sd
option mx%e=mtmm.cov
end