data(multitrait)
# impute missing genotype data
multitrait <- fill.geno(multitrait)
if (FALSE) # Calculate the thresholds
result <- mqmscanfdr(multitrait, threshold=10.0, n.perm=1000)
# \dontshow{
multitrait <- subset(multitrait, chr=1, ind=!apply(multitrait$pheno, 1, function(a) any(is.na(a))))
multitrait$pheno <- multitrait$pheno[,1:3]
result <- mqmscanfdr(multitrait, threshold=10.0, n.perm=2)
# }
Run the code above in your browser using DataLab