Calculate log bayes factors for each hypothesis (Single causal variant assumption)
Examples
library(cophescan)
data(cophe_multi_trait_data)
query_trait_1 <- cophe_multi_trait_data$summ_stat[['Trait_1']]
querysnpid <- cophe_multi_trait_data$querysnpid
res.single.lbf <- cophe.single.lbf(query_trait_1, querysnpid = querysnpid, querytrait='Trait_1')
res.single.lbf
#> lBF.Ha lBF.Hc nsnps querysnp querytrait typeBF
#> 1 15.32189 11.96576 1000 chr19-11182353 Trait_1 ABF
