Simulate Propensity Scores

sim_arguments <- list(
    propensity = list(
        formula = trt ~ 1 + age + ses,
        fixed = list(age = list(var_type = 'ordinal', 
                         levels = -7:7),
                     ses = list(var_type = 'continuous', 
                         mean = 0, sd = 5)),
        sample_size = 1000,
        error = list(variance = 5),
        reg_weights = c(2, 0.3, 0.15),
        outcome_type = 'binary'
    )
)

prop_data <- simulate_propensity(sim_arguments)

# prop_data |> ggformula::gf_density(~ ses, fill = ~factor(trt))
# prop_data |> ggformula::gf_density(~ age, fill = ~factor(trt))

Include Propensity scores into broader process

sim_arguments <- list(
    formula = achievement ~ 1 + motivation + trt + age + ses,
    fixed = list(
        motivation = list(var_type = 'continuous',
                           mean = 0, sd = 20)
    ),
    sample_size = 1000,
    error = list(variance = 10),
    reg_weights = c(50, 0.4, 1.2, 0.1, 0.25),
    propensity = list(
        formula = trt ~ 1 + age + ses,
        fixed = list(age = list(var_type = 'ordinal', 
                         levels = -7:7),
                     ses = list(var_type = 'continuous', 
                         mean = 0, sd = 5)),
        sample_size = 1000,
        error = list(variance = 5),
        reg_weights = c(2, 0.3, -0.5),
        outcome_type = 'binary'
    )
)
simulate_fixed(data = NULL, sim_args = sim_arguments) |>
    simulate_error(sim_args = sim_arguments) |>
    generate_response(sim_arguments)
