data(multitrait) # Load dataset
multitrait <- fill.geno(multitrait) # impute missing genotype data
multitrait <- subset(multitrait, chr=1)
multitrait$pheno <- multitrait$pheno[,1:3]
result <- mqmscanall(multitrait, logtransform=TRUE) # Analyse all 24 traits
Run the code above in your browser using DataLab