Foreign Instruments · archival source criticism meets uncertainty propagation

The Uncertainty the Record Cannot Identify

Cavendish's 1798 table supports several numerical uncertainty statements, including his 1/14 judgment, but the archive does not select one unambiguous component-wise distribution and dependence model for a modern physical-result budget.

Verdict: UNDERDETERMINEDThis is the result, not a loading state.

Calibration · native field first

Make the gauge block come back

Before the JCGM procedure touches an archive, it has to reproduce the core outputs of its published end-gauge example. Press once. The page evaluates the model, checks analytic sensitivities against finite differences, then sends 1,260,000 seeded trials through all nine input distributions.

Published GUM
838 ± 93 nm
99 percent expanded result
Published Monte Carlo
[745, 932] nm
shortest 99 percent interval
Trials in JCGM 101
1.26 million
published adaptive run

Ready. The foreign target remains locked.

Published end-gauge values beside live browser calculations
checkpublishedthis browser
GUM estimate838 nmpending
GUM standard uncertainty32 nmpending
component variancereported 1002 nm²pending
effective degrees of freedom16.7, truncated to 16pending
99 percent expanded uncertainty93 nmpending
relative standard uncertainty6.4 × 10-7pending from reconstructed source-precision components
reported length(50.000 838 ± 0.000 093) mmpending
second-order standard uncertainty34 nmpending
analytic versus finite difference sensitivitiesmust agreepending
pendingGUM frameworkfirst-order model
Monte Carlo estimate838 nmpending
Monte Carlo standard uncertainty36 nmpending
GUM shortest interval[745, 931] nmpending
Monte Carlo shortest interval[745, 932] nmpending
pendingpending trialsseed pending

A precision wrinkle: the four one-decimal contributions shown in Annex H.1 are 25, 9.7, 2.9, and 16.6 nm, whose squares sum to 1003.06 nm². Reconstructing from the source equations gives 1002.25 nm², 16.7 effective degrees of freedom, and 6.3 × 10-7 relative standard uncertainty. The source prints 1002 nm² and 6.4 × 10-7. The core 32 nm result is unchanged, so the page labels rather than hides the remaining display mismatch. The 93 nm expanded value uses a stated conservative ceiling to a whole nanometre.

Open the nine-input home model
Nine input distributions used in the home benchmark
inputassigned distributionpinned information
reference length lSscaled Student t, 18 df50,000,623 nm, scale 25 nm
mean difference Dscaled Student t, 24 df215 nm, scale 6 nm
comparator random d1scaled Student t, 5 df0 nm, scale 4 nm
comparator systematic d2scaled Student t, 8 df0 nm, scale 7 nm
alphaSrectangular[9.5, 13.5] × 10-6 per degree C
theta0normal-0.1 ± 0.2 degree C
cyclic Deltaarcsine[-0.5, 0.5] degree C
delta alphacurvilinear trapezoida=-1.0, b=1.0, d=0.1 × 10-6 per degree C
delta thetacurvilinear trapezoida=-0.050, b=0.050, d=0.025 degree C

l = lS + D + d1 + d2 - lS(delta alpha × (theta0 + Delta) + alphaS × delta theta)

Layer one · ten-second archival audit

Close every envelope, or stop

The ordinary spread of repeated values is one envelope. Cavendish also judged a difference of 1/14 of the whole very unlikely. That is numerical information, but it does not uniquely say whether 0.38 is a 95 percent expanded uncertainty, a bound, or a component separate from variability already in the table. This page therefore applies a deliberately conservative adapter rule: stop unless the component-wise distribution and dependence treatment are explicit. JCGM 100 permits scientific judgment and does not require this exact rule.

Later-wire repeatabilityLOCKEDWaiting for calibration
Two-decimal printingLOCKEDWaiting for calibration
Possible common air-current effectLOCKEDWaiting for calibration
Thermometer calibration and gradientsLOCKEDWaiting for calibration
Clock calibration and serial timing dependenceLOCKEDWaiting for calibration
Shared mass and geometry calibrationLOCKEDWaiting for calibration

Rival baselines · different measurands made explicit

The smaller question still has an answer

The primary historical rival follows Cavendish's own judgment: 5.48 ± 1/14 of the whole, approximately 5.48 ± 0.38 or [5.10, 5.86]. Merkatas and colleagues interpreted 0.38 as an approximate 95 percent expanded uncertainty. That interpretation is informative but not forced by Cavendish's prose. The other intervals below answer narrower process or sampling questions and do not include an unknown common apparatus effect.

