GET
  FILE='C:\rs.sav'.

SAVE OUTFILE='C:\rs1.sav'
 /COMPRESSED.

set mxloop = 999999.
cluster1  y  FILE = 'C:\body.sav'.



define cluster1(!positional !TOKENS(1) /FILE=!CMDEND).
GET   FILE=!File.
QUICK CLUSTER
  	!1
  	/MISSING=LISTWISE
  	/CRITERIA= CLUSTER(2) MXITER(20) CONVERGE(0)
 	 /METHOD=KMEANS(UPDATE )
  	/SAVE CLUSTER
  	/PRINT NONE.



matrix.
get s /variable = qcl_1 /file = *.
compute s = s -1.
compute nr = nrow(s).
compute j = make(nr,1,1).

compute splus = t(j)*s.
print nr.
print splus.
compute splus = splus/nr.

get rs /variables = all /file ="c:\rs1.sav".
compute rs = {rs;splus}.
print rs.

save rs/ outfile ="c:\rs1.sav".
end matrix.


DELETE VARIABLES  qcl_1.


!do  !i = 1 !to  6.

GET   FILE=!File.

matrix.

get x /variables  = !1 /file=*.
compute nr = nrow(x).
compute ystar = make(nr,1,-1).
loop j = 1 to nr.
compute id1 = trunc(nr*uniform(1,1)+1).
compute ystar(j,1) = x(id1,1).
end loop.
save ystar /outfile = *.
end matrix. 

	QUICK CLUSTER
  	col1
  	/MISSING=LISTWISE
  	/CRITERIA= CLUSTER(2) MXITER(20) CONVERGE(0)
 	 /METHOD=KMEANS(UPDATE )
  	/SAVE CLUSTER
  	/PRINT NONE.



matrix.
get s /variable = qcl_1 /file = *.
compute s = s -1.
compute nr = nrow(s).
compute j = make(nr,1,1).

compute splus = t(j)*s.
print nr.
print splus.
compute splus = splus/nr.
get rs /variables = all/file ="c:\rs1.sav".
compute rs = {rs;splus}.
print rs.

save rs/ outfile ="c:\rs1.sav".
end matrix.


DELETE VARIABLES  qcl_1.

!doend.

!enddefine.