##      X.Intercept.    motivation trt age           ses level1_id        error
## 1               1 -20.226141182   1   3   2.954775206         1  0.750026529
## 2               1  -7.490679620   1  -3  -3.528203718         2  0.118249670
## 3               1 -16.125839573   1   6   4.933441999         3 -1.188047073
## 4               1   3.781436764   1   5   0.161795225         4  7.815044107
## 5               1  26.181660690   1  -7  -2.480217866         5 -0.088093375
## 6               1  -2.642434864   0  -7   8.366245297         6  1.712901619
## 7               1  -5.828875327   1   4  -5.256333071         7 -4.205412827
## 8               1  36.588154665   1  -1   0.140079107         8  2.338189630
## 9               1  10.338440849   0   2   2.634450455         9 -2.132025342
## 10              1  40.155655867   1  -4   6.353934600        10  2.403102384
## 11              1  46.297221425   1   4  -2.937962229        11  1.799508193
## 12              1  12.030590784   1  -2  -2.152357221        12 -0.490826976
## 13              1 -16.316030889   1   4  -0.114831313        13 -0.632547275
## 14              1  13.822094047   1   0   0.374369613        14  5.545567904
## 15              1  20.895671968   1   2 -10.449978079        15  1.542430354
## 16              1  -8.192980448   1   1   5.119485865        16 -1.737663422
## 17              1  22.034451410   1   6   3.116200658        17  5.927179127
## 18              1  33.593509262   1  -4  -6.836650042        18 -1.438759904
## 19              1 -17.018234492   0  -4   0.222947613        19  1.257481317
## 20              1   6.042164381   0  -6   1.831368995        20 -1.489699036
## 21              1  18.938440349   1  -6  -0.892348487        21 -2.324748255
## 22              1  36.585653028   1   7  10.508387967        22 -4.132211365
## 23              1   6.013838484   0  -7  -1.212048289        23 -1.953736020
## 24              1 -17.581887031   1   2   0.058586076        24 -2.482515542
## 25              1  10.080013505   1  -5 -13.892168443        25 -0.003340862
## 26              1  -3.902597980   0  -5   1.372456974        26 -4.478585028
## 27              1  16.069000547   1  -5  -3.113659639        27  0.339640501
## 28              1 -16.393086359   1  -4  -5.139158587        28 -3.051612910
## 29              1 -48.608374804   0   2   3.644950546        29 -4.584101012
## 30              1  19.826140837   1   7  -1.962307183        30  0.979951740
## 31              1 -11.277585257   0  -5   1.789227877        31 -0.222967432
## 32              1  28.794787699   1   7  -3.775949300        32  0.802932812
## 33              1  -3.275533676   1  -5   0.333448564        33  3.713799476
## 34              1 -35.596624537   1   1   1.286778252        34 -2.151358230
## 35              1  26.363637498   1   1  -0.454815282        35  0.816803500
## 36              1   9.840401825   1   0  -3.652540133        36  1.315272253
## 37              1  -7.517528891   1  -5  -1.979206232        37 -6.402510286
## 38              1  27.086908281   1   2  -7.079042577        38 -3.219038550
## 39              1  34.397574968   1  -5  -0.544694590        39 -0.496349275
## 40              1  14.589265797   1  -4   0.371499550        40 -3.652505589
## 41              1  17.411507604   1   3   2.250950726        41  1.683850319
## 42              1  21.419894504   1   3   1.936729091        42 -0.700419305
## 43              1 -17.191349622   1   2  -1.291475068        43 -1.684717135
## 44              1  55.619417201   1  -5   7.617453013        44 -0.799285155
## 45              1 -15.834627711   1   1   2.756717191        45  1.675224396
## 46              1 -48.216645129   0  -7   1.649508867        46 -1.135302578
## 47              1  -5.658817814   1  -5   0.864613477        47  5.809993624
## 48              1  22.620101677   0  -5   2.133207681        48 -1.634578147
## 49              1 -12.000164599   1  -3  -4.068497924        49  1.340200373
## 50              1  -0.264383692   1  -2   3.256775470        50 -4.564393000
## 51              1   9.941663793   0  -7   4.599908980        51 -4.879307114
## 52              1  40.066375493   0  -7   1.950508939        52 -0.272184506
## 53              1 -26.336910404   0   1   4.783620490        53 -0.089251954
## 54              1  26.039632657   1   3  -0.849867817        54 -3.571454040
## 55              1 -38.793398003   1   1  -1.899889581        55 -2.372290390
## 56              1 -12.642586688   1   0   0.721039785        56 -0.973305522
## 57              1  -3.865034645   1   7  -3.591420445        57 -2.017293288
## 58              1 -38.381799814   0  -5   1.030543240        58  2.447287783
## 59              1   6.840027644   0   0  -1.836671806        59  0.363596721
## 60              1  -0.432732119   0   2   3.942214385        60  1.465797716
## 61              1  -8.025461949   1   1  -1.554763977        61 -4.209326757
## 62              1  -1.398401682   1   4  -4.295700496        62  1.811654820
## 63              1   3.029043977   1   3   5.426398314        63  0.379193330
## 64              1 -38.578895384   1   1  -3.923896620        64  0.414277708
## 65              1  -8.071686009   0   0  -6.104801168        65 -1.719475229
## 66              1  11.425770271   1   7  -5.283457160        66  1.156773435
## 67              1  17.789329786   1   1  -4.984452770        67  4.798803187
## 68              1  28.043807755   0   1  -3.363025425        68 -4.796199240
## 69              1  23.373781278   0   5   3.326341836        69  1.920440883
## 70              1 -37.867961956   0  -2  11.180388415        70  1.270611912
## 71              1   1.943141337   1   1  -3.008990846        71 -1.951057968
## 72              1   5.135308701   1   4  -1.458088710        72 -0.289790057
## 73              1 -12.135103843   1  -3  -4.204874310        73  2.215131285
## 74              1  35.012491717   1   6  10.158497576        74  2.384127877
## 75              1 -20.799395860   1  -5   2.146860626        75 -0.384662358
## 76              1  18.737823220   0  -1   4.055087606        76  2.563700756
## 77              1 -28.618008686   1  -2  -0.871499743        77  0.414606807
## 78              1   1.317347114   1   6  -1.581857416        78  1.037977841
## 79              1   2.577216044   1   6  -2.465009474        79 -0.521589095
## 80              1  18.919682842   0   0   0.953892371        80  1.285247997
## 81              1  -7.940087692   1   2  -5.518724919        81  2.444418407
## 82              1 -26.334167189   1  -4  -0.532066568        82  1.670246457
## 83              1  -0.426852705   1  -2   7.294963564        83 -0.590128179
## 84              1  28.420255411   1   1  -1.537371491        84  4.911893714
## 85              1  -3.607639143   1   1   3.785907057        85  2.358244144
## 86              1  -1.254163015   1   2  -8.554126950        86 -2.605166287
## 87              1 -12.615604696   1   7   2.976548039        87 -3.190370851
## 88              1   3.291452127   1   0 -10.339833252        88  0.681601838
## 89              1  -1.219427637   1  -7  -0.428698780        89 -2.027925786
## 90              1   5.061885082   1   7   0.161432749        90  2.333466711
## 91              1  18.523670875   1  -4   1.448431987        91 -0.629668080
## 92              1  11.257257990   1  -2   4.044690344        92  2.251839225
## 93              1 -29.451168410   1   7   1.631312357        93  4.845736683
## 94              1  24.432043582   1   3  -8.519287396        94 -4.025908558
## 95              1  16.181078860   0   2   2.627522946        95 -0.141142173
## 96              1  32.426275529   1   7  -4.015169249        96 -2.716030685
## 97              1  -6.917215345   1  -6  -1.188278604        97  0.743261281
## 98              1   6.972488065   1  -3   4.276745293        98  3.119655303
## 99              1  -7.370479797   1  -2  -2.347435224        99  1.627187588
## 100             1  11.090257967   0  -6   0.590063920       100  1.264474003
## 101             1   6.672095471   1   2   6.237207774       101  0.534891919
## 102             1   8.405713663   1  -3  -3.441986493       102 -1.265116413
## 103             1 -12.587260258   1  -3  -6.917847875       103 -3.328225204
## 104             1  -2.441933615   1  -4  -0.214662076       104  3.475620042
## 105             1 -11.008466021   1  -6  -1.046171149       105  1.205449635
## 106             1   9.710972763   1  -1  -7.334967272       106  2.887886682
## 107             1  26.640220616   0  -6  11.178891900       107 -2.272561896
## 108             1   1.312845398   1   1  -3.125724817       108  3.047043597
## 109             1  45.838754849   1  -7  -5.839568927       109  4.692583971
## 110             1 -15.425100671   0  -5   0.747644157       110 -4.207700202
## 111             1  18.383923492   1   5   4.306599502       111  2.261552855
## 112             1  -6.431093599   0  -4  -5.231017035       112 -1.449080464
## 113             1  13.493738784   1  -2   0.861300056       113 -8.700070330
## 114             1 -29.337919296   1   3  -1.535103755       114 -0.671175866
## 115             1  11.319996638   0   3   4.516728049       115 -0.078507145
## 116             1   6.528605880   0   7   4.926243468       116  1.744719158
## 117             1   6.932575528   0   2   8.219142418       117 -3.352145984
## 118             1   0.386333975   0  -6   4.072798622       118 -1.312174991
## 119             1  13.777893890   0   4   4.125565803       119 -1.296239769
## 120             1 -32.600633152   1  -2  -1.728869611       120 -0.304126667
## 121             1   3.268232007   1   7   2.055736674       121 -0.463065012
## 122             1  32.282683507   1  -7  -3.216071915       122  1.319474328
## 123             1   8.954615688   1  -3   4.710297384       123  3.121244453
## 124             1   9.675958641   1   2   2.347958919       124  0.556582527
## 125             1 -36.046231724   1  -5  -3.310561190       125  3.249583844
## 126             1   0.217996750   1  -1  -1.547355511       126 -4.679046580
## 127             1  -6.951356075   1  -3   3.972871318       127 -0.043414673
## 128             1 -11.420423842   1   5  -4.574254782       128 -0.315868863
## 129             1  -2.636162643   1  -7  -2.327601547       129 -1.881740134
## 130             1  38.206366222   1  -2   4.098346815       130 -0.126073105
## 131             1   4.067631967   0   2   7.605178690       131 -2.766454514
## 132             1 -14.895593752   1   5   1.326773952       132  4.366329165
## 133             1  -9.395918296   1   6   3.095334699       133 -0.495104382
## 134             1 -20.190045688   0  -2   2.355404183       134 -0.767287635
## 135             1   6.555414430   0  -7   2.315721067       135 -0.524118053
## 136             1   8.948429395   0   3  -0.800806814       136 -2.022108051
## 137             1  29.935026978   1  -5  -0.178427360       137 -1.616016541
## 138             1   5.299618581   0  -7  -2.873424079       138  2.207440188
## 139             1  -0.288802662   0   6   4.226750084       139  4.457754735
## 140             1 -31.498777022   1  -3  -2.908830613       140 -1.633647205
## 141             1   0.763019685   1   6  -2.511785528       141 -1.145237607
## 142             1 -59.443528118   0  -7   5.440381323       142 -0.680971572
## 143             1 -17.538747773   1  -4  -1.964765097       143 -2.782401365
## 144             1  12.025787332   1   6  -3.302305596       144 -0.414359985
## 145             1  20.936710342   0  -2  -0.745125634       145 -5.746871874
## 146             1   4.734085396   1   6  -4.377426488       146  2.665506883
## 147             1  21.009446201   0   3   8.495007854       147 -2.435973198
## 148             1 -14.402858498   1  -7   0.340951308       148 -1.853056098
## 149             1  -1.967190588   0   1   2.665937901       149 -1.717219964
## 150             1   5.409705212   1   0  -0.404744966       150 -1.718797047
## 151             1  -1.484247345   0   0  12.233907139       151 -1.171699237
## 152             1 -21.831098815   1  -7   1.919495168       152  3.531919140
## 153             1  11.530222700   1  -4  -5.260857890       153  0.031175067
## 154             1  21.020758143   1   3  -8.611169134       154  6.995930232
## 155             1  -8.396890358   0   2  -3.154096248       155 -5.902784453
## 156             1  12.118392958   1   1  -4.767979897       156 -7.213414180
## 157             1  20.721247545   0  -4  -0.399912326       157 -3.238120020
## 158             1 -16.798730745   0  -6   5.313479511       158  2.422403521
## 159             1  30.408828916   0   0   6.755284179       159 -1.221610925
## 160             1  41.021601781   1   6   1.166451102       160 -2.987443816
## 161             1  18.931410247   1  -7   1.669493957       161 -0.387968360
## 162             1 -15.669127659   1  -4  -3.823851290       162  2.142075871
## 163             1  21.955410061   1  -1   3.846922301       163 -6.504228147
## 164             1  42.598631217   0   1  -1.075064807       164  4.451078060
## 165             1  -9.914791683   1   7   1.916902926       165  1.965800010
## 166             1   2.675691584   1   4  -2.712476806       166  4.905213389
## 167             1 -45.048913760   1   1 -10.472891099       167  2.152847637
## 168             1   6.595583352   1   4   1.890250293       168  1.903840981
## 169             1  12.340285768   1   7  -3.529671141       169  2.251984226
## 170             1 -10.840188764   0  -3  -2.935522649       170 -1.516425603
## 171             1   5.998924567   1  -6   0.500431796       171  3.818610940
## 172             1  20.096889641   0  -7   6.010331794       172 -3.863115241
## 173             1   6.790873093   1   1   3.299219980       173 -2.128857312
## 174             1 -22.567561406   1   5  -3.133726778       174  5.741514480
## 175             1   9.222023913   1   0  -2.648880530       175 -0.721643563
## 176             1  15.972577812   1  -3   4.382515935       176  3.369247912
## 177             1  32.459264383   1  -5  -0.250691059       177 -1.916768421
## 178             1   5.657295692   1   7  11.401951764       178 -3.223726828
## 179             1 -29.245752673   1   6   6.899052819       179  0.154324852
## 180             1 -24.353935132   1  -1  -2.868180171       180 -4.068790320
## 181             1   5.036823087   1  -2  -7.084187140       181 -6.041616726
## 182             1   8.265203477   1  -6  -6.224537295       182  0.275729049
## 183             1   6.257315062   0   0   4.500412400       183 -5.458582393
## 184             1 -22.773427764   0  -7   9.947151283       184  0.189414271
## 185             1  -5.373255833   1   3 -13.934876432       185 -0.064392736
## 186             1  40.293651507   1  -2  -0.115676521       186  0.697334114
## 187             1  68.189097548   1  -2   1.020804197       187 -5.648188460
## 188             1  14.517978582   1  -5   0.221142623       188 -3.162949686
## 189             1  21.681233875   1  -2   3.135487439       189  0.515239375
## 190             1  16.238995681   1  -7  -4.736247771       190  7.296652636
## 191             1   2.774816670   1   4  -5.614134719       191  2.435552848
## 192             1  -2.029976563   1   2  -0.101618247       192  2.064486726
## 193             1   3.931604907   1   6   6.125127717       193 -2.402031190
## 194             1   7.209003154   1   0   0.415294923       194  0.604231325
## 195             1  14.222196019   1  -1  -9.429575852       195  0.061181973
## 196             1 -17.739163122   1   1   5.284635638       196 -3.431852160
## 197             1 -35.745067837   1  -4   3.064560943       197  1.431907654
## 198             1  -6.437625491   0  -7   2.948250454       198 -2.187540451
## 199             1   5.285612326   0  -7  -0.332256705       199  0.734808851
## 200             1 -11.358184929   1   5  -3.040539404       200 -1.745398643
## 201             1  22.024710162   0  -4   4.736522872       201  0.856600964
## 202             1  21.394188233   1  -2   1.965849788       202  2.550710056
## 203             1  22.351215906   1   5  -1.926029656       203  3.282163412
## 204             1 -17.364638632   1   4  -3.046315609       204  0.823165117
## 205             1 -42.158112554   1  -3  -1.921925730       205 -2.376136458
## 206             1  42.624943626   1  -3  -4.507118346       206 -0.278558660
## 207             1  39.824029100   0   5   4.364729467       207  0.188527828
## 208             1  11.933221511   0  -5   2.809106293       208  2.385801945
## 209             1 -15.019766022   1  -2  -8.786392610       209 -0.925385458
## 210             1 -21.459190273   1  -7  -4.060580475       210 -3.631694173
## 211             1  16.854441797   1   3   1.193452937       211  1.876633370
## 212             1  41.437551252   0  -2  12.530008597       212 -0.772028141
## 213             1  -6.029000048   1  -5  -4.902935825       213 -3.976985042
## 214             1 -38.770015903   0   0   7.830625277       214 -2.123854207
## 215             1  -7.257890778   1   7  -3.794337177       215  1.026139173
## 216             1  20.861345960   0  -7   2.469225522       216  8.052832513
## 217             1  16.704770142   1   1  -1.635752260       217 -0.693185691
## 218             1   9.606549727   1   3  -2.174395199       218  5.843353001
## 219             1 -15.307745442   0  -2   7.598615530       219 -1.673763388
## 220             1  -5.588543422   1   3   0.018389016       220  5.257484506
## 221             1 -32.353904718   1  -2  -2.011699751       221 -3.182951062
## 222             1   4.196134177   1   3  -7.477332166       222  2.650953673
## 223             1  30.539633906   1  -7  -4.729927566       223  0.017722118
## 224             1 -12.293739090   0   2   4.532147541       224  6.777331454
## 225             1 -34.343081659   1   6  -9.069338114       225 -2.374180943
## 226             1  27.225299905   1  -5  -4.758877158       226  3.569956858
## 227             1  -9.447590037   1  -6   5.634260496       227  2.434484630
## 228             1 -12.058335946   1   7   0.699445324       228 -1.471415870
## 229             1  20.007053130   1  -1   0.501699884       229 -5.510127988
## 230             1 -21.997251451   1  -3  -6.022643147       230  4.719019024
## 231             1 -36.029087227   1  -2  -0.417291482       231 -0.038780139
## 232             1  -9.248830669   1  -3  -3.985411234       232 -2.367867082
## 233             1 -25.011893797   0   1   5.507622893       233 -1.339574310
## 234             1  25.680548634   1  -4   0.542052031       234 -2.378307586
## 235             1  30.740715839   1  -5  -5.952533995       235 -2.227937976
## 236             1 -14.177203017   0  -2   1.399606380       236  0.505982214
## 237             1  -7.676631217   0  -7   3.887038866       237 -1.618652208
## 238             1  -2.024515483   1  -5  -1.105448620       238 -3.214252134
## 239             1 -13.275602319   0  -3  -0.688384177       239  5.346589087
## 240             1  48.155975966   1   4   0.576150216       240 -1.564433108
## 241             1  12.044364423   0  -6   7.666347864       241  0.879008543
## 242             1  33.693692191   0  -4  -3.278872301       242 -1.614392968
## 243             1  -5.949181992   0   1   4.291893504       243 -0.072129274
## 244             1   6.043216789   1   2  -5.938352759       244 -4.847003897
## 245             1  -6.441468531   1   0   2.697701902       245  1.361108880
## 246             1  10.329432062   1   4  -1.331810814       246 -1.993633820
## 247             1 -16.960109363   0  -4   3.790861173       247  1.043990535
## 248             1  -4.686035555   1   0  -5.157010146       248 -0.841655505
## 249             1  -9.972295613   1   1   5.124657552       249 -2.245526705
## 250             1  25.053802261   1  -2  -0.726307961       250  3.321355196
## 251             1   8.229287217   1  -4   0.546119151       251 -2.072303724
## 252             1 -20.678834146   1  -5  -7.392993973       252 -2.627937523
## 253             1  -1.772706716   1   0 -11.395956112       253 -0.356005835
## 254             1 -19.426045687   1  -2   3.000384452       254  5.736613397
## 255             1   2.271477119   1  -3  -1.732890997       255 -7.344838242
## 256             1  37.820669320   1   7  -1.425768360       256 -2.134152601
## 257             1  23.356054306   1   0  -3.238260790       257  1.464915310
## 258             1  -8.594627863   0  -6  -6.361073502       258 -3.368787277
## 259             1  -8.301383195   0  -1  -0.875865904       259  1.371476012
## 260             1 -19.017212057   1   7 -11.006597887       260 -0.560305384
## 261             1  26.707859456   1  -4   2.899379897       261  0.839505012
## 262             1 -18.172505532   1  -7   0.313982355       262  1.142168262
## 263             1 -36.540226561   0   2   8.875622590       263 -0.008316246
## 264             1  -2.086635108   0  -4   6.918473499       264 -4.640163875
## 265             1  41.274048669   1   4   8.154801432       265 -4.606649547
## 266             1  29.133390088   1   7  -4.041956491       266 -6.059305239
## 267             1   4.886892636   1  -2   0.477827257       267 -3.702518056
## 268             1   8.584484286   0   7   4.584040951       268 -2.147386488
## 269             1  27.826047885   1   2   2.369217715       269 -0.562280150
## 270             1  -4.407351609   1   7  -1.976228580       270  2.898838580
## 271             1  41.128800126   1  -1   2.146992719       271  3.471474130
## 272             1  -0.798850626   0  -7   4.254653954       272 -0.995402782
## 273             1  13.600903212   1   0  -0.851973089       273 -0.616233325
## 274             1  -6.207545322   1  -7  -0.003957019       274  0.190536113
## 275             1  -5.389936473   0  -3   4.303207477       275 -3.124158100
## 276             1 -16.992499951   1   4   3.201711901       276  0.192717374
## 277             1  -5.719788659   1  -6   0.534381448       277 -3.604381671
## 278             1 -23.284664930   1  -3  -7.302388822       278 -4.367308695
## 279             1 -14.536669730   1   4   2.275856877       279  2.628707149
## 280             1  15.296058281   0   6   6.360423421       280  0.240835009
## 281             1  -0.536363684   1   1  -9.575937903       281 -0.422119734
## 282             1 -20.874448459   1  -6   1.177945473       282  4.315311981
## 283             1   0.345782309   0   1   5.279571548       283  0.678765095
## 284             1   4.190039758   0  -7   2.158386219       284  2.110847965
## 285             1 -19.326620088   1  -4  -8.729716060       285 -1.030829348
## 286             1  35.745059899   0  -2  13.172877381       286 -4.039121043
## 287             1 -11.244162824   1  -7  -0.154017857       287 -0.314768346
## 288             1 -31.774878442   0  -6 -11.280943819       288 -5.856655152
## 289             1  -6.142735352   1   0   9.184112695       289  1.929643821
## 290             1  -8.558320981   1   2   2.492677830       290 -0.706388846
## 291             1 -11.873440702   0  -2   0.781234234       291 -1.877548116
## 292             1   9.380168440   1  -2   2.782508453       292  3.271282809
## 293             1  16.249850494   1   5   3.435534888       293  4.230856287
## 294             1 -12.174386784   1  -7 -16.520911148       294  2.013654953
## 295             1  -3.847973017   1   1   4.391649815       295  3.691621547
## 296             1   1.566186577   1   6   7.305878432       296 -3.382957955
## 297             1  14.336067956   1  -3  -0.547033404       297  2.899587891
## 298             1 -10.574030304   1  -1   6.628031895       298  1.051179093
## 299             1 -16.536539318   1   5   5.703402869       299 -2.948449401
## 300             1  -2.163071995   1   3  -4.651821073       300 -2.520034070
## 301             1   1.121531168   1   6  -0.973116691       301 -0.063580288
## 302             1 -16.436167721   1   0   4.620036346       302  1.125683019
## 303             1  41.005664173   1   7  -1.418856794       303  2.418177732
## 304             1  28.316315127   0  -6   3.219094216       304 -1.324729744
## 305             1  37.350256473   0  -6   6.024649877       305 -0.631585390
## 306             1  13.768356029   1  -5  -4.932167386       306 -4.451565684
## 307             1  -0.059684825   0  -5   6.814997005       307 -3.900073143
## 308             1  25.133622925   1  -1  -3.185683195       308  3.228015491
## 309             1   2.433185855   0   4  -6.263525723       309 -2.154739091
## 310             1 -16.345390278   0  -2  -2.875688652       310 -1.609693258
## 311             1 -10.187149800   1  -1  -8.381504087       311 -2.245909182
## 312             1   8.632617168   1   6   1.314060091       312  1.884543972
## 313             1 -24.449788298   1  -3  -7.104416464       313  2.859375633
## 314             1  44.667167387   1  -3   3.869233581       314  4.043493468
## 315             1  -4.975509545   1   1  -3.150956990       315 -3.920650320
## 316             1  56.010051145   1   1 -12.616114840       316  3.405287711
## 317             1 -49.783499123   0  -2   5.228380618       317  4.690180897
## 318             1   2.678176063   1   7   3.756701165       318 -6.561601605
## 319             1   3.458070781   0   7   3.990804267       319 -1.061010463
## 320             1  12.427155486   1  -6  -7.838209169       320  4.461252975
## 321             1   0.829692229   0  -5  -5.737384318       321  0.363949873
## 322             1   3.618404903   1   5   2.177933610       322  2.577135616
## 323             1  22.731849782   1  -2   5.266002095       323 -0.527836669
## 324             1 -12.979961291   1  -4   4.817533024       324 -1.200104172
## 325             1 -14.997253213   0  -2   3.469654689       325  0.299257741
## 326             1  -4.558356265   1   6  -1.150354758       326  0.759171804
## 327             1  27.283198625   1   3  -4.677793047       327 -2.000737811
## 328             1  12.313034456   1   3  -4.513101889       328 -1.382378575
## 329             1 -12.986293631   1  -4  -0.154300470       329  0.531712123
## 330             1   1.184126853   1   3  -4.343391123       330  3.351523494
## 331             1  12.034290562   1  -1  10.938937520       331 -0.440103033
## 332             1 -17.699641986   1  -4   7.692650773       332 -2.154332176
## 333             1   5.746035905   1  -3  -9.410363355       333  1.940070069
## 334             1  24.588781318   0  -7   4.579375023       334 -2.263088680
## 335             1  -0.030946115   1  -7  -2.008700070       335  4.077814486
## 336             1  -3.261427563   1   5  -6.284370001       336  2.762294037
## 337             1 -24.303774990   0  -3   5.921984346       337  4.314538989
## 338             1   8.066190913   1   2   5.042862318       338  0.994209178
## 339             1 -16.399339204   1  -5  -5.330760578       339 -0.004625981
## 340             1   8.759096113   1   1  -5.401069989       340 -2.661712864
## 341             1   2.274971799   1   6  -7.632951853       341  1.489884335
## 342             1  16.020061293   1   3  -6.852002131       342  0.686437016
## 343             1 -13.665363292   1  -3  -6.467847165       343 -0.901816338
## 344             1 -36.550732075   1  -7 -14.805148566       344 -1.442226738
## 345             1   4.403695223   0   4   5.627978415       345  4.614729509
## 346             1  12.617035526   1   7   0.020137688       346 -4.735149896
## 347             1  11.900529759   1  -3   0.813355770       347  0.722191983
## 348             1   2.095307399   1   4  -7.419611849       348 -2.037558801
## 349             1 -16.550878074   0   5   5.824301592       349 -0.939159553
## 350             1 -45.408474417   1  -1  -6.362096502       350  2.516177170
## 351             1   9.797071630   1   7  -5.568893329       351  2.319886278
## 352             1 -28.976248140   1  -3  -6.534685224       352 -4.700596363
## 353             1 -12.727451498   1  -3  -0.871041151       353 -1.835317770
## 354             1  12.295106457   1  -6  -3.416257617       354 -4.358686797
## 355             1  -2.315029453   1   0  -3.107841430       355  1.960450500
## 356             1  -7.054769100   0  -6  -0.775816160       356 -0.215087237
## 357             1  33.615176454   1   0  -5.432467914       357 -4.527471446
## 358             1   4.145059091   1   7   2.751110034       358  0.933027151
## 359             1  -3.053072907   0  -1   6.973903389       359  2.668549294
## 360             1  -0.441539889   1   4  -7.941596290       360  1.605273957
## 361             1   5.756825620   1   6  10.337917644       361  5.663392078
## 362             1  37.320147906   1   1  -8.485565235       362  2.232923493
## 363             1  -5.636098511   0   5   4.302538156       363  3.842141126
## 364             1   4.678427340   0   4  -1.898217444       364  0.082077900
## 365             1 -39.103261041   0   0   9.238150534       365  1.026198345
## 366             1   1.140242164   0   7   5.188280895       366  0.402589634
## 367             1  -7.956405831   1   5  -7.347873841       367  0.821447125
## 368             1  -1.120106831   1  -2 -13.884741050       368  4.167215385
## 369             1 -15.269158118   1  -4  -0.384898939       369 -0.312469474
## 370             1   4.964761700   1  -3  -4.772460360       370  5.402075821
## 371             1  19.799742453   0   2   3.543397098       371 -6.333411639
## 372             1  34.528913385   1   3   8.885434041       372 -2.928277803
## 373             1 -26.307232848   0  -6   4.331571539       373 -2.134990483
## 374             1  -1.644883103   0  -7   0.765029086       374 -0.581206984
## 375             1 -15.481218355   1  -2  -0.231945397       375  0.036681554
## 376             1 -27.599637736   1  -5  -8.364237074       376 -1.324372988
## 377             1  17.815734698   1  -1  -5.449453059       377 -3.002485780
## 378             1  10.001183963   1   7  -0.380878524       378  2.282943012
## 379             1  10.976963589   1   1  -6.529158579       379  3.113528550
## 380             1 -11.416961980   1   3  -0.311922814       380 -1.468072719
## 381             1  -3.992753016   1   4   1.711857842       381 -3.128277835
## 382             1 -14.892964729   1  -1  -2.744653075       382 -3.183723767
## 383             1  -2.077224376   0  -5   4.575306633       383 -5.807616493
## 384             1   3.183441440   1   3  -0.438019518       384  1.334612267
## 385             1  -8.792751735   1  -1  -0.127897354       385 -0.436858956
## 386             1   6.643251235   1   7  -6.271787032       386 -0.473784471
## 387             1  15.921119166   0  -1   9.400624114       387  2.464058411
## 388             1  12.758563361   0   7  -0.704746843       388 -3.471016996
## 389             1  -4.577127400   0  -7  -0.533455687       389  0.184324580
## 390             1 -48.134231721   1   4   3.622775452       390 -2.325862945
## 391             1  -7.383420149   0   3   2.144090417       391  1.426333145
## 392             1  -1.745370420   1  -1  -7.901078756       392  2.418517586
## 393             1  49.905538942   1   3   3.873779072       393  1.566547138
## 394             1  -5.184173698   1   3  -5.821460622       394  0.010363493
## 395             1   3.109602071   0  -4   0.044271517       395 -6.060686263
## 396             1 -21.850187011   1   7  -2.275390009       396  3.269735741
## 397             1  -8.005846730   1   1  -5.013450238       397 -3.081244766
## 398             1  -7.433461806   0  -6   2.140347859       398 -1.829792462
## 399             1   1.602439803   0  -4   1.012811508       399 -4.684750740
## 400             1  -7.316114691   0   4   0.907916377       400  4.210537067
## 401             1  -1.388196071   0  -1   2.828469338       401 -4.894467028
## 402             1   8.158364834   0   5  -1.388284899       402  5.174761305
## 403             1  -6.887632914   1  -6  -6.821242407       403  0.458109753
## 404             1   0.277276658   0   5   6.496036248       404  0.602484614
## 405             1   4.186086003   1  -5  -9.964232787       405 -3.351094393
## 406             1 -29.185629126   0   4  -0.518580436       406 -5.831632595
## 407             1  -2.059525930   1   7   1.498831286       407 -1.563453027
## 408             1  40.718907553   1  -1  -2.015002047       408  0.261627207
## 409             1 -37.359165301   0   0   1.692454786       409  1.724219987
## 410             1   2.611711527   1  -5  -7.392135197       410  4.210013229
## 411             1  42.125343727   1   3  -0.232574630       411 -2.963190314
## 412             1  11.983328802   1  -1   0.158440930       412 -1.091051893
## 413             1 -33.127121368   0   3   1.407663036       413 -0.533082672
## 414             1  -1.686616000   0  -3   4.015976268       414 -1.992340451
## 415             1 -33.275108413   1   2  -1.044291744       415  0.995331388
## 416             1   9.681914566   1   5  -1.683640881       416  3.418040631
## 417             1   1.209558487   1  -4  -2.281548758       417 -3.974301598
## 418             1 -16.138965790   0   0   5.121301549       418 -6.500924833
## 419             1 -24.643021830   0  -7   7.286975753       419 -1.526813885
## 420             1 -15.929214376   1   6  -1.124652937       420 -6.603389552
## 421             1 -12.059434813   0   6   1.450342549       421  0.271480533
## 422             1  28.213801623   1   3   1.296382719       422  1.313216712
## 423             1  11.161025617   1  -3   0.049461499       423 -7.523385118
## 424             1   4.045615172   0   0  -0.923369627       424  3.101856598
## 425             1  10.745784433   0  -2   8.271232924       425 -2.907559271
## 426             1  21.117125511   0  -6   2.184271779       426  3.061323307
## 427             1  30.633565390   0   3   4.236580559       427  0.975718691
## 428             1  -7.787113009   0  -5   4.981411501       428 -2.382819719
## 429             1   0.513444938   0  -3   5.337262211       429  0.380915962
## 430             1 -12.349868141   1   4  -0.266712892       430 -0.464485018
## 431             1  38.220702987   1   5   4.964631126       431 -6.299306336
## 432             1  21.043874797   1  -4  -2.496763424       432 -5.639487256
## 433             1 -11.489266088   1   4  -5.581966646       433  4.418100300
## 434             1   2.093370087   1  -3  -4.872583823       434  2.180760911
## 435             1  18.607713293   0  -2  13.937017700       435  0.114003216
## 436             1  34.973028985   1  -5  -3.866567344       436  2.449186137
## 437             1   2.613879111   0   3  17.402698742       437 -9.481702337
## 438             1  23.647517240   1  -6  -3.087109868       438  2.512060665
## 439             1 -17.247748682   0   3   2.603912213       439  3.333908083
## 440             1  20.297169179   1   7   2.962701265       440 -0.223928620
## 441             1 -17.379960078   1   5  -4.501734244       441 -2.058818071
## 442             1  18.471619671   0   1   7.371406605       442  0.897512736
## 443             1  37.750110249   0  -2   5.425663013       443 -1.395575165
## 444             1 -14.266943388   1   7  -8.533437193       444  6.171384360
## 445             1  25.913920887   1   4  -5.843481499       445 -2.791484857
## 446             1 -27.179295855   1  -6  -7.097266670       446 -1.190089033
## 447             1  28.567267165   1   0   3.229854747       447  0.152477824
## 448             1  35.218748343   1  -5  -2.135663945       448  0.608250562
## 449             1   4.780455047   0   0   5.949140734       449 -2.454867734
## 450             1 -12.104742480   1   5  -0.699295712       450  2.641599611
## 451             1   2.614372453   0  -2   0.946979882       451 -2.510418487
## 452             1  19.567412193   0  -3   0.287922026       452 -2.319694247
## 453             1  14.587403310   1  -2  -2.936230633       453  2.486442516
## 454             1  29.347264539   1   5  -1.207072899       454  0.038886062
## 455             1 -24.848442382   1   1   3.819001953       455 -1.176635806
## 456             1  30.804188002   1   1  -3.390483666       456 -7.773042172
## 457             1   1.827789372   1   4   4.184970697       457  4.152502271
## 458             1  10.333562703   1   7   0.209697332       458 -4.073858908
## 459             1 -11.604125081   0   5   4.896547993       459 -3.736324421
## 460             1 -22.687207279   0   0   7.480478272       460  3.411627935
## 461             1   2.596605996   0  -6   6.905835783       461 -2.164000668
## 462             1  13.491364382   0  -5   0.885934254       462  4.871446691
## 463             1 -31.990974918   1  -2  -3.022794429       463  0.042688739
## 464             1 -43.311673510   1   3  -0.032062768       464 -0.618368564
## 465             1  10.442890306   0  -6  -5.154588022       465 -0.834980261
## 466             1  11.222448346   1   7  -2.983096506       466  3.931885667
## 467             1  16.073518529   0  -4   6.150026672       467 -1.274108933
## 468             1   8.058424823   1   4 -10.626319704       468  7.864437335
## 469             1  -0.220472141   1   0   2.732706424       469 -3.428780584
## 470             1  29.706499791   1   5   5.057276703       470 -1.836112277
## 471             1 -26.104286303   0   1   6.835873012       471 -1.042948124
## 472             1 -15.055064173   0  -2   3.539189850       472  1.751334147
## 473             1 -10.881167248   0   1  -0.612579270       473 -1.181805633
## 474             1 -13.816818856   1  -5  -3.997045673       474  3.000787589
## 475             1 -21.383329225   0  -6  -0.896800021       475 -5.051187330
## 476             1  -3.857367530   1   6   0.547264285       476  1.211860256
## 477             1   6.249293261   1   7   1.154632764       477 -1.400399719
## 478             1   4.118886946   1   3  -6.348315810       478 -3.805710574
## 479             1  12.656114850   0  -6   3.102629057       479  2.075567473
## 480             1   7.874582351   0  -6  -3.211742602       480  1.181248321
## 481             1 -17.096972091   0  -7   3.311789069       481  1.796175695
## 482             1 -16.338045855   1   6   1.229421386       482 -1.839149201
## 483             1 -23.535958728   1   2  -4.239468436       483  5.349944586
## 484             1  38.713264808   1  -4  -2.173837053       484 -3.446735547
## 485             1   3.819657069   1   0   4.415879745       485  1.528734173
## 486             1 -22.191876796   0   2   2.931593402       486 -3.737552432
## 487             1  -6.725773600   0  -6  -1.283460204       487 -2.493040952
## 488             1 -28.482175530   1  -4   6.100678524       488 -0.195814138
## 489             1   2.381470224   0   1  11.629439837       489  3.929956729
## 490             1  14.450934515   1   7  -6.535211682       490 -3.214377831
## 491             1  -9.922597628   1   6   0.001357341       491 -0.660077169
## 492             1   7.542880571   0   6   5.712557850       492  2.183740777
## 493             1 -18.327037362   1  -3  -1.109700676       493  1.745966025
## 494             1 -15.164196532   1   6  -1.041572963       494 -3.375480547
## 495             1  21.745717706   1  -2   3.749613559       495  0.372485664
## 496             1 -21.883755867   1   3   0.588028089       496  2.771280726
## 497             1  -8.602795790   0   7   0.495106019       497 -3.558429672
## 498             1  -9.452077622   0   4   6.808442344       498  0.507484701
## 499             1 -21.106238101   1   2  -4.194941193       499 -4.180707119
## 500             1 -29.291624326   1   1  -5.324686095       500 -3.941446224
## 501             1  -4.738501085   0   4  10.565840079       501 -1.568185256
## 502             1   0.647031453   0   1   2.553335934       502 -2.966433722
## 503             1  -0.124839276   0  -4   6.001315172       503 -0.086126484
## 504             1 -10.092912397   1  -2   5.227300276       504 -3.728180844
## 505             1   9.419187926   1  -3   3.602218035       505  0.899889810
## 506             1  31.224162794   0   0   2.514842037       506 -1.483169441
## 507             1 -36.179749697   0   1   4.779206068       507  2.666470882
## 508             1  24.951464937   1  -4   2.460640279       508  2.655987690
## 509             1 -13.015371845   0  -6  -0.779747342       509  1.688987751
## 510             1  -0.906437137   1   0  -2.165310888       510  1.680922311
## 511             1  23.005136463   1   4  -6.117101752       511  4.288098368
## 512             1  -2.065541670   1  -4  -1.451716864       512  0.574693674
## 513             1  -0.503009884   1   6   1.701159944       513  4.374811362
## 514             1  20.141792563   1  -1   1.172760278       514  1.219034974
## 515             1   4.466078386   1   2   0.059119685       515  0.704973866
## 516             1  14.934767471   0  -7   3.023664936       516  1.873657860
## 517             1  19.438336945   1   6   3.687548815       517 -1.716907986
## 518             1   0.205773140   1  -2  -4.490517441       518 -5.351409454
## 519             1 -11.958441447   1   2   0.859106781       519  1.162799761
## 520             1  12.418983564   0  -6   4.566173946       520  0.110842293
## 521             1  -2.143999650   1   5  -5.525767273       521 -4.337649307
## 522             1  14.338496219   1  -6   1.967250723       522 -1.999972939
## 523             1 -20.309333005   1   3  -6.849213355       523  1.045509780
## 524             1  15.465469091   0   3   3.697358077       524  2.876934981
## 525             1 -23.003113662   0  -7   3.191854685       525 -3.049271620
## 526             1   1.815587116   1  -7 -11.401769308       526 -0.539844930
## 527             1 -13.699284712   1  -3 -10.056123492       527 -0.531758219
## 528             1  -3.219267067   0  -5  10.343702799       528  0.737304041
## 529             1  -2.199234038   1   1  -6.575262787       529 -2.915315451
## 530             1 -13.204164004   1  -1  -5.548751998       530  4.247344975
## 531             1  31.983887933   1   5   0.274483598       531 -0.139998404
## 532             1  -6.542502400   1   5  -4.312936953       532  3.531396011
## 533             1 -18.983054880   0  -4   4.653877456       533 -4.835188182
## 534             1  22.457113241   1   4   2.855852004       534 -0.168496212
## 535             1 -21.269173259   0  -7   6.174159632       535  1.045205884
## 536             1  13.190500791   1   6   2.310015632       536 -3.099917213
## 537             1   2.132244085   1  -5  -3.671040972       537 -4.704589470
## 538             1 -13.350029747   0  -5   6.616803121       538  0.557000642
## 539             1   7.254940473   0   1   2.286473852       539  2.710201718
## 540             1 -15.814702569   0   4  -3.682403761       540  0.255559778
## 541             1 -35.874010604   0   5  10.718064157       541 -1.268877369
## 542             1 -29.518806912   1  -7  -1.690796840       542  1.569895882
## 543             1   3.913451058   0  -2  -0.354851999       543  1.389715606
## 544             1  -2.545365330   1   0  -0.441050102       544  0.455629762
## 545             1  -5.372239053   1   1   1.775453584       545 -3.091821857
## 546             1  -4.581198223   1   1  -6.024042278       546  4.175428564
## 547             1  38.624954922   0  -5   7.280002519       547 -5.454206123
## 548             1  -8.268491254   1  -1  -0.002328046       548  0.989416475
## 549             1 -15.478051702   0  -3  -0.588118285       549 -4.005813529
## 550             1  32.853708315   1   4   0.920741069       550  1.245589269
## 551             1 -31.797536648   1   6  -4.306003611       551 -3.655651871
## 552             1  18.584860542   1   6   3.896554234       552 -0.804275358
## 553             1   7.118968316   1   7  -2.270320858       553 -2.700983145
## 554             1   5.428815685   1  -7 -10.143035922       554  0.322470029
## 555             1   4.469292067   0   2   0.687490614       555  0.774759701
## 556             1 -20.902399992   1   3  -3.265826247       556  0.699298571
## 557             1  15.515303327   1   6   3.600175490       557  0.586754669
## 558             1   0.552602790   1  -2   4.042182138       558  3.502203355
## 559             1  23.008571091   1   5   4.812260652       559 -0.917288118
## 560             1   5.068444497   1   5   0.515690362       560  1.995453927
## 561             1 -16.497736462   1   4   3.431800757       561  1.883315334
## 562             1   5.735263783   1   5   5.475994668       562 -0.822304823
## 563             1  18.858456238   0  -5   0.909461077       563  6.413233834
## 564             1 -10.621528005   1   3  -0.984377081       564  1.332167567
## 565             1 -25.952284748   1  -7  -4.240306922       565  2.032965331
## 566             1  24.418679851   0  -5  -9.387215382       566 -0.697689487
## 567             1  29.548028094   1   5 -12.930706531       567 -0.989153292
## 568             1 -24.155964855   1  -1   0.540336137       568 -3.493163560
## 569             1 -18.838691133   1   0  -9.191709732       569 -0.207799940
## 570             1 -11.670314415   1  -6  -0.615426290       570  0.075888196
## 571             1  22.681806521   1   2   0.737439394       571  0.069961899
## 572             1 -35.143441156   1   1  -7.153232049       572 -1.964817504
## 573             1   1.046231391   1   3  -0.529940425       573 -1.425449696
## 574             1  -9.133375155   1   5   2.515909546       574  1.270881918
## 575             1 -19.750931472   1   5  12.563402244       575 -5.136759060
## 576             1   0.425324725   1  -1  -2.206906749       576 -2.794209603
## 577             1   7.080159101   1   4 -10.894995647       577  7.684173150
## 578             1 -23.221235789   1   7  -3.390744840       578 -6.604638764
## 579             1  -9.516685650   0  -7   3.644236329       579  5.630227328
## 580             1   2.307130921   0  -6  -6.217163683       580  1.450361735
## 581             1  -7.337737106   1  -4  -2.941359885       581  4.929989945
## 582             1   3.042185473   0  -3   5.542695656       582 -6.119699052
## 583             1 -16.248624061   1   7  -2.511899938       583  0.377180345
## 584             1  15.658643201   1  -7   1.736493432       584  0.962798216
## 585             1  -8.936129638   0   7   8.740733974       585  6.263877947
## 586             1  -4.900607779   1  -6   2.314276707       586 -2.172237100
## 587             1  18.928414241   0  -5   3.839348537       587 -1.577386529
## 588             1  -3.959327412   1   5   2.255081164       588 -2.328796202
## 589             1  33.537733423   1   3  -1.762886124       589  5.762332891
## 590             1  -0.866363714   1  -5   1.832272429       590  2.261888111
## 591             1   4.084852703   1   5  -2.440624359       591 -3.955190168
## 592             1  16.744395366   0  -3   4.366950072       592  8.816867643
## 593             1  31.932924222   0  -2   8.993903110       593 -2.521396018
## 594             1 -26.800732102   0  -2   5.253976071       594 -0.948469701
## 595             1   2.015702512   0  -6  10.876153198       595  4.785676993
## 596             1 -36.281881898   0  -6  -0.917814314       596 -0.761538957
## 597             1  24.774635899   1   4  -3.433965011       597  0.455809501
## 598             1  13.041438377   1  -1   0.126460921       598  3.125236213
## 599             1  15.465301784   0  -6   3.885156112       599 -0.072652893
## 600             1   8.400850524   0  -4   6.790230334       600 -1.023160516
## 601             1  -6.279683756   0   3   1.559727122       601  1.244384498
## 602             1 -11.819624726   0  -4   4.955988425       602 -2.268506109
## 603             1 -18.994375685   1   6  -1.909493231       603 -3.579096260
## 604             1  -6.297398260   1   3  -4.401307576       604  4.909724872
## 605             1  -3.011687383   1   0   0.115341607       605 -0.406836929
## 606             1 -29.485211531   1   0   3.541983978       606 -3.266316741
## 607             1 -11.107173069   0  -5   7.355355837       607  2.209837640
## 608             1  36.791068621   0   1   7.115925535       608 -1.610641247
## 609             1  -3.264221923   1   6  -3.658786263       609 -0.617789300
## 610             1   3.055846647   1   6  -6.139311590       610  1.132051143
## 611             1   8.570869326   1   2  -4.232786132       611  1.664661207
## 612             1 -43.364986826   0  -3   8.207391213       612 -7.438056511
## 613             1 -31.896551437   1  -4 -12.684913634       613  3.073596296
## 614             1   1.065252778   1   7   7.574834093       614  0.726789290
## 615             1 -15.175888617   0  -4   6.722326437       615 -0.682243855
## 616             1   3.863583059   0  -6   1.711089976       616  1.278652586
## 617             1  20.673811077   1  -7  -1.623352406       617  1.443232087
## 618             1   4.403104834   1   6  -4.754606301       618 -5.426747198
## 619             1  14.971947838   1   5   0.612981617       619  2.438397202
## 620             1  32.849364825   0   5  11.047987729       620 -8.688908020
## 621             1   4.208196693   1   6   6.524640182       621 -2.875985209
## 622             1  -7.806908397   1   6  -2.627570331       622 -2.754809863
## 623             1  29.141125223   1  -6   5.758179613       623 -5.106535657
## 624             1   1.352593717   1   5   4.130834483       624  5.956217603
## 625             1   7.270533842   0  -6   5.781224160       625 -4.157845925
## 626             1  10.230502919   1  -2   2.054863472       626  1.507521603
## 627             1  14.149918373   1   4  -0.280243218       627 -0.944134248
## 628             1  15.546317262   0  -2   3.900885674       628  7.904792291
## 629             1  17.235276314   1   4  -9.615888919       629 -0.998891258
## 630             1  -9.224844747   1  -5  -1.368247441       630 -6.307521687
## 631             1  26.068403544   1   5   3.018588658       631 -1.137699114
## 632             1  25.747779354   1   4  -6.586741567       632  3.020988214
## 633             1  24.443052743   1   7   4.777943150       633  2.383514713
## 634             1  -9.724938401   1   6   0.064267691       634 -0.091394116
## 635             1 -11.537234855   1   1   1.457180470       635 -2.129936390
## 636             1 -16.970549372   1  -1 -10.084027793       636  0.991817641
## 637             1  22.454136857   1  -1  -7.982442018       637  0.314312806
## 638             1  -8.411123764   0  -3   3.952566512       638  1.805984662
## 639             1  -5.984993452   0  -1  -0.279768856       639 -1.559962383
## 640             1 -13.591289536   0   0  13.923282690       640  2.340151229
## 641             1 -30.639279570   0  -2   1.977058522       641 -3.655141348
## 642             1  -4.044362268   0   0  13.489812094       642  2.084691151
## 643             1  34.266691288   0   6   5.283050665       643 -0.773087055
## 644             1   6.191550938   0   0   8.909346560       644 -3.401554731
## 645             1 -18.497402289   0   0   6.593440194       645  0.750516360
## 646             1  -0.585789599   1   2  -3.620415769       646 -1.697897458
## 647             1  -2.297789078   0   1  -1.316455829       647  0.950722541
## 648             1   2.373149280   1  -5  -1.742953775       648 -2.158584263
## 649             1 -25.069849122   1  -1  -1.066259502       649 -3.468056481
## 650             1   5.162001629   1   7   3.743332331       650  0.341706064
## 651             1   4.361125699   1   3 -15.022749173       651  2.087725546
## 652             1 -11.063839732   1   7  -1.028550975       652 -0.504428015
## 653             1   4.597284195   0   2  -0.717583592       653  2.923735450
## 654             1  -8.957412631   1  -4   6.554785984       654  2.680452373
## 655             1 -29.213445070   1   6   3.822866304       655  0.761080210
## 656             1   4.248350079   1  -1   2.822342861       656  2.092055599
## 657             1  19.909549509   1   6  -0.876606925       657 -2.129275800
## 658             1   8.293445934   1   5  -0.012873428       658  0.518038313
## 659             1 -11.382975781   0   1   2.943749690       659 -2.303247334
## 660             1  18.513681876   1   7   0.876585790       660  2.415427452
## 661             1  26.814618206   0   2   6.984614855       661 -4.329210427
## 662             1 -19.913949859   0  -3   4.127083239       662  2.115847972
## 663             1 -13.570663935   1   3   2.773652553       663  0.759712554
## 664             1 -19.474369399   1   5   2.726889913       664  0.703951120
## 665             1   9.525590359   0   0  -2.681922520       665 -5.056259089
## 666             1  20.758198097   1  -7  -9.196337545       666  0.234993866
## 667             1 -40.853088266   1   0  -3.561710969       667 -0.896823949
## 668             1   8.242141157   1  -5  -0.472832279       668  1.658425320
## 669             1   3.263507794   0   2   1.215079844       669  1.858769438
## 670             1  -6.685541461   1  -1   5.102230529       670 -1.038065822
## 671             1  59.850525008   0   4   4.616309957       671  0.008700561
## 672             1 -27.870669696   1   2   1.315899702       672 -1.877559504
## 673             1 -29.357438814   1  -4  -3.437930942       673 -3.226042176
## 674             1  -0.844412832   0   5   1.104305900       674  1.107513091
## 675             1 -17.475088574   1   0  -2.770518887       675  3.274182301
## 676             1   6.204143519   1  -3   2.759231940       676  4.690399173
## 677             1  -7.934523313   1   5  -0.809164926       677  2.432101552
## 678             1 -26.169593175   0  -7   6.757017673       678  0.167977058
## 679             1  20.650528726   1  -7  -0.026665711       679 -3.210270912
## 680             1  -7.777829765   1   2 -12.025667203       680  5.978316277
## 681             1   1.995534297   1   6   4.949046794       681  3.722651172
## 682             1  -8.144763130   1   0  -6.165158399       682  2.201594458
## 683             1  -7.919147770   0  -7   0.137052812       683 -0.519286383
## 684             1  34.793968833   0  -3  -1.955026167       684  0.012304553
## 685             1  -9.268298680   0  -3  12.009050887       685  3.999535350
## 686             1  -8.351226302   0  -5   2.770979994       686  1.695290583
## 687             1 -46.512903550   1   1   2.803879726       687  5.120205666
## 688             1  -9.897295873   1   2  -7.153944041       688  5.255429005
## 689             1  -6.707278271   0  -2   1.282654684       689 -0.446424903
## 690             1  41.731786867   1   6   2.237564921       690 -0.393518693
## 691             1 -27.593686192   0   7   4.545772751       691  4.215871188
## 692             1   5.929213532   1   0  -3.288058960       692  3.760564551
## 693             1   7.546184955   0   5   1.273906463       693  5.898982941
## 694             1  -8.634846672   1   2  -6.814797120       694 -4.463601875
## 695             1  15.212872662   1  -3  -5.758290220       695  0.539295292
## 696             1  14.555633828   1  -7   0.461780927       696 -1.092427318
## 697             1  21.061393831   1   7  -2.630545176       697  0.605607332
## 698             1  -9.664213665   1  -6  -4.706412490       698 -0.344651545
## 699             1 -32.426511374   1  -5  -5.801925593       699 -6.800818960
## 700             1  19.040014220   1   4   7.316083942       700  1.679541693
## 701             1 -15.541414358   1   2  -6.135178685       701  5.377585652
## 702             1 -13.634563713   1   7   6.445611528       702 -0.318412694
## 703             1  15.807201087   0  -3   2.609917072       703  3.137381835
## 704             1 -15.455142563   0  -7   0.193234559       704 -0.057976859
## 705             1 -50.756857573   1   1   0.027404501       705  1.419275379
## 706             1 -17.223198364   0   5  10.254813557       706 -3.007920271
## 707             1  -0.592599802   0  -1   2.920017440       707 -0.786281440
## 708             1  11.109530316   1   6  -3.598201883       708  1.070900067
## 709             1   1.881688678   1  -3   1.586459872       709 -1.537680116
## 710             1  21.468817114   0  -4   0.833578270       710 -2.778867554
## 711             1 -36.490328010   0  -3  12.513392663       711  3.094307766
## 712             1  -4.391830640   0  -3   2.002043547       712  2.125015563
## 713             1  25.115266685   0  -2   5.360559927       713  0.746303585
## 714             1  14.513112708   1  -4   0.275004815       714  6.916181568
## 715             1  13.607011548   1   6  -0.306819301       715 -3.760390831
## 716             1   6.324970956   1   6  -5.382217954       716  3.412740070
## 717             1 -16.608736110   1   6   2.521521465       717  3.065230579
## 718             1  30.873757832   1   5   2.681122749       718  5.481024227
## 719             1  -5.316256542   1   1  -6.774886458       719 -1.499259653
## 720             1  27.390927793   1  -2   1.137531536       720  2.988636462
## 721             1  -2.082229876   1  -4 -12.647093541       721  2.672474974
## 722             1  -1.238750266   1   1  -4.092655717       722 -1.238830481
## 723             1  -2.899649046   1  -3  -7.947649968       723  4.141249477
## 724             1 -11.736083919   1  -7  -5.412466617       724  0.057606328
## 725             1 -57.060529124   1  -3   1.029021799       725  4.129714245
## 726             1  21.188997565   1   3  -5.405911599       726  2.302475735
## 727             1   6.518974574   1   2  -0.943962771       727 -1.103000318
## 728             1  -4.154969069   1   0   1.563782269       728 -0.750132313
## 729             1 -10.375205625   1   3   0.373584993       729  1.545590849
## 730             1  -8.589398628   1  -7  -3.021724271       730 -0.963466478
## 731             1   4.697872016   1   2   0.226233248       731 -7.229608620
## 732             1  24.397640135   0   7   8.242698186       732  0.702418080
## 733             1  -3.526186345   1   7  -3.778414920       733  7.611647058
## 734             1  22.464219412   1   7  -1.440896764       734  0.204961172
## 735             1  31.751462188   1   2  -3.650752530       735  1.788636346
## 736             1   0.268161811   1   0   2.620742601       736  3.063150114
## 737             1 -10.113204337   1  -3  -2.958502646       737 -2.592415444
## 738             1  18.555802352   1  -3  -5.919559931       738 -1.245082049
## 739             1  -1.829610898   1  -3   6.001765659       739 -6.843059480
## 740             1   5.549567335   1   5  -4.540356635       740  3.069461952
## 741             1 -12.740942878   1  -1  -7.744088434       741  0.402581048
## 742             1  36.178378568   1  -4   2.040224764       742 -2.520209970
## 743             1  18.781203097   1   6   4.361714103       743 -1.863037002
## 744             1 -27.712050346   1   2  -9.302388452       744  4.568223552
## 745             1  17.876991143   1   5   1.030588023       745  2.502700415
## 746             1   8.616913864   1   6  -5.327630005       746  1.919339650
## 747             1   0.624672824   1   7   2.221395553       747 -0.003593688
## 748             1 -33.416787740   1   5   1.950318696       748  4.349697517
## 749             1 -39.038230379   1   0  -0.597022551       749  0.870878729
## 750             1  41.614503361   1   4   7.374813242       750 -2.186248924
## 751             1   4.201255348   0  -1  -6.473291305       751  2.546321869
## 752             1   0.463698842   0   2  -7.338063413       752  0.639224532
## 753             1 -37.491744540   1   3  -1.056503652       753  0.391490563
## 754             1  13.568645281   0  -1   2.098093627       754  3.155547957
## 755             1 -13.987392227   0  -2   7.166701779       755 -0.942255799
## 756             1   5.628283881   1   1  -3.687340567       756 -3.397078205
## 757             1 -18.094662744   0  -3  -0.625076714       757  3.483666765
## 758             1 -50.886760111   1   2   1.399120735       758 -1.245552540
## 759             1 -24.506731969   1   4  -2.236860130       759  1.754325339
## 760             1  21.613113628   1   1   0.239788400       760  0.379436984
## 761             1 -29.263363939   0   1   5.464262675       761  1.632584577
## 762             1  -2.243373793   0  -4   2.540048379       762  3.247792463
## 763             1  12.387384748   1   5  -4.349553997       763  2.313721851
## 764             1  12.629182946   1  -7  -3.513382379       764 -0.355599114
## 765             1  10.040816583   1  -7  -6.797036350       765  4.519441927
## 766             1  20.292945201   0  -4   1.728234956       766  2.831189803
## 767             1   1.793729847   1   5   2.109027403       767 -4.744779547
## 768             1  -1.859466434   1   2   1.633558328       768 -2.236807086
## 769             1 -29.457646819   1   5  -4.817875557       769 -0.365431152
## 770             1  25.006632895   1   0   2.393289888       770  0.276885303
## 771             1   7.721891268   1  -4   3.902824583       771  1.246585067
## 772             1  24.798680561   0  -3  -0.927127497       772 -1.925322521
## 773             1 -10.264585378   0  -5  -1.857078810       773 -2.059067240
## 774             1  -1.683744599   1   2  -1.879532076       774  4.692172902
## 775             1 -21.992383813   1   4  -5.248865462       775  0.841855522
## 776             1 -11.083449723   1   4   0.677549433       776 -1.896957914
## 777             1 -17.631015862   1   3  -5.209671246       777  4.643310419
## 778             1  25.206541171   0   1   6.857717307       778 -1.960096036
## 779             1 -33.754812817   1   0   1.058615101       779  1.792524878
## 780             1 -20.142339558   0  -5  -2.099951456       780  4.381125119
## 781             1   9.336058252   1   5  -2.572908713       781  4.422379551
## 782             1  -9.074649573   1   6  -2.405945704       782 -5.569025043
## 783             1  12.316097739   0  -3  -0.114480583       783 -1.558336723
## 784             1  -9.292201384   0  -4  -1.179086875       784  2.701832515
## 785             1  15.462366560   0  -4  -0.214579443       785  1.220403772
## 786             1  24.561383547   1  -5  -2.667360926       786  8.279956147
## 787             1  31.380003484   1   3   1.204618284       787 -0.367716336
## 788             1   4.728441220   1  -6   1.378939997       788  2.770867954
## 789             1  -3.891083468   1   1   4.936689423       789  3.449079695
## 790             1 -22.592984444   1   7  -3.035875092       790  5.982213015
## 791             1 -30.610728341   1   5  -1.215776172       791  1.300576813
## 792             1   0.989411750   1   0  -0.918594481       792 -2.908352950
## 793             1   7.452747530   0   1   2.822745562       793 -2.584161937
## 794             1  23.768449246   1   5  -3.228316527       794  1.514472395
## 795             1  39.505427054   1   2  -9.920987584       795 -1.461049403
## 796             1 -11.991228064   1   2  -4.664053968       796  2.471107025
## 797             1   4.021397347   0   7   9.322798748       797 -0.498152535
## 798             1 -22.952406697   1  -5   1.538341672       798 -2.796144844
## 799             1   8.621854347   1  -7  -9.367384037       799  4.142052124
## 800             1 -14.688343160   1   1  -6.907931863       800 -2.344945348
## 801             1  12.871805728   0  -1   6.071405024       801 -3.708975309
## 802             1  13.309384166   0  -5   0.707409518       802 -0.840674830
## 803             1   1.704396261   1   3   0.315519253       803  3.536460150
## 804             1  47.811435260   0   1   6.177306337       804 -5.555057737
## 805             1  18.252423275   1   5  -5.482731072       805 -4.593770768
## 806             1 -38.234337488   1   6  -3.393918531       806 -3.700395223
## 807             1   2.401389132   0   5  -2.799336189       807 -3.107013067
## 808             1 -16.136497595   1  -2   1.077381781       808 -0.834346511
## 809             1  16.523591389   1   1 -10.277176006       809  2.819007615
## 810             1   5.858196259   1  -6  -2.602043141       810  0.574582956
## 811             1 -42.820244801   0  -6   1.261319528       811 -2.622874545
## 812             1   6.776827931   1  -6   0.765782121       812  1.294946277
## 813             1  27.463759381   0  -5   2.460622122       813 -0.582173801
## 814             1   2.111067696   0  -1   2.669111909       814 -0.820523164
## 815             1 -27.930543947   0  -6  -1.035661810       815 -4.874708943
## 816             1  34.496682895   1   4   0.175599495       816  3.112542455
## 817             1 -13.796633156   1  -1  -5.604880124       817 -0.468865462
## 818             1 -19.941070574   0   1  -1.659829359       818  2.315781472
## 819             1  29.260605796   1  -1  -6.386275635       819 -1.725169418
## 820             1 -19.655585381   0  -5   2.693256638       820 -1.109451644
## 821             1 -18.545632398   0   7   6.744206419       821  4.492248795
## 822             1 -18.312829023   1   5   3.005557511       822 -4.654349072
## 823             1 -24.434039242   1  -7  -0.063597931       823 -1.536626840
## 824             1 -23.406947833   0   4   3.580694978       824 -3.601682779
## 825             1   5.719019321   0  -6   2.028172099       825 -0.071416361
## 826             1 -10.344850678   1  -3  -1.774261469       826 -2.220251127
## 827             1  13.055026054   1  -2  -4.156280408       827 -0.697898377
## 828             1 -44.527992231   0  -1   1.125114510       828 -3.184810862
## 829             1 -22.471932272   1  -4  -3.442059348       829  1.433908276
## 830             1  16.470052953   1  -1  -3.243938072       830  0.701516016
## 831             1  14.822186141   0  -5   2.358656728       831  0.610820505
## 832             1  15.261072080   0  -5   0.873031239       832 -2.091101456
## 833             1 -14.436983618   1   4  -1.966374344       833  0.521933666
## 834             1  10.785867280   1   1   0.616202890       834 -1.296140764
## 835             1 -18.385796828   1   1 -10.033338249       835  1.064837701
## 836             1 -17.582352418   0  -4   6.048917700       836  3.935904657
## 837             1 -16.178681920   0   2   0.701862208       837  3.483019657
## 838             1   5.698693523   1   0  -4.379139370       838  2.248986707
## 839             1  32.249126809   0   5   4.811823675       839 -5.566245888
## 840             1 -27.540028389   0   1   5.951338463       840  0.050827379
## 841             1 -23.098200045   0  -2   0.959889088       841  2.975609486
## 842             1 -35.114569645   1   7  -9.138002655       842 -0.971374285
## 843             1  37.956353850   1  -4  -4.200384099       843 -0.677498883
## 844             1  -6.048295744   0   3   8.116969948       844 -0.270103902
## 845             1 -14.974223666   1  -1   0.800207269       845 -0.547396783
## 846             1 -31.073773801   0   2   8.000756347       846 -0.080624803
## 847             1   0.004788214   1  -2   0.571373540       847 -0.154280002
## 848             1  -2.115682455   1   3  -6.827122747       848 -1.261795801
## 849             1 -22.455178212   0  -5   2.192693672       849  4.226446638
## 850             1 -28.559515180   1   4  -4.807536360       850 -1.111943510
## 851             1 -15.020472601   0   1   6.081744433       851 -0.379142604
## 852             1   2.410524833   0  -6  -3.903559860       852  3.314159401
## 853             1   0.796458040   1  -3  -3.411466049       853  4.117653823
## 854             1   1.755496518   1  -3 -11.217131853       854  1.791930439
## 855             1   3.970950873   0  -2   0.901381965       855  1.102763135
## 856             1   5.246420798   1  -5   3.907533176       856  0.687142840
## 857             1 -25.159060159   1  -1  -4.786156706       857 -3.457132100
## 858             1  40.226492099   1   1   0.457035273       858 -1.180890230
## 859             1  -1.512336924   1   4  -2.078709891       859  0.577957044
## 860             1 -42.090880448   1   5  -3.307117762       860 -0.240422416
## 861             1   9.408360162   0   7   3.419858175       861 -1.687709736
## 862             1  -7.702293998   0   4   4.461921706       862 -3.560555820
## 863             1 -35.407003125   1  -1  -3.634537837       863 -5.983654451
## 864             1  -3.229128463   0  -5   5.078720130       864 -1.455935591
## 865             1 -24.012841053   1   3  -1.977967157       865 -3.163669145
## 866             1   4.853658410   0   1  10.797555340       866 -6.033954910
## 867             1  12.958255724   0  -1  11.077907890       867  1.676583106
## 868             1 -17.357581757   1   1  -1.278585347       868 -0.146687651
## 869             1 -21.922679744   0  -3   1.511458502       869  5.923034883
## 870             1   9.158394940   1   5  -3.329025252       870  3.154762097
## 871             1  20.579719585   1  -1  -3.480306035       871 -0.316517728
## 872             1  14.241152455   0  -6   3.416709668       872  1.036503827
## 873             1  -7.210545030   1   7  -2.946374351       873  1.709296300
## 874             1 -19.665418721   1  -4  -0.370467747       874  1.993851248
## 875             1 -12.208493364   1  -6  -5.561355133       875  1.460482340
## 876             1  -5.502706707   1   0  -2.955311256       876 -2.838330413
## 877             1   3.615310466   0  -6   4.284723048       877  1.695995987
## 878             1   5.895283362   1   7   2.544466949       878 -0.441488356
## 879             1  15.561449321   0   1   1.985519491       879 -1.190258491
## 880             1   9.370462203   0  -5   1.933786257       880 -0.504239054
## 881             1  11.792875055   1   0   1.850680923       881 -2.562035898
## 882             1 -18.165604279   0  -7   5.183615490       882  2.000650161
## 883             1 -15.666534324   0   1  11.793825290       883 -0.795251215
## 884             1   9.916050058   1   0  -4.568010376       884  2.594662101
## 885             1 -19.457915229   1   3   6.377997115       885  1.353816289
## 886             1   9.325950866   1   4  -7.195450459       886 -3.028439455
## 887             1  -2.574640665   0  -6   0.190711035       887 -6.069914164
## 888             1 -14.325383457   1   6  -2.517860975       888 -1.367136916
## 889             1 -12.075177715   1  -4  -8.565887226       889 -3.348603448
## 890             1   4.004885808   0   3   3.762258153       890  0.535363699
## 891             1 -14.695229092   1  -2  -4.417010547       891 -5.639133092
## 892             1   8.045303064   0  -5   0.682593893       892  0.112912808
## 893             1 -11.797065400   1   0   2.131409424       893 -3.097236843
## 894             1 -14.043674640   1  -2  -4.625698753       894  3.966140165
## 895             1  -2.573547300   1   5  -4.845478682       895 -3.578147876
## 896             1  24.539768028   1   4  -1.132913892       896  6.665472959
## 897             1   7.090278186   1  -7  -4.891119818       897 -2.818440342
## 898             1  -3.091110338   1  -4  -2.666770640       898  4.049632890
## 899             1  27.339976395   0  -5   4.331964668       899 -1.356686583
## 900             1  -3.520589393   1  -3  -2.466904707       900 -3.247704769
## 901             1 -10.387096796   1   0  -6.143079194       901  2.480605506
## 902             1 -21.694165200   0  -7   6.422815750       902  0.400672213
## 903             1  -1.667355181   1   6   6.058900861       903  0.203159863
## 904             1 -26.631700420   0  -1   0.089335375       904  2.031541508
## 905             1  10.681919889   1  -6   8.501238880       905  0.290700013
## 906             1  18.911512339   1   6   6.720576765       906  3.377882664
## 907             1 -37.856407789   0   4   1.236247514       907  0.874628601
## 908             1  17.663081662   1  -6  -0.479254973       908  0.989583883
## 909             1  -9.010695937   1   3  -0.529473854       909  3.476805623
## 910             1   5.036178619   1   5  -6.305773568       910  3.460204366
## 911             1  -4.427604463   1   6  -8.801051997       911  6.915436307
## 912             1 -17.088997625   0  -6   8.957154061       912 -3.946420658
## 913             1   1.684973125   1   2  -4.005172892       913  0.155485720
## 914             1   9.848884337   0  -7  -2.789665902       914  1.595510627
## 915             1   7.190214879   0   2   3.780587444       915  5.305466503
## 916             1 -19.666636161   1   6  -0.328271401       916 -0.648982012
## 917             1   0.775218023   0   4   7.019125505       917 -1.380075601
## 918             1   1.253233814   1  -1   4.030615543       918  0.409659823
## 919             1 -24.759664563   1   4  -0.626672191       919  0.693961588
## 920             1  -8.979997317   1   7  -2.822112720       920  4.006254868
## 921             1 -18.172168634   0  -2   4.703889962       921  1.577027198
## 922             1 -25.837131716   0   0   6.349044729       922 -0.793775008
## 923             1  18.476632347   1   5   6.325192452       923 -1.787076500
## 924             1  14.824015832   1   0   3.497652009       924 -2.332953970
## 925             1   3.056214655   1   0  -2.394661437       925 -3.904249257
## 926             1   6.538417148   1   7  -6.571013705       926  5.501517400
## 927             1  -9.665281039   1   2  -4.518229093       927 -1.454378442
## 928             1  -5.475257198   1  -6  -3.413388843       928  3.613363440
## 929             1   7.353031670   1   5   1.872998291       929  2.739355968
## 930             1 -15.723741080   1   4   4.048922108       930  7.433059276
## 931             1 -23.695072603   0  -2  -5.136559471       931 -4.012878816
## 932             1  21.705357843   0   2   1.356023580       932 -3.104227394
## 933             1  -4.564621256   1   2   1.860816928       933 -2.297315484
## 934             1 -19.802108017   1   0  -7.369017571       934 -3.188237854
## 935             1 -42.514659070   1   1  -7.113995746       935  0.304085622
## 936             1 -25.081598894   0   0   3.428573754       936  0.732817667
## 937             1  -2.591969004   1  -3  -7.870950907       937  5.211668715
## 938             1   8.749540264   1  -1  -5.196013186       938 -4.010054012
## 939             1 -15.179083103   0   5   5.314463461       939  2.845758156
## 940             1  15.131325244   1  -1  -6.396904457       940  0.304174088
## 941             1 -12.917248992   1   3  -8.826299801       941 -1.509978822
## 942             1 -20.701509995   1  -6 -11.973424717       942  2.129713732
## 943             1  -0.574319086   1   6   3.309257098       943 -0.754534957
## 944             1 -28.195381234   0  -4   1.396010697       944 -0.759857468
## 945             1   0.565608626   1  -6  -9.878990665       945  1.314963170
## 946             1 -26.882439943   0   0   2.482534665       946  0.501561864
## 947             1  -5.027070486   1   3   6.765158174       947  2.793517851
## 948             1  52.513371065   1   7   3.276741424       948  3.129729337
## 949             1 -15.781306087   1   5   1.043741845       949  0.885175736
## 950             1 -11.049732160   1  -1   2.396522624       950  2.958075159
## 951             1   5.146108399   1   6   5.598609446       951 -4.517475337
## 952             1  -2.813515141   0   6   1.794850294       952  5.622531230
## 953             1 -27.737955292   0  -5   5.992917446       953 -0.761622307
## 954             1  15.313989243   0   0   2.815097324       954 -0.451110937
## 955             1 -13.264300663   1  -4   1.801821553       955  5.349969257
## 956             1 -43.919518515   1   0   4.126819948       956  4.497582793
## 957             1 -48.801006250   1  -4  -0.201962921       957  4.644705056
## 958             1   0.698882555   1  -4  -3.027283681       958  0.051477975
## 959             1  -8.098343135   0  -7  -2.393244191       959  2.900849020
## 960             1   2.823428434   1   4  -8.135184331       960  2.621943368
## 961             1  -3.118576972   0   1   5.479866757       961 -2.536888369
## 962             1   6.710511556   1   2  -5.653078851       962  1.422910666
## 963             1 -10.755604105   1   3   3.747037986       963  0.964838021
## 964             1  -4.220933876   1  -5   2.535204516       964  3.962110080
## 965             1  13.100366202   0  -4   0.746138839       965 -7.404555628
## 966             1  -6.231193978   0   2   3.784574468       966  1.290904516
## 967             1  -7.081779149   1   2  -2.136122648       967 -0.765405508
## 968             1 -21.177923743   1   1  -6.739155605       968  1.685755779
## 969             1 -18.524024599   0   2   0.412977926       969  3.760092725
## 970             1  21.457142284   0   4   5.768640580       970 -0.120764761
## 971             1  -0.656672406   0   4  12.934209501       971  2.842551291
## 972             1  42.290103787   1  -3  -8.206844123       972 -1.653077152
## 973             1   2.903273401   1   7   0.927883767       973 -2.356634166
## 974             1 -21.814428678   1  -6   4.040697526       974 -1.067451042
## 975             1 -12.826256521   1  -2   2.313709508       975  3.638540350
## 976             1  11.894010711   1   3  -0.704314738       976  2.018231174
## 977             1  48.856997196   1   5  -9.150836502       977 -2.993885657
## 978             1   4.951188295   0  -2   2.020439386       978 -0.438623239
## 979             1  21.346868813   0  -1  -0.068753371       979 -2.595854245
## 980             1  -8.211538983   1  -3  -0.511206188       980  5.645519980
## 981             1 -18.177957900   1   5   1.012272223       981 -1.385610932
## 982             1  -9.790658330   0  -2   4.016602740       982 -2.157437830
## 983             1  35.412149255   1  -6  -2.938966227       983 -1.995949788
## 984             1  -9.208231229   0  -3  10.554792720       984  3.963980693
## 985             1   4.930128018   0  -7   0.531595903       985 -2.606651903
## 986             1 -19.994985872   0  -5   1.443191410       986 -7.119843586
## 987             1   5.530327992   0  -6  -4.140645108       987  0.544653781
## 988             1 -11.031243344   1   4   1.296230942       988  3.243390007
## 989             1  22.469278159   1   6  -2.476153692       989  3.268034108
## 990             1  23.275513249   1  -1  -1.447508337       990  3.008368508
## 991             1   3.606519682   0  -7   3.203419766       991  0.475489713
## 992             1 -13.051738850   1  -3  -0.707314307       992  0.190142994
## 993             1  11.007686809   1  -3  -4.399236930       993 -5.839361097
## 994             1 -55.508626717   1  -3  -2.618198057       994  0.076499123
## 995             1   6.153852117   1   7  -3.067182982       995  2.919770493
## 996             1 -32.217410375   1   4   3.747461994       996  0.438158654
## 997             1  26.529532628   0   6   6.840808978       997 -3.161951466
## 998             1  23.885768774   0   4   4.045209623       998  0.439767025
## 999             1  21.301854304   0  -6   5.260219579       999  2.449137045
## 1000            1 -32.249578725   0   5   7.883550864      1000  1.930632169
##      fixed_outcome random_effects achievement
## 1         44.14824              0    44.89826
## 2         47.02168              0    47.13993
## 3         46.58302              0    45.39498
## 4         53.25302              0    61.06807
## 5         60.35261              0    60.26452
## 6         50.33459              0    52.04749
## 7         47.95437              0    43.74895
## 8         65.77028              0    68.10847
## 9         54.99399              0    52.86196
## 10        68.45075              0    70.85385
## 11        69.38440              0    71.18391
## 12        55.27415              0    54.78332
## 13        45.04488              0    44.41233
## 14        56.82243              0    62.36800
## 15        57.14577              0    58.68820
## 16        49.30268              0    47.56502
## 17        61.39283              0    67.32001
## 18        62.52824              0    61.08948
## 19        42.84844              0    44.10592
## 20        52.27471              0    50.78501
## 21        57.95229              0    55.62754
## 22        69.16136              0    65.02915
## 23        51.40252              0    49.44879
## 24        44.38189              0    41.89938
## 25        51.25896              0    51.25562
## 26        48.28208              0    43.80349
## 27        56.34919              0    56.68883
## 28        42.95798              0    39.90636
## 29        31.66789              0    27.08379
## 30        59.33988              0    60.31983
## 31        45.43627              0    45.21331
## 32        62.47393              0    63.27686
## 33        49.47315              0    53.18695
## 34        37.38304              0    35.23169
## 35        61.73175              0    62.54855
## 36        54.22303              0    55.53830
## 37        47.19819              0    40.79568
## 38        60.46500              0    57.24596
## 39        64.32286              0    63.82651
## 40        56.72858              0    53.07608
## 41        59.02734              0    60.71119
## 42        60.55214              0    59.85172
## 43        44.20059              0    42.51587
## 44        74.85213              0    74.05284
## 45        45.65533              0    47.33055
## 46        30.42572              0    29.29042
## 47        48.65263              0    54.46262
## 48        59.08134              0    57.44676
## 49        45.08281              0    46.42301
## 50        51.70844              0    47.14405
## 51        54.42664              0    49.54734
## 52        65.81418              0    65.54199
## 53        40.76114              0    40.67189
## 54        61.70339              0    58.13193
## 55        35.30767              0    32.93538
## 56        46.32323              0    45.34992
## 57        49.45613              0    47.43884
## 58        34.40492              0    36.85220
## 59        52.27684              0    52.64044
## 60        51.01246              0    52.47826
## 61        47.70112              0    43.49180
## 62        49.96671              0    51.77837
## 63        54.06822              0    54.44741
## 64        34.88747              0    35.30175
## 65        45.24513              0    43.52565
## 66        55.14944              0    56.30622
## 67        57.16962              0    61.96842
## 68        60.47677              0    55.68057
## 69        60.68110              0    62.60154
## 70        37.44791              0    38.71852
## 71        51.32501              0    49.37395
## 72        53.28960              0    52.99981
## 73        44.99474              0    47.20987
## 74        68.34462              0    70.72875
## 75        42.91696              0    42.53229
## 76        58.40890              0    60.97260
## 77        39.33492              0    39.74953
## 78        51.93147              0    52.96945
## 79        52.21463              0    51.69304
## 80        57.80635              0    59.09159
## 81        46.84428              0    49.28870
## 82        40.13332              0    41.80356
## 83        52.65300              0    52.06287
## 84        62.28376              0    67.19565
## 85        50.80342              0    53.16167
## 86        48.75980              0    46.15464
## 87        47.59790              0    44.40752
## 88        49.93162              0    50.61322
## 89        49.90505              0    47.87713
## 90        53.96511              0    56.29858
## 91        58.57158              0    57.94191
## 92        56.51408              0    58.76592
## 93        40.52736              0    45.37310
## 94        59.14300              0    55.11709
## 95        57.32931              0    57.18817
## 96        63.86672              0    61.15069
## 97        47.53604              0    48.27931
## 98        54.75818              0    57.87784
## 99        47.46495              0    49.09214
## 100       53.98362              0    55.24809
## 101       55.62814              0    56.16303
## 102       53.40179              0    52.13667
## 103       44.13563              0    40.80741
## 104       49.76956              0    53.24518
## 105       45.93507              0    47.14052
## 106       53.15065              0    56.03853
## 107       62.85081              0    60.57825
## 108       51.04371              0    54.09075
## 109       67.37561              0    72.06819
## 110       43.51687              0    39.30917
## 111       60.13022              0    62.39177
## 112       45.71981              0    44.27073
## 113       56.61282              0    47.91275
## 114       39.38106              0    38.70988
## 115       55.95718              0    55.87867
## 116       54.54300              0    56.28772
## 117       55.02782              0    51.67567
## 118       50.57273              0    49.26056
## 119       56.94255              0    55.64631
## 120       37.52753              0    37.22340
## 121       53.72123              0    53.25816
## 122       62.60906              0    63.92853
## 123       55.65942              0    58.78067
## 124       55.85737              0    56.41396
## 125       35.45387              0    38.70345
## 126       50.80036              0    46.12131
## 127       49.11268              0    49.06926
## 128       45.98827              0    45.67240
## 129       48.86363              0    46.98189
## 130       67.30713              0    67.18106
## 131       53.72835              0    50.96189
## 132       46.07346              0    50.43979
## 133       48.81547              0    48.32036
## 134       42.31283              0    41.54555
## 135       52.50110              0    51.97698
## 136       53.67917              0    51.65706
## 137       62.62940              0    61.01339
## 138       50.70149              0    52.90893
## 139       51.54117              0    55.99892
## 140       37.57328              0    35.93963
## 141       51.47726              0    50.33202
## 142       26.88268              0    26.20171
## 143       43.29331              0    40.51091
## 144       55.78474              0    55.37038
## 145       57.98840              0    52.24153
## 146       52.59928              0    55.26478
## 147       60.82753              0    58.39156
## 148       44.82409              0    42.97104
## 149       49.97961              0    48.26239
## 150       53.26270              0    51.54390
## 151       52.46478              0    51.29308
## 152       42.24743              0    45.77935
## 153       54.09687              0    54.12805
## 154       57.75551              0    64.75144
## 155       46.05272              0    40.14994
## 156       54.95536              0    47.74195
## 157       57.78852              0    54.55040
## 158       44.00888              0    46.43128
## 159       63.85235              0    62.63074
## 160       68.50025              0    65.51281
## 161       58.48994              0    58.10197
## 162       43.57639              0    45.71846
## 163       60.84389              0    54.33967
## 164       66.87069              0    71.32176
## 165       48.41331              0    50.37911
## 166       51.99216              0    56.89737
## 167       30.66221              0    32.81506
## 168       54.71080              0    56.61464
## 169       55.95370              0    58.20568
## 170       44.63004              0    43.11362
## 171       53.12468              0    56.94329
## 172       58.84134              0    54.97822
## 173       54.84115              0    52.71230
## 174       41.88954              0    47.63106
## 175       54.22659              0    53.50495
## 176       58.38466              0    61.75391
## 177       63.62103              0    61.70426
## 178       57.01341              0    53.78968
## 179       41.82646              0    41.98079
## 180       40.64138              0    36.57259
## 181       51.24368              0    45.20207
## 182       52.34995              0    52.62568
## 183       53.62803              0    48.16945
## 184       42.67742              0    42.86683
## 185       45.86698              0    45.80259
## 186       67.08854              0    67.78588
## 187       78.53084              0    72.88265
## 188       56.56248              0    53.39953
## 189       60.45637              0    60.97160
## 190       55.81154              0    63.10819
## 191       51.30639              0    53.74195
## 192       50.56260              0    52.62709
## 193       54.90392              0    52.50189
## 194       54.18742              0    54.79166
## 195       54.43148              0    54.49267
## 196       45.52549              0    42.09364
## 197       37.26811              0    38.70002
## 198       47.46201              0    45.27447
## 199       51.33118              0    52.06599
## 200       46.39659              0    44.65119
## 201       59.59401              0    60.45062
## 202       60.04914              0    62.59985
## 203       60.15898              0    63.44114
## 204       43.89257              0    44.71573
## 205       33.55627              0    31.18014
## 206       66.82320              0    66.54464
## 207       67.52079              0    67.70932
## 208       54.97557              0    57.36137
## 209       42.79550              0    41.87011
## 210       40.90118              0    37.26948
## 211       58.54014              0    60.41677
## 212       69.50752              0    68.73549
## 213       47.06267              0    43.08568
## 214       36.44965              0    34.32580
## 215       48.04826              0    49.07440
## 216       58.26184              0    66.31468
## 217       57.57297              0    56.87978
## 218       54.79902              0    60.64237
## 219       45.57656              0    43.90279
## 220       49.26918              0    54.52666
## 221       37.55551              0    34.37256
## 222       51.30912              0    53.96007
## 223       61.53337              0    61.55109
## 224       46.41554              0    53.19287
## 225       35.79543              0    33.42125
## 226       60.40040              0    63.97036
## 227       48.22953              0    50.66401
## 228       47.25153              0    45.78011
## 229       59.22825              0    53.71812
## 230       40.59544              0    45.31446
## 231       36.48404              0    36.44526
## 232       46.20411              0    43.83625
## 233       41.47215              0    40.13257
## 234       61.20773              0    58.82942
## 235       61.50815              0    59.28021
## 236       44.47902              0    44.98500
## 237       47.20111              0    45.58246
## 238       49.61383              0    46.39958
## 239       44.21766              0    49.56425
## 240       71.00643              0    69.44199
## 241       56.13433              0    57.01334
## 242       62.25776              0    60.64337
## 243       48.79330              0    48.72117
## 244       52.33270              0    47.48569
## 245       49.29784              0    50.65895
## 246       55.39882              0    53.40519
## 247       43.76367              0    44.80766
## 248       48.03633              0    47.19468
## 249       48.59225              0    46.34672
## 250       60.83994              0    64.16130
## 251       54.22824              0    52.15594
## 252       40.58022              0    37.95228
## 253       47.64193              0    47.28592
## 254       43.97968              0    49.71629
## 255       51.37537              0    44.03053
## 256       66.67183              0    64.53767
## 257       59.73286              0    61.19777
## 258       44.37188              0    41.00309
## 259       46.36048              0    47.73196
## 260       41.54147              0    40.98116
## 261       62.20799              0    63.04749
## 262       43.30949              0    44.45166
## 263       37.80282              0    37.79450
## 264       50.49496              0    45.85480
## 265       70.14832              0    65.54167
## 266       62.54287              0    56.48356
## 267       53.07421              0    49.37170
## 268       55.27980              0    53.13242
## 269       63.12272              0    62.56044
## 270       49.64300              0    52.54184
## 271       68.08827              0    71.55974
## 272       50.04412              0    49.04872
## 273       56.42737              0    55.81113
## 274       48.01599              0    48.20653
## 275       48.61983              0    45.49567
## 276       45.60343              0    45.79615
## 277       48.44568              0    44.84130
## 278       39.76054              0    35.39323
## 279       46.35430              0    48.98300
## 280       58.30853              0    58.54936
## 281       48.69147              0    48.26935
## 282       42.54471              0    46.86002
## 283       51.55821              0    52.23697
## 284       51.51561              0    53.62646
## 285       40.88692              0    39.85609
## 286       67.39124              0    63.35212
## 287       45.96383              0    45.64906
## 288       33.86981              0    28.01316
## 289       51.03893              0    52.96858
## 290       48.59984              0    47.89345
## 291       45.24593              0    43.36838
## 292       55.44769              0    58.71898
## 293       59.05882              0    63.28968
## 294       41.50002              0    43.51367
## 295       50.85872              0    54.55034
## 296       54.25294              0    50.86999
## 297       56.49767              0    59.39726
## 298       48.52740              0    49.57857
## 299       46.51123              0    43.56279
## 300       49.47182              0    46.95178
## 301       52.00533              0    51.94175
## 302       45.78054              0    46.90623
## 303       67.94755              0    70.36573
## 304       61.53130              0    60.20657
## 305       65.84627              0    65.21468
## 306       54.97430              0    50.52273
## 307       51.17988              0    47.27980
## 308       60.35703              0    63.58504
## 309       49.80739              0    47.65265
## 310       42.54292              0    40.93323
## 311       44.92976              0    42.68385
## 312       55.58156              0    57.46611
## 313       39.34398              0    42.20336
## 314       69.73418              0    73.77767
## 315       48.52206              0    44.60141
## 316       70.54999              0    73.95528
## 317       31.19370              0    35.88388
## 318       53.91045              0    47.34884
## 319       53.08093              0    52.01992
## 320       53.61131              0    58.07256
## 321       48.39753              0    48.76148
## 322       53.69185              0    56.26898
## 323       61.40924              0    60.88140
## 324       46.81240              0    45.61229
## 325       44.66851              0    44.96777
## 326       49.68907              0    50.44824
## 327       61.24383              0    59.24309
## 328       55.29694              0    53.91456
## 329       45.56691              0    46.09862
## 330       50.88780              0    54.23933
## 331       58.64845              0    58.20835
## 332       45.64331              0    43.48897
## 333       50.84582              0    52.78589
## 334       60.28036              0    58.01727
## 335       49.98545              0    54.06326
## 336       48.82434              0    51.58663
## 337       41.45899              0    45.77353
## 338       55.88719              0    56.88140
## 339       42.80757              0    42.80295
## 340       53.45337              0    50.79166
## 341       50.80175              0    52.29164
## 342       56.19502              0    56.88146
## 343       43.81689              0    42.91508
## 344       32.17842              0    30.73619
## 345       53.56847              0    58.18320
## 346       56.95185              0    52.21670
## 347       55.86355              0    56.58574
## 348       50.58322              0    48.54566
## 349       45.33572              0    44.39656
## 350       31.34609              0    33.86226
## 351       54.42661              0    56.74649
## 352       37.67583              0    32.97523
## 353       45.59126              0    43.75594
## 354       54.66398              0    50.30529
## 355       49.49703              0    51.45748
## 356       46.38414              0    46.16905
## 357       63.28795              0    58.76048
## 358       54.24580              0    55.17883
## 359       50.42225              0    53.09080
## 360       49.43798              0    51.04326
## 361       56.68721              0    62.35060
## 362       64.10667              0    66.33959
## 363       49.32120              0    53.16334
## 364       51.79682              0    51.87889
## 365       36.66823              0    37.69443
## 366       52.45317              0    52.85576
## 367       46.68047              0    47.50192
## 368       47.08077              0    51.24799
## 369       44.59611              0    44.28364
## 370       51.69279              0    57.09487
## 371       59.00575              0    52.67233
## 372       67.53292              0    64.60465
## 373       39.96000              0    37.82501
## 374       48.83330              0    48.25210
## 375       44.74953              0    44.78621
## 376       37.56909              0    36.24471
## 377       56.86393              0    53.86144
## 378       55.80525              0    58.08820
## 379       54.05850              0    57.17202
## 380       46.85523              0    45.38716
## 381       50.43086              0    47.30259
## 382       44.45665              0    41.27293
## 383       49.81294              0    44.00532
## 384       52.66387              0    53.99848
## 385       47.55092              0    47.11407
## 386       52.98935              0    52.51557
## 387       58.61860              0    61.08266
## 388       55.62724              0    52.15622
## 389       47.33579              0    47.52011
## 390       33.25200              0    30.92614
## 391       47.88265              0    49.30899
## 392       48.42658              0    50.84510
## 393       72.43066              0    73.99721
## 394       47.97097              0    47.98133
## 395       50.85491              0    44.79422
## 396       42.59108              0    45.86081
## 397       46.84430              0    43.76305
## 398       46.96170              0    45.13191
## 399       50.49418              0    45.80943
## 400       47.70053              0    51.91107
## 401       50.05184              0    45.15737
## 402       53.41627              0    58.59104
## 403       46.13964              0    46.59775
## 404       52.23492              0    52.83740
## 405       49.88338              0    46.53228
## 406       38.59610              0    32.76447
## 407       51.45090              0    49.88744
## 408       66.88381              0    67.14544
## 409       35.47945              0    37.20367
## 410       49.89665              0    54.10666
## 411       68.29199              0    65.32880
## 412       55.93294              0    54.84189
## 413       37.40107              0    36.86798
## 414       50.02935              0    48.03701
## 415       37.82888              0    38.82422
## 416       55.15186              0    58.56990
## 417       50.71344              0    46.73913
## 418       44.82474              0    38.32381
## 419       41.26454              0    39.73772
## 420       45.14715              0    38.54376
## 421       46.13881              0    46.41029
## 422       63.10962              0    64.42283
## 423       55.37678              0    47.85339
## 424       51.38740              0    54.48926
## 425       56.16612              0    53.25856
## 426       58.39292              0    61.45424
## 427       63.61257              0    64.58829
## 428       47.63051              0    45.24769
## 429       51.23969              0    51.62061
## 430       46.59337              0    46.12889
## 431       68.22944              0    61.93013
## 432       58.59336              0    52.95387
## 433       45.60880              0    50.02690
## 434       50.51920              0    52.69996
## 435       60.72734              0    60.84134
## 436       63.72257              0    66.17176
## 437       55.69623              0    46.21452
## 438       59.28723              0    61.79929
## 439       44.05188              0    47.38579
## 440       60.75954              0    60.53561
## 441       43.62258              0    41.56376
## 442       59.33150              0    60.22901
## 443       66.25646              0    64.86088
## 444       44.05986              0    50.23125
## 445       60.50470              0    57.71321
## 446       37.95396              0    36.76388
## 447       63.43437              0    63.58685
## 448       64.25358              0    64.86183
## 449       53.39947              0    50.94460
## 450       46.68328              0    49.32488
## 451       51.08249              0    48.57208
## 452       57.59895              0    55.27925
## 453       56.10090              0    58.58735
## 454       63.13714              0    63.17602
## 455       42.31537              0    41.13874
## 456       62.77405              0    55.00101
## 457       53.37736              0    57.52986
## 458       56.08585              0    52.01199
## 459       47.08249              0    43.34616
## 460       42.79524              0    46.20686
## 461       52.16510              0    50.00110
## 462       55.11803              0    59.98948
## 463       37.44791              0    37.49060
## 464       34.16731              0    33.54895
## 465       52.28851              0    51.45353
## 466       55.64321              0    59.57509
## 467       57.56691              0    56.29281
## 468       52.16679              0    60.03123
## 469       51.79499              0    48.36621
## 470       64.84692              0    63.01081
## 471       41.36725              0    40.32431
## 472       44.66277              0    46.41411
## 473       45.59439              0    44.41258
## 474       44.17401              0    47.17480
## 475       40.62247              0    35.57128
## 476       50.39387              0    51.60573
## 477       54.68838              0    53.28798
## 478       51.56048              0    47.75477
## 479       55.23810              0    57.31367
## 480       51.74690              0    52.92815
## 481       43.28916              0    45.08533
## 482       45.57214              0    43.73299
## 483       40.92575              0    46.27569
## 484       65.74185              0    62.29511
## 485       53.83183              0    55.36057
## 486       42.05615              0    38.31860
## 487       46.38883              0    43.89578
## 488       40.93230              0    40.73649
## 489       53.95995              0    57.88990
## 490       56.04657              0    52.83219
## 491       47.83130              0    47.17122
## 492       55.04529              0    57.22903
## 493       43.29176              0    45.03773
## 494       45.47393              0    42.09845
## 495       60.63569              0    61.00818
## 496       42.89350              0    45.66479
## 497       47.38266              0    43.82423
## 498       48.32128              0    48.82876
## 499       41.90877              0    37.72806
## 500       38.25218              0    34.31073
## 501       51.14606              0    49.57787
## 502       50.99715              0    48.03071
## 503       51.05039              0    50.96427
## 504       48.26966              0    44.54148
## 505       55.56823              0    56.46812
## 506       63.11838              0    61.63521
## 507       36.82290              0    39.48937
## 508       61.39575              0    64.05173
## 509       43.99891              0    45.68790
## 510       50.29610              0    51.97702
## 511       59.27278              0    63.56088
## 512       49.61085              0    50.18555
## 513       52.02409              0    56.39890
## 514       59.44991              0    60.66894
## 515       53.20121              0    53.90619
## 516       56.02982              0    57.90348
## 517       60.49722              0    58.78031
## 518       49.95968              0    44.60827
## 519       46.83140              0    47.99420
## 520       55.50914              0    55.61998
## 521       49.46096              0    45.12331
## 522       56.82721              0    54.82724
## 523       41.66396              0    42.70947
## 524       57.41053              0    60.28746
## 525       40.89672              0    37.84745
## 526       48.37579              0    47.83595
## 527       42.90626              0    42.37450
## 528       50.79822              0    51.53552
## 529       48.77649              0    45.86118
## 530       44.43115              0    48.67849
## 531       64.56218              0    64.42218
## 532       48.00476              0    51.53616
## 533       43.17025              0    38.33506
## 534       61.29681              0    61.12831
## 535       42.33587              0    43.38108
## 536       57.65370              0    54.55379
## 537       50.63514              0    45.93055
## 538       45.81419              0    46.37119
## 539       53.57359              0    56.28380
## 540       43.15352              0    43.40908
## 541       38.82991              0    37.56103
## 542       38.26978              0    39.83967
## 543       51.27667              0    52.66638
## 544       50.07159              0    50.52722
## 545       49.59497              0    46.50315
## 546       47.96151              0    52.13694
## 547       66.76998              0    61.31578
## 548       47.79202              0    48.78144
## 549       43.36175              0    39.35594
## 550       64.97167              0    66.21726
## 551       38.00448              0    34.34883
## 552       60.20808              0    59.40381
## 553       54.18001              0    51.47902
## 554       50.13577              0    50.45824
## 555       52.15959              0    52.93435
## 556       42.32258              0    43.02188
## 557       58.90617              0    59.49292
## 558       52.23159              0    55.73379
## 559       62.10649              0    61.18921
## 560       53.85630              0    55.85175
## 561       45.85886              0    47.74217
## 562       55.36310              0    54.54080
## 563       57.27075              0    63.68398
## 564       47.00529              0    48.33746
## 565       39.05901              0    41.09197
## 566       56.92067              0    56.22298
## 567       60.28653              0    59.29738
## 568       41.57270              0    38.07953
## 569       41.36660              0    41.15880
## 570       45.77802              0    45.85391
## 571       60.65708              0    60.72704
## 572       35.45432              0    33.48950
## 573       51.78601              0    50.36056
## 574       48.67563              0    49.94651
## 575       46.94048              0    41.80372
## 576       50.71840              0    47.92419
## 577       51.70831              0    59.39249
## 578       41.76382              0    35.15918
## 579       46.40438              0    52.03461
## 580       48.76856              0    50.21892
## 581       47.12957              0    52.05956
## 582       52.30255              0    46.18285
## 583       44.77258              0    45.14976
## 584       57.19758              0    58.16038
## 585       49.31073              0    55.57461
## 586       49.21833              0    47.04609
## 587       58.03120              0    56.45382
## 588       50.68004              0    48.35124
## 589       64.47437              0    70.23670
## 590       50.81152              0    53.07341
## 591       52.72378              0    48.76859
## 592       57.48950              0    66.30636
## 593       64.82165              0    62.30025
## 594       40.39320              0    39.44473
## 595       52.92532              0    57.71100
## 596       34.65779              0    33.89625
## 597       60.65136              0    61.10717
## 598       56.34819              0    59.47343
## 599       56.55741              0    56.48476
## 600       54.65790              0    53.63474
## 601       48.17806              0    49.42244
## 602       46.11115              0    43.84264
## 603       43.72488              0    40.14578
## 604       47.88071              0    52.79044
## 605       50.02416              0    49.61732
## 606       40.29141              0    37.02509
## 607       46.89597              0    49.10581
## 608       66.59541              0    64.98477
## 609       49.57961              0    48.96183
## 610       51.48751              0    52.61956
## 611       53.77015              0    55.43481
## 612       34.40585              0    26.96780
## 613       34.87015              0    37.94375
## 614       54.21981              0    54.94660
## 615       45.21023              0    44.52798
## 616       51.37321              0    52.65186
## 617       58.36369              0    59.80692
## 618       52.37259              0    46.94584
## 619       57.84202              0    60.28042
## 620       66.40174              0    57.71283
## 621       55.11444              0    52.23845
## 622       48.02034              0    45.26553
## 623       63.69599              0    58.58946
## 624       53.27375              0    59.22996
## 625       53.75352              0    49.59567
## 626       55.60592              0    57.11344
## 627       57.18991              0    56.24577
## 628       56.99375              0    64.89854
## 629       56.09014              0    55.09125
## 630       46.66800              0    40.36048
## 631       62.88201              0    61.74431
## 632       60.25243              0    63.27341
## 633       62.87171              0    65.25522
## 634       47.92609              0    47.83470
## 635       47.04940              0    44.91946
## 636       41.79077              0    42.78259
## 637       58.08604              0    58.40036
## 638       47.32369              0    49.12968
## 639       47.43606              0    45.87610
## 640       48.04430              0    50.38446
## 641       38.03855              0    34.38341
## 642       51.75471              0    53.83940
## 643       65.62744              0    64.85435
## 644       54.70396              0    51.30240
## 645       44.24940              0    44.99992
## 646       50.26058              0    48.56268
## 647       48.85177              0    49.80249
## 648       51.21352              0    49.05494
## 649       40.80550              0    37.33744
## 650       54.90063              0    55.24234
## 651       49.48876              0    51.57649
## 652       47.21733              0    46.71290
## 653       51.85952              0    54.78325
## 654       48.85573              0    51.53618
## 655       41.07034              0    41.83142
## 656       53.50493              0    55.59698
## 657       59.54467              0    57.41539
## 658       55.01416              0    55.53220
## 659       46.28275              0    43.97950
## 660       59.52462              0    61.94005
## 661       62.67200              0    58.34279
## 662       42.76619              0    44.88204
## 663       46.76515              0    47.52486
## 664       44.59197              0    45.29593
## 665       53.13976              0    48.08350
## 666       56.50419              0    56.73919
## 667       33.96834              0    33.07151
## 668       53.87865              0    55.53707
## 669       51.80917              0    53.66794
## 670       49.70134              0    48.66328
## 671       75.49429              0    75.50299
## 672       40.58071              0    38.70315
## 673       38.19754              0    34.97150
## 674       50.43831              0    51.54582
## 675       43.51733              0    46.79152
## 676       54.07147              0    58.76186
## 677       48.32390              0    50.75600
## 678       40.52142              0    40.68939
## 679       58.75355              0    55.54327
## 680       45.28245              0    51.26077
## 681       53.83548              0    57.55813
## 682       46.40081              0    48.60240
## 683       46.16660              0    45.64732
## 684       63.12883              0    63.14114
## 685       48.99494              0    52.99448
## 686       46.85225              0    48.54755
## 687       33.39581              0    38.51601
## 688       45.65260              0    50.90802
## 689       47.43775              0    46.99133
## 690       69.05211              0    68.65859
## 691       40.79897              0    45.01484
## 692       52.74967              0    56.51024
## 693       53.83695              0    59.73593
## 694       46.24236              0    41.77876
## 695       55.54558              0    56.08487
## 696       56.43770              0    55.34527
## 697       59.66692              0    60.27253
## 698       45.55771              0    45.21306
## 699       36.27891              0    29.47810
## 700       61.04503              0    62.72457
## 701       43.64964              0    49.02723
## 702       48.05758              0    47.73916
## 703       56.67536              0    59.81274
## 704       43.16625              0    43.10827
## 705       31.00411              0    32.42338
## 706       46.17442              0    43.16650
## 707       50.39296              0    49.60668
## 708       55.34426              0    56.41516
## 709       52.04929              0    50.51161
## 710       58.39592              0    55.61705
## 711       38.23222              0    41.32652
## 712       48.44378              0    50.56879
## 713       61.18625              0    61.93255
## 714       56.67400              0    63.59018
## 715       57.16610              0    53.40571
## 716       52.98443              0    56.39717
## 717       45.78689              0    48.85212
## 718       64.71978              0    70.20081
## 719       47.47978              0    45.98052
## 720       62.24075              0    65.22939
## 721       46.80533              0    49.47781
## 722       49.78134              0    48.54251
## 723       47.75323              0    51.89448
## 724       44.45245              0    44.51006
## 725       28.33304              0    32.46276
## 726       58.62412              0    60.92660
## 727       53.77160              0    52.66860
## 728       49.92896              0    49.17883
## 729       47.44331              0    48.98890
## 730       46.30881              0    45.34534
## 731       53.33571              0    46.10610
## 732       62.51973              0    63.22215
## 733       49.54492              0    57.15657
## 734       60.52546              0    60.73042
## 735       63.18790              0    64.97653
## 736       51.96245              0    55.02560
## 737       46.11509              0    43.52268
## 738       56.84243              0    55.59735
## 739       51.66860              0    44.82554
## 740       52.78474              0    55.85420
## 741       44.06760              0    44.47018
## 742       65.78141              0    63.26120
## 743       60.40291              0    58.53987
## 744       37.98958              0    42.55781
## 745       59.10844              0    61.61114
## 746       53.91486              0    55.83420
## 747       52.70522              0    52.70162
## 748       38.82086              0    43.17056
## 749       35.43545              0    36.30633
## 750       70.08950              0    67.90326
## 751       49.96218              0    52.50850
## 752       48.55096              0    49.19019
## 753       36.23918              0    36.63067
## 754       55.85198              0    59.00753
## 755       45.99672              0    45.05446
## 756       52.62948              0    49.23240
## 757       42.30587              0    45.78953
## 758       31.39508              0    30.14952
## 759       41.23809              0    42.99242
## 760       60.00519              0    60.38463
## 761       39.76072              0    41.39330
## 762       49.33766              0    52.58546
## 763       55.56757              0    57.88129
## 764       54.67333              0    54.31773
## 765       52.81707              0    57.33651
## 766       58.14924              0    60.98043
## 767       52.94475              0    48.19997
## 768       51.06460              0    48.82780
## 769       38.71247              0    38.34704
## 770       61.80098              0    62.07786
## 771       54.86446              0    56.11105
## 772       59.38769              0    57.46237
## 773       44.92990              0    42.87083
## 774       50.25662              0    54.94879
## 775       41.49083              0    42.33269
## 776       47.33601              0    45.43905
## 777       43.14518              0    47.78849
## 778       61.89705              0    59.93695
## 779       37.96273              0    39.75525
## 780       40.91808              0    45.29920
## 781       54.79120              0    59.21358
## 782       47.56865              0    41.99963
## 783       54.59782              0    53.03948
## 784       45.58835              0    48.29018
## 785       55.73130              0    56.95171
## 786       59.85771              0    68.13767
## 787       64.35316              0    63.98544
## 788       52.83611              0    55.60698
## 789       50.97774              0    54.42682
## 790       42.10384              0    48.08605
## 791       39.15176              0    40.45234
## 792       51.36612              0    48.45776
## 793       53.78679              0    51.20262
## 794       60.40030              0    61.91477
## 795       64.72192              0    63.26087
## 796       45.43750              0    47.90860
## 797       54.63926              0    54.14111
## 798       41.90362              0    39.10748
## 799       51.60690              0    55.74895
## 800       43.69768              0    41.35273
## 801       56.56657              0    52.85760
## 802       55.00061              0    54.15993
## 803       52.26064              0    55.79710
## 804       70.76890              0    65.21384
## 805       57.63029              0    53.03652
## 806       35.65779              0    31.95739
## 807       50.76072              0    47.65371
## 808       44.81475              0    43.98040
## 809       55.34014              0    58.15915
## 810       52.29277              0    52.86735
## 811       32.58723              0    29.96436
## 812       53.50218              0    54.79712
## 813       61.10066              0    60.51849
## 814       51.41171              0    50.59118
## 815       37.96887              0    33.09416
## 816       65.44257              0    68.55512
## 817       44.18013              0    43.71126
## 818       41.70861              0    44.02440
## 819       61.20767              0    59.48250
## 820       42.31108              0    41.20163
## 821       44.96780              0    49.46005
## 822       45.12626              0    40.47191
## 823       40.71048              0    39.17386
## 824       41.93239              0    38.33071
## 825       52.19465              0    52.12323
## 826       46.31849              0    44.09824
## 827       55.18294              0    54.48504
## 828       32.37008              0    29.18527
## 829       40.95071              0    42.38462
## 830       56.87704              0    57.57855
## 831       56.01854              0    56.62936
## 832       55.82269              0    53.73159
## 833       45.33361              0    45.85555
## 834       55.76840              0    54.47226
## 835       41.43735              0    42.50218
## 836       44.07929              0    48.01519
## 837       43.90399              0    47.38701
## 838       52.38469              0    54.63368
## 839       64.60261              0    59.03636
## 840       40.57182              0    40.62265
## 841       40.80069              0    43.77630
## 842       35.56967              0    34.59830
## 843       64.93245              0    64.25495
## 844       49.90992              0    49.63982
## 845       45.31036              0    44.76297
## 846       39.77068              0    39.69005
## 847       51.14476              0    50.99048
## 848       48.94695              0    47.68515
## 849       41.06610              0    45.29255
## 850       38.97431              0    37.86237
## 851       45.61225              0    45.23310
## 852       49.38832              0    52.70248
## 853       50.36572              0    54.48337
## 854       48.79792              0    50.58985
## 855       51.61373              0    52.71649
## 856       53.77545              0    54.46259
## 857       39.83984              0    36.38270
## 858       67.50486              0    66.32397
## 859       50.47539              0    51.05334
## 860       34.03687              0    33.79645
## 861       55.31831              0    53.63060
## 862       48.43456              0    44.87401
## 863       36.02856              0    30.04491
## 864       49.47803              0    48.02209
## 865       41.40037              0    38.23670
## 866       54.74085              0    48.70690
## 867       57.85278              0    59.52936
## 868       44.03732              0    43.89063
## 869       41.30879              0    47.23183
## 870       54.53110              0    57.68586
## 871       58.46181              0    58.14529
## 872       55.95064              0    56.98714
## 873       48.27919              0    49.98848
## 874       42.84122              0    44.83507
## 875       44.32626              0    45.78675
## 876       48.26009              0    45.42176
## 877       51.91730              0    53.61330
## 878       54.89423              0    54.45274
## 879       56.82096              0    55.63070
## 880       53.73163              0    53.22739
## 881       56.37982              0    53.81778
## 882       43.32966              0    45.33031
## 883       46.78184              0    45.98659
## 884       54.02442              0    56.61908
## 885       45.31133              0    46.66515
## 886       53.53152              0    50.50308
## 887       48.41782              0    42.34791
## 888       45.44038              0    44.07324
## 889       43.82846              0    40.47985
## 890       52.84252              0    53.37788
## 891       44.01766              0    38.37852
## 892       52.88877              0    53.00168
## 893       47.01403              0    43.91679
## 894       44.22611              0    48.19225
## 895       49.45921              0    45.88106
## 896       61.13268              0    67.79815
## 897       52.11333              0    49.29489
## 898       48.89686              0    52.94650
## 899       61.51898              0    60.16230
## 900       48.87504              0    45.62733
## 901       45.50939              0    47.99000
## 902       42.22804              0    42.62871
## 903       52.64778              0    52.85094
## 904       39.26965              0    41.30120
## 905       56.99808              0    57.28878
## 906       61.04475              0    64.42263
## 907       35.56650              0    36.44113
## 908       57.54542              0    58.53500
## 909       47.76335              0    51.24016
## 910       52.13803              0    55.59823
## 911       47.82870              0    54.74413
## 912       44.80369              0    40.85727
## 913       51.07270              0    51.22818
## 914       52.54214              0    54.13765
## 915       54.02123              0    59.32670
## 916       43.85128              0    43.20230
## 917       52.46487              0    51.08479
## 918       52.60895              0    53.01861
## 919       41.53947              0    42.23343
## 920       47.60247              0    51.60873
## 921       43.70711              0    45.28413
## 922       41.25241              0    40.45863
## 923       60.67195              0    58.88487
## 924       58.00402              0    55.67107
## 925       51.82382              0    47.91957
## 926       52.87261              0    58.37413
## 927       46.40433              0    44.94995
## 928       47.55655              0    51.16991
## 929       55.10946              0    57.84882
## 930       46.32273              0    53.75579
## 931       39.03783              0    35.02495
## 932       59.22115              0    56.11692
## 933       50.03936              0    47.74204
## 934       41.43690              0    38.24866
## 935       32.51564              0    32.81972
## 936       40.82450              0    41.55732
## 937       47.89547              0    53.10714
## 938       53.30081              0    49.29076
## 939       45.75698              0    48.60274
## 940       55.55330              0    55.85748
## 941       44.12653              0    42.61655
## 942       39.32604              0    41.45575
## 943       52.39759              0    51.64305
## 944       38.67085              0    37.91099
## 945       48.35650              0    49.67146
## 946       39.86766              0    40.36922
## 947       51.18046              0    53.97398
## 948       73.72453              0    76.85426
## 949       45.64841              0    46.53359
## 950       47.27924              0    50.23731
## 951       55.25810              0    50.74062
## 952       49.92331              0    55.54584
## 953       39.90305              0    39.14142
## 954       56.82937              0    56.37826
## 955       45.94474              0    51.29470
## 956       34.66390              0    39.16148
## 957       31.22911              0    35.87381
## 958       50.32273              0    50.37421
## 959       45.46235              0    48.36320
## 960       50.69558              0    53.31752
## 961       50.22254              0    47.68565
## 962       52.67093              0    54.09385
## 963       48.13452              0    49.09936
## 964       49.64543              0    53.60754
## 965       55.02668              0    47.62213
## 966       48.65367              0    49.94457
## 967       48.03326              0    47.26785
## 968       41.14404              0    42.82980
## 969       42.89363              0    46.65373
## 970       60.42502              0    60.30425
## 971       53.37088              0    56.21343
## 972       65.76433              0    64.11125
## 973       53.29328              0    50.93665
## 974       42.88440              0    41.81695
## 975       46.44792              0    50.08647
## 976       56.08153              0    58.09976
## 977       68.95509              0    65.96120
## 978       52.28559              0    51.84696
## 979       58.42156              0    55.82570
## 980       47.48758              0    53.13310
## 981       44.68188              0    43.29627
## 982       46.88789              0    44.73045
## 983       64.03012              0    62.03417
## 984       48.65541              0    52.61939
## 985       51.40495              0    48.79830
## 986       41.86280              0    34.74296
## 987       50.57697              0    51.12162
## 988       47.51156              0    50.75495
## 989       60.16867              0    63.43671
## 990       60.04833              0    63.05670
## 991       51.54346              0    52.01895
## 992       45.50248              0    45.69262
## 993       54.20327              0    48.36390
## 994       28.04200              0    28.11850
## 995       53.59475              0    56.51452
## 996       39.64990              0    40.08806
## 997       62.92202              0    59.76006
## 998       60.96561              0    61.40538
## 999       59.23580              0    61.68493
## 1000      39.57106              0    41.50169

