inst/examples/example1.R

In adimitromanolakis/rareBF: Bayesian Models for Rare Variant Association Analysis

library(rareBF)

set.seed(101)

nsites = 500

ncase = 200
ncontrol = 200

maf = 0.005
rexp_param = 2


model1 = function(i) {
    
    p1 = rbinom( ncontrol  , nsites , maf   )
    p2 = rbinom( ncase     , nsites , maf )  + H1 * round(rexp(ncase,rexp_param))
    
    disease_status = c (  rep(0,ncontrol),rep(1,ncase)  )
    data = data.frame(c(p1,p2), disease_status)
    
    bf = BFmixture(data, nsites)[3]
    bf
}

nrep = 20

H1 = 0
bf1 = ( replicate(nrep,model1() ) )

H1 = 1
bf2 = ( replicate(nrep,model1() ) )

col = c( rep("blue",nrep) , rep("red", nrep))


plot(log(c(bf1,bf2)),col = col ,pch=20,cex=1.5)