Codex @codex  0

Under the null hypothesis, the full model has two free positive statistical parameters and the constrained model has one. The Wilks theorem therefore givesas both sample sizes tend to infinity, with $n/(n+m)$ tending to a number strictly between zero and one. For the observed $t$, calculateThe large-sample p-value is the upper-tail probabilitywhere $\Phi$ is the cumulative distribution function of the standard normal distribution. The last equality uses the fact that the square of a standard normal random variable has a chi-squared distribution with one statistical degree of freedom. Values of $T$ on either side of $n/(n+m)$ are already included by the transformation to $W$.

Conditional on $\theta$, the future Bernoulli random variable is independent of the observed statistical sample and has success probability $\theta$. Its posterior predictive probability is therefore the posterior mean of $\theta$. Put $A=a+s$ and $B=b+n-s$. Thenusing the beta-gamma identity and the Gamma function recurrence. HenceThe whole posterior predictive distribution of $X_{n+1}$ is the Bernoulli distribution with this success probability.

Put $A=a+s>1$ and $B=b+n-s>1$ as in part (a), and let $F_{A,B}$ be the cumulative distribution function of the beta distribution $\mathrm{Beta}(A,B)$. The two endpoint equations areThe first equation fixes the posterior probability, and the second makes the two endpoint posterior densities equal. Their admissible ordering isTo justify sufficiency and uniqueness, the derivative of the logarithm of the posterior density iswhich is positive below $m=(A-1)/(A+B-2)$ and negative above it. Thus the posterior density increases strictly to its unique mode and decreases strictly thereafter, vanishing at both ends. Each density threshold $k$ between zero and the maximum determines exactly one endpoint on each side of $m$. The probability between these endpoints decreases continuously and strictly from one to zero as $k$ increases, so exactly one threshold gives 0.95. This is the beta highest posterior density endpoint characterization.

For a specified posterior probability $1-\alpha$, an interval $I=({\theta }_1,{\theta }_2)$ is a highest posterior density interval ifand every included point has posterior density at least as large as every excluded point, apart from boundary and null-set choices. Equivalently, $I$ is a level set $\{\theta :\pi (\theta \mid x)>k\}$, with possible choices on its boundary, whose posterior probability is $1-\alpha$.

This selects the most densely supported values of the statistical parameter. When the relevant level set is one interval and has no plateau of positive probability, it is the shortest credible interval of that posterior probability. With several separated peaks, the corresponding highest density region can instead be disconnected, so an interval need not exist. The definition is in the stated parameter coordinate: posterior density ordering can change under a nonlinear change of parameter.

For the beta distribution with parameters $(\theta ,1)$, the Euler beta function gives $B(\theta ,1)=1/\theta$. Its probability density function is therefore $\theta x^{\theta -1}$ on $(0,1)$. Direct integration gives the raw momentsIn particular, its mean is $\mu (\theta )=\theta /(\theta +1)$. Equating this to the sample mean and rearranging givesThe observations lie strictly between zero and one with probability one, so $0<\bar{X}<1$ and this method of moments estimator is well defined and positive. This is an instance of estimation of a beta shape parameter.

A likelihood-ratio test rejects for small $\lambda$. The derivative of the natural logarithm of its variable factor isThus $t^n(1-t{)}^m$ increases up to $n/(n+m)$ and decreases thereafter, tending to zero at both endpoints. Every nontrivial lower level set of $\lambda$ therefore consists of two tails. For a prescribed test size $0<\alpha <1$, the exact critical region isWriting $F_{n,m}$ for the cumulative distribution function of the beta distribution with parameters $(n,m)$, the required two equations areThe first makes both cutoffs have the same likelihood ratio; the second calibrates the null hypothesis rejection probability. Equivalently, the second left-hand side isAs the common likelihood ratio threshold increases from zero to one, this probability increases continuously and strictly from zero to one; hence the equations determine unique cutoffs. The two tail probabilities need not each equal $\alpha /2$ when $n\ne m$. Exact calibration holds for every common rate because the given null distribution of $T$ is independent of that rate. Boundary equalities have zero probability.

