search

tmatta / lsasim

Simulate large scale assessment data
6 stars 5 forks source link

mean & sigma #16

Closed wleoncio closed 3 years ago

wleoncio commented 3 years ago

0. Setup

I've tested most values below. Not all testings are shown in this report. I only included the testings that are showing errors/warnings or inconsistent results.

cluster_gen_2 <- function(...) {
  cluster_gen(..., verbose = FALSE, calc_weights = FALSE)
}
set.seed(12334)
n1 <- c(3, 6)
n2 <- c(groups = 4, people = 2)
n3 <- c(school = 3, class = 2, student = 5)
n4 <- c(20, 50)
n5 <- list(school = 3, class = c(2, 1, 3), student = c(20, 20, 10, 30, 30, 30))
n5a <- list(school = 3, class = c(2, 3, 3), student = c(20, 20, 10, 30, 30, 30))
n6 <- list(school = 3, class = c(2, 1, 3), student = ranges(10, 50))
n6a <- list(school = 3, class = c(2, 3, 3), student = ranges(10, 50))
n7 <- list(school = 10, student = ranges(10, 50))
n8 <- list(school = 3, student = c(20, 20, 10))
n8a <- list(school = 3, class = c(2, 2, 2),student = c(20, 20, 10))
n8b <- list(school = 3, class = c(2, 3, 3),student = c(20, 20, 10, 5))
n8c <- list(school = 3, class = c(2, 1, 3),student = c(20, 20, 10))
n9 <- list(school = 10, class = c(2,1,3,1,1,1,2,1,2,1), student = ranges(10, 50))
n10 <- list(country = 2, school = 10, class = c(2,1,3,1,1,1,2,1,2,1), student = ranges(10, 50))
n11 <- list(culture = 2, country = 2, school = 10, class = c(2,1,3,1,1,1,2,1,2,1), student = ranges(10, 50))
n12 <- list(culture = 2, country = 2, district = 3, school = 10, class = c(2,1,3,1,1,1,2,1,2,1), student = ranges(10, 50))
N1 <- c(100, 20)

6. mean & sigma

set.seed(12334)

r1 <- cluster_gen_2(n4, n_X = 5, c_mean = c(0.1, 0.5, 0.001, 234, 701), sigma = c(0.111, 0.113, 0.115, 0.117, 0.119))
r2 <- cluster_gen_2(n4, n_X = 5, c_mean = c(10, 55, 0.21, 2.34, 5000), sigma = c(40, 100, 0.11, 3, 1500))
r3 <- cluster_gen_2(n4, n_X = 5, c_mean = c(10, 55, 0.21, 2.34, 5000), sigma = c(40, 100, 0.11, 3, 1500))
r4 <- cluster_gen_2(n4, n_X = 5, c_mean = c(0.1, 0.5, 0.001, 234, 701), sigma = c(0.111, 0.113, 0.115, 0.117, 0.119))
r5 <- cluster_gen_2(n6, n_X = 3, c_mean = c(0.3, 0.35, 0.4), sigma = c(0.11, 0.22, 0.33))
r6 <- cluster_gen_2(n6, n_X = 3, c_mean = c(0.7, 100, 355), sigma = c(0.2, 21, 0.7))
r7 <- cluster_gen_2(n10, n_X = 4, c_mean = c(1, 20, 0.25, 50.54), sigma = c(0.3, 3.5, 0.1, 3))
r8 <- cluster_gen_2(n10, n_X = 4, c_mean = c(0.001, 0.005, 213, 234), sigma = c(0.0001, 0.001, 11, 6))