Model fitting with propensity scores

This is a two-step process, estimate propensity scores, then estimate final model.

sim_arguments <- list(
    formula = achievement ~ 1 + motivation + trt + age + ses,
    fixed = list(
        motivation = list(var_type = 'continuous',
                           mean = 0, sd = 20)
    ),
    sample_size = 1000,
    error = list(variance = 10),
    reg_weights = c(50, 0.4, 1.2, 0.1, 0.25),
    propensity = list(
        formula = trt ~ 1 + age + ses,
        fixed = list(age = list(var_type = 'ordinal', 
                         levels = -7:7),
                     ses = list(var_type = 'continuous', 
                         mean = 0, sd = 5)),
        sample_size = 1000,
        error = list(variance = 5),
        reg_weights = c(2, 0.3, -0.5),
        outcome_type = 'binary'
    ),
    model_fit = list(formula = achievement ~ 1 + motivation + trt + age + ses,
                   model_function = 'lm'),
    propensity_model = list(
        formula = trt ~ 1 + age + ses,
        propensity_type = 'covariate'
    )
)
simulate_fixed(data = NULL, sim_args = sim_arguments) |>
  simulate_error(sim_args = sim_arguments) |>
  generate_response(sim_arguments)|> 
  fit_propensity(sim_arguments) |>
  model_fit(sim_arguments) |>
  extract_coefficients()
