R/HWIlrPlot.R

Defines functions HWIlrPlot

Documented in HWIlrPlot

HWIlrPlot <- function(X,zeroadj=0.5,...) {
  # Plot the rows of X in isometric logratio coordinates.
  nullvalue <- -sqrt(2/3)*log(2)
  Xilr <- HWIlr(X,zeroadj)
  opar <- par(mar=c(5,5,2,1))
  plot(Xilr[,1],Xilr[,2],pch=19,xlab=expression((1/sqrt(2))*ln(f[AA]/f[BB])),
       ylab=expression((1/sqrt(6))*ln(f[AA]*f[BB]/f[AB]^2)),...)
  vx <- c(par("usr")[1],par("usr")[2])
  vy <- c(nullvalue,nullvalue)
  lines(vx,vy)
  par(opar)
}

Try the HardyWeinberg package in your browser

Any scripts or data that you put into this service are public.

HardyWeinberg documentation built on May 7, 2022, 5:05 p.m.