summarize_clusters(r1) #mean q3: -0.004566. sigma: 0.12 0.11  0.12  0.12  0.12
## ───────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## Summary statistics for all schools
##        q1                q2               q3                  q4              q5        q6      q7      q8
##  Min.   :-0.2690   Min.   :0.1565   Min.   :-0.427666   Min.   :233.6   Min.   :700.6   1:312   1:281   1:276
##  Mean   : 0.1003   Mean   :0.5021   Mean   :-0.004566   Mean   :234.0   Mean   :701.0   2:204   2:210   2:214
##  Max.   : 0.4536   Max.   :0.9000   Max.   : 0.387608   Max.   :234.4   Max.   :701.4   3:213   3:178   3:200
##                                                                                         4:271   4:331   4:310
##  Stddev.: 0.12     Stddev.: 0.11    Stddev.: 0.12       Stddev.: 0.12   Stddev.: 0.12
##                                                                                         Prop.   Prop.   Prop.
##                                                                                         1:0.312 1:0.281 1:0.276
##                                                                                         2:0.204 2:0.21  2:0.214
##                                                                                         3:0.213 3:0.178 3:0.2
##                                                                                         4:0.271 4:0.331 4:0.31
##
##
##  q9      q10
##  1:533   1:287
##  2:467   2:169
##          3:222
##  Prop.   4:322
##  1:0.533
##  2:0.467 Prop.
##          1:0.287
##          2:0.169
##          3:0.222
##          4:0.322
##
##
##
##  Heterogeneous correlation matrix
##              q1           q2           q3           q4           q5            q6            q7           q8
## q1   1.00000000 -0.033709555 -0.057207893 -0.050285770 -0.007572790 -0.0487243064 -0.0588806936 -0.093471415
## q2  -0.03370955  1.000000000 -0.005242536  0.060533401  0.041963884 -0.0378406336  0.0107370558  0.169881756
## q3  -0.05720789 -0.005242536  1.000000000  0.059446888  0.036217118  0.0358582554  0.0721081711  0.106717966
## q4  -0.05028577  0.060533401  0.059446888  1.000000000  0.042094746 -0.0550764080  0.0070833221  0.096670634
## q5  -0.00757279  0.041963884  0.036217118  0.042094746  1.000000000 -0.0356164071  0.0188430981  0.014768568
## q6  -0.04872431 -0.037840634  0.035858255 -0.055076408 -0.035616407  1.0000000000  0.0002080099  0.004573422
## q7  -0.05888069  0.010737056  0.072108171  0.007083322  0.018843098  0.0002080099  1.0000000000  0.007002481
## q8  -0.09347141  0.169881756  0.106717966  0.096670634  0.014768568  0.0045734219  0.0070024807  1.000000000
## q9  -0.05046500 -0.022076031  0.122363471  0.020842870 -0.057918876  0.1207264938 -0.0585946370  0.102631478
## q10 -0.02233170 -0.077004252 -0.049903458 -0.047753203 -0.006770408  0.1430814902  0.0416718087 -0.006270074
##              q9          q10
## q1  -0.05046500 -0.022331696
## q2  -0.02207603 -0.077004252
## q3   0.12236347 -0.049903458
## q4   0.02084287 -0.047753203
## q5  -0.05791888 -0.006770408
## q6   0.12072649  0.143081490
## q7  -0.05859464  0.041671809
## q8   0.10263148 -0.006270074
## q9   1.00000000  0.116081191
## q10  0.11608119  1.000000000
summarize_clusters(r2) #mean:  9.332   61.241   0.2120   2.3793   5021.3
## ───────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## Summary statistics for all schools
##        q1                 q2                 q3                q4                q5         q6      q7
##  Min.   :-127.679   Min.   :-227.542   Min.   :-0.1362   Min.   :-6.3438   Min.   : 191.1   1:241   1:299
##  Mean   :   9.332   Mean   :  61.241   Mean   : 0.2120   Mean   : 2.3793   Mean   :5021.3   2:163   2:234
##  Max.   : 133.724   Max.   : 384.065   Max.   : 0.4929   Max.   :10.6954   Max.   :9768.7   3:151   3:199
##                                                                                             4:152   4:268
##  Stddev.: 40.17     Stddev.: 98.13     Stddev.: 0.11     Stddev.: 3.06     Stddev.: 1531.13 5:293
##                                                                                                     Prop.
##                                                                                             Prop.   1:0.299
##                                                                                             1:0.241 2:0.234
##                                                                                             2:0.163 3:0.199
##                                                                                             3:0.151 4:0.268
##                                                                                             4:0.152
##                                                                                             5:0.293
##
##
##
##  Heterogeneous correlation matrix
##             q1          q2           q3          q4          q5           q6          q7
## q1  1.00000000  0.02685254 -0.130647344 -0.21441788 -0.03642083  0.036481204  0.12242260
## q2  0.02685254  1.00000000 -0.110373373 -0.05152180  0.07547174  0.024813427 -0.17689516
## q3 -0.13064734 -0.11037337  1.000000000 -0.03188441 -0.12737462  0.004195498 -0.03845626
## q4 -0.21441788 -0.05152180 -0.031884407  1.00000000 -0.01710220  0.119152277 -0.11637760
## q5 -0.03642083  0.07547174 -0.127374616 -0.01710220  1.00000000  0.077255918 -0.05225078
## q6  0.03648120  0.02481343  0.004195498  0.11915228  0.07725592  1.000000000 -0.11011212
## q7  0.12242260 -0.17689516 -0.038456262 -0.11637760 -0.05225078 -0.110112117  1.00000000
                       #sigma: 40.17 98.13   0.11    3.06    1531.13

summarize_clusters(r3) #mean: 7.109 61.404  0.20981   2.2936 4931.7
## ───────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## Summary statistics for all schools
##        q1                 q2                 q3                 q4                q5          q6      q7
##  Min.   :-114.400   Min.   :-321.481   Min.   :-0.09432   Min.   :-7.1558   Min.   : -278.7   1:320   1:316
##  Mean   :   7.109   Mean   :  61.404   Mean   : 0.20981   Mean   : 2.2936   Mean   : 4931.7   2:180   2:310
##  Max.   : 116.017   Max.   : 344.757   Max.   : 0.56660   Max.   :12.4447   Max.   :10056.8   3:193   3:374
##                                                                                               4:307
##  Stddev.: 37.27     Stddev.: 98.8      Stddev.: 0.11      Stddev.: 2.95     Stddev.: 1486.12          Prop.
##                                                                                               Prop.   1:0.316
##                                                                                               1:0.32  2:0.31
##                                                                                               2:0.18  3:0.374
##                                                                                               3:0.193
##                                                                                               4:0.307
##
##
##
##
##  q8      q9      q10     q11     q12
##  1:231   1:298   1:392   1:294   1:292
##  2:164   2:287   2:237   2:265   2:165
##  3:128   3:415   3:371   3:441   3:177
##  4:198                           4:366
##  5:279   Prop.   Prop.   Prop.
##          1:0.298 1:0.392 1:0.294 Prop.
##  Prop.   2:0.287 2:0.237 2:0.265 1:0.292
##  1:0.231 3:0.415 3:0.371 3:0.441 2:0.165
##  2:0.164                         3:0.177
##  3:0.128                         4:0.366
##  4:0.198
##  5:0.279
##
##
##
##  Heterogeneous correlation matrix
## Warning in log(P): NaNs produced

## Warning in log(P): NaNs produced

