Functions for inspecting population substructure

Description

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

Usage

1
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

A plot

See Also

GenomicControl, WGassociation

Examples

1
2
3
4
5
6
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))

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.