## # A tibble: 6 × 5
##   term        estimate std.error statistic   p.value
##   <chr>          <dbl>     <dbl>     <dbl>     <dbl>
## 1 (Intercept)   51.1     0.848       60.2  0        
## 2 motivation     0.404   0.00523     77.2  0        
## 3 trt            1.12    0.270        4.15 0.0000359
## 4 age            0.158   0.0427       3.71 0.000222 
## 5 ses            0.159   0.0602       2.63 0.00860  
## 6 propensity    -1.83    1.27        -1.44 0.151
sim_arguments <- list(
    formula = achievement ~ 1 + motivation + trt + age + ses,
    fixed = list(
        motivation = list(var_type = 'continuous',
                           mean = 0, sd = 20)
    ),
    sample_size = 1000,
    error = list(variance = 10),
    reg_weights = c(50, 0.4, 1.2, 0.1, 0.25),
    propensity = list(
        formula = trt ~ 1 + age + ses,
        fixed = list(age = list(var_type = 'ordinal', 
                         levels = -7:7),
                     ses = list(var_type = 'continuous', 
                         mean = 0, sd = 5)),
        sample_size = 1000,
        error = list(variance = 5),
        reg_weights = c(2, 0.3, -0.5),
        outcome_type = 'binary'
    ),
    propensity_model = list(
        formula = trt ~ 1 + age + ses,
        propensity_type = 'ipw'
    )
)
simulate_fixed(data = NULL, sim_args = sim_arguments) |>
  simulate_error(sim_args = sim_arguments) |>
  generate_response(sim_arguments)|> 
  fit_propensity(sim_arguments) |>
  model_fit(sim_arguments) |>
  extract_coefficients()
