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))
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
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>