1. Descriptive Analysis

##            vars   n    mean      sd median trimmed     mad min   max range
## SEX*          1 681    1.73    0.45      2    1.78    0.00   1     2     1
## AGE*          2 681    1.99    1.10      2    1.86    1.48   1     4     3
## EDU*          3 681    3.17    1.36      4    3.22    1.48   1     5     4
## CHILD         4  14    1.43    0.51      1    1.42    0.00   1     2     1
## MARRIED       5 681    1.93    0.65      2    1.91    0.00   1     5     4
## EMP           6 681    1.00    0.00      1    1.00    0.00   1     1     0
## PARTFULL*     7 673    1.17    0.38      1    1.09    0.00   1     2     1
## EMPSIZE       8 299    7.32    3.03      8    7.56    2.97   1    11    10
## LYBTINCOME    9 645 2196.18 1366.19   2000 2071.73 1008.17  30 20000 19970
## LYATINCOME   10 651 1973.50 1037.58   1800 1889.99  889.56  30  7440  7410
## OVERALLSAT   11 678    2.80    0.65      3    2.76    0.00   1     5     4
##             skew kurtosis    se
## SEX*       -1.01    -0.99  0.02
## AGE*        0.79    -0.74  0.04
## EDU*       -0.29    -1.32  0.05
## CHILD       0.26    -2.07  0.14
## MARRIED     2.19     9.45  0.03
## EMP          NaN      NaN  0.00
## PARTFULL*   1.73     1.00  0.01
## EMPSIZE    -0.54    -1.12  0.18
## LYBTINCOME  4.05    44.40 53.79
## LYATINCOME  1.07     2.46 40.67
## OVERALLSAT  0.13     0.22  0.03
## 
## Female   Male 
##    187    494
## 
## H30L40 H40L50    H50    L30 
##    298    208     59    116
## 
##    G2yr    G4yr      GG   GHIGH UMIDDLE 
##     102     166      33     272     108
## 
## FULL PART 
##  557  116
## 
##    30    50    80   155   180   200   220   240   260   270   275   288 
##     1     1     2     1     1     1     1     1     2     1     1     1 
##   300   330   360   370   400   420   450   480   500   520   580   600 
##     2     1     1     1     2     2     1     1     2     1     1     8 
##   624   648   650   660   672   680   690   700   720   732   800   840 
##     1     1     1     3     1     1     1     3    10     1     6     7 
##   850   860   880   900   960  1000  1050  1080  1100  1140  1150  1180 
##     2     1     1     6     4    18     1     2     2     1     2     1 
##  1200  1280  1300  1310  1318  1320  1340  1400  1440  1450  1500  1560 
##    21     1    16     2     1    11     1     7    11     2    15    12 
##  1584  1600  1620  1660  1668  1680  1700  1740  1750  1780  1800  1840 
##     1     8     1     1     1     7     7     2     1     2    36     1 
##  1860  1884  1900  1920  1940  1950  2000  2040  2070  2088  2100  2160 
##     2     1     5     9     1     1    33     5     1     1     2     7 
##  2200  2230  2250  2280  2300  2340  2350  2370  2400  2450  2470  2500 
##    17     1     1     3    11     2     1     1    46     1     1    20 
##  2520  2550  2600  2640  2700  2760  2800  2820  2840  2880  2900  2940 
##     3     2    14     5     8     4     8     1     1     2     2     2 
##  2980  3000  3060  3120  3200  3240  3300  3360  3400  3440  3480  3500 
##     2    23     1     3     3     2     8     1     5     1     1    12 
##  3600  3700  3800  3900  3960  4000  4012  4100  4200  4280  4300  4440 
##    10     2     5     3     1     9     1     2     7     1     2     1 
##  4500  4700  4800  4900  5000  5260  5300  5400  5600  6000  6200  6750 
##    11     1     3     1     2     2     2     1     1     2     1     1 
##  7200  8000 20000 
##     2     1     1
## 
##   1   2   3   4   5 
##   7 201 395  71   4
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         58        284         48          8        122        161
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         38        352         42          4         96        149
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         27        382         44          9        117        102
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         26        358         72          6         95        124
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         17        380         87          6         96         95
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         41        403         61          4         69        103
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         44        396         65          8         59        109
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         54        405         63          5         43        111
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         61        400         59          9         53         99
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         63        381         71          6         71         89
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         56        399         64          6         66         90
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         64        389         80          4         41        103
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         61        394         57          5         76         88
## 
## 사회적인정       안정 임금및수입 자율성보장   장래발전   적성취향 
##         72        406         52          8         62         81

2. Sequence

2-1 클러스터 개수 - gap stat 부트스트랩 1000번 7개 클러스터

## 
##   1   2   3   4   5   6   7 
## 145  86  78  67 201  39  65

스퀀스 플롯

세부분석

##      Sequence                                                             
## [1]  (ST,1)-(IN,1)-(ST,1)-(IN,1)-(FU,1)-(ST,4)-(FU,1)-(ST,4)              
## [2]  (ST,1)-(FW,1)-(ST,2)-(FU,1)-(ST,4)-(IN,3)-(SA,1)-(IN,1)              
## [3]  (IN,9)-(FU,1)-(ST,2)-(SA,1)-(ST,1)                                   
## [4]  (FU,1)-(ST,1)-(AU,1)-(FW,1)-(IN,2)-(ST,5)-(IN,1)-(ST,2)              
## [5]  (IN,1)-(FU,1)-(ST,1)-(IN,2)-(ST,2)-(IN,1)-(ST,1)-(FU,2)-(IN,1)-(ST,2)
## [6]  (SA,1)-(ST,1)-(FW,1)-(ST,1)-(FW,1)-(ST,6)-(AU,1)-(ST,2)              
## [7]  (FU,2)-(ST,1)-(FU,1)-(AU,1)-(ST,8)-(FW,1)                            
## [8]  (IN,1)-(ST,1)-(FU,1)-(ST,2)-(AU,1)-(ST,8)                            
## [9]  (FU,1)-(SA,1)-(FU,1)-(SA,1)-(FU,2)-(SA,2)-(FU,1)-(SA,3)-(FW,1)-(SA,1)
## [10] (SA,1)-(ST,9)-(IN,1)-(ST,2)-(FW,1)

##                          Freq Percent
## ST/14                       7    1.03
## FW/1-ST/13                  3    0.44
## FU/1-ST/13                  2    0.29
## IN/1-ST/13                  2    0.29
## IN/1-ST/3-IN/1-ST/9         2    0.29
## ST/1-FU/1-ST/12             2    0.29
## ST/1-IN/1-ST/12             2    0.29
## ST/13-IN/1                  2    0.29
## ST/2-FU/1-ST/1-FW/1-ST/9    2    0.29
## ST/4-FU/1-ST/9              2    0.29

