Stems and cones

Two readings from one stem

Three phase plans in a row have recorded that the two statistics of a divergence sequence want opposite plants — one quiet, one disturbed. Measured on the same stems they do not. The conflict was in the interpretation of a sign, and the window in which both are readable is wide.

Worth reading first: The sequence has a memory · Where the noise gets in · The survey this site cannot do.

There is a worry that has been written into three consecutive phase plans on this site, in almost the same words each time, and it has never been measured.

The earlier statistic reads the correlation between one divergence and the next. It is decisive below about 0.9° of recorded scatter and gone by 1.2°: it wants a quiet plant. The spectrum readout reads the periodicity of the whole thing. It has nothing to autocorrelate at 0.05° of scatter and is clean at 1.1°: it wants a disturbed one. Worse, the first was measured on a stem whose rise was falling and the second needs a rise held fixed, because a band over which the parastichy pair is constant is forty to eighty nodes at the site’s usual rate and forty nodes cannot show a period of eight.

Two statistics, apparently two plants, and a mixture problem that cannot be closed by having two numbers if the two numbers cannot be taken off one specimen.

This essay takes them off one specimen.

Both, on the same sequences

The experiment is the obvious one and its only virtue is that it was actually run. Hold the rise fixed, sweep the disturbance from 0.02 to 1.0 in azimuth degrees per node, five seeded stems at each amplitude, and compute both statistics from the same autocorrelation of the same list of angles.

Both statistics, on the same stems, at a rise of 0.005Five seeded stems at each disturbance, held at a fixed rise. Bars are how many returned the pair the position counter finds; open portions are refusals. The pair comes out from 0.1 to 0.25, and across that whole range the lag-one correlation of the *same* sequences is -0.33, -0.58, -0.59 — decisive, negative and flat. There is no trade between the two: one stem supplies both. Below the window the sequence has locked onto the sampling grid and is a cycle rather than a sample; above it there is no lattice left, at 116° of scatter.012345-1.70-1.30-1-0.824-0.602-0.398-0.2220disturbance amplitude, degrees of azimuth per nodestems out of five returning the counted pair0.020.050.10.150.250.40.61lag-one correlation = 0lag one, on the same sequencesrise 0.005 · 5 stems per point · bars are the pair, line is lag onefilled where the pair agrees with the position counter
Fig. 1 Bars are the pair readout — how many of five stems returned the pair the position counter finds — and the line across them is the lag-one correlation of the very same sequences. The window in which the pair comes out runs from a disturbance of 0.1 to 0.4, and across all of it the lag-one correlation is between −0.50 and −0.65: decisive, negative and flat.

There is no trade. Across the whole range where the pair is readable, the lag-one correlation sits between −0.50 and −0.65 with no trend in it, which is four to six times the sampling band. One stem supplies both, comfortably, at every amplitude where either works.

Where the conflict came from

The sign.

The measurement phase reported a lag-one correlation of +0.54 with no noise at all and called it the sequence’s own memory — a lattice correcting itself, a node placed to one side leaving a gap the next falls into. The previous phase found that explanation is about a lattice, so it ought to hold of one that is not going anywhere, and tested it: at a fixed rise the noiseless rule converges exactly, every divergence identical, and under a disturbance the lag-one correlation is −0.68.

The +0.54 belongs to the rise declining. It is the pattern chasing an equilibrium that is moving under it, and it appears as a threshold in the rate — negative to about fifty-two nodes per rung, positive from fifty-eight.

The memory belongs to the rise, not to the latticeThe lag-one correlation of a noiseless rising stem, against how fast it climbs the ladder. Below about sixty nodes per rung it is negative; above it, 0.54, 0.74, 0.58, 0.58, 0.55 — flat across a fivefold change in rate. The horizontal line is the same rule with the rise held FIXED, where the correlation is -0.68. So the +0.74 the previous phase called the sequence's own memory is the pattern chasing an equilibrium that is moving under it.-0.50000.500100200300nodes per rung of the ladder — how slowly the shoot climbscorrelation between one divergence and the nextthe rise held fixed — -0.68a fast shoot has noneand a slow one saturatesno noise · 254, 323, 438, 553, 922, 1474 nodesthe threshold is in the rate, not the rule
Fig. 2 The correction the previous phase made, which is where the conflict was hiding: the lag-one correlation’s sign is set by how fast the shoot climbs the ladder, not by how quiet the plant is. Everything measured about noise discrimination was measured on a rising stem, so the statistic in that thread is +0.54; everything measured about the comb is measured at a fixed rise, where the same statistic is −0.6.

So the two statistics were never asking for different plants. They were being read on stems that differed in the rise, and the quantity that changed sign was mistaken for a quantity that changed availability. The worry in three phase plans was a worry about a sign convention.

