Dear Russell.
For some reason, my emails don't seem to be getting through to the whole
list -- so my apologies if this ends up being redundant with another reply.
In short, if you supply a vector of tip labels that are identical (e.g.,
c("Cricetomys_emini","Cricetomys_emini")), ape::getMRCA returns the
parent node of the corresponding taxon -- rather than the tip index. I
didn't know that, but now I do.
You can switch to using phytools::fastMRCA. fastMRCA takes two tip label
arguments instead of one -- so in your case the function call would be
fastMRCA(tree,"Cricetomys_emini","Cricetomys_emini").
You can also just use the base R function which. In this case, you would
do which(tree$tip.label=="Cricetomys_emini"). This should have the same
effect.
You do seem to be using the arguments 'where' and 'position' correctly
in phytools::bind.tip. 'where' should point to the node immediately
*above* (i.e., tipward) of the place you want to attach the new tip;
while 'position' should be the depth *below* that node.
Here's an update to your mapply call using fastMRCA instead of getMRCA.
newtree<-mapply(bind.tip,tree=tree,where=lapply(tree,fastMRCA,
sp1="Cercopithecus_mitis",sp2="Cercopithecus_mitis"),
MoreArgs=list(tip.label="Cercopithecus_albogularis",
position=0.59,edge.length=0.59),SIMPLIFY=FALSE)
class(newtree)<-"multiPhylo"
Please let me (& the list) know if this solution works as expected.
All the best, Liam
Liam J. Revell
University of Massachusetts Boston [Assoc. Prof.]
Universidad Católica de la Ssma Concepción [Adj. Res.]
Web & phytools:
http://faculty.umb.edu/liam.revell/, http://www.phytools.org,
http://blog.phytools.org
Academic Director UMass Boston Chile Abroad:
https://www.umb.edu/academics/caps/international/biology_chile
U.S. COVID-19 explorer web application:
https://covid19-explorer.org/
On 6/4/2021 1:17 AM, Russell Engelman wrote:
EXTERNAL SENDER
Dear Everyone,
I think these are all really good suggestions. I tried this with my data
(specifically the mapply with getMRCA) and it sort of worked. Sort of in
that the mapply did apply the transformation to all 100 trees, but I
noticed that after I applied it I got some funky results. Here is
another example with /Cricetomys ansorgei /and some associated code...
hbltree <-read.nexus(file=" hbltree .nex")
class(hbltree)<-"multiPhylo"
newtrees<-mapply(bind.tip,tree=hbltree,where=lapply(hbltree,getMRCA,
tip=c("Cricetomys_emini","Cricetomys_emini")),
MoreArgs=list(tip.label="Cricetomys_ansorgei",
position=3.282,edge.length=3.282),SIMPLIFY=FALSE)
However I noticed that when I printed the phylogeny the function for
some reason does not place /C. ansorgei /as the sister taxon to /C.
emini/, instead it places it as sister to /C. emini/ + /C. gambianus.
/That part I am not sure why this is because getMRCA should find the
terminal tip of /C. emini/, given that it is finding the MRCA of the
same taxon.
It also does not place /C. ansorgei /at the tip of the tree. This might
be a side effect of the function selecting the wrong node to insert the
taxon on, but I think part of it as well is that I'm not familiar with
how to handle the position and edge.length arguments in this case. I
know from consulting the primary literature that /C. ansorgei /is
recoded to have diverged from its nearest sister taxon, /C. emini/,
about 3.282 Ma. So I know the 'edge.length' argument should be set to
3.282. I've often struggled a lot with how to efficiently set the
'position' argument afterwards.
Is it possible I am calling the wrong node or something to insert the
taxon at?
Sincerely,
Russell
On Thu, Jun 3, 2021 at 10:47 AM Liam J. Revell <liamjrev...@gmail.com
<mailto:liamjrev...@gmail.com>> wrote:
Dear Russell et al. (sorry if this is a re-send: my message yesterday
seems to have been blocked by the listserv)
Using a for loop is a great idea! Highly underrated in R, IMO.
However, for future reference, the reason that your code didn't work
with lapply is because the list you're 'applying' over (tree) also
appears among the arguments!
If you want to use apply-family functions instead of a for loop (just,
say, for fun) then you have two basic options: you can write a custom
function; or you can use mapply.
Here's some (untested) code to do it.
## first, using a custom function & lapply:
foo<-function(tree) bind.tip(tree,
tip.label="Cercopithecus_albogularis",
position=0.59,edge.length=0.59,
where=getMRCA(tree,tip=c("Cercopithecus_mitis",
"Cercopithecus_mitis")))
newtree<-lapply(tree,foo)
class(newtree)<-"multiPhylo"
## now, using mapply:
newtree<-mapply(bind.tip,tree=tree,where=lapply(tree,getMRCA,
tip=c("Cercopithecus_mitis","Cercopithecus_mitis")),
MoreArgs=list(tip.label="Cercopithecus_albogularis",
position=0.59,edge.length=0.59),SIMPLIFY=FALSE)
class(newtree)<-"multiPhylo"
(Code is not guaranteed! I don't have the data file, so I didn't
actually test it -- but something like this ought to work.)
Regardless, I recommend using ape::getMRCA (or phytools::fastMRCA)
because otherwise you're computing an N x N matrix in each iteration of
your function call just to get one node index.
Good luck! All the best, Liam
--
Liam J. Revell
University of Massachusetts Boston
Universidad Católica de la Ssma Concepción
Web & phytools:
http://faculty.umb.edu/liam.revell/
<http://faculty.umb.edu/liam.revell/>, http://www.phytools.org
<https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.phytools.org%2F&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154173973%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=gre%2F0WliIwUlZ%2Bir%2B8718MENatjj5HMH1wvGoBWCJMg%3D&reserved=0>,
http://blog.phytools.org
<https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fblog.phytools.org%2F&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154183969%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=GcRchX%2FnKifwnyLCgNB6DSr0qO7nMF%2FpAdJw23AWIT0%3D&reserved=0>
Academic Director UMass Boston Chile Abroad:
https://www.umb.edu/academics/caps/international/biology_chile
<https://www.umb.edu/academics/caps/international/biology_chile>
U.S. COVID-19 explorer web application:
https://covid19-explorer.org/
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fcovid19-explorer.org%2F&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154183969%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=C4N8vAhdHVUy72r8DmrZrXBoZW6xUwuK7ZlBI%2B13Dxs%3D&reserved=0>
On 6/2/2021 8:18 PM, Nathan Upham wrote:
> Hi Russell:
>
> Glad to hear you’re using the VertLife mammal trees — they are
built on a taxonomy of 5,911 species of which only 4,098 are sampled
for DNA, so there is already a ~30% chunk that is placed using
taxonomic constraints and birth-death branch lengths as sampled
during the estimation of 28 Bayesian patch clades.
>
> Adding additional species described since the 2015 cutoff of that
VertLife taxonomy makes sense (e.g., up to ~6,500 species on
mammaldiversity.org
<https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmammaldiversity.org%2F&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154183969%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=1Ua4UnPAguLQi%2BRv7AOP00PBg3Hbuzovhqk%2FGow6uAI%3D&reserved=0>).
However, keep in mind that they will not have birth-death estimated
branch lengths, but rather more likely be added as a polygamy to
given clade and then randomly resolved.
>
> Given the sample code you provided, the key thing you’ll want to
do is run a *loop* rather than using lapply, so that you can specify
a given tree each time, e.g.:
>
> newtrees<-vector(“list”,length(trees))
> for(j in 1:length(trees)){
> newtrees[[j]] <- bind.tip(tree=trees[[j]],
tip.label="Cercopithecus_albogularis”, position=0.59,edge.length =
0.59,
where=mrca(tree[[j]])["Cercopithecus_mitis","Cercopithecus_mitis"])
> }
>
> I also wrote some code to prune mammal trees and add extinct
Caribbean species, which uses a similar approach of making
polytomies and randomly resolving them — here is the repo:
>
https://github.com/n8upham/CaribbeanExtinctions-WTWTW/tree/master/mamPhy_pruningCode
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fn8upham%2FCaribbeanExtinctions-WTWTW%2Ftree%2Fmaster%2FmamPhy_pruningCode&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154193960%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=tbsBqp83Ylw6Icw9%2FIo4FH6gb7YwF8wU%2FdDTsSHGgCg%3D&reserved=0>
> And here is the code file:
>
https://github.com/n8upham/CaribbeanExtinctions-WTWTW/blob/master/mamPhy_pruningCode/pruningCode_MamPhy-to-CaribbeanTaxa.R
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fn8upham%2FCaribbeanExtinctions-WTWTW%2Fblob%2Fmaster%2FmamPhy_pruningCode%2FpruningCode_MamPhy-to-CaribbeanTaxa.R&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154193960%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=iRwABMz4pk73VtYoqL0hipTVO%2FOabNsazgeqF%2FpieR0%3D&reserved=0>
>
> Hope that helps,
> —nate
>
>
>
>
========================================================================
> Nathan S. Upham, Ph.D. (he/him)
> Assistant Research Professor & Associate Curator of Mammals
> Arizona State University, School of Life Sciences
> ~> Check out the new Mammal Tree of Life
<http://vertlife.org/data/mammals/
<https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fvertlife.org%2Fdata%2Fmammals%2F&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154203959%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=l%2F%2FF6%2FBXi53usTw1IVj54Tw%2BHspBK6qiCn3uJS9wles%3D&reserved=0>>
and the Mammal Diversity Database <https://mammaldiversity.org/
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmammaldiversity.org%2F&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154203959%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=SWND%2FEwc2w7FIR2pKk6Zu%2FtB0LpPAS0ICgesFdV4pDU%3D&reserved=0>>
>
> Research Associate, Yale University (Ecology and Evolutionary
Biology)
> Research Associate, Field Museum of Natural History (Negaunee
Integrative Research Center)
> Chair, Biodiversity Committee, American Society of Mammalogists
> Taxonomy Advisor, IUCN/SSC Small Mammal Specialist Group
>
> personal web: n8u.org
<https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fn8u.org%2F&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154213948%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=Pj%2Fq2Gm%2Bw1NMMhiB3rxKoJaG564aftM0VTg9I0QQq5Y%3D&reserved=0>
| Google Scholar
<https://scholar.google.com/citations?hl=en&user=zIn4NoUAAAAJ&view_op=list_works&gmla=AJsN-F6ybkfthmTdjTpow6sgMhWKn1EKcfNtmIF_wzZcev7yeHuEu5_aolFS85rWiVRHpiQgbwg43i6eS6kArrabLdFL4bntzUSRmlRP2CW4lbZqeEcColw
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fscholar.google.com%2Fcitations%3Fhl%3Den%26user%3DzIn4NoUAAAAJ%26view_op%3Dlist_works%26gmla%3DAJsN-F6ybkfthmTdjTpow6sgMhWKn1EKcfNtmIF_wzZcev7yeHuEu5_aolFS85rWiVRHpiQgbwg43i6eS6kArrabLdFL4bntzUSRmlRP2CW4lbZqeEcColw&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154213948%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=wxQgrOr0l0O9Rqfaz8WMRFw4%2FG5MvFwS4okRDtQFbn4%3D&reserved=0>>
| ASU profile <https://isearch.asu.edu/profile/3682356
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fisearch.asu.edu%2Fprofile%2F3682356&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154223944%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=YWuFvYVeBOEFFeq3QB3epfrIT5FBBDBPC%2F1tWZcJ2r0%3D&reserved=0>>
> e: nathan.up...@asu.edu <mailto:nathan.up...@asu.edu> | Skype:
nate_upham | Twitter: @n8_upham <https://twitter.com/n8_upham
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Ftwitter.com%2Fn8_upham&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154223944%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=0A0pE%2B2jkJFzNjV2lxs8AkGRDALmEbl4dzyBpgGAM7A%3D&reserved=0>>
>
========================================================================
>
>
>
>> On Jun 2, 2021, at 4:19 PM, Eliot Miller <eliot.is...@gmail.com
<mailto:eliot.is...@gmail.com>> wrote:
>>
>> Hi Russell,
>>
>> A package I wrote a while back should be able to do that fairly
easily.
>>
https://urldefense.com/v3/__https://github.com/eliotmiller/addTaxa__;!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQrzxl2aIw$
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.com%2Fv3%2F__https%3A%2F%2Fgithub.com%2Feliotmiller%2FaddTaxa__%3B!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQrzxl2aIw%24&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154233943%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=fIBGFZSNi6mvysStlfTCXakm%2F2vpqwgx6zxkTBCf7jg%3D&reserved=0>
The only paper it's described in
>> remains
https://urldefense.com/v3/__https://bsapubs.onlinelibrary.wiley.com/doi/full/10.3732/ajb.1500195__;!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQp7PnRRHg$
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.com%2Fv3%2F__https%3A%2F%2Fbsapubs.onlinelibrary.wiley.com%2Fdoi%2Ffull%2F10.3732%2Fajb.1500195__%3B!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQp7PnRRHg%24&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154233943%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=%2BUOm4n9HRuG7eu7EWDyPLC3%2FYzN7GJcMRdGOuDqfQ5Y%3D&reserved=0>
>> It's a wrapper for bind.tip, with some additional stuff. You
basically
>> would give it a taxonomic file where you identify the clades you're
>> interested in (e.g. both of those Cercopithecus species could be
named some
>> unique clade name and off you go, it'd add the missing one to
the other),
>> then lapply that whole addTaxa command over the list of trees in
>> multiPhylo. At some point I made laser a dependency, and it's
possible I
>> left it in that state. If that's the case, you can still get
laser from old
>> CRAN mirrors I believe. Let me know if you want more help.
>>
>> Best,
>> Eliot
>>
>> On Wed, Jun 2, 2021 at 7:05 PM Russell Engelman
<neovenatori...@gmail.com <mailto:neovenatori...@gmail.com>>
>> wrote:
>>
>>> Dear R-sig-phylo,
>>>
>>> I have been working with a mammalian phylogeny I recently
downloaded from
>>> VertLife
(https://urldefense.com/v3/__http://vertlife.org/phylosubsets/__;!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQoEdEnMAg$
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.com%2Fv3%2F__http%3A%2F%2Fvertlife.org%2Fphylosubsets%2F__%3B!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQoEdEnMAg%24&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154243932%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=IK3uwLncKazEkIwiRuFIOi%2BKgTuBvgHxtg5yFxsPDhY%3D&reserved=0>
). Unfortunately, the phylogeny
>>> is missing a large number of species, so I am trying to
manually add these
>>> taxa to the phylogeny. I have a series of 100 trees that I am
using to do
>>> things such as test for phylogenetic signal. I know how to use
bind.tip to
>>> add new taxa to a single tree, but I am having more trouble with a
>>> multiPhylo object. I am primarily adding these taxa by placing
them as
>>> sister to their nearest included relative (since most of them
are elevated
>>> former subspecies), but the issue here is that in the 100 trees
in the
>>> multiPhylo object the node representing the taxon to bind these
taxa to is
>>> not the same across all trees due to shifting topologies.
>>>
>>> This is an example of the code I have been using, in which
"tree" is the
>>> tree object. This works for a single 'phylo' tree but not
'multiphylo'.
>>>
>>> ```
>>>
newtree<-lapply(tree,bind.tip,tip.label="Cercopithecus_albogularis",
>>> position=0.59,edge.length = 0.59,
>>>
>>> where=mrca(tree)["Cercopithecus_mitis","Cercopithecus_mitis"])
>>> ```
>>>
>>> Now, this code will not work, but I know exactly why: 'tree' is a
>>> multiPhylo object and so the 'where' argument cannot find the
node for the
>>> terminal taxon. However, the issue is how can I tell R to
repeat this
>>> 'where' argument for each of the 100 trees, since the node in
question is
>>> not identical across these trees? Is there an easier way to do
this than
>>> using the 'mrca' call for each terminal taxon? I've noticed
adding a 'mrca'
>>> argument also increases computation time and if I am
reinventing the wheel
>>> it would be nice to know if I am overthinking things.
>>>
>>> Sincerely,
>>> Russell
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> _______________________________________________
>>> R-sig-phylo mailing list - R-sig-phylo@r-project.org
<mailto:R-sig-phylo@r-project.org>
>>>
https://urldefense.com/v3/__https://stat.ethz.ch/mailman/listinfo/r-sig-phylo__;!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQqcmU5jCQ$
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.com%2Fv3%2F__https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-phylo__%3B!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQqcmU5jCQ%24&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154243932%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=%2FeIxhowW77%2BY6ZrPI30ro3VS7UevBvBwgPf1HIZ5bZ8%3D&reserved=0>
>>> Searchable archive at
>>>
https://urldefense.com/v3/__http://www.mail-archive.com/r-sig-phylo@r-project.org/__;!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQq_N5pSwg$
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.com%2Fv3%2F__http%3A%2F%2Fwww.mail-archive.com%2Fr-sig-phylo%40r-project.org%2F__%3B!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQq_N5pSwg%24&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154253924%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=TPxhOEKba7VjjSUKZVBfz831%2FsjficzJ9gi6fHXnix0%3D&reserved=0>
>>>
>>
>> [[alternative HTML version deleted]]
>>
>> _______________________________________________
>> R-sig-phylo mailing list - R-sig-phylo@r-project.org
<mailto:R-sig-phylo@r-project.org>
>>
https://urldefense.com/v3/__https://stat.ethz.ch/mailman/listinfo/r-sig-phylo__;!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQqcmU5jCQ$
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.com%2Fv3%2F__https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-phylo__%3B!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQqcmU5jCQ%24&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154253924%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=gYZwrt8Hfa1eMNZ1grsnmfbUiM%2BeHtq44BpE3khsBKE%3D&reserved=0>
>> Searchable archive at
https://urldefense.com/v3/__http://www.mail-archive.com/r-sig-phylo@r-project.org/__;!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQq_N5pSwg$
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.com%2Fv3%2F__http%3A%2F%2Fwww.mail-archive.com%2Fr-sig-phylo%40r-project.org%2F__%3B!!IKRxdwAv5BmarQ!OZj7-dFRbxvUothKjSj6hr9B0eXscAO6LVWi1-a0w3J_PxlDqvsFDNb0lQq_N5pSwg%24&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154263920%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=%2FN3RohP0qJDPRdbGHAjrc%2FxHDWdO2%2FzGr8gdHrnYFvg%3D&reserved=0>
>
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-phylo mailing list - R-sig-phylo@r-project.org
<mailto:R-sig-phylo@r-project.org>
> https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
<https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-phylo&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154263920%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=JXt6aXIQtEI5axP1sHQegMkCwoSKUnk8HecALpgYgsM%3D&reserved=0>
> Searchable archive at
http://www.mail-archive.com/r-sig-phylo@r-project.org/
<https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.mail-archive.com%2Fr-sig-phylo%40r-project.org%2F&data=04%7C01%7Cliam.revell%40umb.edu%7C31aef1ba99c742b14f6b08d927182bbc%7Cb97188711ee94425953c1ace1373eb38%7C0%7C1%7C637583807154273919%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=GcQZ%2FkOAd7MotgR1MKB53D3mPpv8m4OdG1uFGQuARxo%3D&reserved=0>
>
_______________________________________________
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/