## Warning in log(P): NaNs produced
##               q1          q2            q3           q4          q5            q6          q7           q8
## q1   1.000000000 -0.08864636  0.0218341464 -0.058905622  0.06486101 -0.0289198074 -0.08480181  0.009046702
## q2  -0.088646361  1.00000000  0.0765745493  0.052046387 -0.14230246 -0.1015130659  0.03734397 -0.067643472
## q3   0.021834146  0.07657455  1.0000000000 -0.072484993  0.05144694  0.0006097222  0.10143102  0.024166654
## q4  -0.058905622  0.05204639 -0.0724849934  1.000000000 -0.01652480 -0.0178596621  0.04324850  0.105479602
## q5   0.064861006 -0.14230246  0.0514469412 -0.016524804  1.00000000  0.1328678115  0.07814149  0.078148404
## q6  -0.028919807 -0.10151307  0.0006097222 -0.017859662  0.13286781  1.0000000000 -0.02145899 -0.103971288
## q7  -0.084801805  0.03734397  0.1014310232  0.043248501  0.07814149 -0.0214589879  1.00000000  0.043683957
## q8   0.009046702 -0.06764347  0.0241666538  0.105479602  0.07814840 -0.1039712878  0.04368396  1.000000000
## q9  -0.109399902 -0.04715158 -0.0663456174  0.084293038 -0.13180438  0.0004820534  0.03421766 -0.045339385
## q10  0.038166021 -0.06920142 -0.0480908181  0.002481838  0.05997579  0.0248628631 -0.02228583  0.089648392
## q11 -0.066316089  0.07208785  0.0892913631 -0.064050259  0.13313303  0.1555895444  0.07193066 -0.066730059
## q12  0.154679626 -0.09820972  0.0753307604  0.010423188  0.03983716 -0.1209274603  0.01037961 -0.103200410
##                q9          q10         q11          q12
## q1  -0.1093999021  0.038166021 -0.06631609  0.154679626
## q2  -0.0471515816 -0.069201419  0.07208785 -0.098209721
## q3  -0.0663456174 -0.048090818  0.08929136  0.075330760
## q4   0.0842930376  0.002481838 -0.06405026  0.010423188
## q5  -0.1318043772  0.059975789  0.13313303  0.039837160
## q6   0.0004820534  0.024862863  0.15558954 -0.120927460
## q7   0.0342176631 -0.022285826  0.07193066  0.010379609
## q8  -0.0453393849  0.089648392 -0.06673006 -0.103200410
## q9   1.0000000000 -0.140481162  0.08009391 -0.003382839
## q10 -0.1404811625  1.000000000  0.01950217 -0.005711982
## q11  0.0800939104  0.019502170  1.00000000  0.098195969
## q12 -0.0033828386 -0.005711982  0.09819597  1.000000000
                       #sigma:37.27     98.8  0.11   2.95   1486.12