The memory of a divergence sequence, at 0.75° of scatterWith no noise at all the lag-one correlation is 0.54: the rule corrects itself, so a lattice arrives with a memory in it. Matched at the same recorded scatter, placement noise leaves -0.10, jostle noise leaves 0.66, field noise leaves 0.47. The band is ±0.13, which is what an uncorrelated sequence of this length gives.-0.25000.2500.500123456lag, in nodescorrelation between a divergence and the one that many nodes laterno noiseplacement noisejostle noisefield noisesampling band4 runs each · 243 divergences per runmatched at 0.75° of scatter
Fig. 3 The original measurement, on rising stems, where the lag-one correlation is large and positive with no noise and is erased by placement noise. Nothing in that thread’s conclusions changes. What changes is which stem the number describes, and therefore whether it can be read off the same stem as the comb.

The window’s two ends, and they fail differently

The pair is readable from a disturbance of about 0.1 to about 0.4, which is a recorded divergence scatter of roughly 0.4° to 1.0°. Both ends are hard edges and they fail in ways that are not equivalent.

The loud end is safe. At 0.6 and above the lattice is gone: the divergences wander over the whole circle, the recorded scatter is above a hundred degrees, and every run refuses. A botanist looking at such a stem would not call it patterned either, and the instrument’s coherence test is the same judgement made arithmetically.

The quiet end is dangerous, and this is the finding of the essay that was not looked for. At a disturbance of 0.02 the rule locks onto the grid of azimuths it is sampled on. Every divergence comes out at 137.8125° — which is exactly 147 of 384 grid steps — on every seed, to the last digit. What is left to autocorrelate is a short deterministic repeat rather than a sample of a process, its main comb is 0.98, and the pair it names is 8 and 11 where the counter says 8 and 13.

It reports, and it is wrong. That is a worse failure than a refusal and it took a sweep to find, because at the amplitude the earlier essays use it does not happen.

Two combs, at a rise of 0.005The autocorrelation of 760 divergence angles from one stem held at a rise of 0.005. The filled teeth are the lags at multiples of 8; the open teeth are the second comb, at the same spacing offset by 5. Reading the spacing off the first and the offset off the second gives the pair 8 and 13, which is what the position counter reports for the same stem — from angles alone, with no coordinate anywhere in the calculation.-0.50000.500125810131621242629lag, in internodescorrelation between a divergence and the one that many internodes later816245132129spacing 8 · offset 5pair 8/13 — counter says 8/13sampling bandone stem · 760 divergences · disturbance 0.25generated from a stated rule, not drawn to look right
Fig. 4 The spectrum at the disturbance the working essays use, where the main comb is 0.64 and the second 0.50. At a fiftieth of that disturbance the same figure has a main comb of 0.98 — a correlation that high is not a strongly patterned plant, it is a sequence with almost nothing random in it, and that is what the cap on the readout now refuses.

The cap, and what it does not fix

A correlation above 0.95 means the sequence is very nearly deterministic, which no real measurement of a real plant will be. So the readout refuses above it.

That removes the failures at the quietest amplitude — all five runs — and three of five at the next. One wrong answer survives at 0.05 and one at 0.15, and both are left in the record rather than tuned out, because moving the cap far enough to catch them also refuses runs that are right.

The rest of the requirement is a condition on the specimen and not something the angles can enforce: the plant has to be disturbed. That is an odd thing to have to ask for and it is worth saying why it is not a defect. The comb is a measurement of how a disturbance propagates. A pattern with no disturbance in it has nothing propagating, so there is nothing to measure — in the same way that a perfectly still fluid reports nothing about its viscosity.

A lattice or a wreck, with nothing in betweenEvery run in the sweep, at both kinds of noise. The line at 6° is where a run stops being counted as a lattice, and nothing lands near it: the intact runs reach 1.97° and the destroyed ones start at 8.06°, a factor of 4.1 away.00.50011.5001234567amplitude, by stepscatter, log₁₀ degreesintactno latticeplacementfield48 runs · both kindsan empty factor of 4.1 at the cut
Fig. 5 What a ruler reads against what was put in, from the phase that separated the kinds of noise. The recorded scatter is the quantity a botanist can check before attempting either statistic, and the window this essay measures is stated in it: below about 0.4° there is not enough sequence to read, above about 1.2° there is not enough lattice.

The grid artefact, and what it says about the model

The quiet-end failure is a property of the simulation rather than of plants, and it is worth being clear about which.

The rule samples 384 azimuths and takes the minimum, so its output is quantised to 0.9375°. At large disturbances the noise moves the chosen sample around and the quantisation is invisible. At small ones it does not, and the run settles onto a single grid point and stays there. The “sequence” is then a repeat with a period set by the seeded stretch.