## data$SEX: Female
##    Mean
## SA 0.90
## ST 8.28
## FW 1.21
## AU 0.17
## FU 1.21
## IN 2.22
## -------------------------------------------------------- 
## data$SEX: Male
##    Mean
## SA 1.04
## ST 7.65
## FW 1.29
## AU 0.11
## FU 1.70
## IN 2.20
##         [-> SA] [-> ST] [-> FW] [-> AU] [-> FU] [-> IN]
## [SA ->]    0.23    0.43    0.06    0.00    0.15    0.13
## [ST ->]    0.06    0.68    0.07    0.01    0.08    0.11
## [FW ->]    0.05    0.44    0.25    0.01    0.10    0.14
## [AU ->]    0.06    0.42    0.05    0.08    0.19    0.20
## [FU ->]    0.08    0.43    0.11    0.01    0.20    0.17
## [IN ->]    0.05    0.45    0.07    0.01    0.12    0.29
##     [State frequencies]
##     DJ06   DJ07  DJ08   DJ09   DJ10   DJ11  DJ12   DJ13  DJ14   DJ15
## SA 0.085 0.0558 0.040 0.0382 0.0250 0.0602 0.065 0.0793 0.090 0.0925
## ST 0.417 0.5169 0.561 0.5257 0.5580 0.5918 0.581 0.5947 0.587 0.5595
## FW 0.070 0.0617 0.065 0.1057 0.1278 0.0896 0.095 0.0925 0.087 0.1043
## AU 0.012 0.0059 0.013 0.0088 0.0088 0.0059 0.012 0.0073 0.013 0.0088
## FU 0.179 0.1410 0.172 0.1395 0.1410 0.1013 0.087 0.0631 0.078 0.1043
## IN 0.236 0.2188 0.150 0.1821 0.1395 0.1512 0.160 0.1630 0.145 0.1307
##      DJ16   DJ17   DJ18  DJ19
## SA 0.0822 0.0940 0.0896 0.106
## ST 0.5859 0.5712 0.5786 0.596
## FW 0.0940 0.1175 0.0837 0.076
## AU 0.0088 0.0059 0.0073 0.012
## FU 0.0969 0.0602 0.1116 0.091
## IN 0.1322 0.1512 0.1292 0.119
## 
##      [Valid states]
##    DJ06 DJ07 DJ08 DJ09 DJ10 DJ11 DJ12 DJ13 DJ14 DJ15 DJ16 DJ17 DJ18 DJ19
## N   681  681  681  681  681  681  681  681  681  681  681  681  681  681
## 
##      [Entropy index]
##    DJ06 DJ07 DJ08 DJ09 DJ10 DJ11 DJ12 DJ13 DJ14 DJ15 DJ16 DJ17 DJ18 DJ19
## H  0.82 0.73 0.71 0.74 0.71 0.69 0.71 0.69 0.71 0.74 0.71 0.71 0.72 0.71

##   Trans.
## 1      7
## 2      7
## 3      4
## 4      7
## 5      9
## 6      7
##   SA ST FW AU FU IN
## 1  0 10  0  0  2  2
## 2  1  7  1  0  1  4
## 3  1  3  0  0  1  9
## 4  0  8  1  1  1  3
##   Subseq.
## 1     128
## 2     168
## 3      28
## 4     184

카이제곱 분석 & ANOVA

