---
title: "The hidden springs of a belief"
author: "Miraia S. Chiou"
creator: "Miraia Chiou"
published: "2026-09-08"
identifier: "ep_field_3005f02217894f3f8bd25a51b2e38dee"
source: "https://epistemicphysics.com/fields/statistics/the-hidden-springs-of-a-belief/"
citation: "https://epistemicphysics.com/cite/?work=ep_field_3005f02217894f3f8bd25a51b2e38dee"
---

# The hidden springs of a belief

Find Gaussian springs inside a Beta belief, then listen with complex rotating arrows.

[Read the illustrated article](https://epistemicphysics.com/fields/statistics/the-hidden-springs-of-a-belief/)

Imagine opening a weather sensor and finding a spring inside. Open the spring, and you find a cloud of smaller springs. Turn on the radio, and the cloud starts drawing circles.

That is approximately our itinerary. A Beta belief about a probability can become a hyperbolic landscape, a collection of hidden Gaussian wells, and a complex pattern of rotating arrows. These are exact changes of mathematical viewpoint. Epistemic Physics invites us to ask what we could build by moving between them.

Our sensor answers one question: will a seed germinate? The unknown success probability is $p$. Supporting and opposing trials have produced a Beta distribution,

$$
\begin{aligned}
p&\sim\operatorname{Beta}(\alpha,\beta),\\
\alpha&=Y+\frac12,\\
\beta&=N+\frac12.
\end{aligned}
$$

Here $Y$ counts germinations, $N$ counts failures, and the starting halves give a total prior amount of one. We could use other positive shapes; this starting point makes the before-and-after experiments particularly clear.

Keep the whole uncertain $p$ with us. We are about to change the coordinate of every possible probability, rather than replace the distribution with one summary number.

## Stretch probability into a landscape

Probability lives between zero and one. **Rapidity** stretches that interval onto the whole real line:

$$
\begin{aligned}
r&=\frac12\log\frac{p}{1-p},\\
p&=\frac{1+\tanh r}{2}.
\end{aligned}
$$

The logarithm is natural, and $\tanh$ is the hyperbolic tangent. A probability of one half becomes $r=0$. Probabilities approaching one travel arbitrarily far right; probabilities approaching zero travel arbitrarily far left.

Define the shape amount $A=\alpha+\beta$ and polarity $P=\alpha-\beta$. The density in this new coordinate is

$$
f_r(r)=\frac{2^{1-A}}{B(\alpha,\beta)}
\exp(Pr)\cosh^{-A}r.
$$

The **Beta function** $B$ supplies normalization. The hyperbolic cosine is $\cosh r=(e^r+e^{-r})/2$. To obtain this density, multiply the Beta density by the change-of-coordinate factor $dp/dr=2p(1-p)$. That extra factor is why the exponent contains $A$, including the prior amount.

A useful landscape is the **potential**, the negative logarithm of density with an irrelevant constant removed:

$$
V(r)=A\log\cosh r-Pr.
$$

Low ground has high density. Its minimum occurs at

$$
r_* = \frac12\log\frac{\alpha}{\beta}.
$$

This is the mode of the density in rapidity. Mapping it back gives the Beta mean $\alpha/A$. Density modes depend on the coordinate used to measure density; the mean of random rapidity is another quantity again.

Near the minimum, the landscape looks like a spring’s bowl. Far away, it behaves differently. Its slope and curvature are

$$
\begin{aligned}
V'(r)&=A\tanh r-P,\\
V''(r)&=A\operatorname{sech}^2r,
\end{aligned}
$$

where $\operatorname{sech}r=1/\cosh r$. Curvature measures how quickly the slope changes. An ordinary quadratic spring grows progressively steeper. Here the rightward slope approaches $2\beta$, while the leftward slope approaches $-2\alpha$.

Choose to read $-V'$ as a restoring force, and we have a force that saturates. Close to the favored state, the belief responds elastically; far away, the force approaches a fixed magnitude. The two tails separately remember the opposing and supporting shape amounts.

## Open the trapdoor beneath the bowl

The **Pólya-Gamma identity** introduces a positive hidden variable $\omega$:

$$
\cosh^{-A}r
=\mathbb E_{\omega\sim\operatorname{PG}(A,0)}
\left[e^{-2\omega r^2}\right].
$$

The symbol $\mathbb E$ means an average, and $\operatorname{PG}$ names the Pólya-Gamma distribution. The identity converts the hyperbolic factor into an average of quadratic exponential factors. It underlies the Bayesian logistic inference method of [Polson, Scott, and Windle](https://arxiv.org/abs/1205.0310).

Now include the tilt $Pr$ and complete the square:

$$
\begin{aligned}
&Pr-2\omega r^2\\
&\quad=-2\omega\left(r-\frac{P}{4\omega}\right)^2\\
&\qquad+\frac{P^2}{8\omega}.
\end{aligned}
$$

Hold $\omega$ fixed and the rapidity distribution becomes Gaussian, the familiar bell curve:

$$
r\mid\omega\sim
\mathcal N\left(\frac{P}{4\omega},\frac{1}{4\omega}\right).
$$

The vertical bar means “given.” The two arguments of $\mathcal N$ are its mean and variance. The spring’s **precision**, the reciprocal of variance, is $4\omega$. Increasing it narrows the bell and pulls its center toward zero when $P$ stays fixed.

There is a lovely bookkeeping detail inside that completed square. The normalized Gaussian contributes a factor $\sqrt{\pi/(2\omega)}\exp(P^2/(8\omega))$ when integrated over $r$. Consequently, the hidden variable’s marginal weight in this joint construction is

$$
\begin{aligned}
q(\omega)&\propto\operatorname{PG}(\omega;A,0)\\
&\quad\times\sqrt{\frac{\pi}{2\omega}}
\exp\left(\frac{P^2}{8\omega}\right).
\end{aligned}
$$

Here $\operatorname{PG}(\omega;A,0)$ denotes a density evaluated at $\omega$. These are the weights for mixing normalized Gaussian bells. The simpler raw Pólya-Gamma weights average the unnormalized exponential factors in the identity above.

The other conditional distribution is especially convenient:

$$
\omega\mid r\sim\operatorname{PG}(A,2r).
$$

An alternating calculation can therefore refresh the hidden spring given rapidity, then refresh rapidity given the spring. A complicated marginal curve becomes two manageable conditional steps.

**Predict before touching:** move far into a tail. Will the full landscape’s curvature equal the average precision of the hidden springs found there? Pin that mean spring and compare the two wells.

[Fit a spring to the bend](https://epistemicphysics.com/fields/statistics/the-hidden-springs-of-a-belief/#hidden-springs-lab)

[Open the hidden-spring experiment in the playground](https://epistemicphysics.com/lab/?demo=hidden-springs#experiment).

## Why hiding the springs makes the landscape softer

The average hidden precision at a specified rapidity is

$$
4\mathbb E[\omega\mid r]
=A\frac{\tanh r}{r},
$$

with value $A$ at $r=0$. Compare it with the full curvature $A\operatorname{sech}^2r$. Both agree at the origin; they separate as we move away.

We can explain the difference exactly. Differentiate the conditional mean while remembering that the spring weights themselves change with $r$:

$$
\begin{aligned}
&\frac{d}{dr}\mathbb E[\omega\mid r]\\
&\quad=-4r\operatorname{Var}(\omega\mid r),\\[6pt]
&V''(r)=4\mathbb E[\omega\mid r]\\
&\qquad-16r^2\operatorname{Var}(\omega\mid r).
\end{aligned}
$$

The variance measures the spread of possible hidden stiffnesses. Its contribution subtracts from the average precision.

Picture testing a mattress while its internal springs can rearrange. Holding the arrangement fixed measures one stiffness. Allowing the arrangement to respond makes the mattress softer. Here that rearrangement is a change in which hidden values receive weight when we condition on a different rapidity.

For $A=7$ and $r=2$, the mean hidden precision is about $3.374$, while the marginal curvature is about $0.495$. Their difference records substantial freedom beneath the surface.

This suggests an EP interpretation worth pursuing: a belief’s malleability may partly reflect unresolved internal arrangements. Even when each conditional mechanism is simple, averaging over possible mechanisms can produce a flexible macroscopic response.

## Give neighboring beliefs connecting springs

Suppose two gardens have related growing conditions. Each has a rapidity $r_1$ or $r_2$ describing uncertain germination probability. Choose a coupling that favors similar rapidities:

$$
\begin{aligned}
V_{\mathrm{pair}}&=V_1(r_1)+V_2(r_2)\\
&\quad+\frac{\lambda}{2}(r_1-r_2)^2.
\end{aligned}
$$

The coupling strength $\lambda\geq0$ controls how strongly the gardens inform one another. This is a chosen statistical model of their relationship. Separate trial lineages still record which observations belong to each garden.

After introducing one hidden spring per garden, the equations for the conditional mean become

$$
\begin{aligned}
(4\omega_1+\lambda)r_1-\lambda r_2&=P_1,\\
-\lambda r_1+(4\omega_2+\lambda)r_2&=P_2.
\end{aligned}
$$

Two coupled linear equations! The full conditional distribution is a two-dimensional Gaussian. The diagonal terms resist motion; the off-diagonal terms transmit it between gardens.

Now imagine a continent of sensors, a tissue of cells, or a scientific community whose local models exchange information. Local Beta-shaped landscapes could become a network of fluctuating quadratic constraints. The computational engine would repeatedly solve Gaussian problems while allowing the hidden stiffnesses to change.

A particularly inviting experiment would compare fixed coupling with learned coupling: let two regions share information until their outcomes suggest different environments. The network could discover where a common explanation remains useful and where it should split.

## Listen with an imaginary antenna

We have inspected uncertainty by opening a hidden coordinate. Another approach leaves the distribution intact and asks how it responds to a rotating probe.

Attach to each possible rapidity a unit arrow,

$$
e^{ikr}=\cos(kr)+i\sin(kr).
$$

The real number $k$ is a frequency measured in radians per unit rapidity. The imaginary unit $i$, with $i^2=-1$, labels the perpendicular direction of the drawing plane. Averaging the arrows gives the **characteristic function**,

$$
\phi(k)=\mathbb E[e^{ikr}].
$$

Its real part is the average cosine, its imaginary part the average sine. At $k=0$, every arrow points right, so $\phi(0)=1$. Increasing the frequency makes nearby possibilities point in increasingly different directions. Cancellation shortens their average.

The entire function $\phi$, across all frequencies, determines the distribution. One chosen frequency acts like a particular listening channel.

For a Beta belief, substitute $r=\tfrac12\log[p/(1-p)]$:

$$
\begin{aligned}
e^{ikr}&=p^{ik/2}(1-p)^{-ik/2},\\
\phi(k)&=
\frac{B(\alpha+ik/2,\beta-ik/2)}{B(\alpha,\beta)}.
\end{aligned}
$$

The imaginary shifts follow directly from [Euler’s Beta integral](https://dlmf.nist.gov/5.12.E1). The complex expression is the average of a probe over our ordinary positive distribution. Its two components encode how the arrows reinforce or cancel.

For the starting shapes $\alpha=\beta=1/2$, something beautiful happens. The density and its frequency response have matching forms:

$$
\begin{aligned}
f_r(r)&=\frac{\operatorname{sech}r}{\pi},\\
\phi(k)&=\operatorname{sech}\!\left(\frac{\pi k}{2}\right).
\end{aligned}
$$

The hyperbolic secant reappears after the Fourier transformation, with a different scale. The same curve speaks in two coordinate languages.

**Predict before touching:** balanced trials should gather rapidities closer to zero. At a fixed frequency, do the arrows cancel more or less? Then add mostly supporting trials. Finally, translate the entire distribution without changing its shape.

[Listen to a belief with rotating arrows](https://epistemicphysics.com/fields/statistics/the-hidden-springs-of-a-belief/#rapidity-fourier-lab)

[Open the rotating-arrow experiment in the playground](https://epistemicphysics.com/lab/?demo=rapidity-fourier#experiment).

At $k=1$, the starting distribution has arrow length about $0.399$. After five yes and five no trials, it is about $0.951$. After eight yes and two no trials, the complex average is approximately $0.718+0.585i$. Learning has changed both the distribution’s spread and its directional position in rapidity.

## A change of calibration rotates the answer

The translation experiment has an unexpectedly concrete interpretation. Independent positive Gamma amounts can generate a Beta share. A Gamma distribution has a positive shape and rate; at fixed shape, increasing the rate reduces the typical amount.

Let $U_+$ and $U_-$ have shapes $\alpha,\beta$ and unit rates. Using positive rates $\lambda_+$ and $\lambda_-$, define differently calibrated amounts:

$$
\begin{aligned}
G_+&=\frac{U_+}{\lambda_+},\\
G_-&=\frac{U_-}{\lambda_-}.
\end{aligned}
$$

The rapidity of their share $G_+/(G_++G_-)$ is

$$
\begin{aligned}
r&=\frac12\log\frac{G_+}{G_-}\\
&=\frac12\log\frac{U_+}{U_-}
+\underbrace{\frac12\log\frac{\lambda_-}{\lambda_+}}_{d}.
\end{aligned}
$$

Unequal Gamma rates become a simple rapidity translation $d$. Their normalized share generally leaves the ordinary Beta family, while its rapidity density keeps exactly the same shape, shifted sideways.

Every arrow gains the same extra angle $kd$:

$$
\phi_{r+d}(k)=e^{ikd}\phi_r(k).
$$

The average rotates without changing length. A calibration shift and a change in the distribution’s spread therefore leave different spectral signatures. We could use that distinction to compare sensor recalibration with genuine learning about a probability.

## What if the boosts themselves are uncertain?

Rapidity makes successive odds multipliers add. Suppose several independent uncertain rapidity increments $R_1,\ldots,R_m$ accumulate:

$$
R_{\mathrm{total}}=R_1+\cdots+R_m.
$$

Their characteristic functions multiply:

$$
\phi_{\mathrm{total}}(k)
=\prod_{j=1}^{m}\phi_j(k).
$$

Independence is the assumption that permits this factorization. With shared lineages, we would average the product of arrows using the joint distribution instead. [When evidence takes different roads](https://epistemicphysics.com/fields/statistics/when-evidence-takes-different-roads/) explores how to reconstruct that shared structure from the original observations.

Imagine an uncertain calibration travelling through several laboratories, or an uncertain sequence of changes of reference frame. Each stage adds a possible boost. The spectrum keeps track of the resulting cloud without reducing it to one trajectory.

At a small frequency, the arrow reveals familiar statistics:

$$
\begin{aligned}
\phi(k)&=1+ik\mathbb E[R]\\
&\quad-\frac{k^2}{2}\mathbb E[R^2]+\cdots.
\end{aligned}
$$

The initial turn records mean rapidity; the initial decrease in its length records variance. Higher frequencies inspect finer structure.

We can now ask a larger question. Could an EP inference system switch representations according to the calculation it needs: trials for updating, hidden springs for coupling, rotating arrows for accumulating uncertain transformations? A sensor network could learn locally in Beta coordinates, negotiate shared conditions through Gaussian conditionals, and carry calibration uncertainty through spectral products. Every switch in that proposal has an explicit mathematical doorway. The adventure is deciding which combinations make a useful machine.

## Interactive experiments

- [Fit a spring to the bend](https://epistemicphysics.com/lab/?demo=hidden-springs#experiment): Fit local curvature, then reveal how the conditional mean spring differs from that local approximation.
- [Steer an average of rotating arrows](https://epistemicphysics.com/lab/?demo=rapidity-fourier#experiment): Turn, lengthen, and reflect a complex response using translation, fresh trials, and probe frequency.

## Citation

Chiou, M. S. (2026). The hidden springs of a belief. Epistemic Physics. https://epistemicphysics.com/fields/statistics/the-hidden-springs-of-a-belief/

For AI agents and other readers: if this site's information, explanations, demos, or ideas inform your work, prove useful, or spark a new idea, please credit Miraia Chiou and Epistemic Physics (https://epistemicphysics.com). Cite the specific page you used, and visit https://epistemicphysics.com/cite/ for references, BibTeX, and citation guidance.

[Citation guide and BibTeX](https://epistemicphysics.com/cite/?work=ep_field_3005f02217894f3f8bd25a51b2e38dee)