9.19 Plotting PBS trees
Another useful way to visualize PBS is by comparing trees for the top PBS outliers to the genome-wide average tree. Run the code blocks below to plot these trees:
# create average tree
<- rtree(n = 3,
tr_mean rooted = FALSE,
br = c(mean(pbs$T_rps_png),
mean(pbs$T_rps_chb),
mean(pbs$T_png_chb)))
# plot average tree
ggtree(tr_mean, layout = "daylight") +
geom_treescale(width = 0.1) +
geom_tiplab(label = c("RPS", "PNG", "CHB"))
## Average angle change [1] 0.407407407407407
## Average angle change [2] 0
# create tree for top snp
<- rtree(n = 3,
tr_top rooted = FALSE,
br = c(pbs[1,]$T_rps_png,
1,]$T_rps_chb,
pbs[1,]$T_png_chb))
pbs[
# plot top snp tree
ggtree(tr_top, layout = "daylight") +
geom_treescale(width = 0.1) +
geom_tiplab(label = c("RPS", "PNG", "CHB"))
## Average angle change [1] 0.407407407407407
## Average angle change [2] 0