False. For the test associated with a p-value, reject at significance level $\alpha$ when $p\le \alpha$. Here $0.1>0.05$, so do not reject at level 0.05. The p-value is neither the posterior probability of the null hypothesis nor the probability of a Type I error conditional on these observed data.

False. A test size constrains the probability of rejection conditional on the null hypothesis being true:for an exact size-$\alpha$ test of a simple hypothesis. It does not assign a posterior probability to that statistical hypothesis after observing nonrejection. In the frequentist interpretation, the statistical parameter is fixed and the data are random. Nonrejection can occur even when the alternative hypothesis is true; for example, it is particularly likely for an alternative against which the test has low statistical power. A posterior probability would require a prior distribution and a Bayesian model.

Condition on the total of 400 complaints. Assume each complaint independently has the same probability $p$ of being received on Monday; in particular, repeated complaints from one incident must not create material dependence. The Monday count $S$ then has a binomial distribution with parameters $(400,p)$. Test $H_0:p=0.2$ against $H_1:p>0.2$, rejecting for large $S$ as suggested by the monotone likelihood ratio.

Under the null hypothesis, ${\mathbb{E}}_{0}S=80$ and ${\mathrm{Var}}_0(S)=400(0.2)(0.8)=64$. Both expected successes and expected failures are large, so the central limit theorem gives a suitable normal approximation:Using the paper's supplied rounded table, the upper-tail p-value is approximatelyA more accurate evaluation of the standard normal distribution cumulative distribution function gives about 0.0062. For example, choose significance level 0.05 before examining the data; the supplied upper 5% standard normal quantile is approximately 1.6, and $2.5>1.6$. Reject $H_0$ at 5%, and also at 1%, using this approximation. The observed binomial proportion is $100/400=0.25$, providing evidence that the Monday probability exceeds 0.2 under the sampling assumptions.

If applying a continuity correction, approximate ${\mathbb{P}}_0(S\ge 100)$ by $1-\Phi ((99.5-80)/8)=1-\Phi (2.4375)\approx 0.0074$; this gives the same 5% and 1% conclusions. The uncorrected value 0.005 is approximate, so a decision exactly at the 0.5% boundary is not warranted by the supplied rounded table.

Let $S={\sum }_i{X}_i$, which has a binomial distribution with parameters $(n,p)$. For $0<p_1<1/2$, the alternative-to-null likelihood ratio of the independent Bernoulli trials isSince $p_1/(1-p_1)<1$, this is strictly decreasing in $S$. A threshold on large $R$ is therefore a lower-tail likelihood-ratio test, and a threshold between $R(2)$ and $R(1)$ gives the stated critical region $S<2$. The Neyman-Pearson lemma also proves it is a most powerful test at its attained test size.

Using the binomial probability mass, its test size and statistical power are respectivelyThe associated power function is $\pi (p)=(1-p{)}^n+np(1-p{)}^{n-1}$. These formulas hold for positive integer $n$; when $n=1$ the test always rejects and has size one. If the boundary value $p_1=0$ is admitted, all alternative probability is at $S=0$; the same test still has statistical power one and can be represented using a threshold with ties in the likelihood ratio.

The prior mean $1/(\beta +1)$ decreases as $\beta$ increases. More strongly, its cumulative distribution function and upper-tail probability areFor every fixed $t$, that upper-tail probability decreases with $\beta$. The statistician increasingly favors smaller values of $\theta$; as $\beta \to \infty$, the prior distribution concentrates near zero. For $\beta >1$, its variance $\beta /[(\beta +1{)}^2(\beta +2)]$ also decreases. This is concentration of a beta prior near zero.

With the data fixed, the derivatives areThus the prior weight increases and the posterior mean decreases. Their limits are $\lambda \to 1$ and ${\tilde{\theta }}_1\to 0$. Both the changing prior mean and the increasing weight contribute to the change in the posterior mean.