## 
##  Pearson's Chi-squared test
## 
## data:  tabSEX
## X-squared = 12.114, df = 6, p-value = 0.05948
## 
##  
##    Cell Contents
## |-------------------------|
## |                       N |
## | Chi-square contribution |
## |           N / Row Total |
## |           N / Col Total |
## |         N / Table Total |
## |-------------------------|
## 
##  
## Total Observations in Table:  681 
## 
##  
##              | data$cluster 
##     data$SEX |         1 |         2 |         3 |         4 |         5 |         6 |         7 | Row Total | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##       Female |        34 |        20 |        22 |        12 |        68 |         8 |        23 |       187 | 
##              |     0.850 |     0.553 |     0.016 |     2.225 |     2.971 |     0.685 |     1.487 |           | 
##              |     0.182 |     0.107 |     0.118 |     0.064 |     0.364 |     0.043 |     0.123 |     0.275 | 
##              |     0.234 |     0.233 |     0.282 |     0.179 |     0.338 |     0.205 |     0.354 |           | 
##              |     0.050 |     0.029 |     0.032 |     0.018 |     0.100 |     0.012 |     0.034 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##         Male |       111 |        66 |        56 |        55 |       133 |        31 |        42 |       494 | 
##              |     0.322 |     0.210 |     0.006 |     0.842 |     1.125 |     0.259 |     0.563 |           | 
##              |     0.225 |     0.134 |     0.113 |     0.111 |     0.269 |     0.063 |     0.085 |     0.725 | 
##              |     0.766 |     0.767 |     0.718 |     0.821 |     0.662 |     0.795 |     0.646 |           | 
##              |     0.163 |     0.097 |     0.082 |     0.081 |     0.195 |     0.046 |     0.062 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
## Column Total |       145 |        86 |        78 |        67 |       201 |        39 |        65 |       681 | 
##              |     0.213 |     0.126 |     0.115 |     0.098 |     0.295 |     0.057 |     0.095 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
## 
## 
## 
##  Pearson's Chi-squared test
## 
## data:  tabAGE
## X-squared = 53.953, df = 18, p-value = 1.866e-05
## 
##  
##    Cell Contents
## |-------------------------|
## |                       N |
## | Chi-square contribution |
## |           N / Row Total |
## |           N / Col Total |
## |         N / Table Total |
## |-------------------------|
## 
##  
## Total Observations in Table:  681 
## 
##  
##              | data$cluster 
##     data$AGE |         1 |         2 |         3 |         4 |         5 |         6 |         7 | Row Total | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##       H30L40 |        65 |        31 |        34 |        33 |        88 |        18 |        29 |       298 | 
##              |     0.038 |     1.169 |     0.001 |     0.462 |     0.000 |     0.051 |     0.011 |           | 
##              |     0.218 |     0.104 |     0.114 |     0.111 |     0.295 |     0.060 |     0.097 |     0.438 | 
##              |     0.448 |     0.360 |     0.436 |     0.493 |     0.438 |     0.462 |     0.446 |           | 
##              |     0.095 |     0.046 |     0.050 |     0.048 |     0.129 |     0.026 |     0.043 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##       H40L50 |        36 |        35 |        19 |        13 |        67 |         7 |        31 |       208 | 
##              |     1.551 |     2.903 |     0.977 |     2.722 |     0.512 |     2.025 |     6.259 |           | 
##              |     0.173 |     0.168 |     0.091 |     0.062 |     0.322 |     0.034 |     0.149 |     0.305 | 
##              |     0.248 |     0.407 |     0.244 |     0.194 |     0.333 |     0.179 |     0.477 |           | 
##              |     0.053 |     0.051 |     0.028 |     0.019 |     0.098 |     0.010 |     0.046 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##          H50 |        12 |         9 |         3 |         3 |        25 |         4 |         3 |        59 | 
##              |     0.025 |     0.322 |     2.090 |     1.355 |     3.305 |     0.114 |     1.230 |           | 
##              |     0.203 |     0.153 |     0.051 |     0.051 |     0.424 |     0.068 |     0.051 |     0.087 | 
##              |     0.083 |     0.105 |     0.038 |     0.045 |     0.124 |     0.103 |     0.046 |           | 
##              |     0.018 |     0.013 |     0.004 |     0.004 |     0.037 |     0.006 |     0.004 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##          L30 |        32 |        11 |        22 |        18 |        21 |        10 |         2 |       116 | 
##              |     2.158 |     0.909 |     5.715 |     3.802 |     5.118 |     1.696 |     7.433 |           | 
##              |     0.276 |     0.095 |     0.190 |     0.155 |     0.181 |     0.086 |     0.017 |     0.170 | 
##              |     0.221 |     0.128 |     0.282 |     0.269 |     0.104 |     0.256 |     0.031 |           | 
##              |     0.047 |     0.016 |     0.032 |     0.026 |     0.031 |     0.015 |     0.003 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
## Column Total |       145 |        86 |        78 |        67 |       201 |        39 |        65 |       681 | 
##              |     0.213 |     0.126 |     0.115 |     0.098 |     0.295 |     0.057 |     0.095 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
## 
## 
## 
##  Pearson's Chi-squared test
## 
## data:  tabEDU
## X-squared = 67.203, df = 24, p-value = 5.751e-06
## 
##  
##    Cell Contents
## |-------------------------|
## |                       N |
## | Chi-square contribution |
## |           N / Row Total |
## |           N / Col Total |
## |         N / Table Total |
## |-------------------------|
## 
##  
## Total Observations in Table:  681 
## 
##  
##              | data$cluster 
##     data$EDU |         1 |         2 |         3 |         4 |         5 |         6 |         7 | Row Total | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##         G2yr |        30 |        15 |         8 |        11 |        25 |         8 |         5 |       102 | 
##              |     3.158 |     0.349 |     1.161 |     0.093 |     0.866 |     0.798 |     2.304 |           | 
##              |     0.294 |     0.147 |     0.078 |     0.108 |     0.245 |     0.078 |     0.049 |     0.150 | 
##              |     0.207 |     0.174 |     0.103 |     0.164 |     0.124 |     0.205 |     0.077 |           | 
##              |     0.044 |     0.022 |     0.012 |     0.016 |     0.037 |     0.012 |     0.007 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##         G4yr |        39 |        13 |        26 |        27 |        43 |         6 |        12 |       166 | 
##              |     0.378 |     3.025 |     2.567 |     6.969 |     0.734 |     1.293 |     0.933 |           | 
##              |     0.235 |     0.078 |     0.157 |     0.163 |     0.259 |     0.036 |     0.072 |     0.244 | 
##              |     0.269 |     0.151 |     0.333 |     0.403 |     0.214 |     0.154 |     0.185 |           | 
##              |     0.057 |     0.019 |     0.038 |     0.040 |     0.063 |     0.009 |     0.018 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##           GG |         5 |         3 |        11 |         4 |         5 |         4 |         1 |        33 | 
##              |     0.584 |     0.327 |    13.793 |     0.175 |     2.307 |     2.356 |     1.467 |           | 
##              |     0.152 |     0.091 |     0.333 |     0.121 |     0.152 |     0.121 |     0.030 |     0.048 | 
##              |     0.034 |     0.035 |     0.141 |     0.060 |     0.025 |     0.103 |     0.015 |           | 
##              |     0.007 |     0.004 |     0.016 |     0.006 |     0.007 |     0.006 |     0.001 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##        GHIGH |        58 |        41 |        25 |        18 |        83 |        15 |        32 |       272 | 
##              |     0.000 |     1.288 |     1.216 |     2.868 |     0.092 |     0.021 |     1.404 |           | 
##              |     0.213 |     0.151 |     0.092 |     0.066 |     0.305 |     0.055 |     0.118 |     0.399 | 
##              |     0.400 |     0.477 |     0.321 |     0.269 |     0.413 |     0.385 |     0.492 |           | 
##              |     0.085 |     0.060 |     0.037 |     0.026 |     0.122 |     0.022 |     0.047 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##      UMIDDLE |        13 |        14 |         8 |         7 |        45 |         6 |        15 |       108 | 
##              |     4.345 |     0.010 |     1.544 |     1.237 |     5.403 |     0.006 |     2.135 |           | 
##              |     0.120 |     0.130 |     0.074 |     0.065 |     0.417 |     0.056 |     0.139 |     0.159 | 
##              |     0.090 |     0.163 |     0.103 |     0.104 |     0.224 |     0.154 |     0.231 |           | 
##              |     0.019 |     0.021 |     0.012 |     0.010 |     0.066 |     0.009 |     0.022 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
## Column Total |       145 |        86 |        78 |        67 |       201 |        39 |        65 |       681 | 
##              |     0.213 |     0.126 |     0.115 |     0.098 |     0.295 |     0.057 |     0.095 |           | 
## -------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
## 
## 
## 
##  Pearson's Chi-squared test
## 
## data:  tabPARTFULL
## X-squared = 12.087, df = 6, p-value = 0.06005
## 
##  
##    Cell Contents
## |-------------------------|
## |                       N |
## | Chi-square contribution |
## |           N / Row Total |
## |           N / Col Total |
## |         N / Table Total |
## |-------------------------|
## 
##  
## Total Observations in Table:  673 
## 
##  
##               | data$cluster 
## data$PARTFULL |         1 |         2 |         3 |         4 |         5 |         6 |         7 | Row Total | 
## --------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##          FULL |       120 |        71 |        69 |        61 |       159 |        30 |        47 |       557 | 
##               |     0.093 |     0.006 |     0.306 |     0.555 |     0.197 |     0.067 |     0.859 |           | 
##               |     0.215 |     0.127 |     0.124 |     0.110 |     0.285 |     0.054 |     0.084 |     0.828 | 
##               |     0.851 |     0.835 |     0.885 |     0.910 |     0.799 |     0.789 |     0.723 |           | 
##               |     0.178 |     0.105 |     0.103 |     0.091 |     0.236 |     0.045 |     0.070 |           | 
## --------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##          PART |        21 |        14 |         9 |         6 |        40 |         8 |        18 |       116 | 
##               |     0.449 |     0.029 |     1.469 |     2.666 |     0.947 |     0.321 |     4.123 |           | 
##               |     0.181 |     0.121 |     0.078 |     0.052 |     0.345 |     0.069 |     0.155 |     0.172 | 
##               |     0.149 |     0.165 |     0.115 |     0.090 |     0.201 |     0.211 |     0.277 |           | 
##               |     0.031 |     0.021 |     0.013 |     0.009 |     0.059 |     0.012 |     0.027 |           | 
## --------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
##  Column Total |       141 |        85 |        78 |        67 |       199 |        38 |        65 |       673 | 
##               |     0.210 |     0.126 |     0.116 |     0.100 |     0.296 |     0.056 |     0.097 |           | 
## --------------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|-----------|
## 
## 
##              Df    Sum Sq Mean Sq F value Pr(>F)  
## cluster       1 5.901e+06 5901412   3.172 0.0754 .
## Residuals   643 1.196e+09 1860198                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 36 observations deleted due to missingness
##              Df Sum Sq Mean Sq F value  Pr(>F)   
## cluster       1   3.92   3.917   9.298 0.00238 **
## Residuals   676 284.80   0.421                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 3 observations deleted due to missingness

