data mixed; 

do row = 1 to 24; do col=1 to 24; 
d=abs(row-col); 
ar1 = .8**d; 
ar1_ran = .6*ar1 + .4; 
CS = 1 - .8*(d>0); 
toep = (.95**d)*cos(d/2); 
banded = toep*(d<4); 
vc = (d=0); 
output; end; end; 

proc g3d;  scatter row*col=ar1/shape ="prism" zmin=0; title h=1.5 "AR(1)"; 
proc g3d; scatter row*col=ar1_ran/shape ="prism" zmin=0; title h=1.5 "Repeated AR(1)& Random Subject"; 
proc g3d; scatter row*col= toep/shape ="prism" zmin=-1; title h=1.5 "Toeplitz"; 
proc g3d; scatter row*col=banded/shape="prism" zmin=0; title h=1.5 "Banded Toeplitz"; 
proc g3d; scatter row*col=VC/shape ="prism" zmin=0; title h=1.5 "Variance Components"; 
proc g3d; scatter row*col=CS/shape ="prism" zmin=0; 
title h=1.5 "Compound Symmetry"; title2 "or Repeated Variance Components & Random Subject";
run;