Branching and transport

The band decides the answer

A fit over a whole tree's junctions returns the exponent the tree was built at, even though most of its junctions are from bands that on their own return 1.6. Least squares is already weighting by leverage. The dangerous sample is not the mixed one — it is the one a person can reach.

Worth reading first: Fitting the exponent · The cube law · How many plants would it take.

The previous phase found that a sample of fifty lopsided junctions, taken from a tree built at an exponent of exactly three, returns 1.69 with an interval of [1.56, 1.79] — excluding three, excluding two, and only a fifth wider than a good sample’s. The cause is a selection nobody applies on purpose: at a daughter ratio of a twentieth the parent is genuinely thicker than its larger daughter by 0.11%, against a measurement noise of 2.8%, so about half the measured junctions come out with the parent thinner than a branch it carries. Those have no exponent at all, every estimator drops them silently, and the half that survives is the half where noise inflated the parent.

The obvious next worry follows immediately and it is the one this essay tests.

A real tree is not a band of daughter ratios. It is all of them at once, and the lopsided ones vastly outnumber the even ones — a tree has one trunk fork and hundreds of twigs. So a single fit over a whole tree ought to be dragged towards the twigs’ 1.7, and the amount it is dragged ought to depend on how many twigs happened to be in the sample.

It is not dragged

Two hundred junctions spanning daughter ratios from 0.05 to 1, from trees built at an exponent of exactly three with the standard two per cent measurement error, fit 2.85 to 2.99 across five seeds. The twig band alone, from the same trees with the same fitter, gives 1.60 to 1.73.

A tree built at an exponent of 3, measured band by bandSix bands of daughter ratio, 300 junctions in each, all from trees built at an exponent of exactly 3 with the same 2% measurement error. The median implied exponent falls from 3.00 at an even fork to 1.58 at a twig, and the share of junctions that have any exponent at all falls from 100% to 53% over the same range. The rule across the top is a single fit over 200 junctions spanning the whole range: 2.85. A mixed sample is safe because least squares already weights by leverage — 42% of it sits in the most symmetric band and 0.06% in the twigs. The sample that is not safe is the one a person can reach.11.5022.5033.50012345daughter ratio of the junctions measuredexponent implied, and the share of junctions that imply one0.9–1100%0.7–0.9100%0.5–0.799%0.3–0.578%0.15–0.356%0.05–0.1553%all 200 junctions together: 2.85the tree was built at 3bars are the share with any exponent · whiskers are the tenth and ninetieth percentilesmeasurement error 2%
Fig. 1 Six bands of daughter ratio, three hundred junctions in each, from trees built at an exponent of three. The dots are the median exponent each band implies on its own and they fall from 3.00 at an even fork to 1.58 at a twig; the bars behind them are the share of junctions that have any exponent at all, falling from 100% to 53%. The rule across the top is a single fit over all two hundred junctions of a whole tree, and it is on the true value.

The prediction is refuted, and the refutation is more useful than the prediction would have been, because the reason is a property of least squares rather than a happy accident.

The leverage does it

The residual a junction contributes to the fit moves with the exponent at a rate set by how far apart the candidate exponents’ predictions are at that junction — which is exactly the quantity the previous essay measured as the separation, and which collapses as the cube of the daughter ratio. At an even fork the two rules predict parent radii 0.154 apart; at a daughter ratio of a twentieth, 0.0012 apart.

A least-squares fit weights each observation by how much moving the parameter changes its residual, which is the square of that separation. So the weights across a tree are in the ratio 0.154² to 0.0012², which is more than sixteen thousand to one.

Computed band by band, the shares are:

  • daughter ratio 0.9–1.0 — 41.8% of the tree’s leverage
  • 0.7–0.9 — 32.7%
  • 0.5–0.7 — 18.3%
  • 0.3–0.5 — 6.2%
  • 0.15–0.3 — 1.0%
  • 0.05–0.15 — 0.06%

Three quarters of a whole tree’s leverage sits in the junctions with a daughter ratio above 0.7, and the twigs carry six hundredths of one per cent of it.

So the estimator is already doing the right thing, without anybody choosing it. Adding a hundred twigs to a sample of ten even forks changes the answer by essentially nothing, because the twigs contribute essentially nothing.