아노바 디테일 06, 19 및 사후검증

06

## # A tibble: 7 x 2
##   cluster LYBTINCOME
##     <int>      <dbl>
## 1       1      2338.
## 2       2      2106.
## 3       3      2443.
## 4       4      2285.
## 5       5      2054.
## 6       6      2449.
## 7       7      1912.
##                     Df    Sum Sq Mean Sq F value Pr(>F)
## as.factor(cluster)   6 1.945e+07 3240919   1.748  0.107
## Residuals          638 1.183e+09 1853548               
## 36 observations deleted due to missingness
## 
##  Pairwise comparisons using t tests with pooled SD 
## 
## data:  data$LYBTINCOME and data$cluster 
## 
##   1    2    3    4    5    6   
## 2 1.00 -    -    -    -    -   
## 3 1.00 1.00 -    -    -    -   
## 4 1.00 1.00 1.00 -    -    -   
## 5 1.00 1.00 0.80 1.00 -    -   
## 6 1.00 1.00 1.00 1.00 1.00 -   
## 7 0.92 1.00 0.54 1.00 1.00 1.00
## 
## P value adjustment method: bonferroni
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = LYBTINCOME ~ as.factor(cluster), data = data)
## 
## $`as.factor(cluster)`
##            diff        lwr       upr     p adj
## 2-1 -231.660294  -788.4411  325.1206 0.8820901
## 3-1  104.930197  -479.3341  689.1945 0.9983857
## 4-1  -52.837150  -673.3715  567.6972 0.9999789
## 5-1 -283.644506  -734.4788  167.1897 0.5072592
## 6-1  111.540342  -635.1144  858.1950 0.9994348
## 7-1 -425.819118 -1049.9131  198.2749 0.4041329
## 3-2  336.590492  -305.9858  979.1668 0.7146502
## 4-2  178.823144  -496.9005  854.5468 0.9865749
## 5-2  -51.984212  -576.1894  472.2210 0.9999481
## 6-2  343.200636  -449.9153 1136.3165 0.8610867
## 7-2 -194.158824  -873.1529  484.8352 0.9799487
## 4-3 -157.767348  -856.3103  540.7756 0.9942442
## 5-3 -388.574704  -941.8838  164.7344 0.3675573
## 6-3    6.610144  -806.0353  819.2556 1.0000000
## 7-3 -530.749315 -1232.4564  170.9577 0.2770950
## 5-4 -230.807356  -822.2877  360.6730 0.9107210
## 6-4  164.377492  -674.7237 1003.4787 0.9973720
## 7-4 -372.981967 -1105.1642  359.2003 0.7408319
## 6-5  395.184848  -327.5042 1117.8739 0.6711383
## 7-5 -142.174611  -737.3885  453.0393 0.9922081
## 7-6 -537.359459 -1379.0965  304.3776 0.4888922

