data(vignette)
gene <- "YNL037C"
variables <- which(rownames(brem.data) != gene)
control <- iBMAcontrolLM(OR = 50, nbest = 20, thresProbne0 = 5)
iBMAmodel.YNL037C <- iterateBMAlm( x = t(brem.data[variables,]),
y = unlist(brem.data[gene,]), control = control)
Run the code above in your browser using DataLab