单体型软件LBLGXE用法

#file.rda文件中的变量名就是LBL函数的第一个参数LBL.data

#运行结果out.LBL有7个指标（BF，OR，CI.OR，freq，CI.freq，CI.lambda，CI.D），将它们都写入res.txt中

load("file.rda")

library(hapassoc)

library(dummies)

library(LBLGXE)

out.LBL <- LBL(LBL.data, numSNPs = 5, burn.in = 100, num.it = 1000)

write.table("$BF\n", file = "res.txt", quote = F, row.names = F, col.names = F, append = T)

write.table(out.LBL$BF, file = "res.txt", sep = "\t", quote = F, row.names = T, append = T)

write.table("\n\n$OR\n", file = "res.txt", quote = F, row.names = T, col.names = F, append = T)

write.table(out.LBL$OR, file = "res.txt", sep = "\t", quote = F, row.names = T, append = T)

write.table("\n\n$CI.OR\n", file = "res.txt", quote = F, row.names = T, col.names = F,  append = T)

write.table(out.LBL$CI.OR, file = "res.txt", sep = "\t", quote = F, row.names = T,  append = T)

write.table("\n\n$freq\n", file = "res.txt",  quote = F, row.names = T, col.names = F, append = T)

write.table(out.LBL$freq, file = "res.txt", sep = "\t", quote = F, row.names = T, append = T)

write.table("\n\n$CI.freq\n", file = "res.txt",  quote = F, row.names = T, col.names = F, append = T)

write.table(out.LBL$CI.freq, file = "res.txt", sep = "\t", quote = F, row.names = T, append = T)

write.table("\n\n$CI.lambda\n", file = "res.txt",  quote = F, row.names = T, col.names = F, append = T)

write.table(out.LBL$CI.lambda, file = "res.txt", sep = "\t", quote = F, row.names = T, append = T)

write.table("\n\n$CI.D\n", file = "res.txt",  quote = F, row.names = T, col.names = F, append = T)

write.table(out.LBL$CI.D, file = "res.txt", sep = "\t", quote = F, row.names = T, append = T)