trees_with_root_edge.R
In nLTT: Calculate the NLTT Statistic

## ---- include = FALSE---------------------------------------------------------
# Set high warning levels
options(warn = 2)

## -----------------------------------------------------------------------------
set.seed(42)
tree1 <- ape::rcoal(2)
tree1$edge.length <- tree1$edge.length / tree1$edge.length[1] # nolint ape variable name

## -----------------------------------------------------------------------------
tree2 <- tree1
tree2$root.edge <- 1 # nolint ape variable name

## ----fig.width = 3, fig.height = 3--------------------------------------------
ape::plot.phylo(tree1, root.edge = TRUE)
ape::add.scale.bar() #nolint

## ----fig.width = 3, fig.height = 3--------------------------------------------
ape::plot.phylo(tree2, root.edge = TRUE)
ape::add.scale.bar() #nolint

## ----fig.width = 7, fig.height = 7--------------------------------------------
nLTT::nltt_plot(tree1, xlim = c(0, 1), ylim = c(0, 1))
nLTT::nltt_lines(tree2, col = "red")

## -----------------------------------------------------------------------------
print(nLTT::nLTTstat_exact(tree1, tree2, ignore_stem = FALSE))

## -----------------------------------------------------------------------------
print(nLTT::nLTTstat_exact(tree1, tree2, ignore_stem = TRUE))