The uninformative sample does not say so — it says something elseAbove: how many junctions of a given asymmetry it takes to distinguish an exponent of 3 from an exponent of 2, at 2% measurement error on each radius. An even fork needs one; a junction whose small daughter is a twentieth of the large one needs 2195. Below: 50 junctions from each end of the range, on a synthetic tree built at exactly 3. The even forks return 3.01; the twigs return 1.69, with an interval no wider — because 47% of them measure as a parent thinner than its own larger daughter, and dropping those keeps only the half where the noise ran the right way.1101001000100000.2000.4000.6000.8001daughter ratio — the smaller branch over the largerjunctions needed to tell Murray's rule from Da Vinci'san even fork — one is enougha twig — 2195the same tree, built at an exponent of 3, measured twicethe truth — 350 even forks3.0150 twigs1.692% error per radius · 47% of twigs impossiblemargin 0.11% against 2.8% noise
Fig. 2 The separation curve the leverage is computed from, and the two samples the previous phase measured. The top panel is the arithmetic: how many junctions of a given asymmetry it takes to tell the two rules apart, running from one to two thousand across the range. The leverage is the square of that separation, which is why the range in weight is a factor of sixteen thousand rather than two thousand.

Which relocates the danger

The result does not make the previous phase’s finding harmless. It relocates it, and the new location is worse.

The bias appears when a sample is restricted to the uninformative band, and the question is when that happens in practice. It happens whenever the sample is selected by something correlated with the daughter ratio — and the most obvious such thing is whether a person can reach the junction and put callipers on it.

On a standing tree, twigs are at eye level and in the hand. The trunk fork is twenty metres up. A sample gathered by walking through a wood with callipers is not a mixed sample; it is a twig sample, and a twig sample returns 1.7 with a tight interval from data generated at 3.

So the practical advice inverts. It is not “take a big sample and let the mixture sort itself out” — that works, but it is not what anybody does. It is “measure the few junctions that are worth measuring”, which on this arithmetic means the symmetric ones, and one of them settles the question at the top of the range.

A branching tree in which every junction obeys the cube lawr₀³ = r₁³ + r₂³ at all 63 junctions, to 2e-16. The widths in the drawing are the radii the law gives, not widths chosen to look right.branch angle 32° · daughter ratio 1.60127 segments · 63 junctions checked
Fig. 3 A tree built to Murray’s law at a daughter ratio of 0.7, which is inside the informative band. Almost all of the leverage in a fit sits in junctions like these near the base, and almost all of the junctions on a real tree are the lopsided ones further out.

The retention rate, which is the diagnostic

The previous phase ended with a recommendation: report attempted and retained, not retained alone. This essay measures what that number does across the range, and it tracks the bias exactly.

  • daughter ratio 0.9–1.0 — median exponent 3.00 — 100% of junctions have one
  • 0.7–0.9 — 2.98 — 100%
  • 0.5–0.7 — 2.95 — 99%
  • 0.3–0.5 — 2.73 — 78%
  • 0.15–0.3 — 2.11 — 56%
  • 0.05–0.15 — 1.58 — 53%

Two columns of the same fact. Where the parent’s genuine margin over its larger daughter exceeds the measurement noise, every junction has an exponent and the median is the true one. Where it does not, half the junctions are impossible, the survivors are the ones noise inflated, and the median collapses.

The value of the retention column is that it needs no model. A person reporting a branching exponent knows how many junctions they measured and how many the fitter could use; the ratio of those two is one integer over another and it says whether the answer is a measurement or a selection. A sample that retained 53% of its junctions is not a sample with a slightly larger error bar. It is a sample whose answer is wrong by more than the distance between the hypotheses being tested.

The uninformative sample does not say so — it says something elseAbove: how many junctions of a given asymmetry it takes to distinguish an exponent of 3 from an exponent of 2, at 2% measurement error on each radius. An even fork needs one; a junction whose small daughter is a twentieth of the large one needs 2195. Below: 60 junctions from each end of the range, on a synthetic tree built at exactly 3. The even forks return 2.99; the twigs return 1.66, with an interval no wider — because 47% of them measure as a parent thinner than its own larger daughter, and dropping those keeps only the half where the noise ran the right way.1101001000100000.2000.4000.6000.8001daughter ratio — the smaller branch over the largerjunctions needed to tell Murray's rule from Da Vinci'san even fork — one is enougha twig — 2195the same tree, built at an exponent of 3, measured twicethe truth — 360 even forks2.9960 twigs1.662% error per radius · 47% of twigs impossiblemargin 0.11% against 2.8% noise
Fig. 4 The previous phase’s version of the same finding: fifty twig junctions from a tree built at exactly three, returning 1.69 with an interval that excludes both hypotheses. The retention rate is the number that would have shown it, and nothing in the output reports it unless somebody chooses to.

Why the median and not the fit, band by band

The table above reports a median of per-junction exponents rather than a fitted exponent per band, and the choice is not cosmetic.

