library(SNPassoc)
library(haplo.stats)
data(SNPs)
datSNP<-setupSNP(SNPs, 6:40, sep = "")
mod <- haplo.interaction(casco~int(sex)+blood.pre, data = datSNP,
SNPs.sel = c("snp10001","snp10004","snp10005"))
Run the code above in your browser using DataLab