Historical physical-result judgment [5.10, 5.86] Cavendish's 1/14 judgment, interpreted by Merkatas et al. as approximate 95 percent expanded uncertainty

center 5.48
expanded amount 0.38
component decomposition and coverage meaning not uniquely specified

BESIDE
Experiment-cluster sensitivity locked 95 percent t interval for the later-wire row mean, clustered by 14 printed experiments

cluster SE locked
equal experiment-weighted mean locked
within-experiment correlation is tested as a sensitivity, not asserted as fact

Narrower sampling diagnostics

iid Student t interval locked
row mean locked; sample SD locked; SEM locked
median locked; 10 percent trimmed mean locked
sampling plus optional independent rounding u locked
Huber estimate locked; row-wise Huber BCa locked
locked resamples, seed locked

The baseline runs when the archive audit is unlocked.

These are distinct outputs, not interchangeable estimates of the same thing. The finite printed mean is exact for those 23 numbers. Iid and cluster intervals concern a repeatable process. A partially identified statement is mu_physical = mu_process + b_common, with the sampling part numerical and the common correction left symbolic. Cavendish's 1/14 judgment concerns the physical result. The conservative completeness gate asks whether a modern component-wise physical-result model is uniquely specified.

Layer two · source against adapter

Look at what the reduction discarded

Public-domain facsimile crop of Cavendish's page 520 conclusion table, showing experimental columns and 29 printed density values
Conclusion table, original page 520. Public-domain scan. The reduced computation uses the density column and row identity. Every other visible numerical column is dropped, not silently absorbed.
rowcohortdensity
Adapter audit · loss first

Four mappings enter the collision

Source-to-formal-object mappings and their losses
sourceformal objectwhat is lost
Printed final densityRepeated output DiUpstream covariance and the input responsible for each residual.
Rows 7 to 29 after the wire changePrimary wire_2 cohortConditioning, creep, ageing, and gaps inside the cohort.
Two-decimal printingOptional independent rectangular rounding inputsTruncation, recomputation, and editorial intervention cannot be recovered.
Air-current warningPossible common Type B effect or model discrepancyThe 1/14 judgment survives, but its component-wise distribution, coverage meaning, and overlap with observed scatter do not.

Dropped on purpose

  • the first six first-wire densities from the primary estimate
  • the unresolved 4.88 to 5.88 hypothesis
  • today's accepted density as a tuning target
  • qualitative weather notes as invented numbers
  • weight motion, arm motion, corrected arm, vibration time, and corrected vibration time

Conventions supplied here

  • the physical result, finite printed table, and later-wire process are kept as different measurands
  • blank means missing, never zero
  • target coverage is 95 percent and intervals are shortest unless labeled otherwise
  • rounding is one independent rectangular input per printed value only when its switch is on
  • absent covariance is not renamed zero covariance
  • the stop rule is this page's conservative adapter choice, not a command quoted from JCGM

Disable each mapping

Each button clones the current adapter data, disables one mapping, and reruns the same auditTarget, propagation, and endpoint-validation path. An ablation survives only if the verdict stays the same and the mean, trimmed, and Huber centers each move by less than 0.01 density unit. The rounding ablation explicitly starts from the enabled rounding branch. If a mapping removal makes the calculation undefined, the row says so.

not run
not run
not run
not run

Secondary reproduction · sealed until reader-state freeze

Reproduce Lauginie's narrower calculation

Lauginie's 2007 reanalysis helped motivate the later-wire cohort, so it is not independent validation of that choice. After freezing the reader's current options, this panel deliberately recomputes the immutable later-wire comparison: mean 5.48 and twice the standard error about 1.50 percent. It is a reproduction of a secondary calculation, not evidence that the adapter was untuned.

Fault bench · one catch, one miss

Break what the instrument can see

The nonlinear branch is explicitly synthetic. It doctors the model input consumed by the production audit to use a positive denominator-like quantity with 50 percent relative uncertainty. The same unmodified completeness, propagation, and endpoint-validation path then checks a first-order interval against 200,000 Monte Carlo draws. The frozen target record is untouched.

Agreement between two propagation methods is not truth. A common bias can move every result together while leaving scatter, interval width, and method agreement intact.

The check

What is fixed, free, and still missing

