noi di _dup(72) "-" _n "cluster"
noi which cluster
drop _all
label drop _all
prog drop _all
discard
set seed 1001  /* make deterministic */
input case cbr cdr imr
  1.       42        10         81  
  2.       38         9         93  
  3.       39         8         74  
  4.       36        10         90  
  5.       45        16        112  
  6.       31         9         71  
  7.       49        23          .  
  8.       51        20        115  
  9.       48        19        145  
 10.       35         8         77  
 11.       46        15        105  
 12.       49        28        169  
 13.       42        11         72  
 14.       47        23        153  
 15.       44        20        132  
 16.       48        16        127  
 17.       50        22        175  
 18.       50        20        132  
 19.       51        22        141  
 20.       46        17        122  
 21.       46        20        137  
 22.       47        29        175  
 23.       47        14        117  
 24.       47        18        119  
 25.       47        14         96  
 26.       43        18        132  
 27.       46        15        118  
 28.       54        13         76  
 29.       44        16         63  
 30.       53        21        157  
 31.       19         7       27.2  
 32.       45        19        147  
 33.       24         6         12  
 34.       53        16        122  
 35.       26         8       17.4  
 36.       48        17        147  
 37.       50        15        111  
 38.       50        16        108  
 39.       50        13         87  
 40.       47        12         76  
 41.       47        21        143  
 42.       43        16        126  
 43.       44        19        148  
 44.       43        23        143  
 45.       47        13        112  
 46.       38        20        130  
 47.       34        18        112  
 48.       36         9       61.7  
 49.       45        15        103  
 50.       48        14         67  
 51.       41        16        106  
 52.       45        17        111  
 53.       32        10         66  
 54.       47        17        124  
 55.       32         5         32  
 56.       20         8         12  
 57.       48         7         58  
 58.       45        10         86  
 59.       23         7       11.2  
 60.       42         6         56  
 61.       32         3       18.4  
 62.       28         7         51  
 63.       47        14        117  
 64.       30         2         42  
 65.       42         9         85  
 66.       47         9         59  
 67.       31         9         95  
 68.       30         4         38  
 69.       55        22        175  
 70.       48        15        116  
 71.       48        24        183  
 72.       43        17        135  
 73.       38        18        139  
 74.       33        13        104  
 75.       45        13        113  
 76.       47        10         83  
 77.       42        17        112  
 78.       43        15        121  
 79.       24         6       30.5  
 80.       31         3         11  
 81.       34        13        103  
 82.       44        22        166  
 83.       27        10         88  
 84.       40        17        134  
 85.       41        16        122  
 86.       31         7         30  
 87.       35         7         51  
 88.       15         5        9.4  
 89.       29         8         52  
 90.       34         8         53  
 91.       21         7         44  
 92.       13         5        7.7  
 93.       11         6        5.2  
 94.       31         6         33  
 95.       19         6         30  
 96.       23         6         12  
 97.       37        11         53  
 98.       16         5        6.9  
 99.       15         7        7.9  
100.       16         9         10  
101.       36         7         36  
102.       34         6       18.9  
103.       35         5         65  
104.       41         7         65  
105.       39         9         69  
106.       30         8         50  
107.       43         6         69  
108.       27         8         25  
109.       15         5         10  
110.       24         5       26.4  
111.       16         8       10.8  
112.       16         6       13.6  
113.       21         5       23.9  
114.       32         8         70  
115.       26         7         22  
116.       20         7         17  
117.       41        13        117  
118.       23         6         20  
119.       18         6         13  
120.       19         6         10  
121.       19         7       14.9  
122.       24        11       27.8  
123.       28         6       21.5  
124.       26         7       26.5  
125.       29         7       12.7  
126.       40        14        110  
127.       28         8         63  
128.       28         7         48  
129.       36         8         66  
130.       26         6         36  
131.       36         7         45  
132.       34         9         88  
133.       27         7         33  
134.       29         5         36  
135.       24         9       35.3  
136.       22         6       19.5  
137.       18        10       29.5  
138.       11        11        8.4  
139.       12        10        5.8  
140.       16         7        5.4  
141.       17        10        8.7  
142.       13        11        8.5  
143.       12        11        5.9  
144.       13        12        9.5  
145.       12        12       10.3  
146.       12        11        9.7  
147.       14        10          8  
148.       10        11        8.6  
149.       12        11        7.9  
150.       13         9        7.7  
151.       12         9        6.8  
152.       13        11       14.5  
153.       14        12       13.9  
154.       13        13        9.2  
155.       12        14         19  
156.       17        10       17.5  
157.       16        11       25.6  
158.       27         6         43  
159.       11         9       12.3  
160.       10        10       10.1  
161.       15         8         13  
162.       12         9       15.8  
163.       12         8          9  
164.       15         9       27.1  
165.       20        10         25  
166.       15         7        9.8  
167.       28         5         21  
168.       28         4         23  
169.       24         6         36  
170.       16         8       10.8  
171.       34         5         33  
172.       36        12        100  
173.       41         6         42  
174.       38         5         38  
175.       37         8         50
end
label data "1988 World Population Data"
sort case
cluster cbr cdr imr, groups(3) gen(clust)
lab var clust "3 GROUPS, 1/K START"
sort case
cluster cbr cdr imr, groups(3) gen(clust2) random
lab var clust2 "3 GROUPS, RANDOM START"
sort case
cluster cbr cdr imr, groups(3) gen(clust3) first
lab var clust3 "3 GROUPS, FIRST K START"
tab2 clust clust2 clust3