22.01.2015 Views

Refined Buneman Trees

Refined Buneman Trees

Refined Buneman Trees

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

2.5 Diagonal quartets<br />

Looking at the definition of a quartet and its associated <strong>Buneman</strong> score, we<br />

might observe that the quartet q = ab|cd canbeviewedastwodiagonal quartets,<br />

denoted ab||cd and ab||dc. This intuition is not simple, but we can make the<br />

following rewrite to illustrate this point:<br />

β q = 1 2 (min{δ ac + δ bd ,δ ad + δ bc }−(δ ab + δ cd )) (2.2)<br />

{ 1<br />

= min 2 (−δ }<br />

ab + δ bc − δ cd + δ da )<br />

1<br />

2 (−δ (2.3)<br />

ab + δ bd − δ dc + δ ca )<br />

= min{η ab||cd ,η ab||dc } (2.4)<br />

where η ab||cd = 1 2 (δ bc − δ ab + δ ad − δ cd ) is the score of the diagonal quartet<br />

ab||cd.<br />

Looking at Figure 2.5 we can see that the terms η ab||cd and η ab||dc correspond<br />

to a “tour” in one of the two diagonal quartets. In words, we go through either<br />

a → b → c → d → a or a → b → d → c → a, adding or subtracting chords<br />

as appropriate. We might also notice that diagonal quartets are symmetric<br />

on either side of their central edge, i.e. starting out with the diagonal quartet<br />

ab||cd we can swap a ↔ b and c ↔ d or ab ↔ cd to obtain ba||dc or cd||ab, where<br />

η ab||cd = η ba||dc = η cd||ab and all three diagonal quartets still identify the same<br />

quartet ab|cd.<br />

a<br />

c<br />

b<br />

ab|cd<br />

d<br />

a<br />

c<br />

a<br />

c<br />

b<br />

ab||cd<br />

d<br />

b<br />

ab||dc<br />

d<br />

Figure 2.5: The symmetric quartet and its associated asymmetric diagonal quartets.<br />

To keep an ordering on diagonal quartets we shall define the minimum diagonal<br />

of ab|cd: Leta be the smallest (by index in X) of the species a, b, c, d ∈ X.<br />

ab||cd is the minimum diagonal if η ab||cd

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!