## # A tibble: 5 × 5
##   term        estimate std.error statistic  p.value
##   <chr>          <dbl>     <dbl>     <dbl>    <dbl>
## 1 (Intercept)   49.8     0.204      245.   0       
## 2 motivation     0.396   0.00498     79.5  0       
## 3 trt            1.49    0.260        5.74 1.24e- 8
## 4 age            0.109   0.0239       4.55 6.11e- 6
## 5 ses            0.245   0.0226      10.8  7.57e-26
sim_arguments <- list(
    formula = achievement ~ 1 + motivation + trt + age + ses,
    fixed = list(
        motivation = list(var_type = 'continuous',
                           mean = 0, sd = 20)
    ),
    sample_size = 1000,
    error = list(variance = 10),
    reg_weights = c(50, 0.4, 1.2, 0.1, 0.25),
    propensity = list(
        formula = trt ~ 1 + age + ses,
        fixed = list(age = list(var_type = 'ordinal', 
                         levels = -7:7),
                     ses = list(var_type = 'continuous', 
                         mean = 0, sd = 5)),
        sample_size = 1000,
        error = list(variance = 5),
        reg_weights = c(2, 0.3, -0.5),
        outcome_type = 'binary'
    ),
    propensity_model = list(
        formula = trt ~ 1 + age + ses,
        propensity_type = 'sbw'
    )
)
simulate_fixed(data = NULL, sim_args = sim_arguments) |>
  simulate_error(sim_args = sim_arguments) |>
  generate_response(sim_arguments)|> 
  fit_propensity(sim_arguments) |>
  model_fit(sim_arguments) |>
  extract_coefficients()