summarize_clusters(r4) #q3 mean: 0.002425
## ───────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## Summary statistics for all schools
##        q1                 q2               q3                  q4              q5        q6      q7      q8
##  Min.   :-0.27537   Min.   :0.1656   Min.   :-0.374855   Min.   :233.6   Min.   :700.7   1:268   1:511   1:305
##  Mean   : 0.10507   Mean   :0.5028   Mean   : 0.002425   Mean   :234.0   Mean   :701.0   2:162   2:489   2:199
##  Max.   : 0.41851   Max.   :0.8674   Max.   : 0.342597   Max.   :234.3   Max.   :701.3   3:155           3:183
##                                                                                          4:166   Prop.   4:313
##  Stddev.: 0.11      Stddev.: 0.11    Stddev.: 0.11       Stddev.: 0.11   Stddev.: 0.12   5:249   1:0.511
##                                                                                                  2:0.489 Prop.
##                                                                                          Prop.           1:0.305
##                                                                                          1:0.268         2:0.199
##                                                                                          2:0.162         3:0.183
##                                                                                          3:0.155         4:0.313
##                                                                                          4:0.166
##                                                                                          5:0.249
##
##
##  q9      q10     q11     q12
##  1:257   1:290   1:265   1:258
##  2:196   2:189   2:172   2:216
##  3:219   3:197   3:165   3:208
##  4:328   4:324   4:169   4:318
##                  5:229
##  Prop.   Prop.           Prop.
##  1:0.257 1:0.29  Prop.   1:0.258
##  2:0.196 2:0.189 1:0.265 2:0.216
##  3:0.219 3:0.197 2:0.172 3:0.208
##  4:0.328 4:0.324 3:0.165 4:0.318
##                  4:0.169
##                  5:0.229
##
##
##
##  Heterogeneous correlation matrix
##               q1           q2          q3          q4           q5          q6           q7           q8
## q1   1.000000000 -0.019719497  0.04466830 -0.03741994  0.022338641 -0.02845867 -0.109125008  0.004814079
## q2  -0.019719497  1.000000000 -0.03269803  0.01423826  0.007502822  0.01929364  0.090217716  0.129328507
## q3   0.044668304 -0.032698029  1.00000000 -0.10895263 -0.061704820 -0.12422195 -0.115374051 -0.052316370
## q4  -0.037419939  0.014238263 -0.10895263  1.00000000 -0.073782195  0.04064608  0.084671220  0.082966734
## q5   0.022338641  0.007502822 -0.06170482 -0.07378220  1.000000000 -0.02064799 -0.002729842  0.026669219
## q6  -0.028458665  0.019293637 -0.12422195  0.04064608 -0.020647994  1.00000000  0.134760711 -0.017948907
## q7  -0.109125008  0.090217716 -0.11537405  0.08467122 -0.002729842  0.13476071  1.000000000 -0.007172670
## q8   0.004814079  0.129328507 -0.05231637  0.08296673  0.026669219 -0.01794891 -0.007172670  1.000000000
## q9  -0.061023239  0.031632846 -0.05504252 -0.05950947 -0.066292607  0.13449174  0.060820982  0.094248221
## q10 -0.052132032 -0.028416020  0.08242385  0.01533040 -0.042201133  0.07216644 -0.073334408 -0.073337215
## q11  0.017282983  0.043484013 -0.07078090  0.09062775  0.003010196 -0.01070286 -0.029458268  0.014836451
## q12  0.060922375  0.016960325 -0.00651158 -0.12534302  0.168403801 -0.01811740  0.004779470 -0.020126254
##              q9         q10          q11         q12
## q1  -0.06102324 -0.05213203  0.017282983  0.06092238
## q2   0.03163285 -0.02841602  0.043484013  0.01696032
## q3  -0.05504252  0.08242385 -0.070780902 -0.00651158
## q4  -0.05950947  0.01533040  0.090627753 -0.12534302
## q5  -0.06629261 -0.04220113  0.003010196  0.16840380
## q6   0.13449174  0.07216644 -0.010702863 -0.01811740
## q7   0.06082098 -0.07333441 -0.029458268  0.00477947
## q8   0.09424822 -0.07333721  0.014836451 -0.02012625
## q9   1.00000000 -0.05393121 -0.064194877 -0.18618997
## q10 -0.05393121  1.00000000 -0.024018830 -0.16256718
## q11 -0.06419488 -0.02401883  1.000000000  0.11830563
## q12 -0.18618997 -0.16256718  0.118305633  1.00000000
summarize_clusters(r5) #level1 mean: 0.29996   0.3262 0.4013
## ───────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## Summary statistics for all schools
##        q1               q2               q3          q4      q5      q6
##  Min.   :0.1507   Min.   :0.0673   Min.   :-0.1348   1:2     1:3     2:2
##  Mean   :0.2953   Mean   :0.3793   Mean   : 0.4047   3:2     3:2     4:3
##  Max.   :0.4669   Max.   :0.6312   Max.   : 0.9826   2:2     2:1     1:1
##
##  Stddev.: 0.13    Stddev.: 0.25    Stddev.: 0.41     Prop.   Prop.   Prop.
##                                                      1:0.333 1:0.5   2:0.333
##                                                      3:0.333 3:0.333 4:0.5
##                                                      2:0.333 2:0.167 1:0.167
##
##
##
##  Heterogeneous correlation matrix
## Warning in hetcor.data.frame(df): the correlation matrix has been adjusted to make it positive-definite
##             q1         q2         q3         q4          q5          q6
## q1  1.00000000 -0.2793414 -0.5500657  0.1168326  0.04093963 -0.01335278
## q2 -0.27934136  1.0000000  0.5101079 -0.1901446 -0.15447229 -0.60722141
## q3 -0.55006567  0.5101079  1.0000000 -0.3372862 -0.49335512  0.18728918
## q4  0.11683255 -0.1901446 -0.3372862  1.0000000 -0.57679236  0.34395316
## q5  0.04093963 -0.1544723 -0.4933551 -0.5767924  1.00000000 -0.58266650
## q6 -0.01335278 -0.6072214  0.1872892  0.3439532 -0.58266650  1.00000000
## Summary statistics for all classes
##        q1                 q2                q3          q4      q5      q6      q7      q8      q9      q10
##  Min.   :-0.04591   Min.   :-0.1006   Min.   :-0.4625   1:48    1:38    1:36    1:47    1:84    1:55    1:47
##  Mean   : 0.29996   Mean   : 0.3262   Mean   : 0.4013   2:85    2:38    2:28    2:19    2:49    3:23    2:44
##  Max.   : 0.53427   Max.   : 0.9026   Max.   : 1.2129           4:35    3:29    3:13            4:37    3:42
##                                                         Prop.   3:22    4:40    4:16    Prop.   2:18
##  Stddev.: 0.1       Stddev.: 0.23     Stddev.: 0.34     1:0.361                 5:38    1:0.632         Prop.
##                                                         2:0.639 Prop.   Prop.           2:0.368 Prop.   1:0.353
##                                                                 1:0.286 1:0.271 Prop.           1:0.414 2:0.331
##                                                                 2:0.286 2:0.211 1:0.353         3:0.173 3:0.316
##                                                                 4:0.263 3:0.218 2:0.143         4:0.278
##                                                                 3:0.165 4:0.301 3:0.098         2:0.135
##                                                                                 4:0.12
##                                                                                 5:0.286
##
##
##
##  Heterogeneous correlation matrix
##              q1          q2           q3          q4          q5           q6          q7          q8          q9
## q1   1.00000000  0.04780571  0.128963825 -0.01825320  0.14623416  0.203421352  0.10230697 -0.16528312 -0.04956281
## q2   0.04780571  1.00000000  0.077266212 -0.05374943  0.13480197  0.086272477 -0.09625784 -0.14374188  0.10943339
## q3   0.12896383  0.07726621  1.000000000 -0.20827985  0.12076343 -0.076654839 -0.07297245  0.14233118  0.05226638
## q4  -0.01825320 -0.05374943 -0.208279850  1.00000000  0.04799061  0.216402440 -0.09012485  0.19752335  0.03426720
## q5   0.14623416  0.13480197  0.120763432  0.04799061  1.00000000  0.086957174 -0.34533010 -0.02019908  0.02501169
## q6   0.20342135  0.08627248 -0.076654839  0.21640244  0.08695717  1.000000000 -0.17824535  0.05263700 -0.37243935
## q7   0.10230697 -0.09625784 -0.072972450 -0.09012485 -0.34533010 -0.178245346  1.00000000  0.10395805  0.08307646
## q8  -0.16528312 -0.14374188  0.142331180  0.19752335 -0.02019908  0.052637000  0.10395805  1.00000000 -0.21863038
## q9  -0.04956281  0.10943339  0.052266381  0.03426720  0.02501169 -0.372439353  0.08307646 -0.21863038  1.00000000
## q10 -0.39066739 -0.30383465  0.004060791  0.11463426 -0.07641638  0.001061375 -0.16283615  0.07634474  0.25058932
##              q10
## q1  -0.390667394
## q2  -0.303834654
## q3   0.004060791
## q4   0.114634255
## q5  -0.076416376
## q6   0.001061375
## q7  -0.162836152
## q8   0.076344739
## q9   0.250589320
## q10  1.000000000
summarize_clusters(r6)
## ───────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## Summary statistics for all schools
##        q1               q2               q3        q4      q5      q6
##  Min.   :0.5282   Min.   : 66.16   Min.   :354.0   1:5     1:2     3:5
##  Mean   :0.7042   Mean   :104.93   Mean   :354.6   2:1     4:3     1:1
##  Max.   :0.9691   Max.   :138.14   Max.   :355.9           3:1
##                                                    Prop.           Prop.
##  Stddev.: 0.18    Stddev.: 25.32   Stddev.: 0.69   1:0.833 Prop.   3:0.833
##                                                    2:0.167 1:0.333 1:0.167
##                                                            4:0.5
##                                                            3:0.167
##
##
##
##  Heterogeneous correlation matrix
## Warning in polyserial(y, x, ML = ML, std.err = std.err, bins = bins): initial correlation inadmissible,
## -1.11894475763533, set to -0.9999
## Warning in polyserial(y, x, ML = ML, std.err = std.err, bins = bins): inadmissible correlation set to 0.9999
## Warning in polyserial(y, x, ML = ML, std.err = std.err, bins = bins): initial correlation inadmissible,
## 1.39394961212503, set to 0.9999
## Warning in hetcor.data.frame(df): the correlation matrix has been adjusted to make it positive-definite
##            q1         q2         q3         q4         q5         q6
## q1  1.0000000  0.7464256  0.7043458 -0.3526262 -0.3958877  0.7057745
## q2  0.7464256  1.0000000  0.8500206 -0.8421612  0.2280144  0.9585307
## q3  0.7043458  0.8500206  1.0000000 -0.7708158  0.1701872  0.9586432
## q4 -0.3526262 -0.8421612 -0.7708158  1.0000000 -0.7000367 -0.8553525
## q5 -0.3958877  0.2280144  0.1701872 -0.7000367  1.0000000  0.2482209
## q6  0.7057745  0.9585307  0.9586432 -0.8553525  0.2482209  1.0000000
## Summary statistics for all classes
##        q1               q2               q3        q4      q5     q6      q7
##  Min.   :0.3088   Min.   : 52.75   Min.   :353.6   1:42    1:51   1:40    1:39
##  Mean   :0.7080   Mean   : 99.14   Mean   :355.1   3:37    2:51   2:30    3:12
##  Max.   :1.2033   Max.   :151.86   Max.   :357.2   2:23           3:32    4:32
##                                                            Prop.          2:19
##  Stddev.: 0.19    Stddev.: 19.11   Stddev.: 0.69   Prop.   1:0.5  Prop.
##                                                    1:0.412 2:0.5  1:0.392 Prop.
##                                                    3:0.363        2:0.294 1:0.382
##                                                    2:0.225        3:0.314 3:0.118
##                                                                           4:0.314
##                                                                           2:0.186
##
##
##
##  Heterogeneous correlation matrix
##             q1          q2          q3          q4          q5          q6          q7
## q1  1.00000000  0.06804983  0.02884317  0.05488177 -0.16294487 -0.11188934 -0.09254774
## q2  0.06804983  1.00000000 -0.30234075 -0.13869510 -0.22826624  0.02920063  0.15485272
## q3  0.02884317 -0.30234075  1.00000000  0.04230270 -0.04015550  0.12322271 -0.13459225
## q4  0.05488177 -0.13869510  0.04230270  1.00000000 -0.01691490 -0.04635040  0.17748941
## q5 -0.16294487 -0.22826624 -0.04015550 -0.01691490  1.00000000  0.09716052 -0.27616800
## q6 -0.11188934  0.02920063  0.12322271 -0.04635040  0.09716052  1.00000000  0.08689299
## q7 -0.09254774  0.15485272 -0.13459225  0.17748941 -0.27616800  0.08689299  1.00000000
summarize_clusters(r7) #sigma: l0.31    3.54    0.1  3
## ───────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## Summary statistics for all countries
##        q1               q2              q3                 q4        q5     q6     q7     q8     q9
##  Min.   :0.3102   Min.   :15.41   Min.   :0.007886   Min.   :44.70   1: 5   1:9    1:9    1:3    1:5
##  Mean   :0.9938   Mean   :21.01   Mean   :0.223429   Mean   :51.21   3:13   2:3    2:4    2:8    2:9
##  Max.   :1.4725   Max.   :29.54   Max.   :0.346146   Max.   :57.02   2: 2   3:5    3:7    3:2    3:6
##                                                                             4:3           4:7
##  Stddev.: 0.32    Stddev.: 4.37   Stddev.: 0.08      Stddev.: 2.91   Prop.         Prop.         Prop.
##                                                                      1:0.25 Prop.  1:0.45 Prop.  1:0.25
##                                                                      3:0.65 1:0.45 2:0.2  1:0.15 2:0.45
##                                                                      2:0.1  2:0.15 3:0.35 2:0.4  3:0.3
##                                                                             3:0.25        3:0.1
##                                                                             4:0.15        4:0.35
##
##
##
##  Heterogeneous correlation matrix
## Warning in log(P): NaNs produced
## Warning in log(P): NaNs produced

