approx.alpha.mle: #' @keywords internal TIMBR.scan.plot <-...

In wesleycrouse/TIMBR: Tree-based Inference of Multiallelism via Bayesian Regression

Description

#' @keywords internal TIMBR.scan.plot <- function(TIMBR.scan.output, col="black", add=F, ylim=NULL) plot.data <- c(TIMBR.scan.output$intervals[,1], TIMBR.scan.output$intervals[,2]-1) plot.data <- cbind(plot.data, c(TIMBR.scan.output$ln.BFs, TIMBR.scan.output$ln.BFs)) plot.data <- plot.data[order(plot.data[,1]),]

if (add) lines(plot.data[,1], plot.data[,2], col=col) else if (is.null(ylim)) plot(plot.data[,1], plot.data[,2], type="l", col=col, ylab="lnBF", xlab="Location") else plot(plot.data[,1], plot.data[,2], type="l", col=col, ylab="lnBF", xlab="Location", ylim=ylim)

Usage

approx.alpha.mle(TIMBR.output, ln.BF.full = NULL)

wesleycrouse/TIMBR documentation built on Feb. 19, 2021, 7:31 a.m.