## # A tibble: 5 × 5
##   term        estimate std.error statistic  p.value
##   <chr>          <dbl>     <dbl>     <dbl>    <dbl>
## 1 (Intercept)   49.8     0.207      241.   0       
## 2 motivation     0.396   0.00509     77.9  0       
## 3 trt            1.25    0.262        4.77 2.08e- 6
## 4 age            0.106   0.0239       4.44 1.00e- 5
## 5 ses            0.244   0.0237      10.3  1.24e-23
sim_arguments <- list(
    formula = achievement ~ 1 + motivation + trt + age + ses,
    fixed = list(
        motivation = list(var_type = 'continuous',
                           mean = 0, sd = 20)
    ),
    sample_size = 1000,
    error = list(variance = 10),
    reg_weights = c(50, 0.4, 1.2, 0.1, 0.25),
    propensity = list(
        formula = trt ~ 1 + age + ses,
        fixed = list(age = list(var_type = 'ordinal', 
                         levels = -7:7),
                     ses = list(var_type = 'continuous', 
                         mean = 0, sd = 5)),
        sample_size = 1000,
        error = list(variance = 5),
        reg_weights = c(2, 0.3, -0.5),
        outcome_type = 'binary'
    ),
    model_fit = list(formula = achievement ~ 1 + motivation + trt + age + ses,
                   model_function = 'lm'),
    propensity_model = list(
        formula = trt ~ 1 + age + ses,
        propensity_type = 'covariate'
    ),
    replications = 20,
    extract_coefficients = TRUE
)
replicate_simulation(sim_arguments) |>
  compute_statistics(sim_arguments, alternative_power = FALSE, 
                     type_1_error = FALSE)
