powered by
qqpval(p, pch=16, col=4, ...)
GenomicControl
WGassociation
data(SNPs) datSNP<-setupSNP(SNPs,6:40,sep="") res<-scanWGassociation(casco,datSNP,model=c("do","re","log-add")) # observed vs expected p values for recessive model qqpval(recessive(res))
Run the code above in your browser using DataLab