HWGenotypePlot {HardyWeinberg}

Scatter plot of the genotype frequencies

Description

HWGenotypePlot makes a scatterplots of the AB or BB frequency versus the AA frequency and represents a blue curve indicating the Hardy-Weinberg equilibrium condition.

Usage

HWGenotypePlot(X, plottype = 1, xlab = expression(f[AA]), ylab =
ifelse(plottype == 1, expression(f[AB]), expression(f[BB])), asp = 1,
pch = 19, xlim = c(0, 1), ylim = c(0, 1), cex = 1, cex.axis = 2, cex.lab = 2, ...)

Arguments

X

A matrix of genotype counts or frequencies with three columns (AA, AB, BB)

plottype

plottype=1 produces a plot of AB versus AA, plottype=2 produced a plot of BB versus AA.

xlab

A label for the x axis

ylab

A label for the y axis

asp

Aspec ratio (1 by default)

pch

Plotting charachter (19 by default)

xlim

Limits for the x axis (0-1 by default)

ylim

Limits for the y axis (0-1 by default)

cex

Character expansion factor (1 by default)

cex.axis

Character expansion factor for the axes (2 by default)

cex.lab

Character expansion factor for labels of axis (2 by default)

...

Additional arguments for the plot function

Value

NULL

Author(s)

Jan Graffelman [e-mail]

See Also

HWTernaryPlot

Examples

n <- 100 # sample size
m <- 100 # number of markers
Xc <- HWClo(HWData(n,m))
HWGenotypePlot(Xc,plottype=1,main="Heterozygote-homozygote scatterplot")