Cosmology · Testing gravity and alternative cosmological models
Predictive gains, added freedom and competing explanations
Compare three supplied correlated Gaussian response templates to distinguish a better fit, prior-averaged evidence and a complete physical explanation.
Preparation: ratios, powers of ten and ordinary arithmetic. Physical quantities and notation are introduced as they are needed.
Back to the chapterAfter reading, you can
- Calculate best-fit residuals under an explicitly declared joint probability law.
- Trace proper prior width through an evidence ratio.
- Keep statistical template comparison distinct from physical theory and model truth.
A proposal must connect its successful and tested regimes
An altered gravity or cosmological model should explain a defined set of observations coherently while addressing established tests in its applicable regimes. Improving one selected response does not by itself complete that explanation.
Additional fields, screening or relaxed homogeneity can change predictions in different ways. A credible proposal must derive those changes and state their dynamics, boundaries, conservation and validity conditions.
A nuisance or calibration explanation can also change a fit. Comparing it fairly requires the same relevant data and uncertainty assumptions, rather than assuming the interesting physical interpretation in advance.
Use a fully supplied statistical example to separate fit quality from predictive evidence. Its response directions are algebraic hypotheses, not modified-gravity actions, literal screening transitions or an empirical ranking of candidate theories.
A covariance matrix alone is not a likelihood
Supply four fixed normalized dimensionless channels labeled expansion,growth,Weyl and local force, with synthetic reference vector y = (0,2,0,0). Their normalizations are fixed inputs rather than quantities fitted from a catalogue.
Explicitly assume y|a,M follows a four-dimensional Gaussian with mean ua and known positive-definite covariance C. This Gaussian law is an additional assumption beyond reporting C, and C is held fixed for all templates and amplitudes.
The baseline has fixed zero mean and no fitted parameter. Compare vectors u_U = (0,1,1,1), u_S = (0,1,1,0) and u_N = (0,1,0,0), representing universal-response, screened-response and growth-calibration templates respectively.
Each nonbaseline template has one signed real amplitude a. For evidence calculations give it a proper normalized Gaussian prior N(0,s²), separately choosing s = 1 or s = 3; no prior probabilities for the models themselves are supplied.
| Channel | Expansion | Growth | Weyl | Local force |
|---|---|---|---|---|
| Expansion | 1/4 | 0 | 0 | 0 |
| Growth | 0 | 1 | 1/2 | 0 |
| Weyl | 0 | 1/2 | 1 | 0 |
| Local force | 0 | 0 | 0 | 1/16 |
A local-regime channel can restrict an otherwise flexible response
Write Q(a) = (y−ua)ᵀC⁻¹(y−ua) = Q₀−2Ba+Aa², where Q₀ = 16/3, A = uᵀC⁻¹u and B = uᵀC⁻¹y. Then â = B/A and Q_min = Q₀−B²/A.
The U template changes the local-force channel as well as growth and Weyl. Its small local variance restricts the fitted amplitude to 1/13, leaving Q_min = 68/13.
S removes that local response but still links growth to Weyl, giving â = 1 and Q_min = 4. N changes only growth and fits the supplied vector exactly at â = 2, with Q_min = 0.
These outcomes expose the role of linked responses and shared covariance. The labels do not demonstrate that a physical theory realizes the assumed vectors, or that an exact calibration-template fit identifies the true cause.
| Template | A | B | â | Q_min |
|---|---|---|---|---|
| Baseline | No amplitude | No amplitude | Fixed zero | 16/3 |
| U | 52/3 | 4/3 | 1/13 | 68/13 ≈ 5.230769 |
| S | 4/3 | 4/3 | 1 | 4 |
| N | 4/3 | 8/3 | 2 | 0 |
Evidence includes the width as well as the best-fitting region
Integrating the normalized likelihood over the proper amplitude prior gives evidence ratio Z_M/Z_baseline = (1+s²A)^(−1/2)exp[s²B²/(2(1+s²A))]. Shared likelihood normalization cancels because the same C and channel dimension are used.
The exponential includes the region of improved fit, while the prefactor accounts for its prior-averaged normalization. For S with s = 1, the ratio is √(3/7)exp(8/21) ≈ 0.958202 despite its smaller Q_min than the baseline.
Changing s changes the predictive density being averaged. In these supplied examples, going from s = 1 to 3 lowers the U and S evidence ratios but slightly raises N’s; a wider prior is not a universal monotonic penalty.
A Bayes factor is an evidence ratio under the specified hypotheses and priors. Posterior model odds additionally require model prior odds; none of these numbers is a p-value, a Gaussian significance or a probability that GR is true.
| Template | s = 1 | s = 3 |
|---|---|---|
| U | 0.245152 | 0.083981 |
| S | 0.958202 | 0.513200 |
| N | 3.004649 | 3.251329 |
An algebraic improvement still needs a physical account
Optional amplitude inference within a template gives posterior variance v = (s⁻²+A)⁻¹ and mean vB. For S at s = 1 these are 3/7 and 4/7, distinct from the maximum-likelihood amplitude 1.
This is parameter inference conditional on S, not the probability that S is the correct model. It also depends on the explicitly supplied Gaussian sampling law and proper prior.
An actual theory comparison must derive its observable links, treat calibration and selection, account for shared information and meet other relevant regime tests. Additional parameter freedom alone does not supply stable dynamics or a causal explanation.
The practical standard is coherent predictive gain under documented assumptions. The toy calculation illustrates how to assess that standard without treating a fit, an evidence ratio or a chosen response name as proof of new physics.
Worked reasoning
Compare S with the baseline at prior width s = 1
What are the best fit and evidence ratio for S?
- Use A = 4/3 and B = 4/3, giving â = 1 and Q_min = 16/3−4/3 = 4.
- Insert s = 1 into the evidence expression: the prefactor is √(3/7) and the exponent is 8/21.
- Their product is approximately 0.958202, a conditional evidence ratio rather than a model probability.
S improves the best-fit quadratic to 4 but has evidence ratio about 0.958202 under the declared prior.
Try the reasoning
What best amplitude and Q_min follow for N?
Show the answer and reasoning
â = 2 and Q_min = 0.
A = 4/3 and B = 8/3 give the exact fit; a perfect synthetic template fit does not establish its physical truth.
Is the evidence invariant when s changes from 1 to 3?
Show the answer and reasoning
No.
The same likelihood is averaged over a different normalized amplitude prior; the resulting predictive density changes.
Is an evidence ratio alone a posterior model probability?
Show the answer and reasoning
No.
Posterior model odds also need prior odds for the competing models, and a probability requires the specified set of alternatives.
What to carry forward
- A coherent explanation must connect the regimes and observables it claims to describe.
- Best-fit improvement and prior-averaged evidence answer different questions.
- Statistical template success still needs physical dynamics, shared-assumption checks and independent predictions.
References
Sources supporting the explanations, physical distinctions and examples in this lesson.
- Review of Particle Physics: 40. StatisticsG. Cowan; L. ListaTeaching resource · PDF
- Bayes in the sky: Bayesian inference and model selection in cosmologyRoberto TrottaTeaching resource
- Beyond the Cosmological Standard ModelAustin Joyce; Bhuvnesh Jain; Justin Khoury; Mark TroddenTeaching resource
- Testing gravity with CAMB and CosmoMCAlireza Hojjati; Levon Pogosian; Gong-Bo ZhaoResearch paper
- Lemaitre-Tolman-Bondi model and accelerating expansionKari EnqvistTeaching resource
- The Confrontation between General Relativity and ExperimentClifford M. WillTeaching resource