SSCP matrix function
SSCP(A) Sums of squares and cross-products. Returns a square matrix having as many rows (and columns) as the argument has columns. The diagonal values are sums of squares and the off-diagonal values are sums of cross products. SSCP(A) equals A′A (for A′ see TRANS function above).