## 
##   Posthoc multiple comparisons of means : Fisher LSD 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##            diff      lwr.ci    upr.ci   pval    
## 2-1 -231.660294  -601.31127 137.99068 0.2189    
## 3-1  104.930197  -282.96723 492.82762 0.5955    
## 4-1  -52.837150  -464.81454 359.14024 0.8012    
## 5-1 -283.644506  -582.95675  15.66773 0.0632 .  
## 6-1  111.540342  -384.16931 607.24999 0.6587    
## 7-1 -425.819118  -840.15982 -11.47842 0.0440 *  
## 3-2  336.590492   -90.02074 763.20173 0.1218    
## 4-2  178.823144  -269.79485 627.44113 0.4341    
## 5-2  -51.984212  -400.00799 296.03957 0.7694    
## 6-2  343.200636  -183.35495 869.75623 0.2010    
## 7-2 -194.158824  -644.94807 256.63043 0.3980    
## 4-3 -157.767348  -621.53525 306.00056 0.5044    
## 5-3 -388.574704  -755.92075 -21.22866 0.0382 *  
## 6-3    6.610144  -532.91127 546.13156 0.9808    
## 7-3 -530.749315  -996.61788 -64.88075 0.0256 *  
## 5-4 -230.807356  -623.49562 161.88091 0.2489    
## 6-4  164.377492  -392.70806 721.46305 0.5625    
## 7-4 -372.981967  -859.08326 113.11933 0.1324    
## 6-5  395.184848   -84.61382 874.98352 0.1063    
## 7-5 -142.174611  -537.34156 252.99234 0.4801    
## 7-6 -537.359459 -1096.19500  21.47608 0.0594 .  
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##   Posthoc multiple comparisons of means : Tukey HSD 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##            diff     lwr.ci    upr.ci   pval    
## 2-1 -231.660294  -788.4411  325.1206 0.8821    
## 3-1  104.930197  -479.3341  689.1945 0.9984    
## 4-1  -52.837150  -673.3715  567.6972 1.0000    
## 5-1 -283.644506  -734.4788  167.1897 0.5073    
## 6-1  111.540342  -635.1144  858.1950 0.9994    
## 7-1 -425.819118 -1049.9131  198.2749 0.4041    
## 3-2  336.590492  -305.9858  979.1668 0.7147    
## 4-2  178.823144  -496.9005  854.5468 0.9866    
## 5-2  -51.984212  -576.1894  472.2210 0.9999    
## 6-2  343.200636  -449.9153 1136.3165 0.8611    
## 7-2 -194.158824  -873.1529  484.8352 0.9799    
## 4-3 -157.767348  -856.3103  540.7756 0.9942    
## 5-3 -388.574704  -941.8838  164.7344 0.3676    
## 6-3    6.610144  -806.0353  819.2556 1.0000    
## 7-3 -530.749315 -1232.4564  170.9577 0.2771    
## 5-4 -230.807356  -822.2877  360.6730 0.9107    
## 6-4  164.377492  -674.7237 1003.4787 0.9974    
## 7-4 -372.981967 -1105.1642  359.2003 0.7408    
## 6-5  395.184848  -327.5042 1117.8739 0.6711    
## 7-5 -142.174611  -737.3885  453.0393 0.9922    
## 7-6 -537.359459 -1379.0965  304.3776 0.4889    
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##   Posthoc multiple comparisons of means : Scheffe Test 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##            diff     lwr.ci    upr.ci   pval    
## 2-1 -231.660294  -901.8861  438.5655 0.9583    
## 3-1  104.930197  -598.3788  808.2392 0.9996    
## 4-1  -52.837150  -799.8063  694.1320 1.0000    
## 5-1 -283.644506  -826.3369  259.0479 0.7487    
## 6-1  111.540342  -787.2464 1010.3271 0.9999    
## 7-1 -425.819118 -1177.0733  325.4350 0.6669    
## 3-2  336.590492  -436.9118 1110.0927 0.8791    
## 4-2  178.823144  -634.5802  992.2265 0.9962    
## 5-2  -51.984212  -682.9971  579.0286 1.0000    
## 6-2  343.200636  -611.5139 1297.9151 0.9496    
## 7-2 -194.158824 -1011.4990  623.1813 0.9941    
## 4-3 -157.767348  -998.6395  683.1048 0.9984    
## 5-3 -388.574704 -1054.6214  277.4720 0.6343    
## 6-3    6.610144  -971.6131  984.8334 1.0000    
## 7-3 -530.749315 -1375.4303  313.9316 0.5437    
## 5-4 -230.807356  -942.8028  481.1881 0.9697    
## 6-4  164.377492  -845.6918 1174.4468 0.9993    
## 7-4 -372.981967 -1254.3475  508.3836 0.8929    
## 6-5  395.184848  -474.7532 1265.1229 0.8549    
## 7-5 -142.174611  -858.6642  574.3150 0.9978    
## 7-6 -537.359459 -1550.6018  475.8828 0.7351    
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## # A tibble: 7 x 2
##   cluster OVERALLSAT
##     <int>      <dbl>
## 1       1       2.73
## 2       2       2.80
## 3       3       2.59
## 4       4       2.76
## 5       5       2.86
## 6       6       2.92
## 7       7       2.98
##                     Df Sum Sq Mean Sq F value  Pr(>F)   
## as.factor(cluster)   6   7.58   1.263   3.015 0.00648 **
## Residuals          671 281.14   0.419                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 3 observations deleted due to missingness
## 
##  Pairwise comparisons using t tests with pooled SD 
## 
## data:  data$OVERALLSAT and as.factor(data$cluster) 
## 
##   1      2      3      4      5      6     
## 2 1.0000 -      -      -      -      -     
## 3 1.0000 0.7574 -      -      -      -     
## 4 1.0000 1.0000 1.0000 -      -      -     
## 5 1.0000 1.0000 0.0468 1.0000 -      -     
## 6 1.0000 1.0000 0.1857 1.0000 1.0000 -     
## 7 0.2086 1.0000 0.0063 1.0000 1.0000 1.0000
## 
## P value adjustment method: bonferroni
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = OVERALLSAT ~ as.factor(cluster), data = data)
## 
## $`as.factor(cluster)`
##            diff          lwr        upr     p adj
## 2-1  0.06805985 -0.193155332 0.32927503 0.9876427
## 3-1 -0.14452214 -0.413972481 0.12492819 0.6913167
## 4-1  0.02692830 -0.256473202 0.31032979 0.9999597
## 5-1  0.12073427 -0.088899848 0.33036838 0.6141481
## 6-1  0.18881119 -0.156995417 0.53461779 0.6728879
## 7-1  0.25034965 -0.036005424 0.53670472 0.1320271
## 3-2 -0.21258199 -0.511893427 0.08672944 0.3535434
## 4-2 -0.04113155 -0.353061448 0.27079835 0.9997256
## 5-2  0.05267442 -0.194166245 0.29951508 0.9957847
## 6-2  0.12075134 -0.248797113 0.49029980 0.9609897
## 7-2  0.18228980 -0.132325963 0.49690557 0.6072166
## 4-3  0.17145044 -0.147407482 0.49030836 0.6887692
## 5-3  0.26525641  0.009716927 0.52079589 0.0360146
## 6-3  0.33333333 -0.042081336 0.70874800 0.1198668
## 7-3  0.39487179  0.073385879 0.71635771 0.0055800
## 5-4  0.09380597 -0.176403860 0.36401580 0.9478225
## 6-4  0.16188289 -0.223667491 0.54743328 0.8775636
## 7-4  0.22342135 -0.109844484 0.55668719 0.4267366
## 6-5  0.06807692 -0.267003891 0.40315774 0.9967860
## 7-5  0.12961538 -0.143690620 0.40292139 0.8005375
## 7-6  0.06153846 -0.326188142 0.44926506 0.9992011