## Warning in log(P): NaNs produced

## Warning in log(P): NaNs produced

## Warning in log(P): NaNs produced

## Warning in log(P): NaNs produced
## Warning in hetcor.data.frame(df): the correlation matrix has been adjusted to make it positive-definite
##              q1         q2          q3          q4           q5          q6         q7          q8          q9
## q1  1.000000000 -0.7759417  0.40464045  0.42614527 -0.001097924 -0.04057159 -0.6087368 -0.16223638 -0.31238365
## q2 -0.775941653  1.0000000 -0.68723066 -0.34758159 -0.158353860  0.21168268  0.6980966  0.37004914  0.58529503
## q3  0.404640447 -0.6872307  1.00000000  0.32269732  0.099193353 -0.22992829 -0.4307596 -0.51270019 -0.43214067
## q4  0.426145273 -0.3475816  0.32269732  1.00000000 -0.136343401 -0.13551096 -0.1815548 -0.03910875 -0.07596138
## q5 -0.001097924 -0.1583539  0.09919335 -0.13634340  1.000000000  0.34879554 -0.2971719  0.56970924 -0.12571132
## q6 -0.040571590  0.2116827 -0.22992829 -0.13551096  0.348795541  1.00000000  0.3408737  0.79427747  0.39984354
## q7 -0.608736823  0.6980966 -0.43075962 -0.18155481 -0.297171910  0.34087369  1.0000000  0.26627561  0.39432577
## q8 -0.162236384  0.3700491 -0.51270019 -0.03910875  0.569709240  0.79427747  0.2662756  1.00000000  0.52844732
## q9 -0.312383655  0.5852950 -0.43214067 -0.07596138 -0.125711318  0.39984354  0.3943258  0.52844732  1.00000000
## Summary statistics for all schools
##        q1               q2              q3               q4        q5      q6      q7      q8
##  Min.   :0.4793   Min.   :10.88   Min.   :0.1051   Min.   :46.04   1:13    1:11    2: 9    5: 8
##  Mean   :0.9669   Mean   :19.47   Mean   :0.2383   Mean   :51.35   3:11    4: 7    3:10    1:10
##  Max.   :1.6280   Max.   :25.64   Max.   :0.3865   Max.   :56.05   2: 6    3: 9    1:11    2: 4
##                                                                            2: 3            3: 3
##  Stddev.: 0.28    Stddev.: 3.53   Stddev.: 0.09    Stddev.: 2.53   Prop.           Prop.   4: 5
##                                                                    1:0.433 Prop.   2:0.3
##                                                                    3:0.367 1:0.367 3:0.333 Prop.
##                                                                    2:0.2   4:0.233 1:0.367 5:0.267
##                                                                            3:0.3           1:0.333
##                                                                            2:0.1           2:0.133
##                                                                                            3:0.1
##                                                                                            4:0.167
##
##
##
##  Heterogeneous correlation matrix
##             q1          q2          q3          q4          q5          q6          q7          q8
## q1  1.00000000  0.29816944 -0.44009270  0.27961929  0.45086436  0.07307666  0.14646562 -0.06457049
## q2  0.29816944  1.00000000 -0.26727752 -0.25904855 -0.02146690  0.15757947  0.04409052 -0.07729463
## q3 -0.44009270 -0.26727752  1.00000000  0.04628507  0.03303019 -0.18474098 -0.16280506  0.19042179
## q4  0.27961929 -0.25904855  0.04628507  1.00000000 -0.02030120 -0.09730632  0.11533347 -0.08043289
## q5  0.45086436 -0.02146690  0.03303019 -0.02030120  1.00000000  0.41141681  0.11413634  0.38156363
## q6  0.07307666  0.15757947 -0.18474098 -0.09730632  0.41141681  1.00000000  0.08722410  0.11664850
## q7  0.14646562  0.04409052 -0.16280506  0.11533347  0.11413634  0.08722410  1.00000000  0.08148543
## q8 -0.06457049 -0.07729463  0.19042179 -0.08043289  0.38156363  0.11664850  0.08148543  1.00000000
## Summary statistics for all classes
##        q1               q2               q3                 q4        q5      q6      q7
##  Min.   :0.1396   Min.   : 9.202   Min.   :-0.05129   Min.   :39.80   1:197   1:234   1:211
##  Mean   :0.9669   Mean   :20.080   Mean   : 0.24889   Mean   :50.27   2:131   2:177   2:209
##  Max.   :1.8995   Max.   :30.233   Max.   : 0.53463   Max.   :60.40   3:137   3:246   3:237
##                                                                       4:192
##  Stddev.: 0.31    Stddev.: 3.54    Stddev.: 0.1       Stddev.: 3              Prop.   Prop.
##                                                                       Prop.   1:0.356 1:0.321
##                                                                       1:0.3   2:0.269 2:0.318
##                                                                       2:0.199 3:0.374 3:0.361
##                                                                       3:0.209
##                                                                       4:0.292
##
##
##
##  Heterogeneous correlation matrix
##             q1           q2          q3           q4           q5           q6          q7
## q1  1.00000000  0.122069526 0.104612626 -0.069616594  0.034610049 -0.015444333  0.02639731
## q2  0.12206953  1.000000000 0.006427058 -0.002347094 -0.036437619 -0.150208033  0.07851958
## q3  0.10461263  0.006427058 1.000000000  0.126636062  0.141740425  0.072373650  0.12136895
## q4 -0.06961659 -0.002347094 0.126636062  1.000000000  0.114398690 -0.048082905  0.08368643
## q5  0.03461005 -0.036437619 0.141740425  0.114398690  1.000000000  0.001140093 -0.09704255
## q6 -0.01544433 -0.150208033 0.072373650 -0.048082905  0.001140093  1.000000000 -0.07485062
## q7  0.02639731  0.078519580 0.121368951  0.083686434 -0.097042547 -0.074850623  1.00000000
summarize_clusters(r8) #sigma: q1: 0 q2: 0
## ───────────────────────────────────────────────────────────────────────────────────────────────────────────────────
## Summary statistics for all countries
##        q1                  q2                 q3              q4        q5     q6     q7     q8     q9     q10
##  Min.   :0.0008541   Min.   :0.002601   Min.   :191.1   Min.   :228.8   1:8    1:8    1:6    1:7    1:5    1: 9
##  Mean   :0.0010028   Mean   :0.004634   Mean   :213.0   Mean   :237.6   2:5    3:3    2:9    2:5    2:2    2:11
##  Max.   :0.0011256   Max.   :0.007706   Max.   :245.1   Max.   :246.1   3:3    4:1    3:5    3:8    3:7
##                                                                         4:4    5:5                  4:6    Prop.
##  Stddev.: 0          Stddev.: 0         Stddev.: 14.07  Stddev.: 5.08          2:3    Prop.  Prop.         1:0.45
##                                                                         Prop.         1:0.3  1:0.35 Prop.  2:0.55
##                                                                         1:0.4  Prop.  2:0.45 2:0.25 1:0.25
##                                                                         2:0.25 1:0.4  3:0.25 3:0.4  2:0.1
##                                                                         3:0.15 3:0.15               3:0.35
##                                                                         4:0.2  4:0.05               4:0.3
##                                                                                5:0.25
##                                                                                2:0.15
##
##
##
##  Heterogeneous correlation matrix
## Warning in log(P): NaNs produced

