Return the Bayes Factor for each effect
get_log10bf(m)the mash result (from joint or 1by1 analysis); must have been computed using usepointmass=TRUE
if m was fitted using usepointmass=TRUE then returns a vector of the log10(bf) values for each effect. That is, the jth element lbf[j] is log10(Pr(Bj | g=ghat-nonnull)/Pr(Bj | g = 0)) where ghat-nonnull is the non-null part of ghat. Otherwise returns NULL.
simdata = simple_sims(50,5,1)
data = mash_set_data(simdata$Bhat, simdata$Shat)
m = mash(data, cov_canonical(data))
#> - Computing 200 x 141 likelihood matrix.
#> - Likelihood calculations took 0.00 seconds.
#> - Fitting model with 141 mixture components.
#> - Model fitting took 0.15 seconds.
#> - Computing posterior matrices.
#> - Computation allocated took 0.00 seconds.
get_log10bf(m)
#> [,1]
#> [1,] -0.273530427
#> [2,] 0.037220895
#> [3,] -0.190549578
#> [4,] -0.310071537
#> [5,] -0.049194803
#> [6,] -0.179990976
#> [7,] -0.338519846
#> [8,] -0.292484477
#> [9,] -0.376218214
#> [10,] -0.417569562
#> [11,] -0.241382394
#> [12,] -0.442630801
#> [13,] 0.289602358
#> [14,] -0.457239027
#> [15,] -0.528192409
#> [16,] 0.574263994
#> [17,] -0.051771650
#> [18,] -0.467483267
#> [19,] -0.121408194
#> [20,] -0.425396713
#> [21,] -0.194776450
#> [22,] 0.223089615
#> [23,] -0.355094271
#> [24,] -0.332012148
#> [25,] -0.442640890
#> [26,] -0.258305954
#> [27,] 0.128111787
#> [28,] -0.517651655
#> [29,] 0.310052087
#> [30,] 0.316442447
#> [31,] 0.078427468
#> [32,] -0.543752660
#> [33,] -0.456489891
#> [34,] -0.554861984
#> [35,] 0.326322121
#> [36,] -0.355048058
#> [37,] 0.079290313
#> [38,] -0.339115235
#> [39,] -0.432013099
#> [40,] 0.203249275
#> [41,] 0.124357072
#> [42,] 0.511506553
#> [43,] 0.115247358
#> [44,] -0.477146309
#> [45,] -0.276398031
#> [46,] -0.009884408
#> [47,] 0.045542480
#> [48,] 0.229794817
#> [49,] -0.308104694
#> [50,] -0.252475086
#> [51,] 0.085567729
#> [52,] -0.414749814
#> [53,] 0.295979652
#> [54,] 0.331257258
#> [55,] 0.286958612
#> [56,] 0.019322361
#> [57,] 0.415839160
#> [58,] 0.033940772
#> [59,] 0.253712505
#> [60,] 0.071300252
#> [61,] 0.358433574
#> [62,] -0.179039286
#> [63,] -0.526604053
#> [64,] -0.192525441
#> [65,] -0.300542521
#> [66,] 0.201920289
#> [67,] -0.007620657
#> [68,] -0.043540688
#> [69,] -0.381958268
#> [70,] 1.433308358
#> [71,] 1.686133078
#> [72,] 0.264015470
#> [73,] -0.070430153
#> [74,] 1.577734564
#> [75,] -0.383924184
#> [76,] -0.328470742
#> [77,] 1.496768509
#> [78,] -0.008000826
#> [79,] -0.390466020
#> [80,] -0.044133508
#> [81,] 0.467426156
#> [82,] -0.144045996
#> [83,] -0.472585414
#> [84,] 0.553273114
#> [85,] 0.570425838
#> [86,] 0.204376794
#> [87,] 0.766469256
#> [88,] 1.054604697
#> [89,] 1.161535907
#> [90,] -0.132042726
#> [91,] 0.107840181
#> [92,] -0.443488562
#> [93,] 0.101283455
#> [94,] 0.911214187
#> [95,] 1.830934389
#> [96,] 0.372441478
#> [97,] 0.802164804
#> [98,] -0.520666371
#> [99,] 0.046137284
#> [100,] 0.008193687
#> [101,] 0.140456759
#> [102,] -0.274531416
#> [103,] -0.207017393
#> [104,] 0.939820416
#> [105,] -0.541232330
#> [106,] -0.463274994
#> [107,] 0.182837807
#> [108,] 0.416530091
#> [109,] 0.218271490
#> [110,] 0.059898948
#> [111,] -0.307860261
#> [112,] -0.010881738
#> [113,] -0.511564291
#> [114,] -0.091273781
#> [115,] -0.123735380
#> [116,] -0.459544538
#> [117,] -0.467600192
#> [118,] -0.185403118
#> [119,] -0.443262684
#> [120,] 0.271472218
#> [121,] -0.399003839
#> [122,] -0.034014979
#> [123,] -0.332332175
#> [124,] 0.395604542
#> [125,] 0.049534418
#> [126,] -0.270746523
#> [127,] -0.282688456
#> [128,] -0.440209288
#> [129,] 0.403840512
#> [130,] -0.216429544
#> [131,] 1.115180494
#> [132,] -0.441712310
#> [133,] 0.061556507
#> [134,] -0.067509601
#> [135,] -0.162351525
#> [136,] -0.094082348
#> [137,] -0.406242850
#> [138,] -0.510985094
#> [139,] -0.338824829
#> [140,] -0.179022318
#> [141,] 0.226304368
#> [142,] -0.303006937
#> [143,] 0.482386187
#> [144,] -0.366087753
#> [145,] -0.255818570
#> [146,] 0.339416791
#> [147,] -0.096294523
#> [148,] -0.106480231
#> [149,] -0.351172225
#> [150,] 0.327356978
#> [151,] -0.385969867
#> [152,] 0.092273506
#> [153,] -0.175808683
#> [154,] 9.339558641
#> [155,] -0.534889905
#> [156,] 0.016359558
#> [157,] 1.948875484
#> [158,] -0.108609428
#> [159,] 2.619508888
#> [160,] 3.185976976
#> [161,] 0.685112929
#> [162,] -0.250307623
#> [163,] -0.248476226
#> [164,] 0.300865593
#> [165,] -0.191988100
#> [166,] 0.279151614
#> [167,] -0.242068166
#> [168,] -0.032141085
#> [169,] -0.437453354
#> [170,] -0.503292048
#> [171,] 0.223210148
#> [172,] 0.066814564
#> [173,] -0.519906990
#> [174,] 0.961974989
#> [175,] -0.297690017
#> [176,] -0.540703451
#> [177,] -0.265396325
#> [178,] 0.217142893
#> [179,] -0.443222129
#> [180,] -0.089148771
#> [181,] 1.138716349
#> [182,] 1.148240230
#> [183,] -0.438798794
#> [184,] -0.403012328
#> [185,] -0.227453353
#> [186,] 1.079178968
#> [187,] 0.540687959
#> [188,] -0.446377189
#> [189,] -0.078673635
#> [190,] -0.476267382
#> [191,] -0.344689185
#> [192,] 3.952185497
#> [193,] -0.215141758
#> [194,] 0.780828289
#> [195,] 0.156065551
#> [196,] 0.017711061
#> [197,] -0.265198340
#> [198,] -0.429192465
#> [199,] 0.066538109
#> [200,] -0.398115257