## 
##   Posthoc multiple comparisons of means : Fisher LSD 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##            diff       lwr.ci      upr.ci   pval    
## 2-1  0.06805985 -0.105373971  0.24149367 0.4413    
## 3-1 -0.14452214 -0.323423695  0.03437941 0.1132    
## 4-1  0.02692830 -0.161236129  0.21509272 0.7788    
## 5-1  0.12073427 -0.018452304  0.25992084 0.0890 .  
## 6-1  0.18881119 -0.040787107  0.41840948 0.1068    
## 7-1  0.25034965  0.060224199  0.44047510 0.0099 ** 
## 3-2 -0.21258199 -0.411309819 -0.01385416 0.0361 *  
## 4-2 -0.04113155 -0.248237407  0.16597430 0.6967    
## 5-2  0.05267442 -0.111215440  0.21656428 0.5282    
## 6-2  0.12075134 -0.124610355  0.36611304 0.3342    
## 7-2  0.18228980 -0.026599335  0.39117894 0.0871 .  
## 4-3  0.17145044 -0.040255277  0.38315616 0.1123    
## 5-3  0.26525641  0.095590970  0.43492185 0.0022 ** 
## 6-3  0.33333333  0.084076763  0.58258990 0.0088 ** 
## 7-3  0.39487179  0.181421221  0.60832237 0.0003 ***
## 5-4  0.09380597 -0.085599847  0.27321179 0.3050    
## 6-4  0.16188289 -0.094103285  0.41786907 0.2148    
## 7-4  0.22342135  0.002149501  0.44469321 0.0478 *  
## 6-5  0.06807692 -0.154399984  0.29055383 0.5482    
## 7-5  0.12961538 -0.051846137  0.31107691 0.1612    
## 7-6  0.06153846 -0.195892617  0.31896954 0.6390    
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##   Posthoc multiple comparisons of means : Tukey HSD 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##            diff       lwr.ci     upr.ci   pval    
## 2-1  0.06805985 -0.193155332 0.32927503 0.9876    
## 3-1 -0.14452214 -0.413972481 0.12492819 0.6913    
## 4-1  0.02692830 -0.256473202 0.31032979 1.0000    
## 5-1  0.12073427 -0.088899848 0.33036838 0.6141    
## 6-1  0.18881119 -0.156995417 0.53461779 0.6729    
## 7-1  0.25034965 -0.036005424 0.53670472 0.1320    
## 3-2 -0.21258199 -0.511893427 0.08672944 0.3535    
## 4-2 -0.04113155 -0.353061448 0.27079835 0.9997    
## 5-2  0.05267442 -0.194166245 0.29951508 0.9958    
## 6-2  0.12075134 -0.248797113 0.49029980 0.9610    
## 7-2  0.18228980 -0.132325963 0.49690557 0.6072    
## 4-3  0.17145044 -0.147407482 0.49030836 0.6888    
## 5-3  0.26525641  0.009716927 0.52079589 0.0360 *  
## 6-3  0.33333333 -0.042081336 0.70874800 0.1199    
## 7-3  0.39487179  0.073385879 0.71635771 0.0056 ** 
## 5-4  0.09380597 -0.176403860 0.36401580 0.9478    
## 6-4  0.16188289 -0.223667491 0.54743328 0.8776    
## 7-4  0.22342135 -0.109844484 0.55668719 0.4267    
## 6-5  0.06807692 -0.267003891 0.40315774 0.9968    
## 7-5  0.12961538 -0.143690620 0.40292139 0.8005    
## 7-6  0.06153846 -0.326188142 0.44926506 0.9992    
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##   Posthoc multiple comparisons of means : Scheffe Test 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##            diff       lwr.ci    upr.ci   pval    
## 2-1  0.06805985 -0.246375811 0.3824955 0.9965    
## 3-1 -0.14452214 -0.468870806 0.1798265 0.8664    
## 4-1  0.02692830 -0.314213963 0.3680706 1.0000    
## 5-1  0.12073427 -0.131611102 0.3730796 0.8209    
## 6-1  0.18881119 -0.227450714 0.6050731 0.8560    
## 7-1  0.25034965 -0.094347953 0.5950473 0.3523    
## 3-2 -0.21258199 -0.572875710 0.1477117 0.6213    
## 4-2 -0.04113155 -0.416614640 0.3343515 0.9999    
## 5-2  0.05267442 -0.244458033 0.3498069 0.9989    
## 6-2  0.12075134 -0.324089620 0.5655923 0.9879    
## 7-2  0.18228980 -0.196426380 0.5610060 0.8165    
## 4-3  0.17145044 -0.212372203 0.5552731 0.8649    
## 5-3  0.26525641 -0.042347174 0.5728600 0.1530    
## 6-3  0.33333333 -0.118569037 0.7852357 0.3321    
## 7-3  0.39487179  0.007885725 0.7818579 0.0414 *  
## 5-4  0.09380597 -0.231456926 0.4190689 0.9834    
## 6-4  0.16188289 -0.302220262 0.6259860 0.9565    
## 7-4  0.22342135 -0.177744702 0.6245874 0.6860    
## 6-5  0.06807692 -0.335273895 0.4714277 0.9991    
## 7-5  0.12961538 -0.199374507 0.4586053 0.9224    
## 7-6  0.06153846 -0.405184300 0.5282612 0.9998    
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

19

## # A tibble: 7 x 2
##   cluster LYBTINCOME19
##     <int>        <dbl>
## 1       1        4435.
## 2       2        4149.
## 3       3        5081.
## 4       4        5050.
## 5       5        4116.
## 6       6        3964.
## 7       7        3552.
##                     Df    Sum Sq  Mean Sq F value  Pr(>F)   
## as.factor(cluster)   6 1.363e+08 22721690   3.361 0.00285 **
## Residuals          673 4.550e+09  6760943                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 1 observation deleted due to missingness
## 
##  Pairwise comparisons using t tests with pooled SD 
## 
## data:  data$LYBTINCOME19 and data$cluster 
## 
##   1     2     3     4     5     6    
## 2 1.000 -     -     -     -     -    
## 3 1.000 0.468 -     -     -     -    
## 4 1.000 0.729 1.000 -     -     -    
## 5 1.000 1.000 0.117 0.243 -     -    
## 6 1.000 1.000 0.606 0.820 1.000 -    
## 7 0.489 1.000 0.010 0.022 1.000 1.000
## 
## P value adjustment method: bonferroni
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = LYBTINCOME19 ~ as.factor(cluster), data = data)
## 
## $`as.factor(cluster)`
##            diff        lwr        upr     p adj
## 2-1  -285.29864 -1331.8710  761.27372 0.9843477
## 3-1   646.02635  -433.7098 1725.76245 0.5693044
## 4-1   615.13124  -526.6470 1756.90946 0.6867879
## 5-1  -319.04560 -1156.8691  518.77786 0.9200188
## 6-1  -470.76852 -1857.8089  916.27185 0.9531430
## 7-1  -882.50186 -2030.2999  265.29614 0.2583297
## 3-2   931.32499  -271.0001 2133.65003 0.2499341
## 4-2   900.42988  -357.9075 2158.76725 0.3442209
## 5-2   -33.74696 -1024.5565  957.06260 0.9999999
## 6-2  -185.46989 -1669.9349 1298.99515 0.9997994
## 7-2  -597.20322 -1861.0053  666.59886 0.8032551
## 4-3   -30.89510 -1316.9469 1255.15673 1.0000000
## 5-3  -965.07195 -1990.8497   60.70578 0.0807448
## 6-3 -1116.79487 -2624.8243  391.23457 0.3022142
## 7-3 -1528.52821 -2819.9275 -237.12894 0.0089262
## 5-4  -934.17684 -2025.0698  156.71609 0.1493978
## 6-4 -1085.89977 -2638.9549  467.15537 0.3732290
## 7-4 -1497.63310 -2841.3368 -153.92945 0.0177691
## 6-5  -151.72292 -1497.1860 1193.74015 0.9998898
## 7-5  -563.45626 -1660.6482  533.73568 0.7335734
## 7-6  -411.73333 -1969.2194 1145.75278 0.9866655