A per-junction exponent is the exact solution of r₀^k = r₁^k + r₂^k for that junction, which exists whenever the parent is thicker than its larger daughter and does not otherwise. So the median of the per-junction exponents in a band is a statement about that band’s junctions individually, and the count of junctions that have one is available as a by-product — which is the retention column.

A fitted exponent per band would hide both. It returns a number whether or not the junctions in the band are individually solvable, because the fit minimises a residual rather than solving an equation, and a residual can be minimised over data containing impossible junctions without any indication that they were there.

That is the same distinction the site keeps running into in other subjects: an instrument that always returns a value is not the same as an instrument that measures. The angle recovery refuses when no divergence makes a pair the closest. The comb readout refuses when no spacing wins by a band. The per-junction exponent refuses when a branch is thinner than what it carries — and the count of refusals is the diagnostic.

The exponent fitted from the junctions, rather than assumedSweeping k and asking where r₀ᵏ = Σ rᵢᵏ holds best gives 3.000 — Murray's 3, recovered rather than imposed.02468234exponent khow badly r₀ᵏ = Σ rᵢᵏ fails at that kk = 3127 junctionsfitted k = 3.000
Fig. 5 The site’s fitter recovering an exponent from a synthetic tree, which is the control this whole thread depends on: built at 2.4 it returns 2.4, so a departure from three is a finding rather than a bug. The band table above is that same fitter applied six times to six populations of the same tree.

What a mixed sample does to the interval

The point estimate survives a mixed sample. The interval does not, and it moves in the direction that matters.

Bootstrapping the whole-tree fit gives an interval of [2.80, 2.91] on two hundred junctions. The even-fork band alone, on the same number, gives [2.88, 2.99]. So adding a hundred and fifty uninformative junctions to fifty informative ones does not widen the interval much — it barely narrows it either, which is the leverage argument again — but it does shift it, by about a twentieth, downwards.

That shift is small and it is not nothing. It is the twigs contributing their six hundredths of a per cent, biased, and it says the correct summary is not “a mixed sample is safe” but “a mixed sample is safe in proportion to how much of it is informative”. A tree sampled with ten even forks and a thousand twigs would be dragged more than one sampled with fifty and a hundred and fifty, because the weight is a total rather than a share.

The version of the advice that survives all of this is therefore about a ratio rather than about a count: what matters is the leverage of the informative junctions in the sample against the leverage of the biased ones, and both are computable from the daughter ratios alone before any exponent is fitted. A person with a list of measured junctions can work out, in advance, whether their sample can answer the question — which is the same service the internode counts and specimen counts elsewhere in this collection provide, in a subject where nobody usually asks.

The two costs, and where their sum is leastPumping cost falls as r⁻⁴ and upkeep rises as r². Minimising the sum gives flow proportional to r³ — checked here across flows of 1.2, 2.4, where Q/r³ holds to 0.2%.010200.50011.502vessel radiuscost per unit lengthleast total costpumping falls, upkeep risesthe sum has one minimum
Fig. 6 The derivation the exponent comes from, drawn at two flow exponents. What this essay measures is not whether three is right but whether a given set of junctions could tell three from two — and the answer depends entirely on which junctions, and not much on how many.

The limit of the claim

All of this is arithmetic about a model. It says what a junction could report if the model is right, and what a sample of junctions does to an estimator under a stated measurement error. It does not say Murray’s derivation applies to a tree, and this site’s own essay on fitting the exponent is careful that a fitted value departing from three is the interesting case rather than an error.

Two specific assumptions carry the numbers. The measurement error is two per cent of a radius, stated rather than measured, and every count in the thread scales as its square — a reader who can do one per cent should divide the sample sizes by four and will find the retention rates rise correspondingly. And the daughter ratios are drawn evenly across the band, which no real tree does; a real distribution concentrated at the lopsided end has less leverage in total than the one measured here, though its shares by band are unchanged.

What does not depend on either is the shape: leverage falls as the sixth power of the daughter ratio, retention falls with it, and a sample selected by convenience is selected against information. That is a statement about a class of measurement rather than about trees, and it is why this essay sits in a collection about phyllotaxis at all.

The two costs, and where their sum is leastPumping cost falls as r⁻⁴ and upkeep rises as r². Minimising the sum gives flow proportional to r³ — checked here across flows of 1.5, 2.4, where Q/r³ holds to 0.7%.010200.50011.502vessel radiuscost per unit lengthleast total costpumping falls, upkeep risesthe sum has one minimum
Fig. 7 Where the exponent of three comes from — the trade between the cost of pumping and the cost of maintaining the vessel, which is Murray’s derivation and is the reason three is the hypothesis rather than an arbitrary number. Everything in this essay is about whether a measurement could distinguish it from two, not about whether it is right.

