powered by
This function plots ranked observed p values against the corresponding expected p values in -log scale.
qqpval(p, pch=16, col=4, ...)
No return value, just the plot
a vector of p values
symbol to use for points
color for points
other plot arguments
GenomicControl, WGassociation
GenomicControl
WGassociation
data(SNPs) datSNP<-setupSNP(SNPs,6:40,sep="") res<-WGassociation(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