## # A tibble: 6 × 7
##   term   avg_estimate power param_estimate_sd avg_standard_error precision_ratio
##   <chr>         <dbl> <dbl>             <dbl>              <dbl>           <dbl>
## 1 (Inte…     50.0      1              0.883              0.837             1.06 
## 2 age         0.0961   0.5            0.0461             0.0396            1.16 
## 3 motiv…      0.400    1              0.00515            0.00502           1.03 
## 4 prope…      0.00924  0              1.40               1.23              1.14 
## 5 ses         0.256    0.95           0.0592             0.0568            1.04 
## 6 trt         1.28     1              0.217              0.263             0.824
## # ℹ 1 more variable: replications <dbl>
sim_arguments <- list(
    formula = achievement ~ 1 + motivation + trt + age + ses,
    fixed = list(
        motivation = list(var_type = 'continuous',
                           mean = 0, sd = 20)
    ),
    sample_size = 1000,
    error = list(variance = 10),
    reg_weights = c(50, 0.4, 1.2, 0.1, 0.25),
    propensity = list(
        formula = trt ~ 1 + age + ses,
        fixed = list(age = list(var_type = 'ordinal', 
                         levels = -7:7),
                     ses = list(var_type = 'continuous', 
                         mean = 0, sd = 5)),
        sample_size = 1000,
        error = list(variance = 5),
        reg_weights = c(2, 0.3, -0.5),
        outcome_type = 'binary'
    ),
    model_fit = list(formula = achievement ~ 1 + motivation + trt + age + ses,
                   model_function = 'lm'),
    propensity_model = list(
        formula = trt ~ 1 + age + ses,
        propensity_type = 'ipw'
    ),
    replications = 500,
    extract_coefficients = TRUE
)
replicate_simulation(sim_arguments) |>
  compute_statistics(sim_arguments, alternative_power = FALSE, 
                     type_1_error = FALSE)