## 
##   Posthoc multiple comparisons of means : Fisher LSD 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##            diff      lwr.ci     upr.ci   pval    
## 2-1  -285.29864  -980.17291  409.57564 0.4204    
## 3-1   646.02635   -70.86708 1362.91977 0.0773 .  
## 4-1   615.13124  -142.95519 1373.21768 0.1116    
## 5-1  -319.04560  -875.32053  237.22933 0.2605    
## 6-1  -470.76852 -1391.69737  450.16032 0.3159    
## 7-1  -882.50186 -1644.58514 -120.41857 0.0233 *  
## 3-2   931.32499   133.03834 1729.61163 0.0223 *  
## 4-2   900.42988    64.95371 1735.90605 0.0347 *  
## 5-2   -33.74696  -691.59739  624.10347 0.9198    
## 6-2  -185.46989 -1171.08407  800.14430 0.7119    
## 7-2  -597.20322 -1436.30770  241.90126 0.1627    
## 4-3   -30.89510  -884.77236  822.98215 0.9434    
## 5-3  -965.07195 -1646.13958 -284.00432 0.0055 ** 
## 6-3 -1116.79487 -2118.05471 -115.53504 0.0289 *  
## 7-3 -1528.52821 -2385.95590 -671.10051 0.0005 ***
## 5-4  -934.17684 -1658.47787 -209.87581 0.0116 *  
## 6-4 -1085.89977 -2117.05452  -54.74502 0.0390 *  
## 7-4 -1497.63310 -2389.78842 -605.47778 0.0010 ** 
## 6-5  -151.72292 -1045.04642  741.60057 0.7389    
## 7-5  -563.45626 -1291.93952  165.02701 0.1293    
## 7-6  -411.73333 -1445.83004  622.36338 0.4346    
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##   Posthoc multiple comparisons of means : Tukey HSD 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##            diff     lwr.ci     upr.ci   pval    
## 2-1  -285.29864 -1331.8710  761.27372 0.9843    
## 3-1   646.02635  -433.7098 1725.76245 0.5693    
## 4-1   615.13124  -526.6470 1756.90946 0.6868    
## 5-1  -319.04560 -1156.8691  518.77786 0.9200    
## 6-1  -470.76852 -1857.8089  916.27185 0.9531    
## 7-1  -882.50186 -2030.2999  265.29614 0.2583    
## 3-2   931.32499  -271.0001 2133.65003 0.2499    
## 4-2   900.42988  -357.9075 2158.76725 0.3442    
## 5-2   -33.74696 -1024.5565  957.06260 1.0000    
## 6-2  -185.46989 -1669.9349 1298.99515 0.9998    
## 7-2  -597.20322 -1861.0053  666.59886 0.8033    
## 4-3   -30.89510 -1316.9469 1255.15673 1.0000    
## 5-3  -965.07195 -1990.8497   60.70578 0.0807 .  
## 6-3 -1116.79487 -2624.8243  391.23457 0.3022    
## 7-3 -1528.52821 -2819.9275 -237.12894 0.0089 ** 
## 5-4  -934.17684 -2025.0698  156.71609 0.1494    
## 6-4 -1085.89977 -2638.9549  467.15537 0.3732    
## 7-4 -1497.63310 -2841.3368 -153.92945 0.0178 *  
## 6-5  -151.72292 -1497.1860 1193.74015 0.9999    
## 7-5  -563.45626 -1660.6482  533.73568 0.7336    
## 7-6  -411.73333 -1969.2194 1145.75278 0.9867    
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##   Posthoc multiple comparisons of means : Scheffe Test 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##            diff     lwr.ci     upr.ci   pval    
## 2-1  -285.29864 -1545.1011  974.50379 0.9955    
## 3-1   646.02635  -653.6967 1945.74935 0.7920    
## 4-1   615.13124  -759.2744 1989.53690 0.8638    
## 5-1  -319.04560 -1327.5683  689.47714 0.9732    
## 6-1  -470.76852 -2140.4064 1198.86933 0.9852    
## 7-1  -882.50186 -2264.1538  499.15006 0.5228    
## 3-2   931.32499  -515.9634 2378.61336 0.5131    
## 4-2   900.42988  -614.2829 2415.14261 0.6125    
## 5-2   -33.74696 -1226.4254 1158.93148 1.0000    
## 6-2  -185.46989 -1972.3818 1601.44206 0.9999    
## 7-2  -597.20322 -2118.4940  924.08760 0.9237    
## 4-3   -30.89510 -1578.9689 1517.17867 1.0000    
## 5-3  -965.07195 -2199.8430  269.69913 0.2594    
## 6-3 -1116.79487 -2932.0723  698.48253 0.5706    
## 7-3 -1528.52821 -3083.0389   25.98249 0.0581 .  
## 5-4  -934.17684 -2247.3298  378.97610 0.3797    
## 6-4 -1085.89977 -2955.3765  783.57692 0.6396    
## 7-4 -1497.63310 -3115.1048  119.83855 0.0945 .  
## 6-5  -151.72292 -1771.3125 1467.86662 1.0000    
## 7-5  -563.45626 -1884.1916  757.27905 0.8892    
## 7-6  -411.73333 -2286.5438 1463.07711 0.9962    
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## # A tibble: 7 x 2
##   cluster OVERALLSAT19
##     <int>        <dbl>
## 1       1         2.56
## 2       2         2.52
## 3       3         2.40
## 4       4         2.64
## 5       5         2.56
## 6       6         2.46
## 7       7         2.65
##                     Df Sum Sq Mean Sq F value Pr(>F)  
## as.factor(cluster)   6   3.36  0.5599   1.873 0.0829 .
## Residuals          674 201.43  0.2989                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##  Pairwise comparisons using t tests with pooled SD 
## 
## data:  data$OVERALLSAT19 and as.factor(data$cluster) 
## 
##   1    2    3    4    5    6   
## 2 1.00 -    -    -    -    -   
## 3 0.76 1.00 -    -    -    -   
## 4 1.00 1.00 0.16 -    -    -   
## 5 1.00 1.00 0.60 1.00 -    -   
## 6 1.00 1.00 1.00 1.00 1.00 -   
## 7 1.00 1.00 0.15 1.00 1.00 1.00
## 
## P value adjustment method: bonferroni
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = OVERALLSAT19 ~ as.factor(cluster), data = data)
## 
## $`as.factor(cluster)`
##             diff         lwr        upr     p adj
## 2-1 -0.035364876 -0.25540387 0.18467411 0.9991418
## 3-1 -0.161184792 -0.38819637 0.06582679 0.3539207
## 4-1  0.083170355 -0.15565114 0.32199185 0.9470451
## 5-1 -0.001406759 -0.17755686 0.17474334 1.0000000
## 6-1 -0.097082228 -0.38870369 0.19453924 0.9573025
## 7-1  0.087533156 -0.15378825 0.32885456 0.9359401
## 3-2 -0.125819917 -0.37860549 0.12696565 0.7615065
## 4-2  0.118535231 -0.14490735 0.38197781 0.8374389
## 5-2  0.033958116 -0.17435690 0.24227313 0.9990695
## 6-2 -0.061717352 -0.37382209 0.25038738 0.9972346
## 7-2  0.122898032 -0.14281292 0.38860898 0.8188193
## 4-3  0.244355147 -0.02493854 0.51364884 0.1040264
## 5-3  0.159778033 -0.05588894 0.37544501 0.3017480
## 6-3  0.064102564 -0.25295652 0.38116165 0.9968719
## 7-3  0.248717949 -0.02279523 0.52023113 0.0975439
## 5-4 -0.084577114 -0.31264222 0.14348799 0.9290741
## 6-4 -0.180252583 -0.50587186 0.14536669 0.6582472
## 7-4  0.004362801 -0.27709920 0.28582480 1.0000000
## 6-5 -0.095675469 -0.37855542 0.18720448 0.9539230
## 7-5  0.088939916 -0.14174170 0.31962153 0.9154517
## 7-6  0.184615385 -0.14284183 0.51207260 0.6382300