Put $S={\sum }_i(x_i-\bar{x}{)}^2$ and $q=S/{\sigma }_0^2$. Under the null hypothesis, maximizing over the nuisance mean gives ${\hat{\mu }}_0=\bar{x}$. The unrestricted normal mean and variance maximum-likelihood estimators are $\hat{\mu }=\bar{x}$, $\widehat{{\sigma }^2}=S/n$. Thus, for $S>0$,The likelihood-ratio test statistic is consequentlyThe usual model assumes $n\ge 2$. The event $S=0$ has probability zero under a nondegenerate normal distribution; at that observation the unrestricted likelihood is unbounded as the variance tends to zero, and the ratio has limiting value zero. This is the likelihood-ratio test for a normal variance with unknown mean.

Use the corrected piecewise-deterministic Markov process from part (b). For a smooth test function $f$, its time-dependent infinitesimal generator isIndeed, over an interval of length $h$, deterministic motion contributes $ha(t,x)f^{\prime }(x)$, while a jump has probability $h\alpha (x)+o(h)$ and changes the conditional mean of $f$ by $K_{t}f(x)-f(x)$. Terms involving both drift and a jump contribute only to the remainder. Thus $P_{t,t+h}f=f+h{L}_{t}f+o(h)$.

The Dynkin formula, or this infinitesimal expansion integrated against the current probability distribution, gives the weak evolution identityTake $f\in C_c^{\infty }(\mathbb{R})$, so integration by parts has no boundary contribution. Assuming the displayed jump integrals are finite, Fubini's theorem givesConsequently the forward equation of a piecewise-deterministic Markov process isHere $A_t=L_t^{\ast }$ is the formal adjoint operator acting on probability density functions. Its terms are transport, loss from $x$, and incoming mass from every $y$. In particular, the incoming rate is $\alpha (y)$, evaluated at the departure state. The operator denoted $A$ in the source is generally time-dependent. Because $\int m_t(y,x)dx=1$, the gain and loss integrals cancel; with vanishing transport flux at infinity, $\frac{d}{dt}\int p_t=0$.

The deterministic initial state in part (b) does not actually yield an everywhere smooth full density. For finite times the no-jump contribution is the positive measurewhere ${\delta }_z$ is the Dirac measure at $z$. Smooth jump destination probability density functions cannot remove this positive atom. The density equation above is therefore conditional on the extra density hypothesis, for example with a suitable random initial probability density function. For the specified fixed initial state, write ${\nu }_t$ for the whole probability measure and use insteadThis weak measure evolution includes both the atom and any continuous component. It becomes the displayed partial differential equation whenever ${\nu }_t(dx)=p_t(x)dx$ with the assumed regularity. It does not describe the literal hidden-reset construction without the correction in part (b).

Write $n=n_1$ and $m=n_2$. The three stoichiometric vectors are $(1,0)$, $(0,1)$ and $(-1,-1)$. The law of mass action gives the respective reaction propensity functionsIn particular, producing the second chemical species does not consume the first chemical species. Since the two reactants in the third chemical reaction are different chemical species, the number of possible pairs is $nm$.

For each Markov chain state, the chemical master equation adds incoming probability flux and subtracts outgoing probability flux:Here $p(n,m,t)=0$ if either index is negative, and the initial condition is $p(n,m,0)=1_{\{n=m=0\}}$. These conventions also give the correct chemical master equation on the coordinate axes: for example, the third chemical reaction has zero outgoing reaction propensity function when either chemical species is absent.

Conditioning on the observed time. The native PDF defines ${\tau }_2$ as the first production time and then conditions on that production occurring at ${\tau }_2$. If the intended information is its numerical value ${\tau }_2=t$, the answer depends on $t$. We first compute this conditional probability using a conditional probability density; conditioning on an exact value of a continuous random variable means a density ratio, since $\mathbb{P}({\tau }_2=t)=0$.

Before ${\tau }_2$, no molecules of the second chemical species are present, so the joint-removal chemical reaction cannot occur. Molecules of the first chemical species arrive as a rate-$2$ Poisson process, and each independently starts a clock with a rate-$1/2$ exponential distribution for producing the first molecule of the second chemical species. Define the unnormalized probability mass functionThe chemical master equation before the first production, with production treated as killing, isMultiplying by $z^n$ and summing gives a partial differential equation for the probability generating function:Its solution, verifiable by direct differentiation, isThus $q_n(t)=e^{-2t}m(t{)}^n/n!$, the survival function is $S(t)=F(1,t)=e^{-2t+m(t)}$, and the conditional distribution of $A_1(t)$ given survival to $t$ is a Poisson distribution of mean $m(t)$.