Result branchUNDERDETERMINED under this page's conservative rule because the source does not select a component-wise distribution and dependence model for its named possible common air effect. Three further apparatus categories are analyst-inferred cautions and are not assigned invented non-zero sensitivities.
Home acceptanceCore checks: GUM 838 nm, u=32 nm, conservatively reported U99=93 nm; Monte Carlo 838 nm and u=36 nm. The fixed-seed shortest interval is diagnostic because this page does not implement the full JCGM adaptive endpoint-stability test.
Transcription29 density values read from the page 520 crop and checked against the independently published HistData CSV. Default primary cohort is the final 23.
Free choicesmeasurand, cohort policy, printed 4.88 or hypothesized 5.88, optional independent round-to-nearest inputs, 95 percent target coverage, shortest interval, Huber c=1.345, and bootstrap seed 17981345.
UncertaintiesCavendish's 1/14 judgment is retained. What remains unclear is its distribution, coverage meaning, dependence, and overlap with the variability in the table. Thermometer, timing, and geometry categories are analyst-inferred cautions, not effects claimed as source-named.
Secondary evidenceLauginie's rounded mean and twice-SEM result are reproduced, not presented as independent validation. Manuscript images are not shipped; provisional manuscript values are not treated as print corrections.
Source acquisition and file identities
Fetched source identities, retrieval results, and reuse decisions
sourceretrievedresponse and hashlicence decision
JCGM 100:20082026-09-10HTTP 200, 1,888,806 bytes
41bbf068fbc0d7986c98691b2d1af6680cb3044f6a1a89b3560933ed9ef9626c
JCGM 2008, all rights reserved; linked, not redistributed
JCGM 101:20082026-09-10HTTP 200, 1,489,478 bytes
6d8548af875df112dfc5cf14eb974f5544341f4a31cbbdfc19fc5ed155d8fa20
JCGM 2008, all rights reserved; linked, not redistributed
JCGM 100 Amendment 1:20262026-09-10HTTP 200, 620,150 bytes
ce50a6ba8b5a338d81082494e81988f137c89d951c8dcb647a1fc856b7b5f809
JCGM 2026, all rights reserved; linked, not redistributed
Cavendish printed paper2026-09-10HTTP 200, 14,601,468 bytes
cc1054b89a9494e8977daef89b873519b990648c4134d56ac84882a502d5cdf9
Public Domain Mark 1.0
HistData CSV2026-09-10HTTP 200, 516 bytes
df2d7b063f3bae5fd518e4cc313280d37e14f833dc26a43148ce397bae797c3d
GPL-2 or GPL-3
Royal Society manuscript object2026-09-10HTTP 200, 18,255,161 bytes
1654026949ba2cbce5fe055cbcbed997b7d52da23cbc0ea17df3b3d4ba4c6ef0
Research access; permission required for web reproduction, so no image is shipped

We searched the general web, Crossref Works, and OpenAlex Works on 10 September 2026 for an application of JCGM 100 and JCGM 101 to Cavendish's 1798 published observation tables, and did not find a source that constructs a source-audited Type A and Type B budget and validates its linear propagation against Monte Carlo. This reports the search, not universal priority.

The check that stands behind this page recomputes the page's numerical results, imports the browser engine, compares the independent frozen transcription, checks the manifest ledger, reads every static displayed figure back out of this HTML, validates the share surface, and carries mutation controls that must turn those assertions red.

Sources

The documents, not substitutes for them

  1. Henry Cavendish, Experiments to Determine the Density of the Earth, Philosophical Transactions 88 (1798), 469 to 526. Conclusion table at original page 520.
  2. JCGM, JCGM 100:2008 and JCGM 101:2008. The 50 mm end-gauge example is Annex H.1 and Supplement 1 section 9.5.
  3. Pierre Lauginie, Weighing the Earth, weighing the Worlds (2007), reproduced as a secondary rounded calculation.
  4. Merkatas, Toman, Possolo, and Schlamminger, Shades of dark uncertainty and consensus value for the Newtonian constant of gravitation, Metrologia 56 (2019) 054001. They interpret Cavendish's 0.38 as an approximate 95 percent expanded uncertainty.
  5. JCGM, JCGM GUM-6:2020. Its hierarchical-model guidance motivates the experiment-cluster sensitivity without asserting that within-experiment correlation is proven.
  6. Stephen M. Stigler, Do Robust Estimators Work with Real Data?, Annals of Statistics 5(6), 1977, 1055 to 1098. Huber and trimmed centers are retained as secondary robustness diagnostics, not privileged as the strongest rival.
  7. Michael Friendly, HistData Cavendish documentation, machine-readable comparison path for the density column.