A real apex has no grid. So the specific failure — locking to 137.8125° on every seed — will not happen to a plant. What will happen to a plant is the general version of it: a measurement quantised more coarsely than the pattern’s own variation produces a sequence with structure that is the quantisation’s rather than the plant’s. A protractor read to the nearest degree, on a plant whose divergences vary by half a degree, is exactly this failure with a different grid in it.

So the cap earns its place twice. It guards the model against its own sample grid, and it guards a real measurement against a coarse protractor — and in the second case it is the only guard there is, because the dilution law that prices reading error assumes the error is fine and random, which a rounding is not.

Two stems at 0.75° of scatter, one angle at a timeThe divergence of each node, with the slow climb of the ladder removed. Both stems scatter by about 25.97° and a botanist would report them identically. The jostled one wanders in runs — its lag-one correlation is 0.70 — and the one with placement noise alternates about the mean at 0.21, because a displacement of one node enters two consecutive divergences with opposite signs.jostle noise — correlation 0.70placement noise — correlation 0.21243 nodes each, both at 25.97° of scattercorrelations 0.70 and 0.21
Fig. 6 Two stems with the same recorded scatter and different processes behind them, angle by angle. The quantity a ruler reports is the same for both; what differs is the order, and the same is true of a quantised sequence against a genuine one — the summary statistics agree and the structure does not.

What is now available from one stem

Setting out the readings, in the order of what they cost:

  • the recorded scatter, from a ruler, which says whether the other two are worth attempting;
  • the lag-one correlation, from sixty internodes, which on a fixed rise is the rule’s self-correction at about −0.6 and on a rising stem is the pattern chasing the rise at about +0.55 — and whose sign therefore reports which of those the specimen is;
  • the parastichy pair, from two hundred and fifty internodes at a quarter of a degree.

All three from one list of angles up one stem. That is the thing three phase plans doubted, and it is settled — for a stem held at a fixed rise, which no plant is.

The remaining question is whether it survives a stem that climbs, and it does, on a condition that turns out to be about the shoot rather than about the statistics: the window read has to fit inside one rung of the ladder. That is the next essay, and it is where the mixture problem finally gets an answer with a shape.

What the sequence sees that the scatter cannotEach point is an ensemble at one amplitude, placed at the scatter it produces. A stem at three quarters of a degree of scatter has a lag-one correlation near zero if its noise arrived after the primordium was placed, and near 0.7 if it arrived before — and no measurement of the scatter can tell those apart. The separation closes above about a degree, because what the other two kinds preserve is the correlation of a lattice.-0.20000.2000.4000.6000.5000.75011.251.50divergence scatter, in degrees — the one quantity a plant offerscorrelation between one divergence and the nextplacement noisejostle noisefield noise4 runs per point · band ±0.13every point is a lattice
Fig. 7 The lag-one statistic against the recorded scatter for three kinds of noise, from the measurement phase. Its useful range and the comb’s overlap over most of a decade, which is the overlap this essay’s sweep measures directly rather than inferring from two separate experiments.

The window is wide, and that is the surprise

It is worth dwelling on how much room there is, because the worry that prompted the experiment implied there might be none.

The pair is readable from a disturbance of 0.1 to 0.4, which is a factor of four in amplitude. The recorded scatter across that range runs from about 0.4° to about 1.0°, which is most of the range in which the site’s noise work has found a lattice to exist at all. And the lag-one correlation is not merely present across it — it is between −0.50 and −0.65, which is four to six sampling bands, at every amplitude.

Compare that with the two statistics’ individual ranges as the earlier phases measured them: the lag-one statistic decisive below 0.9° and gone by 1.2°, the comb absent at 0.05° and clean at 1.1°. Those two ranges overlap over most of a decade, and the overlap is where both work. The phase plans that recorded the conflict were reading the two endpoints of one range against the two endpoints of the other and concluding that the ranges were disjoint, when what they had was two overlapping intervals described from opposite ends.

That is a small error of reading and it survived three restatements, which is worth recording as a caution about how a worry propagates. Each restatement was a paraphrase of the one before, none of them was a measurement, and the paraphrase got sharper each time — “two statistics with opposite requirements” in one phase became “they want opposite specimens” in the next and “one a disturbed plant on a steady rung, one a quiet plant on a rising stem” in the third. The last of those is the most specific and the most wrong.

The tolerance falls from 2.9° to 0.8° up the ladderEach dot is one configuration: a stem carried down to the stated rise, swept over the whole amplitude range, with the largest divergence scatter any surviving run showed. It ends on 8/13 at the coarse end and 8/13 at the fine one. The open marks are the width of the band of divergence angles that gives each pair at all — a quantity from a different calculation entirely, moving the same way.00.2000.4002.402.502.602.702.80final rise, as −log₁₀degrees, log₁₀8/138/13scatter survivedthe pair's band5 runs per amplitude · placement noise2.89° down to 0.83°
Fig. 8 Why the window has a loud end at all: the amount of scatter a lattice tolerates tightens as the pattern gets finer, which the measurement phase established and which sets where the coherence test cuts. The instrument’s window is bounded above by the pattern’s own survival rather than by anything about the statistic.

