# Use the multitrait dataset
data(multitrait)
multitrait <- subset(multitrait, chr=1:2,
ind=!apply(multitrait$pheno, 1, function(a) any(is.na(a))))
multitrait$pheno <- log(multitrait$pheno)
multitrait <- calc.genoprob(multitrait)
result <- mqmpermutation(multitrait,pheno.col=7, n.perm=2, batchsize=2)
if (FALSE) #Set 50 cofactors
cof <- mqmautocofactors(multitrait,50)
cof <- mqmautocofactors(multitrait,5)
multitrait <- fill.geno(multitrait)
result <- mqmpermutation(multitrait,scanfunction=mqmscan,cofactors=cof,
pheno.col=7, n.perm=2,batchsize=2,verbose=FALSE)
#Create a permutation object
f2perm <- mqmprocesspermutation(result)
#Get Significant LOD thresholds
summary(f2perm)
Run the code above in your browser using DataLab