## # A tibble: 5 × 7
##   term   avg_estimate power param_estimate_sd avg_standard_error precision_ratio
##   <chr>         <dbl> <dbl>             <dbl>              <dbl>           <dbl>
## 1 (Inte…       50.0   1               0.307              0.216              1.42
## 2 age           0.101 0.88            0.0420             0.0252             1.66
## 3 motiv…        0.400 1               0.00975            0.00502            1.94
## 4 ses           0.245 0.996           0.0544             0.0205             2.65
## 5 trt           1.17  0.938           0.374              0.263              1.42
## # ℹ 1 more variable: replications <dbl>
sim_arguments <- list(
    formula = achievement ~ 1 + motivation + trt + age + ses,
    fixed = list(
        motivation = list(var_type = 'continuous',
                           mean = 0, sd = 20)
    ),
    sample_size = 1000,
    error = list(variance = 10),
    reg_weights = c(50, 0.4, 1.2, 0.1, 0.25),
    propensity = list(
        formula = trt ~ 1 + age + ses,
        fixed = list(age = list(var_type = 'ordinal', 
                         levels = -7:7),
                     ses = list(var_type = 'continuous', 
                         mean = 0, sd = 5)),
        sample_size = 1000,
        error = list(variance = 5),
        reg_weights = c(2, 0.3, -0.5),
        outcome_type = 'binary'
    ),
    model_fit = list(formula = achievement ~ 1 + motivation + trt + age + ses,
                   model_function = 'lm'),
    propensity_model = list(
        formula = trt ~ 1 + age + ses,
        propensity_type = 'sbw'
    ),
    replications = 500,
    extract_coefficients = TRUE
)
replicate_simulation(sim_arguments) |>
  compute_statistics(sim_arguments, alternative_power = FALSE, 
                     type_1_error = FALSE)
## # A tibble: 5 × 7
##   term   avg_estimate power param_estimate_sd avg_standard_error precision_ratio
##   <chr>         <dbl> <dbl>             <dbl>              <dbl>           <dbl>
## 1 (Inte…       50.0   1               0.301              0.216              1.39
## 2 age           0.100 0.864           0.0401             0.0252             1.59
## 3 motiv…        0.400 1               0.00927            0.00502            1.85
## 4 ses           0.250 1               0.0524             0.0205             2.55
## 5 trt           1.19  0.936           0.381              0.262              1.45
## # ℹ 1 more variable: replications <dbl>