Compute LD scores from an LD reference, optionally stratified by
annotations. Each r-squared is ldsc's small-sample unbiased estimate,
\(r^2 - (1 - r^2)/(n - 2)\) for a reference of \(n\) samples, as
buildLdScore stores, so the two agree.
Usage
computeLdScores(ldRef, annotations = NULL, ...)
# S4 method for class 'LdEigen'
computeLdScores(ldRef, annotations = NULL)
# S4 method for class 'LdScore'
computeLdScores(ldRef, annotations = NULL)Examples
data(ldEigenExample)
computeLdScores(ldEigenExample)
#> base_l2
#> [1,] 3.146379
#> [2,] 3.436336
#> [3,] 2.655487
#> [4,] 3.313009
#> [5,] 2.464155
#> [6,] 3.369796
#> [7,] 4.598756
#> [8,] 3.076745
#> [9,] 4.206343
#> [10,] 2.655669
#> [11,] 2.680233
#> [12,] 3.880546
#> [13,] 2.880485
#> [14,] 5.734343
#> [15,] 6.193682
#> [16,] 6.883334
#> [17,] 6.883334
#> [18,] 6.883334
#> [19,] 3.851683
#> [20,] 3.963069
#> [21,] 3.005638
#> [22,] 2.950587
#> [23,] 3.108698
#> [24,] 2.850877
#> [25,] 3.345546
#> [26,] 3.279622
#> [27,] 3.437041
#> [28,] 3.278083
#> [29,] 4.413120
#> [30,] 4.600659
#> [31,] 4.371835
#> [32,] 6.782334
#> [33,] 3.463930
#> [34,] 6.411946
#> [35,] 7.080253
#> [36,] 4.610600
#> [37,] 7.164481
#> [38,] 7.510986
#> [39,] 7.237913
#> [40,] 4.626483
#> [41,] 7.660934
#> [42,] 8.791799
#> [43,] 9.236815
#> [44,] 7.741542
#> [45,] 11.216092
#> [46,] 6.742422
#> [47,] 9.275131
#> [48,] 11.216092
#> [49,] 9.282903
#> [50,] 7.833519
#> [51,] 5.804487
#> [52,] 11.216092
#> [53,] 11.216092
#> [54,] 11.216092
#> [55,] 10.781312
#> [56,] 10.262660
#> [57,] 3.620936
#> [58,] 3.079627
#> [59,] 5.016769
#> [60,] 4.605945
#> [61,] 3.688850
#> [62,] 3.755102
#> [63,] 3.301777
#> [64,] 2.899539
#> [65,] 4.100630
#> [66,] 3.898339
#> [67,] 3.129712
#> [68,] 3.534907
#> [69,] 3.163733
#> [70,] 2.829869
#> [71,] 5.049307
#> [72,] 4.627447
#> [73,] 5.066096
#> [74,] 5.066096
#> [75,] 2.654630
#> [76,] 3.183732
#> [77,] 3.476311
#> [78,] 3.726583
#> [79,] 2.615500
#> [80,] 3.133914
#> [81,] 4.212416
#> [82,] 4.429248
#> [83,] 3.905235
#> [84,] 3.305307
#> [85,] 4.546984
#> [86,] 3.298820
#> [87,] 3.247811
#> [88,] 2.769390
#> [89,] 3.501962
#> [90,] 4.019066
#> [91,] 4.392086
#> [92,] 3.735797
#> [93,] 3.786163
#> [94,] 3.033887
#> [95,] 3.874965
#> [96,] 4.727651
#> [97,] 4.466406
#> [98,] 3.092609
#> [99,] 4.581036
#> [100,] 4.217729
#> [101,] 4.674471
#> [102,] 5.230218
#> [103,] 5.385014
#> [104,] 3.592222
#> [105,] 3.377900
#> [106,] 3.542912
#> [107,] 2.937860
#> [108,] 1.915661
#> [109,] 3.197017
#> [110,] 2.344438
#> [111,] 3.275075
#> [112,] 2.988499
#> [113,] 3.463272
#> [114,] 3.225642
#> [115,] 2.681422
#> [116,] 3.128038
#> [117,] 2.799897
#> [118,] 2.508748
#> [119,] 3.364890
#> [120,] 4.487167
#> [121,] 3.299198
#> [122,] 3.295066
#> [123,] 4.728542
#> [124,] 4.728542
#> [125,] 3.677150
#> [126,] 2.971147
#> [127,] 3.892594
#> [128,] 4.658054
#> [129,] 4.376736
#> [130,] 2.954840
#> [131,] 4.132142
#> [132,] 4.353946
#> [133,] 4.414126
#> [134,] 4.129864
#> [135,] 2.772231
#> [136,] 4.160859
#> [137,] 2.285214
#> [138,] 4.018283
#> [139,] 3.856979
#> [140,] 2.711219
#> [141,] 3.880577
#> [142,] 4.282457
#> [143,] 3.431143
#> [144,] 7.844940
#> [145,] 7.844940
#> [146,] 7.485647
#> [147,] 8.053997
#> [148,] 8.053997
#> [149,] 8.053997
#> [150,] 7.504634
#> [151,] 4.220820
#> [152,] 3.930637
#> [153,] 3.469913
#> [154,] 3.617238
#> [155,] 3.455787
#> [156,] 6.198782
#> [157,] 4.853993
#> [158,] 3.183240
#> [159,] 7.081749
#> [160,] 3.029728
#> [161,] 2.852346
#> [162,] 7.081749
#> [163,] 7.081749
#> [164,] 6.982345
#> [165,] 6.545356
#> [166,] 4.434405
#> [167,] 2.919684
#> [168,] 3.727815
#> [169,] 3.884053
#> [170,] 3.475057
#> [171,] 3.636684
#> [172,] 2.413500
#> [173,] 2.183060
#> [174,] 2.923933
#> [175,] 2.095836
#> [176,] 2.749381
#> [177,] 3.422833
#> [178,] 2.402625
#> [179,] 2.497229
#> [180,] 3.002456
#> [181,] 2.735667
#> [182,] 4.028946
#> [183,] 3.562310
#> [184,] 4.812618
#> [185,] 4.866661
#> [186,] 4.048286
#> [187,] 2.913980
#> [188,] 3.904910
#> [189,] 3.572608
#> [190,] 3.664642
#> [191,] 3.016233
#> [192,] 4.094486
#> [193,] 4.578867
#> [194,] 5.024018
#> [195,] 5.440073
#> [196,] 4.014667
#> [197,] 6.505837
#> [198,] 6.515631
#> [199,] 5.938360
#> [200,] 5.856521
#> [201,] 3.531141
#> [202,] 4.149924
#> [203,] 3.539181
#> [204,] 3.688841
#> [205,] 3.229552
#> [206,] 3.473660
#> [207,] 3.983048
#> [208,] 3.157058
#> [209,] 2.719767
#> [210,] 4.836877
#> [211,] 5.017058
#> [212,] 5.428133
#> [213,] 5.490252
#> [214,] 3.448937
#> [215,] 5.431948
#> [216,] 4.886577
#> [217,] 3.296513
#> [218,] 2.885028
#> [219,] 2.615375
#> [220,] 2.557367
#> [221,] 5.267715
#> [222,] 2.829113
#> [223,] 3.317106
#> [224,] 5.050281
#> [225,] 6.182850
#> [226,] 3.463070
#> [227,] 5.547460
#> [228,] 5.210518
#> [229,] 5.048456
#> [230,] 3.218308
#> [231,] 3.458408
#> [232,] 3.439478
#> [233,] 3.704014
#> [234,] 3.806976
#> [235,] 7.098080
#> [236,] 7.098080
#> [237,] 7.098080
#> [238,] 7.098080
#> [239,] 6.309042
#> [240,] 2.903681
#> [241,] 6.299647
#> [242,] 4.580608
#> [243,] 3.752004
#> [244,] 5.350014
#> [245,] 2.821509
#> [246,] 4.702909
#> [247,] 3.559061
#> [248,] 3.537576
#> [249,] 3.366484
#> [250,] 3.449143
#> [251,] 3.545213
#> [252,] 3.867969
#> [253,] 4.094892
#> [254,] 3.110709
#> [255,] 3.551507
#> [256,] 3.068172
#> [257,] 3.256956
#> [258,] 3.293777
#> [259,] 2.813277
#> [260,] 2.886096
#> [261,] 3.513101
#> [262,] 3.776486
#> [263,] 4.793157
#> [264,] 4.427311
#> [265,] 3.829941
#> [266,] 8.901644
#> [267,] 8.901644
#> [268,] 5.934050
#> [269,] 4.427311
#> [270,] 8.425452
#> [271,] 6.670812
#> [272,] 6.911087
#> [273,] 8.425452
#> [274,] 8.582890
#> [275,] 6.260003
#> [276,] 5.378521
#> [277,] 4.061832
#> [278,] 5.771459
#> [279,] 4.217907
#> [280,] 8.508789
#> [281,] 3.006863
#> [282,] 2.946679
#> [283,] 4.136097
#> [284,] 4.023659
#> [285,] 2.925562
#> [286,] 3.068012
#> [287,] 3.671544
#> [288,] 5.164242
#> [289,] 4.962477
#> [290,] 3.048550
#> [291,] 3.288928
#> [292,] 5.164242
#> [293,] 7.858718
#> [294,] 7.976924
#> [295,] 8.011907
#> [296,] 8.156716
#> [297,] 5.411844
#> [298,] 5.411844
#> [299,] 7.123499
#> [300,] 8.051311
#> [301,] 7.710619
#> [302,] 5.097448
#> [303,] 2.780421
#> [304,] 3.702262
#> [305,] 4.231930
#> [306,] 4.851822
#> [307,] 2.849129
#> [308,] 3.090644
#> [309,] 5.002064
#> [310,] 5.098176
#> [311,] 4.951482
#> [312,] 3.760157
#> [313,] 2.660994
#> [314,] 3.031169
#> [315,] 2.673737
#> [316,] 4.194007
#> [317,] 3.766483
#> [318,] 4.270583
#> [319,] 3.247372
#> [320,] 3.339743
#> [321,] 2.932399
#> [322,] 2.544670
#> [323,] 2.437847
#> [324,] 2.581154
#> [325,] 2.339075
#> [326,] 3.120555
#> [327,] 3.083279
#> [328,] 3.616313
#> [329,] 3.384208
#> [330,] 3.040242
#> [331,] 2.859939
#> [332,] 3.342720
#> [333,] 2.386435
#> [334,] 3.363229
#> [335,] 3.208545
#> [336,] 2.420827
#> [337,] 3.326139
#> [338,] 3.522352
#> [339,] 3.688738
#> [340,] 2.253171
#> [341,] 2.508853
#> [342,] 3.774302
#> [343,] 4.343694
#> [344,] 4.682407
#> [345,] 4.153825
#> [346,] 2.442456
#> [347,] 3.359604
#> [348,] 2.631665
#> [349,] 2.933599
#> [350,] 2.902631
#> [351,] 2.760889
#> [352,] 2.917805
#> [353,] 2.622331
#> [354,] 2.871880
#> [355,] 2.836567
#> [356,] 2.268110
#> [357,] 4.339308
#> [358,] 4.505025
#> [359,] 4.799028
#> [360,] 4.799028
#> [361,] 3.157447
#> [362,] 4.605921
#> [363,] 4.591211
#> [364,] 3.722917
#> [365,] 2.114006
#> [366,] 3.060509
#> [367,] 3.204471
#> [368,] 3.778719
#> [369,] 3.163962
#> [370,] 4.199194
#> [371,] 2.871382
#> [372,] 3.388623
#> [373,] 4.248736
#> [374,] 4.261197
#> [375,] 4.480415
#> [376,] 3.293698
#> [377,] 3.836563
#> [378,] 4.711619
#> [379,] 4.686704
#> [380,] 2.607942
#> [381,] 2.917881
#> [382,] 2.916046
#> [383,] 3.126898
#> [384,] 2.420781
#> [385,] 3.297247
#> [386,] 4.185760
#> [387,] 4.185760
#> [388,] 4.185760
#> [389,] 3.361458
#> [390,] 3.379113
#> [391,] 2.529218
#> [392,] 3.697809
#> [393,] 5.779907
#> [394,] 5.779907
#> [395,] 5.779907
#> [396,] 4.309931
#> [397,] 4.304329
#> [398,] 4.929399
#> [399,] 5.060448
#> [400,] 4.426639
#> [401,] 4.625790
#> [402,] 3.779551
#> [403,] 3.649816
#> [404,] 3.633904
#> [405,] 4.106840
#> [406,] 3.034661
#> [407,] 3.426320
#> [408,] 3.697012
#> [409,] 4.942735
#> [410,] 3.326765
#> [411,] 3.005672
#> [412,] 3.450345
#> [413,] 3.915883
#> [414,] 4.475136
#> [415,] 3.737004
#> [416,] 6.106859
#> [417,] 12.154135
#> [418,] 12.154135
#> [419,] 12.154135
#> [420,] 12.154135
#> [421,] 12.154135
#> [422,] 12.154135
#> [423,] 11.169539
#> [424,] 11.539813
#> [425,] 11.539813
#> [426,] 7.443475
#> [427,] 6.027091
#> [428,] 7.341558
#> [429,] 4.785584
#> [430,] 3.090597
#> [431,] 3.792225
#> [432,] 4.256355
#> [433,] 4.268732
#> [434,] 4.351659
#> [435,] 4.492883
#> [436,] 3.599733
#> [437,] 3.988852
#> [438,] 4.984513
#> [439,] 6.158353
#> [440,] 5.885947
#> [441,] 5.595742
#> [442,] 6.025069
#> [443,] 5.562915
#> [444,] 5.395335
#> [445,] 5.395335
#> [446,] 3.579027
#> [447,] 4.868866
#> [448,] 4.847536
#> [449,] 4.847536
#> [450,] 5.090319
#> [451,] 3.716246
#> [452,] 5.035294
#> [453,] 4.930697
#> [454,] 3.072137
#> [455,] 4.316905
#> [456,] 4.423738
#> [457,] 4.884349
#> [458,] 10.123994
#> [459,] 10.806244
#> [460,] 11.315645
#> [461,] 5.937332
#> [462,] 10.712577
#> [463,] 12.790348
#> [464,] 12.790348
#> [465,] 12.790348
#> [466,] 12.790348
#> [467,] 12.790348
#> [468,] 12.790348
#> [469,] 8.439350
#> [470,] 11.409894
#> [471,] 11.409894
#> [472,] 5.369604
#> [473,] 5.338473
#> [474,] 3.986649
#> [475,] 6.398803
#> [476,] 6.153218
#> [477,] 5.344499
#> [478,] 4.817526
#> [479,] 6.889701
#> [480,] 6.802403
#> [481,] 6.749200
#> [482,] 6.424704
#> [483,] 6.578240
#> [484,] 6.424704
#> [485,] 4.095237
#> [486,] 4.291848
#> [487,] 3.281262
#> [488,] 4.300438
#> [489,] 4.016652
#> [490,] 4.402343
#> [491,] 3.746471
#> [492,] 3.725541
#> [493,] 4.637141
#> [494,] 4.146708
#> [495,] 4.653332
#> [496,] 3.346756
#> [497,] 3.136614
#> [498,] 5.159034
#> [499,] 5.346301
#> [500,] 5.346301
#> [501,] 3.361680
#> [502,] 6.511572
#> [503,] 7.023566
#> [504,] 7.023566
#> [505,] 7.023566
#> [506,] 4.247923
#> [507,] 3.062658
#> [508,] 3.436264
#> [509,] 4.801915
#> [510,] 4.964942
#> [511,] 6.748593
#> [512,] 6.580333
#> [513,] 6.010064
#> [514,] 3.065854
#> [515,] 4.562758
#> [516,] 4.311486
#> [517,] 4.173805
#> [518,] 4.253296
#> [519,] 4.797537
#> [520,] 4.797537
#> [521,] 5.537161
#> [522,] 4.830859
#> [523,] 5.366995
#> [524,] 3.948373
#> [525,] 3.680345
#> [526,] 7.174054
#> [527,] 3.840026
#> [528,] 4.053022
#> [529,] 4.202143
#> [530,] 4.311165
#> [531,] 4.615420
#> [532,] 3.599558
#> [533,] 2.805561
#> [534,] 3.947280
#> [535,] 4.054863
#> [536,] 3.726297
#> [537,] 3.432417
#> [538,] 3.719568
#> [539,] 3.777844
#> [540,] 2.773518
#> [541,] 2.415827
#> [542,] 4.170972
#> [543,] 3.522416
#> [544,] 2.974366
#> [545,] 2.987278
#> [546,] 8.674346
#> [547,] 8.966086
#> [548,] 8.966086
#> [549,] 8.966086
#> [550,] 7.025748
#> [551,] 3.016035
#> [552,] 2.797109
#> [553,] 3.085006
#> [554,] 2.136600
#> [555,] 2.231598
#> [556,] 2.345050
#> [557,] 4.924883
#> [558,] 3.922284
#> [559,] 3.191789
#> [560,] 3.103604
#> [561,] 4.885618
#> [562,] 3.186557
#> [563,] 10.153149
#> [564,] 10.153149
#> [565,] 11.046549
#> [566,] 5.076627
#> [567,] 9.184280
#> [568,] 11.555762
#> [569,] 10.639238
#> [570,] 4.397404
#> [571,] 11.584879
#> [572,] 11.584879
#> [573,] 11.273320
#> [574,] 11.584879
#> [575,] 6.080871
#> [576,] 10.947028
#> [577,] 10.506090
#> [578,] 9.495012
#> [579,] 9.644286
#> [580,] 2.831326
#> [581,] 11.050840
#> [582,] 11.254464
#> [583,] 3.299417
#> [584,] 3.962902
#> [585,] 3.962902
#> [586,] 4.138192
#> [587,] 4.083229
#> [588,] 2.874345
#> [589,] 3.311728
#> [590,] 2.511084
#> [591,] 3.573911
#> [592,] 4.476533
#> [593,] 3.364651
#> [594,] 3.522718
#> [595,] 3.641155
#> [596,] 3.126358
#> [597,] 3.348497
#> [598,] 2.448244
#> [599,] 2.438454
#> [600,] 2.323531
#> [601,] 6.171522
#> [602,] 6.171522
#> [603,] 6.171522
#> [604,] 5.055863
#> [605,] 3.135812
#> [606,] 4.505496
#> [607,] 3.480845
#> [608,] 3.211378
#> [609,] 2.706215
#> [610,] 2.422820
#> [611,] 2.973501
#> [612,] 2.838424
#> [613,] 3.737003
#> [614,] 3.296964
#> [615,] 2.882379
#> [616,] 2.512521
#> [617,] 2.570925
#> [618,] 2.580453
#> [619,] 3.564370
#> [620,] 3.300073
#> [621,] 3.784443
#> [622,] 3.668596
#> [623,] 2.982575
#> [624,] 3.881112
#> [625,] 2.727567
#> [626,] 4.179181
#> [627,] 3.221891
#> [628,] 3.440976
#> [629,] 2.429210
#> [630,] 3.999092
#> [631,] 3.326662
#> [632,] 2.755662
#> [633,] 3.671522
#> [634,] 3.926946
#> [635,] 3.310369
#> [636,] 2.793092
#> [637,] 3.139811
#> [638,] 4.119797
#> [639,] 4.140225
#> [640,] 4.200202
#> [641,] 3.990050
#> [642,] 3.448593
#> [643,] 3.600084
#> [644,] 2.665172
#> [645,] 3.567584
#> [646,] 3.270154
#> [647,] 3.915288
#> [648,] 3.590624
#> [649,] 3.247009
#> [650,] 3.305148
#> [651,] 3.355047
#> [652,] 2.378917
#> [653,] 2.956925
#> [654,] 3.183433
#> [655,] 2.663633
#> [656,] 4.252297
#> [657,] 3.918418
#> [658,] 3.316325
#> [659,] 4.308608
#> [660,] 4.165275
#> [661,] 3.875890
#> [662,] 4.674693
#> [663,] 5.758010
#> [664,] 5.528767
#> [665,] 5.179821
#> [666,] 6.202695
#> [667,] 6.202695
#> [668,] 3.585585
#> [669,] 3.968846
#> [670,] 4.298556
#> [671,] 4.134332
#> [672,] 3.272108
#> [673,] 5.207848
#> [674,] 4.106869
#> [675,] 3.976174
#> [676,] 4.515800
#> [677,] 5.109426
#> [678,] 4.159660
#> [679,] 8.275187
#> [680,] 9.194492
#> [681,] 4.045443
#> [682,] 8.530601
#> [683,] 8.963424
#> [684,] 9.024258
#> [685,] 4.765443
#> [686,] 7.329319
#> [687,] 5.892274
#> [688,] 2.979732
#> [689,] 2.982804
#> [690,] 3.761452
#> [691,] 4.194433
#> [692,] 5.269417
#> [693,] 3.212166
#> [694,] 3.244711
#> [695,] 4.330749
#> [696,] 5.821143
#> [697,] 6.777025
#> [698,] 5.975106
#> [699,] 5.749357
#> [700,] 6.630824
#> [701,] 6.648320
#> [702,] 6.648320
#> [703,] 7.179926
#> [704,] 5.914339
#> [705,] 6.564283
#> [706,] 4.964507
#> [707,] 4.078030
#> [708,] 4.539795
#> [709,] 5.120143
#> [710,] 7.659784
#> [711,] 3.388077
#> [712,] 6.429866
#> [713,] 5.912788
#> [714,] 3.166225
#> [715,] 6.015126
#> [716,] 5.818549
#> [717,] 5.629337
#> [718,] 5.444193
#> [719,] 6.255684
#> [720,] 5.934118
#> [721,] 6.125828
#> [722,] 3.713487
#> [723,] 3.130183
#> [724,] 5.821709
#> [725,] 6.331113
#> [726,] 4.740779
#> [727,] 4.987815
#> [728,] 6.006096
#> [729,] 4.140808
#> [730,] 5.177875
#> [731,] 5.518034
#> [732,] 5.086509
#> [733,] 6.523003
#> [734,] 7.799783
#> [735,] 7.799783
#> [736,] 7.799783
#> [737,] 7.799783
#> [738,] 5.625716
#> [739,] 5.520784
#> [740,] 6.058793
#> [741,] 3.610127
#> [742,] 6.680999
#> [743,] 6.680999
#> [744,] 4.842809
#> [745,] 5.995762
#> [746,] 4.367258
#> [747,] 5.715873
#> [748,] 3.349999
#> [749,] 3.449525
#> [750,] 3.214904
#> [751,] 2.521851
#> [752,] 3.382707
#> [753,] 4.564109
#> [754,] 4.269624
#> [755,] 3.258424
#> [756,] 4.564807
#> [757,] 3.838508
#> [758,] 3.434762
#> [759,] 2.190297
#> [760,] 2.716325
#> [761,] 2.038824
#> [762,] 1.878400
#> [763,] 2.861317
#> [764,] 3.203631
#> [765,] 2.206672
#> [766,] 3.200344
#> [767,] 2.074896
#> [768,] 2.205244
#> [769,] 4.814830
#> [770,] 3.959708
#> [771,] 5.889257
#> [772,] 5.573177
#> [773,] 4.717337
#> [774,] 5.348297
#> [775,] 4.669799
#> [776,] 3.388205
#> [777,] 6.077974
#> [778,] 2.286197
#> [779,] 4.204665
#> [780,] 4.028194
#> [781,] 3.149826
#> [782,] 3.558502
#> [783,] 3.275251
#> [784,] 5.268392
#> [785,] 3.052991
#> [786,] 4.843614
#> [787,] 3.352271
#> [788,] 3.415918
#> [789,] 2.459404
#> [790,] 2.879515
#> [791,] 4.291912
#> [792,] 2.993510
#> [793,] 3.567353
#> [794,] 2.669828
#> [795,] 2.300038
#> [796,] 2.419694
#> [797,] 3.119735
#> [798,] 4.169143
#> [799,] 3.634300
#> [800,] 4.288234
#> [801,] 2.877077
#> [802,] 4.372082
#> [803,] 4.949861
#> [804,] 3.841591
#> [805,] 4.819611
#> [806,] 4.333476
#> [807,] 4.242428
#> [808,] 4.687397
#> [809,] 3.651078
#> [810,] 5.123994
#> [811,] 4.664774
#> [812,] 5.230880
#> [813,] 4.653373
#> [814,] 4.717543
#> [815,] 3.523253
#> [816,] 4.152296
#> [817,] 4.587996
#> [818,] 2.990551
#> [819,] 3.298656
#> [820,] 2.295550
#> [821,] 2.625594
#> [822,] 2.824481
#> [823,] 2.985958
#> [824,] 2.608164
#> [825,] 7.608837
#> [826,] 7.460202
#> [827,] 7.578228
#> [828,] 7.460202
#> [829,] 3.776462
#> [830,] 7.034410
#> [831,] 6.665359
#> [832,] 4.135256
#> [833,] 5.015569
#> [834,] 6.643942
#> [835,] 3.830943
#> [836,] 6.286417
#> [837,] 4.223000
#> [838,] 6.662307
#> [839,] 4.183614
#> [840,] 4.944958
#> [841,] 6.269880
#> [842,] 3.862961
#> [843,] 2.875517
#> [844,] 3.321967
#> [845,] 6.114126
#> [846,] 6.114126
#> [847,] 2.936418
#> [848,] 3.112109
#> [849,] 7.342780
#> [850,] 6.163789
#> [851,] 3.353521
#> [852,] 5.446205
#> [853,] 7.297527
#> [854,] 2.852632
#> [855,] 6.122113
#> [856,] 7.217071
#> [857,] 6.850336
#> [858,] 6.850336
#> [859,] 5.989051
#> [860,] 8.500815
#> [861,] 9.080457
#> [862,] 4.773844
#> [863,] 7.574427
#> [864,] 6.995340
#> [865,] 4.032783
#> [866,] 7.774416
#> [867,] 7.774416
#> [868,] 7.346835
#> [869,] 6.679798
#> [870,] 7.962968
#> [871,] 5.238803
#> [872,] 9.279522
#> [873,] 8.725614
#> [874,] 4.835824
#> [875,] 5.250534
#> [876,] 2.920574
#> [877,] 6.566419
#> [878,] 3.711096
#> [879,] 6.839758
#> [880,] 4.124912
#> [881,] 5.876282
#> [882,] 2.823430
#> [883,] 3.553678
#> [884,] 6.055310
#> [885,] 5.813363
#> [886,] 6.122797
#> [887,] 3.024513
#> [888,] 2.457249
#> [889,] 2.458227
#> [890,] 2.776793
#> [891,] 2.757788
#> [892,] 2.342067
#> [893,] 2.733485
#> [894,] 3.278824
#> [895,] 5.353230
#> [896,] 3.448564
#> [897,] 2.679106
#> [898,] 3.485224
#> [899,] 3.306152
#> [900,] 5.175321
#> [901,] 2.839005
#> [902,] 4.962841
#> [903,] 2.362465
#> [904,] 3.148843
#> [905,] 3.761101
#> [906,] 4.934611
#> [907,] 5.041799
#> [908,] 5.134877
#> [909,] 5.403140
#> [910,] 4.379372
#> [911,] 4.585337
#> [912,] 2.857915
#> [913,] 4.756494
#> [914,] 2.728263
#> [915,] 4.317905
#> [916,] 4.464485
#> [917,] 4.332322
#> [918,] 3.200952
#> [919,] 5.080835
#> [920,] 4.578976
#> [921,] 5.839262
#> [922,] 5.604592
#> [923,] 3.451145
#> [924,] 3.975284
#> [925,] 4.044674
#> [926,] 5.883415
#> [927,] 4.709790
#> [928,] 3.169702
#> [929,] 3.614770
#> [930,] 4.618931
#> [931,] 4.680913
#> [932,] 2.931881
#> [933,] 2.675721
#> [934,] 3.603009
#> [935,] 5.336688
#> [936,] 3.767246
#> [937,] 5.092574
#> [938,] 3.187937
#> [939,] 2.869965
#> [940,] 2.963775
#> [941,] 2.972422
#> [942,] 3.006919
#> [943,] 2.608282
#> [944,] 4.163998
#> [945,] 3.095145
#> [946,] 3.155935
#> [947,] 4.259970
#> [948,] 4.259970
#> [949,] 3.641226
#> [950,] 3.641226
#> [951,] 3.810117
#> [952,] 2.809729
#> [953,] 2.854382
#> [954,] 2.977084
#> [955,] 3.945644
#> [956,] 4.432747
#> [957,] 3.949639
#> [958,] 2.713726
#> [959,] 3.110990
#> [960,] 3.954678
#> [961,] 3.762200
#> [962,] 5.621933
#> [963,] 5.616838
#> [964,] 4.028419
#> [965,] 6.057162
#> [966,] 6.310744
#> [967,] 3.393901
#> [968,] 4.038527
#> [969,] 4.416792
#> [970,] 4.416792
#> [971,] 3.347052
#> [972,] 4.238916
#> [973,] 3.862604
#> [974,] 3.442203
#> [975,] 2.665440
#> [976,] 3.031749
#> [977,] 3.448242
#> [978,] 3.272397
#> [979,] 3.440061
#> [980,] 4.027429
#> [981,] 4.475924
#> [982,] 3.315595
#> [983,] 3.315595
#> [984,] 3.449869
#> [985,] 3.193698
#> [986,] 4.023489
#> [987,] 4.283013
#> [988,] 4.023489
#> [989,] 5.543120
#> [990,] 2.799796
#> [991,] 3.479875
#> [992,] 5.133541
#> [993,] 5.285251
#> [994,] 3.469639
#> [995,] 3.462908
#> [996,] 3.926281
#> [997,] 2.784467
#> [998,] 3.150928
#> [999,] 3.257456
#> [1000,] 3.484738
#> [1001,] 2.967805
#> [1002,] 3.115348
#> [1003,] 2.916994
#> [1004,] 2.703871
#> [1005,] 3.196217
#> [1006,] 3.059642
#> [1007,] 2.695923
#> [1008,] 3.451800
#> [1009,] 3.427350
#> [1010,] 3.152853
#> [1011,] 3.894005
#> [1012,] 3.894005
#> [1013,] 2.188211
#> [1014,] 2.608755
#> [1015,] 2.630573
#> [1016,] 2.725560
#> [1017,] 2.893157
#> [1018,] 2.673463
#> [1019,] 2.893520
#> [1020,] 3.517630
#> [1021,] 2.877488
#> [1022,] 3.493626
#> [1023,] 3.629165
#> [1024,] 2.098475
#> [1025,] 2.434023
#> [1026,] 2.363170
#> [1027,] 2.867590
#> [1028,] 2.717779
#> [1029,] 3.414679
#> [1030,] 4.537510
#> [1031,] 5.811045
#> [1032,] 5.237814
#> [1033,] 4.348837
#> [1034,] 4.403798
#> [1035,] 3.082293
#> [1036,] 3.082293
#> [1037,] 3.358242
#> [1038,] 3.667885
#> [1039,] 5.284512
#> [1040,] 4.532769
#> [1041,] 4.981772
#> [1042,] 5.901085
#> [1043,] 4.579338
#> [1044,] 4.094730
#> [1045,] 2.978670
#> [1046,] 4.976612
#> [1047,] 2.706788
#> [1048,] 2.781508
#> [1049,] 2.916159
#> [1050,] 1.864908
#> [1051,] 3.235899
#> [1052,] 2.706788
#> [1053,] 2.978252
#> [1054,] 2.607164
#> [1055,] 2.348745
#> [1056,] 3.103728
#> [1057,] 2.115094
#> [1058,] 2.173397
#> [1059,] 2.345612
#> [1060,] 2.312233
#> [1061,] 4.103216
#> [1062,] 4.103216
#> [1063,] 4.103216
#> [1064,] 1.939810
#> [1065,] 5.481907
#> [1066,] 5.481907
#> [1067,] 3.122370
#> [1068,] 4.635612
#> [1069,] 4.122984
#> [1070,] 3.681146
#> [1071,] 5.344442
#> [1072,] 3.010078
#> [1073,] 4.961671
#> [1074,] 4.616440
#> [1075,] 3.848914
#> [1076,] 3.740028
#> [1077,] 3.520467
#> [1078,] 3.241367
#> [1079,] 2.972166
#> [1080,] 3.769410
#> [1081,] 2.303825
#> [1082,] 3.209601
#> [1083,] 2.523008
#> [1084,] 4.117213
#> [1085,] 2.598642
#> [1086,] 3.327462
#> [1087,] 2.729228
#> [1088,] 2.985115
#> [1089,] 3.059547
#> [1090,] 2.229978
#> [1091,] 3.594578
#> [1092,] 3.836482
#> [1093,] 2.933885
#> [1094,] 4.022593
#> [1095,] 3.834692
#> [1096,] 3.084951
#> [1097,] 2.628109
#> [1098,] 2.975511
#> [1099,] 2.840565
#> [1100,] 3.763966
#> [1101,] 4.792743
#> [1102,] 4.792743
#> [1103,] 4.447017
#> [1104,] 2.667234
#> [1105,] 4.399855
#> [1106,] 4.090874
#> [1107,] 4.399855
#> [1108,] 3.635846
#> [1109,] 3.435083
#> [1110,] 2.462986
#> [1111,] 2.722248
#> [1112,] 3.395836
#> [1113,] 3.129625
#> [1114,] 3.085057
#> [1115,] 3.700442
#> [1116,] 3.614990
#> [1117,] 3.922137
#> [1118,] 3.381374
#> [1119,] 3.795068
#> [1120,] 2.652668
#> [1121,] 2.711771
#> [1122,] 3.119488
#> [1123,] 2.449373
#> [1124,] 2.805867
#> [1125,] 8.863514
#> [1126,] 8.863514
#> [1127,] 8.863514
#> [1128,] 8.863514
#> [1129,] 8.228015
#> [1130,] 8.863514
#> [1131,] 8.863514
#> [1132,] 8.863514
#> [1133,] 3.551140
#> [1134,] 3.551140
#> [1135,] 2.384519
#> [1136,] 2.552125
#> [1137,] 3.419903
#> [1138,] 2.661156
#> [1139,] 3.509488
#> [1140,] 2.602649
#> [1141,] 2.569596
#> [1142,] 2.139862
#> [1143,] 4.270820
#> [1144,] 3.572920
#> [1145,] 3.988333
#> [1146,] 2.845540
#> [1147,] 3.688991
#> [1148,] 3.682929
#> [1149,] 3.972287
#> [1150,] 3.146689
#> [1151,] 3.162695
#> [1152,] 2.811009
#> [1153,] 3.597380
#> [1154,] 2.872496
#> [1155,] 2.426765
#> [1156,] 3.317585
#> [1157,] 2.340523
#> [1158,] 2.139079
#> [1159,] 3.318773
#> [1160,] 1.928566
#> [1161,] 1.988021
#> [1162,] 2.107818
#> [1163,] 2.898309
#> [1164,] 2.898309
#> [1165,] 2.124907
#> [1166,] 2.170453
#> [1167,] 2.529131
#> [1168,] 2.670861
#> [1169,] 2.848810
#> [1170,] 2.541264
#> [1171,] 2.246855
#> [1172,] 2.985644
#> [1173,] 2.222912
#> [1174,] 3.455101
#> [1175,] 4.333310
#> [1176,] 4.162875
#> [1177,] 1.860659
#> [1178,] 2.839757
#> [1179,] 3.041363
#> [1180,] 4.172152
#> [1181,] 4.243664
#> [1182,] 3.260194
#> [1183,] 3.136197
#> [1184,] 2.965319
#> [1185,] 2.852950
#> [1186,] 2.847636
#> [1187,] 1.928782
#> [1188,] 2.487622
#> [1189,] 2.368516
#> [1190,] 2.109926
#> [1191,] 2.376308
#> [1192,] 2.248121
#> [1193,] 1.776544
#> [1194,] 2.220401
#> [1195,] 3.078377
#> [1196,] 3.937702
#> [1197,] 3.552705
#> [1198,] 2.280772
#> [1199,] 3.129692
#> [1200,] 3.357921
#> [1201,] 3.022261
#> [1202,] 3.937702
#> [1203,] 3.937702
#> [1204,] 3.361281
#> [1205,] 3.818788
#> [1206,] 3.648661
#> [1207,] 3.397792
#> [1208,] 3.248495
#> [1209,] 3.101112
#> [1210,] 4.210122
#> [1211,] 4.220411
#> [1212,] 3.240295
#> [1213,] 2.902305
#> [1214,] 2.391227
#> [1215,] 2.572713
#> [1216,] 2.137377
#> [1217,] 1.828757
#> [1218,] 1.958946
#> [1219,] 2.099362
#> [1220,] 2.159552
#> [1221,] 1.842804
#> [1222,] 1.890814
#> [1223,] 2.051806
#> [1224,] 3.900025
#> [1225,] 3.952731
#> [1226,] 2.130187
#> [1227,] 2.063060
#> [1228,] 2.703662
#> [1229,] 3.164842
#> [1230,] 3.408398
#> [1231,] 3.022988
#> [1232,] 2.431207
#> [1233,] 3.441341
#> [1234,] 3.401706
#> [1235,] 2.265136
#> [1236,] 4.248590
#> [1237,] 4.880526
#> [1238,] 4.740186
#> [1239,] 4.740186
#> [1240,] 4.504482
#> [1241,] 4.346871
#> [1242,] 4.096873
#> [1243,] 3.301283
#> [1244,] 2.152750
#> [1245,] 2.051708
#> [1246,] 2.655088
#> [1247,] 3.712056
#> [1248,] 5.289205
#> [1249,] 2.393780
#> [1250,] 6.535380
#> [1251,] 3.207480
#> [1252,] 6.965069
#> [1253,] 6.510191
#> [1254,] 2.189077
#> [1255,] 4.100766
#> [1256,] 2.511964
#> [1257,] 2.127884
#> [1258,] 5.377878
#> [1259,] 5.696519
#> [1260,] 3.897053
#> [1261,] 4.679343
#> [1262,] 6.922848
#> [1263,] 4.003160
#> [1264,] 8.235795
#> [1265,] 5.850806
#> [1266,] 8.521514
#> [1267,] 7.432843
#> [1268,] 7.717090
#> [1269,] 8.015525
#> [1270,] 6.407930
#> [1271,] 7.868879
#> [1272,] 2.519787
#> [1273,] 2.310635
#> [1274,] 3.458351
#> [1275,] 9.180662
#> [1276,] 3.774051
#> [1277,] 3.012195
#> [1278,] 2.880081
#> [1279,] 2.957019
#> [1280,] 2.121947
#> [1281,] 9.180662
#> [1282,] 9.180662
#> [1283,] 8.223591
#> [1284,] 4.959304
#> [1285,] 8.395862
#> [1286,] 6.428264
#> [1287,] 4.575469
#> [1288,] 7.073966
#> [1289,] 6.089900
#> [1290,] 4.795327
#> [1291,] 6.856658
#> [1292,] 6.856658
#> [1293,] 6.856658
#> [1294,] 5.699272
#> [1295,] 4.538094
#> [1296,] 4.879082
#> [1297,] 3.406341
#> [1298,] 3.406341
#> [1299,] 6.150439
#> [1300,] 5.351724