MESSAGE
DATE | 2017-01-26 |
FROM | Ruben Safir
|
SUBJECT | Subject: [Learn] (fwd) Re: Felsenstein Phylogenies
|
From learn-bounces-at-nylxs.com Thu Jan 26 18:50:19 2017 Return-Path: X-Original-To: archive-at-mrbrklyn.com Delivered-To: archive-at-mrbrklyn.com Received: from www.mrbrklyn.com (www.mrbrklyn.com [96.57.23.82]) by mrbrklyn.com (Postfix) with ESMTP id EC4EC16131B; Thu, 26 Jan 2017 18:50:18 -0500 (EST) X-Original-To: learn-at-nylxs.com Delivered-To: learn-at-nylxs.com Received: from mailbackend.panix.com (mailbackend.panix.com [166.84.1.89]) by mrbrklyn.com (Postfix) with ESMTP id 6A721161315 for ; Thu, 26 Jan 2017 18:49:52 -0500 (EST) Received: from panix2.panix.com (panix2.panix.com [166.84.1.2]) by mailbackend.panix.com (Postfix) with ESMTP id 2C324138B6 for ; Thu, 26 Jan 2017 18:49:50 -0500 (EST) Received: by panix2.panix.com (Postfix, from userid 20529) id 6518033CC7; Thu, 26 Jan 2017 18:49:50 -0500 (EST) From: Ruben Safir To: learn-at-nylxs.com User-Agent: tin/2.2.1-20140504 ("Tober an Righ") (UNIX) (NetBSD/7.0.2 (i386)) Message-Id: <20170126234950.6518033CC7-at-panix2.panix.com> Date: Thu, 26 Jan 2017 18:49:50 -0500 (EST) Subject: [Learn] (fwd) Re: Felsenstein Phylogenies X-BeenThere: learn-at-nylxs.com X-Mailman-Version: 2.1.17 Precedence: list List-Id: List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , MIME-Version: 1.0 Content-Type: text/plain; charset="us-ascii" Content-Transfer-Encoding: 7bit Errors-To: learn-bounces-at-nylxs.com Sender: "Learn"
-- forwarded message -- Path: reader1.panix.com!panix!bloom-beacon.mit.edu!bloom-beacon.mit.edu!168.235.88.217.MISMATCH!feeder.erje.net!2.us.feeder.erje.net!feeder.usenetexpress.com!feeder1.iad1.usenetexpress.com!216.166.98.84.MISMATCH!border1.nntp.dca1.giganews.com!border2.nntp.dca1.giganews.com!nntp.giganews.com!buffer2.nntp.dca1.giganews.com!news.giganews.com.POSTED!not-for-mail NNTP-Posting-Date: Thu, 26 Jan 2017 10:56:45 -0600 Subject: Re: Felsenstein Phylogenies Newsgroups: sci.bio.paleontology References: <429a9d55-1da8-b068-0049-4029944f897c-at-mrbrklyn.com> From: John Harshman Date: Thu, 26 Jan 2017 08:56:45 -0800 User-Agent: Mozilla/5.0 (Macintosh; Intel Mac OS X 10.12; rv:45.0) Gecko/20100101 Thunderbird/45.6.0 MIME-Version: 1.0 In-Reply-To: <429a9d55-1da8-b068-0049-4029944f897c-at-mrbrklyn.com> Content-Type: text/plain; charset=windows-1252; format=flowed Content-Transfer-Encoding: 7bit Message-ID: Lines: 80 X-Usenet-Provider: http://www.giganews.com X-Trace: sv3-6N1BtAi0GtL9xFoAEdQ5dyYBUJwBnQGSBJGFnimSxO082lAFdze6JFfeGPyIp8df4YMkBVzNfAJkwqt!yjZKodV1Kpyb2by4cFGht/qWAXeufhwTOUxj7p+J68vKkyuHPZbcm/ywZIM0j/DKAyEt5Ll8ejA/ X-Complaints-To: abuse-at-giganews.com X-DMCA-Notifications: http://www.giganews.com/info/dmca.html X-Abuse-and-DMCA-Info: Please be sure to forward a copy of ALL headers X-Abuse-and-DMCA-Info: Otherwise we will be unable to process your complaint properly X-Postfilter: 1.3.40 Bytes: 4288 X-Original-Bytes: 4226 Xref: panix sci.bio.paleontology:67128
On 1/26/17 3:46 AM, ruben safir wrote: > On 01/25/2017 09:35 PM, John Harshman wrote: >> On 1/25/17 3:53 PM, Ruben Safir wrote: >>> Does anyone have the above text handy? I think what he wrote with >>> regard to Subtree Programming and Grafting is incorrect. >>> >>> If you have 2 subtrees n1 and n2,the number of neighbors should be (2n1 >>> -4) * (2n2 - 4) --- not addition >>> >>> each spot has 2n-3 - 1 permutations. >>> >>> He doesn't explain what external branches are either. >>> >>> >> The book you're referring to is called Inferring Phylogenies and the >> procedure you're talking about is called subtree *pruning* and >> *regrafting*. The number of rearrangements given a particular subtree >> should be equal to the number of branches on the second subtree, which >> is twice the number of taxa minus 3. >> >> I don't currently have a copy handy. Please explain more clearly what >> Felsenstein says about it and what you think it should say. > > Correct, what it does say is that once you divide the tree there would > be 2n1 - 3 - 1 reassertion points for the tree. Then after that he is > not clear to me. He says > > "In fact considering both subtrees (no having n1 species and the one > having n2 species, there are > (2n1-3-1) + (2n2-3-1 ) = (2n-3-1) = 2n-8 > neighbors generated at each interior branch." > > This assumes n1 + n2 = n.
Which it must.
> I guess that is all the possible combinations assuming the same > attachment locations for the trees, examining one tree at a time.
Not sure what you meant by that.
> Then he states that external nodes (which is not defined) is 2n-6. > Without proof I'll accept that for a moment (and I think it corresponds > to binary tree theory), but I'm not sure that an exterior node is. That > is a node that connects to leafs?
I don't know what 2n-6 is, based on your description, but from the formula below it appears to be the number of subtrees that could be attached to any terminal branch of the tree, i.e. the number of subtrees not containing that branch.
> Finally, the last unclear sentence, at least to me, states: > > "Thus, as there are n exterior branches on an unrooted bifurcating tree > and n-3 interior branches, the total number of neighbors examined by SPR > will be > > n(2n-6)+(n+3)(2n-8)" > > That is where he lost me.
Why? Did his formula have n+3 when it should have had n-3? Otherwise I see no problem.
> Then he follows up > and he says tha there are 288 neighbors for n=11 > and > "Of course, 2(n-3)=16 of them are the same as NNI"
Why is that a problem?
> For TBR he says that there is no general formula for the number of > neighbors that will be examine. That made be stand on the edge of my > seat? Say what? Then what are we doing?
I do not understand your problem there.
For all of these, you need to explain what you think the problem is.
-- end of forwarded message -- _______________________________________________ Learn mailing list Learn-at-nylxs.com http://lists.mrbrklyn.com/mailman/listinfo/learn
|
|