xy <- locator(1)
Ntip <- .last_plot.phylo$Ntip
d <- sqrt((xy$x - .last_plot.phylo$xx)^2 +
(xy$y - .last_plot.phylo$yy)^2)
NODE <- which.min(d)
xy <- locator(1)
Ntip <- .last_plot.phylo$Ntip
d <- sqrt((xy$x - .last_plot.phylo$xx)^2 +
(xy$y - .last_plot.phylo$yy)^2)
NODE <- which.min(d)