NAM (version 1.4.2)

Manhattan: Manhattan plot for Association Studies

Description

Generates a graphical visualization for the output of the function gwas/gwas2.

Usage

## S3 method for class 'NAM':
plot(x,...,alpha=0.05,colA=2,colB=4,find=NULL,FDR=NULL,gtz=FALSE)

Arguments

x

Output of the gwas/gwas2 function.

...

Further arguments passed to or from other methods.

alpha

Numberic. Significance threshold to display in the Manhattan plot.

colA

Color of odd chromosomes in the Manhattan plot.

colB

Color of even chromosomes in the Manhattan plot.

find

Integer. If provided, you can click on the specified number of hits in the Manhattan plot to obtain the name of the markers.

FDR

Null or numeric between zero and one. If provided, it will display the Manhattan plot with Bonferroni threshold by chromosome, adjusted for the specified false discovery rate (FDR). Thus, zero provides the Bonferroni correction.

gtz

Logical. If TRUE, the argument FDR will just take into account markers with p-value Greater Than Zero (GTZ).

Examples

data(tpod)
test=gwas2(y=y,gen=gen,fam=fam,chr=chr)

par(mfrow=c(2,1))

# Example Manhattan 1
plot(x=test,colA=3,colB=1,alpha=0.05/ncol(gen),type="h",main="Example 1: Genome-Wide Bonferroni")

# Example Manhattan 2
Title = "Example 2: FDR 0.25 by Chromosome"
plot(x=test,alpha=0.05,FDR=0.25,gtz=TRUE,pch=20,main=Title)