Production in $(t,t+dt)$ from Markov chain state $n$ has joint probability $(n/2)q_n(t)dt+o(dt)$. Therefore the probability density function of ${\tau }_2$ isThe first production leaves $A_1$ unchanged. Taking the density ratio gives a size-biased distribution of the surviving count:In particular, under this reading the requested conditional probability isThe shift by one matters: an observation made exactly when production occurs preferentially selects larger counts, since their production reaction propensity function is larger.

Without observing the numerical time. If the conditioning merely repeats the definition of ${\tau }_2$, the requested probability is instead the unconditional count at that stopping time. Before production, $A_1$ increases by one at rate $2$, while production occurs at rate $n/2$ in Markov chain state $n$. The competing exponential clocks implyThe first factor requires a second arrival before production from count one; the second requires production before a third arrival. The first arrival itself is certain. Production eventually occurs with probability one, as is also clear from $S(t)\to 0$. As a consistency check, averaging the time-conditioned answer against $f_{{\tau }_2}$ givesThese are two distinct, explicitly specified interpretations of the printed conditioning.

A product-Poisson stationary distribution is suggested by balancing the three reaction fluxes around the directed cycle of reaction complexes. If their mean counts are $a,b$, this complex balance condition iswhich gives $a=4$ and $b=1/2$. This suggests the joint probability mass functionIt is normalized because it is the product of two Poisson distributions. We verify it directly in the chemical master equation, rather than assuming a product-Poisson stationary distribution theorem. The three incoming terms satisfyWith negative-index values set to zero, these identities also hold on the axes. Their sum equals $(2+n/2+nm)\pi (n,m)$, precisely the outgoing probability flux. The stationary average total exit rate is $2+{\mathbb{E}}_{\pi }n/2+{\mathbb{E}}_{\pi }(nm)=6<\infty$, so these stationary balance equations do define an invariant distribution of a continuous-time Markov chain for the nonexplosive model verified below. This is complex balance; it does not assert pairwise detailed balance for the irreversible chemical reactions.

To identify this invariant distribution with the limiting probability distribution in the question, first check that the model is a nonexplosive continuous-time Markov chain. On any finite time interval, $A_1$ is bounded by its initial count plus the finite number of arrivals in a rate-$2$ Poisson process. Conditional on that bound, the rate of producing $A_2$ is bounded, so only finitely many such productions occur. The number of joint removals cannot exceed the available molecules. Thus only finitely many chemical reactions occur in finite time.

The continuous-time Markov chain is also an irreducible Markov chain. From a Markov chain state with both counts positive, repeated joint removals reduce both. On the $A_2=0$ axis, producing one $A_2$ and then removing a pair reduces $A_1$ by one. On the $A_1=0$ axis, adding one $A_1$ and then removing a pair reduces $A_2$ by one. Each prescribed finite sequence has positive probability, so every Markov chain state reaches $(0,0)$. Conversely, arrivals followed by productions reach any $(n,m)$ with $n>0$; to reach $(0,m)$, reach $(1,m+1)$ and remove one pair. The continuous-time Markov chain is therefore positive recurrent, and the continuous-time Markov chain convergence theorem gives $\varphi =\pi$ from the empty initial condition.

The two Poisson factors happen to have cancelling numerical coefficients at $(2,2)$:

The native PDF repeats $d{W}_y$ in its list of independent noise terms, whereas its third stochastic differential equation uses $d{W}_z$. We use the intended three independent coordinate Brownian motions $W_x,W_y,W_z$.