Why an even fork is rare, and what that costs

The arithmetic says one even fork settles the question. It is worth asking how often a tree provides one, because “one is enough” is only useful advice if one is available.

A branching architecture that divides its flow evenly at every junction is a particular kind of tree and most are not. The common case is a trunk that sheds branches — a dominant axis with subordinate ones coming off it — which puts every junction at a small daughter ratio by construction. Species with genuinely dichotomous branching, where a shoot divides into two comparable daughters, are the ones this measurement is cheap on, and they are a minority.

So the practical position is uncomfortable in a specific way. On a dichotomous tree the exponent is settled by a handful of junctions and the measurement is almost free. On an excurrent tree, where the interesting architecture is, the informative junctions are the few near the base and they are the ones a person cannot reach — and the thousands that are reachable carry, between all of them, less weight than one fork at the crown.

That is not a limitation of the estimator and no better estimator removes it. It is the geometry of the question: every exponent predicts nearly the same parent radius at a lopsided junction, so a lopsided junction contains almost no information about the exponent, however carefully it is measured and however many of them there are. The right response is to say so in the sample design rather than to hope that a large n compensates, because it does not — it is the one case where a bigger sample buys almost exactly nothing.

A tree built at an exponent of 2, measured band by bandSix bands of daughter ratio, 300 junctions in each, all from trees built at an exponent of exactly 2 with the same 2% measurement error. The median implied exponent falls from 2.00 at an even fork to 1.50 at a twig, and the share of junctions that have any exponent at all falls from 100% to 58% over the same range. The rule across the top is a single fit over 200 junctions spanning the whole range: 1.96. A mixed sample is safe because least squares already weights by leverage — 42% of it sits in the most symmetric band and 0.06% in the twigs. The sample that is not safe is the one a person can reach.11.5022.5033.50012345daughter ratio of the junctions measuredexponent implied, and the share of junctions that imply one0.9–1100%0.7–0.9100%0.5–0.7100%0.3–0.599%0.15–0.382%0.05–0.1558%all 200 junctions together: 1.96the tree was built at 2bars are the share with any exponent · whiskers are the tenth and ninetieth percentilesmeasurement error 2%
Fig. 8 The same six bands on a tree built at an exponent of two rather than three, which is the control the whole essay needs: the whole-tree fit returns two, the bands fall away from it in the same pattern, and the retention rates collapse at the same place. What the bands report is the geometry of the junctions rather than the exponent they were built at.

What it adds to the branching thread

The thread now has three results and they fit together.

Which junctions could ever say anything: the separation between the two rules collapses as the cube of the daughter ratio, so one even fork settles it and two thousand twigs do not.

What happens when the ones that cannot are measured: a confident wrong answer, because half of them are impossible and the survivors are biased.

And what happens when everything is measured: nothing bad, because least squares is already weighting by the first result — which means the danger is never the mixture and always the selection.

The third is the one that changes what a person should do. The first two invite the reading “be careful with twigs”; together with the third the advice is sharper. Take whatever junctions are available and report the retention rate — or take the few good ones deliberately. What must not be done is to take the ones that are easy to reach and treat the number that comes out as though it came from a tree.

The number to report

If one line comes out of this thread it is: report attempted and retained.

It costs nothing, it needs no model, it is two integers, and it separates a measurement from a selection in a case where nothing else in the output does. A sample that retained every junction has an answer. A sample that retained half has a number produced by which half survived.

That the same diagnostic works on this site’s sequence instrument is a coincidence of shape rather than of subject, and it is the reason the two threads ended in the same phase. A readout with a refusal in it produces two numbers — how many were attempted and how many returned — and the second on its own is never enough.

Shares its objects with

Essays that name at least two of the same things, and that neither author linked.

  • What a quiet plant is worth — both name discrimination, evidence, honest limits, identifiability, measurement, measurement error, sample size, specimen
  • What the pair costs — both name discrimination, honest limits, identifiability, measurement, measurement error, residual, sample size, specimen
  • The test a plant could settle — both name discrimination, evidence, identifiability, measurement, measurement error, sample size, specimen
  • What a refusal does not say — both name discrimination, evidence, honest limits, identifiability, measurement, sample size, specimen
  • A window inside a rung — both name honest limits, identifiability, measurement, specimen
  • Two readings from one stem — both name artefact, identifiability, measurement, specimen

Named objects

A flat tag is an object no other essay names yet.

ArtefactBranchBranching exponentDa Vinci's ruleDiscriminationEvidenceFittingHonest limitsIdentifiabilityMeasurementMeasurement errorMurray's lawResidualSample sizeSpecimen