## Warning in log(P): the correlation matrix has been adjusted to make it positive-definite
##               q1           q2          q3          q4          q5          q6          q7          q8          q9
## q1   1.000000000  0.009780447  0.34088398  0.35507484 -0.09988053 -0.24448188 -0.50019599 -0.45758029 -0.21854124
## q2   0.009780447  1.000000000  0.11576567 -0.01320661  0.42763465 -0.07666214 -0.57366725 -0.08764248 -0.06711720
## q3   0.340883979  0.115765675  1.00000000  0.59303496 -0.38966057 -0.35424101 -0.07493377 -0.32610054 -0.28937639
## q4   0.355074839 -0.013206607  0.59303496  1.00000000 -0.29752435  0.13384336 -0.10769214 -0.02226312 -0.15774298
## q5  -0.099880532  0.427634652 -0.38966057 -0.29752435  1.00000000  0.16249075 -0.25009682  0.47260327 -0.30371971
## q6  -0.244481876 -0.076662143 -0.35424101  0.13384336  0.16249075  1.00000000  0.07256728  0.65354742  0.47304489
## q7  -0.500195994 -0.573667250 -0.07493377 -0.10769214 -0.25009682  0.07256728  1.00000000  0.20207761  0.19728126
## q8  -0.457580295 -0.087642475 -0.32610054 -0.02226312  0.47260327  0.65354742  0.20207761  1.00000000 -0.09657972
## q9  -0.218541244 -0.067117204 -0.28937639 -0.15774298 -0.30371971  0.47304489  0.19728126 -0.09657972  1.00000000
## q10 -0.103765493 -0.621912692  0.05476263 -0.10003128 -0.50936500  0.37284708  0.46584952  0.23963092  0.21168341
##             q10
## q1  -0.10376549
## q2  -0.62191269
## q3   0.05476263
## q4  -0.10003128
## q5  -0.50936500
## q6   0.37284708
## q7   0.46584952
## q8   0.23963092
## q9   0.21168341
## q10  1.00000000
## Summary statistics for all schools
##        q1                  q2                 q3              q4        q5      q6      q7     q8     q9
##  Min.   :0.0008677   Min.   :0.002402   Min.   :192.0   Min.   :222.8   1:13    3:11    1:18   1: 9   1:10
##  Mean   :0.0010027   Mean   :0.005066   Mean   :210.3   Mean   :232.4   2:17    2:13    2: 6   2:12   2: 9
##  Max.   :0.0012104   Max.   :0.006664   Max.   :242.5   Max.   :245.1           1: 6    3: 6   3: 9   3: 4
##                                                                         Prop.                         4: 7
##  Stddev.: 0          Stddev.: 0         Stddev.: 12.2   Stddev.: 6.56   1:0.433 Prop.   Prop.  Prop.
##                                                                         2:0.567 3:0.367 1:0.6  1:0.3  Prop.
##                                                                                 2:0.433 2:0.2  2:0.4  1:0.333
##                                                                                 1:0.2   3:0.2  3:0.3  2:0.3
##                                                                                                       3:0.133
##                                                                                                       4:0.233
##
##
##  q10     q11     q12     q13
##  2: 7    3: 7    1:14    2:10
##  3: 5    4:10    2: 5    3:11
##  1: 8    2: 7    3:11    1: 9
##  4:10    1: 6
##                  Prop.   Prop.
##  Prop.   Prop.   1:0.467 2:0.333
##  2:0.233 3:0.233 2:0.167 3:0.367
##  3:0.167 4:0.333 3:0.367 1:0.3
##  1:0.267 2:0.233
##  4:0.333 1:0.2
##
##
##
##  Heterogeneous correlation matrix
##               q1           q2          q3           q4           q5          q6           q7          q8
## q1   1.000000000  0.133975478  0.01622790 -0.009867378 -0.081550954  0.15864389  0.067322031  0.21378933
## q2   0.133975478  1.000000000  0.09023096 -0.108703527 -0.095881413 -0.51399440 -0.243428508 -0.20560496
## q3   0.016227900  0.090230962  1.00000000 -0.227139740  0.160892665  0.33708826  0.250078089  0.12602260
## q4  -0.009867378 -0.108703527 -0.22713974  1.000000000  0.004666439  0.07382203  0.002347129 -0.16742935
## q5  -0.081550954 -0.095881413  0.16089267  0.004666439  1.000000000 -0.02678469 -0.013606310  0.00000000
## q6   0.158643886 -0.513994402  0.33708826  0.073822034 -0.026784690  1.00000000  0.321774454  0.14794094
## q7   0.067322031 -0.243428508  0.25007809  0.002347129 -0.013606310  0.32177445  1.000000000  0.44899801
## q8   0.213789335 -0.205604963  0.12602260 -0.167429352  0.000000000  0.14794094  0.448998009  1.00000000
## q9   0.141020776  0.155602596  0.27651935  0.040572956  0.011602840  0.44506725 -0.068817003 -0.05146368
## q10  0.071371390 -0.304661235 -0.38668941  0.319679507 -0.238427290  0.08597810  0.450211333  0.08178613
## q11  0.021918933 -0.387752768  0.29109020 -0.123460442 -0.006190976  0.12306411  0.437251505  0.14930461
## q12 -0.014873105 -0.005718078  0.17280370  0.052718213 -0.032039509  0.15252463  0.310527044  0.30238666
## q13  0.212828511  0.041805453 -0.11798894  0.306596948 -0.283699975  0.20264615 -0.176391291 -0.25885454
##               q9         q10          q11          q12         q13
## q1   0.141020776  0.07137139  0.021918933 -0.014873105  0.21282851
## q2   0.155602596 -0.30466124 -0.387752768 -0.005718078  0.04180545
## q3   0.276519348 -0.38668941  0.291090205  0.172803702 -0.11798894
## q4   0.040572956  0.31967951 -0.123460442  0.052718213  0.30659695
## q5   0.011602840 -0.23842729 -0.006190976 -0.032039509 -0.28369997
## q6   0.445067252  0.08597810  0.123064111  0.152524632  0.20264615
## q7  -0.068817003  0.45021133  0.437251505  0.310527044 -0.17639129
## q8  -0.051463685  0.08178613  0.149304607  0.302386657 -0.25885454
## q9   1.000000000  0.17904889  0.005794865  0.004541747  0.32689580
## q10  0.179048894  1.00000000  0.241296037  0.032841496  0.18980039
## q11  0.005794865  0.24129604  1.000000000  0.232951390 -0.26997963
## q12  0.004541747  0.03284150  0.232951390  1.000000000 -0.54364869
## q13  0.326895795  0.18980039 -0.269979634 -0.543648691  1.00000000
## Summary statistics for all classes
##        q1                  q2                 q3              q4        q5      q6      q7      q8      q9
##  Min.   :0.0007098   Min.   :0.001869   Min.   :184.6   Min.   :217.0   1:238   1:159   1:234   1:213   1:261
##  Mean   :0.0009992   Mean   :0.005042   Mean   :213.3   Mean   :233.5   2:194   2:117   2:138   2:139   2:150
##  Max.   :0.0013015   Max.   :0.007609   Max.   :248.1   Max.   :261.6   3:208   3:111   3:268   3:103   3:229
##                                                                                 4:105           4:185
##  Stddev.: 0          Stddev.: 0         Stddev.: 11.04  Stddev.: 5.82   Prop.   5:148   Prop.           Prop.
##                                                                         1:0.372         1:0.366 Prop.   1:0.408
##                                                                         2:0.303 Prop.   2:0.216 1:0.333 2:0.234
##                                                                         3:0.325 1:0.248 3:0.419 2:0.217 3:0.358
##                                                                                 2:0.183         3:0.161
##                                                                                 3:0.173         4:0.289
##                                                                                 4:0.164
##                                                                                 5:0.231
##
##
##  q10     q11
##  1:182   1:215
##  2:100   2:101
##  4: 96   3:125
##  5:175   4:199
##  3: 87
##          Prop.
##  Prop.   1:0.336
##  1:0.284 2:0.158
##  2:0.156 3:0.195
##  4:0.15  4:0.311
##  5:0.273
##  3:0.136
##
##
##
##  Heterogeneous correlation matrix
##               q1          q2           q3            q4            q5          q6           q7           q8
## q1   1.000000000  0.09769012  0.051350380 -0.1072436909  0.0458736585  0.03235084  0.018640128 -0.072942822
## q2   0.097690116  1.00000000  0.175532049 -0.0486588081  0.1584892266 -0.03418544 -0.036849970 -0.018688698
## q3   0.051350380  0.17553205  1.000000000  0.0464788844 -0.0753195698  0.06459476 -0.004405361  0.007288513
## q4  -0.107243691 -0.04865881  0.046478884  1.0000000000  0.0005744827 -0.18203934 -0.017194017  0.085701600
## q5   0.045873658  0.15848923 -0.075319570  0.0005744827  1.0000000000  0.03994391  0.104589942 -0.001320349
## q6   0.032350842 -0.03418544  0.064594756 -0.1820393429  0.0399439139  1.00000000  0.091271817 -0.060375470
## q7   0.018640128 -0.03684997 -0.004405361 -0.0171940172  0.1045899420  0.09127182  1.000000000 -0.022285469
## q8  -0.072942822 -0.01868870  0.007288513  0.0857016001 -0.0013203488 -0.06037547 -0.022285469  1.000000000
## q9  -0.011202132 -0.12422065 -0.080860507 -0.0308860725 -0.0875450659  0.02223704  0.096642682  0.004443342
## q10  0.006436406 -0.01819545  0.039495483 -0.0072452324 -0.0712768727 -0.02228744  0.085970550 -0.046567131
## q11  0.078475378  0.11911179  0.032173028 -0.0022691267 -0.1227604302 -0.14893931 -0.045939746 -0.010390039
##               q9          q10          q11
## q1  -0.011202132  0.006436406  0.078475378
## q2  -0.124220652 -0.018195446  0.119111787
## q3  -0.080860507  0.039495483  0.032173028
## q4  -0.030886072 -0.007245232 -0.002269127
## q5  -0.087545066 -0.071276873 -0.122760430
## q6   0.022237036 -0.022287436 -0.148939307
## q7   0.096642682  0.085970550 -0.045939746
## q8   0.004443342 -0.046567131 -0.010390039
## q9   1.000000000  0.123594626  0.049692506
## q10  0.123594626  1.000000000 -0.040206184
## q11  0.049692506 -0.040206184  1.000000000
wleoncio commented 3 years ago

All comments addressed on the aforementioned commits.