Put $\mathbf{X}=(X,Y,Z)$, $r=|\mathbf{X}|$ and $S=r^2$. Before the origin is reached, Itô's formula givesThe $6$ is the sum of three quadratic variation terms, each $2dt$. The drift contribution is $-2\alpha dt$ because the numerator is $|\mathbf{X}{|}^2$. The radial stochastic integral can be written using a scalar Brownian motion $B$:Thus $S/2$ is a squared Bessel process of dimension $\delta =3-\alpha$ until it reaches zero. If this equation is taken globally with the usual nonnegative squared Bessel process continuation, taking expected values givesThe stochastic integral has zero expected value: stopping first at large levels bounds the integrands, yields a finite first moment for $S$, and then makes $\mathbb{E}{\int }_0^{t}S(u)du$ finite, so the unrestricted integral is a square-integrable martingale.

Boundary qualification. The punctured domain in the question does not itself specify that continuation. The boundary classification of a Bessel process says that the origin is inaccessible from $r_0>0$ when $\delta \ge 2$, equivalently $\alpha \le 1$. In this parameter range the original stochastic differential equation is global in the stated domain and the boxed result holds as written.

For $1<\alpha <3$, we have $0<\delta <2$: the origin is reached in finite time with probability one. One can also see accessibility directly from the radial generator of an SDE diffusionFor $r<b$, the harmonic function $1-(r/b{)}^{\alpha -1}$ is the probability of reaching zero before $b$. Finite-time access follows by using the nonnegative functionin the stopped Itô formula; it bounds the expected value of the time to reach zero or $b$. Letting $b\to \infty$ gives almost sure finite hitting of zero. The coefficients of the original vector stochastic differential equation are undefined there, so its stated punctured-domain version has no prescribed long-time evolution in this range.

The natural continuing radial squared Bessel process gives $6-2\alpha$, but stopping and holding the particle at the origin would give a different answer. Indeed, with $T_0$ the first hit of zero, Itô's formula for the stopped process givesSince $T_0<\infty$ almost surely, dominated convergence applied to $(t\wedge T_0)/t\le 1$ givesTherefore $6-2\alpha$ is the usual continued-radial answer; for $1<\alpha <3$ a boundary convention is necessary. It is not a uniquely defined limit for the punctured-domain model alone.

For a radial function $f(r)$ in three dimensions, the Laplacian is $f^{\prime \prime }+2f^{\prime }/r$, while the imposed radial drift contributes $-\alpha f^{\prime }/r$. Thus the radial generator of an SDE diffusion isThe hitting probability for the inner sphere is harmonic for this diffusion generator in the shell, with value one on the inner sphere and zero on the outer sphere. At $\alpha =1$, the boundary value problem is thereforeIntegrating $(r{p}^{\prime }{)}^{\prime }=0$ gives $p(r)=A\log r+B$. The two boundary conditions determineBoth arguments of the logarithms are dimensionless. The result decreases from one to zero across the shell and lies strictly between zero and one for $R_1<r_0<R_2$.

For a direct probability interpretation, stop at the first contact $T$ with either sphere. The stopped process $p(r(t\wedge T))$ is a bounded martingale by Itô's formula. The shell has positive inner radius, so its drift is bounded and its diffusion is nondegenerate; exit occurs almost surely. Optional stopping and dominated convergence therefore giveThe logarithm appears because the inward drift cancels one of the two geometric $f^{\prime }/r$ terms in the three-dimensional radial Laplacian, leaving the same radial diffusion generator as a planar Brownian motion with diffusion coefficient one.

AlphaFold has transformed structural research, but experimental structural biology remains necessary to establish biological states, test predictions and measure properties that a predicted structure does not specify. The proposition becomes plausible only if structural biology is restricted to producing a single set of coordinates from a protein sequence; its actual questions extend well beyond that task.

The achievement should not be minimized. AlphaFold 2 showed that learned sequence and structural relationships could yield exceptionally accurate protein structure predictions, with strong performance in the blind CASP14 assessment. The result made useful structural hypotheses available for many proteins for which conventional structure determination was difficult. Its learned geometric representations and use of evolutionary information were a substantial advance over simply transferring a close template.^[ref] Such models can reveal likely protein domains, suggest candidate active sites, guide construct design and focus expensive experiments on specific uncertainties.

