#####################################################################           
#                                                                   #           
#                         P1K2DRVR MINITABC                         #           
#                         -----------------                         #           
#                                                                   #           
#      CLUSTER OPTIMIZATION FOR P=1, K=2, WGSS (BGSS) CRITERION     #           
#                  Version 1.4       26-Jan-94                      #           
#                                                                   #           
#                                                                   #           
# INPUT:  ordered data, in a dataset  P1K2CLUS DAT                  #           
#                                                                   #           
#                                                                   #           
#####################################################################           
#                                                                               
#                                                                               
#                                                                               
      OUTFILE 'P1K2CLUS MTBOUTP A'                                              
      READ 'P1K2CLUS DAT' C1                                                    
      PARSUMS C1, PUT INTO C11                                                  
      LET C12 = SUM(C1)-C11                                                     
# Note that the sample size is N(C1).                                           
      let k1 = n(C1)-1                                                          
      COPY C11 TO C11;                                                          
        USE ROWS 1:k1.                                                          
      COPY C12 TO C12;                                                          
         USE ROWS 1:k1.                                                         
      PRINT C11 C12                                                             
      SET C21                                                                   
      1:k1                                                                      
      LET C22 = n(c1) - C21                                                     
      PRINT C22                                                                 
      LET C31 = C11**2/C21                                                      
      LET C32 = C31       +  C12**2/C22                                         
      PRINT C21 C32                                                             
      HEIGHT 40                                                                 
      name c32 'BGSS'                                                           
      let c32  =  c32  -  ( SUM(c1)**2 )/N(C1)                                  
      name c33 'WGSS'                                                           
# Compute SSTot:                                                                
      let k2  =  ( STAN(c1)**2 )*k1                                             
      let c33  =  k2 - c32                                                      
      width = 70                                                                
      PLOT C33 C21;                                                             
         xincr = 2.                                                             
      PRINT C1 'WGSS'