What the two statistics are about, restated

Since the sign is what caused the trouble, it is worth writing down plainly what each number means on which kind of stem.

On a stem at a fixed rise, the pattern has an equilibrium and stays at it. A node displaced to one side leaves a gap the next node falls into, so consecutive divergences alternate: the lag-one correlation is negative, around −0.6, and it is the rule correcting itself. The comb is present, and its spacing is the parastichy number.

On a stem whose rise is falling, the equilibrium is moving. The pattern is always slightly behind where the current rise would put it, so consecutive divergences are displaced in the same direction: the lag-one correlation is positive, around +0.55, and it is the pattern chasing. The comb is absent, because the parastichy number is changing.

So the sign of the lag-one correlation reports which regime the specimen is in, and it does so before either instrument is attempted. That is a use for it nobody had, and it costs nothing — it is the same number, read for its sign rather than its magnitude.

What the divergence does while the pattern climbsThe stem produces a sequence rather than a constant. Over the second half of the run it stays within 5.6° of 137.51°, and the vertical marks are where the counted pair changed — the wander is largest around them.136138140142100200300400nodedivergence from the node before (°)137.51°415 nodes at 87 per rungspread 5.63° over the second half
Fig. 9 A growing stem’s divergence angle against node number, which is the picture the positive lag-one correlation is a summary of: the value drifts as the rise falls, so consecutive readings are correlated in the same direction. On a fixed rise this curve is flat and the correlation flips sign.

The methodological note

This essay’s result is that a worry recorded three times was not a real constraint, and the reason it survived three phases is worth naming: it was never measured because it was never cheap to measure the two together. Each statistic lived in its own library, each library had its own runs, and the comparison was made by reading two numbers out of two docstrings.

The fix was to compute both from one call. bothStatistics takes a list of angles and returns the lag-one correlation, the comb scores and the pair, all off a single autocorrelation, precisely so the two cannot be quoted from different runs by accident. That is not a clever piece of engineering; it is a shape that makes a particular mistake impossible, and the mistake had been made three times.

What is left of the worry

Nothing about the fixed rise, and one thing about the rising one.

At a fixed rise the two statistics come off one stem across a wide window of disturbance, and the conflict recorded three times was a misread sign. That much is settled.

What remains is that a real plant’s rise is not fixed, so the comb has to be read over a window and the lag-one statistic over the whole climb — and those are different lengths of the same stem with different requirements. The next essay measures what that costs, and the answer is a condition on the shoot’s rate rather than on either statistic.

Why the sweep is five stems and not fifty

Every point in the disturbance sweep is five seeded stems, which is small, and it is worth saying why it is not larger.

The quantity being measured at each amplitude is not a mean with an error bar. It is a count — how many of five returned the counted pair — and what the sweep establishes is where that count is zero, where it is five, and roughly where it changes. Five is enough to distinguish those three states and it is not enough to locate the boundary precisely, which the essay does not claim to do: the window is quoted as “about 0.1 to about 0.4” rather than with edges.

Fifty stems per point would locate the edges and would cost fifty times as much, on a sweep whose points are already grown stems of nine hundred nodes each. The edges are not worth that, because a real plant’s disturbance is not a number this model can predict — what matters is that the window is wide, and five stems at seven amplitudes says so.

Shares its objects with

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

  • A counter that sees no positions — both name autocorrelation, discretisation, divergence angle, ensemble, equilibrium, measurement, noise, rise, sampling, summary statistic
  • What the protractor has to be — both name autocorrelation, discretisation, divergence angle, ensemble, equilibrium, measurement, noise, sampling, specimen, summary statistic
  • A shoot too fast to remember — both name autocorrelation, divergence angle, ensemble, equilibrium, measurement, noise, rise, sampling, self correction
  • The order carries the count — both name autocorrelation, discretisation, divergence angle, ensemble, equilibrium, measurement, noise, rise, summary statistic
  • What one angle says about the next — both name autocorrelation, divergence angle, ensemble, identifiability, measurement, noise, self correction, summary statistic, tolerance
  • A harmonic is a step taken twice — both name artefact, autocorrelation, discretisation, divergence angle, equilibrium, measurement, self correction, summary statistic

Named objects

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

ArtefactAutocorrelationDiscretisationDivergence angleEnsembleEquilibriumIdentifiabilityMeasurementNoiseRiseSamplingSelf correctionSpecimenSummary statisticTolerance