qqpval {SNPassoc}

Functions for inspecting population substructure

Description

This function plots ranked observed p values against the corresponding expected p values in -log scale.

Usage

qqpval(p, pch=16, col=4, ...)

Arguments

p

a vector of p values

pch

symbol to use for points

col

color for points

...

other plot arguments

Value

No return value, just the plot

See Also

GenomicControl, WGassociation

Examples

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))