## 
##   Posthoc multiple comparisons of means : Fisher LSD 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##             diff      lwr.ci      upr.ci   pval    
## 2-1 -0.035364876 -0.18146057  0.11073082 0.6347    
## 3-1 -0.161184792 -0.31190996 -0.01045962 0.0361 *  
## 4-1  0.083170355 -0.07539605  0.24173676 0.3034    
## 5-1 -0.001406759 -0.11836226  0.11554875 0.9812    
## 6-1 -0.097082228 -0.29070537  0.09654091 0.3252    
## 7-1  0.087533156 -0.07269307  0.24775939 0.2838    
## 3-2 -0.125819917 -0.29365782  0.04201799 0.1415    
## 4-2  0.118535231 -0.05637843  0.29344889 0.1838    
## 5-2  0.033958116 -0.10435340  0.17226963 0.6299    
## 6-2 -0.061717352 -0.26894043  0.14550573 0.5589    
## 7-2  0.122898032 -0.05352172  0.29931779 0.1718    
## 4-3  0.244355147  0.06555662  0.42315368 0.0075 ** 
## 5-3  0.159778033  0.01658515  0.30297091 0.0288 *  
## 6-3  0.064102564 -0.14640998  0.27461511 0.5501    
## 7-3  0.248717949  0.06844578  0.42899012 0.0069 ** 
## 5-4 -0.084577114 -0.23600178  0.06684755 0.2732    
## 6-4 -0.180252583 -0.39644869  0.03594353 0.1021    
## 7-4  0.004362801 -0.18251492  0.19124052 0.9635    
## 6-5 -0.095675469 -0.28349465  0.09214371 0.3176    
## 7-5  0.088939916 -0.06422199  0.24210182 0.2546    
## 7-6  0.184615385 -0.03280103  0.40203180 0.0959 .  
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##   Posthoc multiple comparisons of means : Tukey HSD 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##             diff      lwr.ci     upr.ci   pval    
## 2-1 -0.035364876 -0.25540387 0.18467411 0.9991    
## 3-1 -0.161184792 -0.38819637 0.06582679 0.3539    
## 4-1  0.083170355 -0.15565114 0.32199185 0.9470    
## 5-1 -0.001406759 -0.17755686 0.17474334 1.0000    
## 6-1 -0.097082228 -0.38870369 0.19453924 0.9573    
## 7-1  0.087533156 -0.15378825 0.32885456 0.9359    
## 3-2 -0.125819917 -0.37860549 0.12696565 0.7615    
## 4-2  0.118535231 -0.14490735 0.38197781 0.8374    
## 5-2  0.033958116 -0.17435690 0.24227313 0.9991    
## 6-2 -0.061717352 -0.37382209 0.25038738 0.9972    
## 7-2  0.122898032 -0.14281292 0.38860898 0.8188    
## 4-3  0.244355147 -0.02493854 0.51364884 0.1040    
## 5-3  0.159778033 -0.05588894 0.37544501 0.3017    
## 6-3  0.064102564 -0.25295652 0.38116165 0.9969    
## 7-3  0.248717949 -0.02279523 0.52023113 0.0975 .  
## 5-4 -0.084577114 -0.31264222 0.14348799 0.9291    
## 6-4 -0.180252583 -0.50587186 0.14536669 0.6582    
## 7-4  0.004362801 -0.27709920 0.28582480 1.0000    
## 6-5 -0.095675469 -0.37855542 0.18720448 0.9539    
## 7-5  0.088939916 -0.14174170 0.31962153 0.9155    
## 7-6  0.184615385 -0.14284183 0.51207260 0.6382    
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
##   Posthoc multiple comparisons of means : Scheffe Test 
##     95% family-wise confidence level
## 
## $`as.factor(cluster)`
##             diff      lwr.ci    upr.ci   pval    
## 2-1 -0.035364876 -0.30023485 0.2295051 0.9998    
## 3-1 -0.161184792 -0.43444795 0.1120784 0.6217    
## 4-1  0.083170355 -0.20430889 0.3706496 0.9831    
## 5-1 -0.001406759 -0.21344587 0.2106324 1.0000    
## 6-1 -0.097082228 -0.44811896 0.2539545 0.9867    
## 7-1  0.087533156 -0.20295533 0.3780216 0.9791    
## 3-2 -0.125819917 -0.43010829 0.1784685 0.9035    
## 4-2  0.118535231 -0.19858142 0.4356519 0.9393    
## 5-2  0.033958116 -0.21679922 0.2847155 0.9998    
## 6-2 -0.061717352 -0.43741064 0.3139759 0.9993    
## 7-2  0.122898032 -0.19694915 0.4427452 0.9309    
## 4-3  0.244355147 -0.07980473 0.5685150 0.3042    
## 5-3  0.159778033 -0.09982917 0.4193852 0.5701    
## 6-3  0.064102564 -0.31755448 0.4457596 0.9992    
## 7-3  0.248717949 -0.07811362 0.5755495 0.2922    
## 5-4 -0.084577114 -0.35910845 0.1899542 0.9766    
## 6-4 -0.180252583 -0.57221388 0.2117087 0.8475    
## 7-4  0.004362801 -0.33444457 0.3431702 1.0000    
## 6-5 -0.095675469 -0.43618969 0.2448387 0.9855    
## 7-5  0.088939916 -0.18874102 0.3666209 0.9715    
## 7-6  0.184615385 -0.20955832 0.5787891 0.8356    
## 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1