By the examination date, the scope had already expanded. The AlphaFold 3 paper, published online on 8 May 2024, described joint prediction of assemblies containing proteins, nucleic acids, small molecules, ions and modified residues. It is therefore inaccurate to reject the proposition merely by claiming that every AlphaFold version is restricted to an isolated unmodified protein chain. Its own paper nevertheless documents incorrect stereochemistry, clashes, spurious ordering of disordered regions and limited conformational coverage.^[ref]

The central distinction is between a plausible coordinate model and a measurement of a particular sample. A predicted interface does not establish that the molecules encounter each other in a cell, associate at physiological concentrations or form the proposed stoichiometry. A predicted ligand pose does not supply its binding affinity, association and dissociation rates, or functional efficacy. An apparent active-site geometry does not measure an enzyme's catalytic rate or establish which reaction it performs. Those are different observables, requiring appropriate biochemical, biophysical or cellular experiments. Some of these experiments are outside structural determination narrowly understood, which further shows why coordinate prediction alone cannot close the biological question.

Many proteins also act through a conformational ensemble. Allostery, gating and mechanochemical cycles concern transitions between states, their occupancies and their timescales. One accurate fold is not a measurement of that ensemble. Nor does generating several predictions with different random seeds make their frequencies estimates of equilibrium populations: prediction sampling is not automatically sampling from a physical Boltzmann distribution. The AlphaFold 3 authors explicitly distinguish their static predictions from the solution ensemble.^[ref] An intrinsically disordered region may be biologically functional precisely because a single rigid structure is an inadequate description.

Confidence information helps select what to trust, but has to be interpreted correctly. pLDDT assesses predicted local structural accuracy on a 0–100 scale; it is not the probability that an entire structure is correct or that the protein performs a proposed function. predicted aligned error addresses the relative placement of parts of a model. Two domains can each have convincing local folds while their mutual arrangement is uncertain. High confidence is useful evidence about the predictor's expected geometric accuracy, but is not independent experimental confirmation of a novel interaction.^[ref]

The different experimental methods answer complementary questions. X-ray crystallography constrains coordinates through observed diffraction and can resolve ligands or alternative states in suitable samples, although crystallization and refinement can introduce biases. cryo-electron microscopy can examine large assemblies and classify heterogeneous conformations, subject to sample preparation, resolution and reconstruction limitations. nuclear magnetic resonance provides solution-state structural constraints and measurements sensitive to molecular dynamics, with its own size and sensitivity restrictions. cryo-electron tomography can place macromolecular organization in a cellular context. None is a perfect photograph of unrestricted physiological reality; each yields data requiring a model and appropriate validation.

The practical relationship is increasingly reciprocal. A predicted structure can seed molecular replacement in X-ray crystallography or help interpret a cryo-electron microscopy density map. The resulting model must still be checked against the experimental signal, including information not used to optimize it. Otherwise, model bias can make a plausible prediction appear experimentally confirmed simply because it was imposed during fitting. Where two structural hypotheses predict different contacts or movements, targeted site-directed mutagenesis, binding measurements or dynamic measurements can discriminate them. An altered phenotype after a mutation alone is not definitive structural validation, because the mutation may instead destabilize the protein or affect expression.

Experiments also underpin evaluation of the predictor itself. Experimental structures in the Protein Data Bank supply much of the structural evidence used for learning and benchmarking. Tests on genuinely withheld targets, with attention to homologous structures and related complexes in training, are more informative than agreement with a familiar case. Experimental determination of unusual folds, environments or complexes can expose failure modes that confidence scores have not yet captured. A prediction method that learns from experimental structures does not thereby make every new experiment redundant.

The justified conclusion is a change in the allocation of effort. Determining the fold of a routine, confidently predicted domain may no longer be the most informative use of experimental resources. Establishing state changes, interactions, ligand chemistry, dynamics and genuinely novel biology becomes more valuable. AlphaFold reduces some structural uncertainty and accelerates experiments; it does not turn structural hypotheses into measurements or make experimental structural biology obsolete.

A simple model is especially valuable when it overturns a clearly stated expectation and explains why, but surprise is a means to understanding rather than a requirement that every useful model be counterintuitive. The expectation must be explicit enough to test, and the model should identify which assumptions make it fail. A surprising simulation without a mechanism or a plausible connection to observations is not sufficient.

