In rdinnager/fibre: Fast Evolutionary Trait Modelling on Phylogenies using Branch Regression Models

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(fibre)
library(phytools)
library(RRphylo)

set.seed(230045)
## simulated tree
phy <- pbtree(n = 1000, scale = 1, d = 0.2)

L1 <- makeL(phy)[ , -1]
L2 <- make_L(phy)

all.equal(L1, L2)

L11 <- makeL1(phy)[ , -1]
L12 <- make_L(phy, return_nodes = "internal")[ , 1:(phy$Nnode - 1)]

all.equal(L11, L12)

t1 <- make_L(phy, return_nodes = "both")
glimpse(t1)
t2 <- make_L(phy, return_nodes = "both", return_type = "list")
glimpse(t2)

rdinnager/fibre documentation built on Dec. 14, 2024, 10:33 a.m.