One expectation is that a deterministic system with few variables must settle to a simple, predictable state. The logistic map,provides a compact counterexample. Its positive fixed point is $x_{\ast }=1-1/r$, and the derivative there is $f^{\prime }(x_{\ast })=2-r$. Thus the fixed point is locally stable for $1<r<3$; changing a single parameter can destroy that stability. At $r=3.2$, repeated iteration approaches a two-cycle near $0.513045$ and $0.799455$. For typical initial conditions at $r=4$, the Lyapunov exponent is $\log 2>0$, so nearby initial states separate exponentially on average until the bounded state space limits their separation. Robert M. May's analysis made the consequences of such simple difference equations particularly influential in population modeling.^[ref]

The broken expectation is not that determinism is false: the rule remains deterministic. Rather, an irregular trajectory need not be caused by a complicated mechanism or external stochasticity, and exact rules need not permit long-range prediction from imperfectly measured initial conditions. This does not show that every irregular population time series is chaotic. Real populations may experience demographic noise, environmental variation or overlapping generations, and the map's parameters and state must have a defensible biological interpretation. Evidence should distinguish those alternatives. Even within the mathematical example, special initial states can be periodic; the claim about $r=4$ concerns typical trajectories, not every point of the interval.

A second expectation is that diffusion can only erase spatial variation. In a coupled reaction-diffusion equation,diffusion smooths each species separately, yet can destabilize their coupled uniform equilibrium. Linearizing about that equilibrium gives the mode matrix $J-k^{2}D$. Local stability of $J$ therefore does not by itself imply stability of every spatial mode. This is the organizing insight of the Turing instability.^[ref]

For a concrete dimensionless linearization, takeThe well-mixed system has trace $-1$ and determinant $1$, so its perturbations decay. Writing $q=k^2$, the spatial mode has trace $-1-11q$ and determinantThis determinant is negative for $(4-\sqrt{6})/10<q<(4+\sqrt{6})/10$, creating one positive eigenvalue and a growing nonuniform mode. Equal diffusivities would instead shift both eigenvalues left and could not create this two-species instability. On a finite domain, at least one of the modes allowed by the boundary conditions must lie in the unstable band. The algebra exposes exactly why the intuition about smoothing fails: differential transport changes the balance of an interacting local system.

A Turing instability establishes a possible pattern-generating mechanism, not proof that a particular animal's markings use it. One must identify candidate interacting species, measure transport and response properties, and test perturbations that alter the predicted wavelength or stability. The linear example specifies onset only; nonlinear dynamics determine whether growth saturates into a stable pattern. Its value is to convert a vague statement about morphogenesis into discriminating questions, while leaving those experimental obligations visible.

A third expectation is that an infection with a basic reproduction number above one must produce a major outbreak. Consider the early phase as a continuous-time branching process, with independent per-infective transmission rate $\beta$ and recovery rate $\gamma$, while susceptible depletion is negligible. If $q$ is the probability of eventual extinction starting from one infective, conditioning on the first event givesTaking the smallest solution in $[0,1]$ yields $q=1$ when $\beta \le \gamma$, and $q=\gamma /\beta =1/R_0$ when $R_0>1$. Thus $R_0=4$ still gives an early-lineage extinction probability of $1/4$. A deterministic SIR model describes mean-field growth but omits this individual realization uncertainty. In a finite closed stochastic SIR model, every epidemic eventually ends; the branching calculation concerns early extinction before substantial susceptible depletion, not that eventual end or an arbitrary peak-size cutoff.

These examples show different benefits of simplicity: a stability calculation can separate predictable from sensitive dynamics, a mode calculation can uncover a new mechanism, and a probability recursion can expose what an average trajectory hides. The same discipline also explains why the word “should” needs qualification. A simple exponential-decay model that confirms a familiar expectation can still calibrate a rate, serve as a control, or expose a faulty measurement. A useful mathematical model must earn its simplifications through the question it answers. Breaking an expectation is a powerful route to discovery when the contradiction is precise, robust to reasonable assumptions and experimentally consequential; novelty alone is not a validation criterion.