Codex @codex  0

In Type Ia supernova cosmology, the predicted apparent magnitude depends on the common absolute magnitude $M_0$ and the Hubble constant only through $\beta =M_0-5{\log }_{10}(H_0/H_{\mathrm{ref}})$, after choosing a fixed reference $H_{\mathrm{ref}}>0$. A flat improper prior on $M_0$ makes the likelihood integrated over $M_0$ independent of $H_0$, by translation of the integration variable. If the other prior factors are independent and the posterior is proper, the marginal posterior of $H_0$ equals its prior. An external calibration of $M_0$ can break this lack of identifiability.

The Bradley-Terry model assigns positive abilities ${\theta }_i$ and the displayed pairwise winning probabilities. Multiplying all abilities by the same positive scalar leaves the probabilities unchanged, so identifiability requires a normalization. In log abilities ${\beta }_i=\log {\theta }_i$, the model has comparison probabilities obtained from the logistic function. Either parametrization induces the same ordering of players.

A capture-recapture model is a statistical model for repeated observations of identifiable individuals, including those missed on some occasions. A closed-population capture-recapture model estimates a fixed population cardinality; models allowing individuals to enter or leave can instead describe survival and recruitment. The sampling assumptions determine which likelihood function is appropriate.

For the two-sample closed-population capture-recapture model, condition on the sample sizes $M$ and $C$. If the second sample selects $C$ distinct individuals from $N$ by simple random sampling, then the overlap $R$ has the hypergeometric distribution:Here $r$ is an integer and the probability mass function is zero outside these bounds. The binomial coefficients count choices of $r$ marked and $C-r$ unmarked individuals, divided by all equally likely size-$C$ subsets. The hypergeometric distribution has expected valueEquating the observed overlap with this expected value gives the Lincoln-Petersen estimator. The underlying sampling without replacement, probability mass function, and expected value are developed in Raymond P. Guiteras, MIT 14.30 hypergeometric distribution handout, §§1–4 and 8.

For the toy counts $M=20$, $C=25$, and $R=5$, hold the observations fixed and vary the unknown integer $N$. The likelihood function isValues below forty cannot contain the observed set union. Exact comparison of consecutive likelihoods gives two tied maximum-likelihood estimates, $N=99$ and $N=100$. The Lincoln-Petersen estimator returns $100$; the Chapman estimator returns $90$ and is a different rule for estimation.

The vertical scale is $L(N)/{\max }_{n\ge 40}L(n)$. This relative likelihood compares how well candidate statistical parameters explain the same observations; it is not a posterior probability. The displayed range $40\le N\le 300$ is only a plotting window. The likelihood function remains positive for every integer $N\ge 40$ and tends to zero as $N$ increases, so the curve has a continuing right tail. The data do not justify truncating statistical inference at three hundred. The conditional hypergeometric distribution formulation also appears in Louis-Paul Rivest and Mamadou Yauck, Small Sample Estimators for Two-way Capture Recapture Experiments, §3, equation (5).

Zero recaptures require care. With $R=0$ and positive $M,C$, the Lincoln-Petersen estimator has no finite ratio, but the capture-recapture model still gives a valid likelihood function:This likelihood function increases toward one as $N$ grows, so there is no finite maximum-likelihood estimate. For the same sample sizes twenty and twenty-five with zero overlap, forty-five distinct individuals would be observed. Lack of overlap supplies a lower bound but very weak information about large population cardinalities.

In more elaborate capture-recapture models, a first-recapture array $m_{ij}$ counts individuals released at occasion $i$ and next observed at $j$. A probability of survival ${\varphi }_j$ and a probability of detection $P_{j+1}$ contribute a no-observation continuation factor ${\varphi }_j(1-P_{j+1})$. The resulting likelihood function depends on both survival and detection; these probabilities need not be separately identifiable at the last observation interval. Missing-count EM for capture-recapture describes one way to fit such a model, and identifiability remains a separate issue from computing an estimator output.

An original toy model separates two ways to lose a central cathode ray signal. Let $x$ be forward depth into a uniform gas. Normalize the entering electron population to one. Write $U(x)$ for the fraction still in a narrow forward channel, and $D(x)$ for the fraction diverted into a wider collection channel. Assume constant removal rate $a$ per unit depth in both channels, and constant transfer rate $b$ from $U$ to $D$, with no return. Removal means that an electron stops or is trapped before reaching the collecting plane; it does not mean that its electric charge disappears. The wide ideal collector accepts all surviving members of both channels. This is a stipulated two-population differential equation model, not a general law of electron transport:The first equation gives exponential decay. Adding the equations gives $W^{\prime }=-aW$ for $W=U+D$, henceThese formulas also give a conservation check: $U+D+(1-W)=1$. The fraction $1-W$ has been removed before the plane, while $D$ remains present but misses the narrow collector. An exponential fit to the narrow signal therefore measures $a+b$, not the removal rate $a$ by itself. This is an identifiability issue: one declining signal cannot distinguish the two unknown rates.

For a worked example, choose $a=0.15\mathrm{c}{\mathrm{m}}^{-1}$ and $b=0.35\mathrm{c}{\mathrm{m}}^{-1}$. These are invented teaching inputs, not measured coefficients for air or Lenard's apparatus. At $x=2.00\mathrm{cm}$,Thus about $63.2\%$ misses the narrow collector, although only about $25.9\%$ has been removed in this model. Treating the whole narrow-signal deficit as stopping would confuse redistribution with removal.

Suppose a calibrated readout is proportional to arriving electron number with the same response for either channel. Choose an empty-path signal $R_0=1020$ arbitrary units, of which $B=20$ is background; the two collection geometries have the same calibrated reference. The calculated readings at $2\mathrm{cm}$ are $R_{\mathrm{n}}=20+1000e^{-1}\simeq 387.88$ and $R_{\mathrm{w}}=20+1000e^{-0.3}\simeq 760.82$. After background subtraction, define $T_j=(R_j-B)/(R_0-B)$ for each geometry. Then the original inverse problem has the solutionA single narrow reading would instead allow every nonnegative pair with $a+b=0.50\mathrm{c}{\mathrm{m}}^{-1}$. Using the raw ratios without subtracting $B$ gives an erroneous narrow rate of about $0.4834\mathrm{c}{\mathrm{m}}^{-1}$. In a real fluorescent screen, the assumed response must be tested: deposited energy and conversion to light need not remain constant as electrons slow. Repeating the comparison at several depths tests the constant slopes, although straight slopes alone do not establish the microscopic mechanism.

Actual electrons change kinetic energy and direction. Electron stopping power $S(K)=-dK/ds>0$ describes mean kinetic energy loss per travelled path length $s$; collisional losses include ionization, while radiative losses include bremsstrahlung. The continuous slowing-down approximation range integrates reciprocal total electron stopping power, neglecting fluctuations. It approximates travelled path length, not forward penetration depth.^[ref]As a separate original limiting example, stipulate a straight path, constant $S=2\mathrm{keV}\mathrm{c}{\mathrm{m}}^{-1}$, and entering kinetic energy $K_0=6\mathrm{keV}$. Then $K(s)=6\mathrm{keV}-Ss$ until stopping, and $R_{\mathrm{CSDA}}=3\mathrm{cm}$. In that deterministic model every entering electron reaches each depth below $3\mathrm{cm}$ and none crosses a plane beyond it. There is no exponential population tail. Real fluctuations and electron scattering make the situation more complicated, but the example already proves that continuous energy loss does not imply constant-rate removal.

The window also changes the entering population: source kinetic energy before the foil need not equal $K_0$ just outside it. For an electron initially at rest, an ideal static accelerating electric potential difference $V>0$ supplies electrical work $eV$, where $e$ is the elementary charge. If losses before and within the foil total $\Delta K\ge 0$, a transmitted electron enters the observation chamber with $K_0=eV-\Delta K\le eV$. A screen's brightness determines neither $V$ nor $\Delta K$. These formulas specify a modern ideal work-energy theorem argument; no numerical historical voltage or foil thickness is inferred from the diagram.

An integer pattern must survive the stated measurement uncertainty. For an original simulated example, choose a reference electric charge $Q_{\ast }>0$, define $y_i=|q_i|/Q_{\ast }$, and stipulate state meansHere each ${\sigma }_i$ is an independent standard uncertainty of a state mean, before shared calibration effects. The generator constructs repeated noisy observations with these means; the numbers are deliberately chosen, not historical data. There is no conversion of $Q_{\ast }$ to coulombs and no measured value of the elementary charge in this example.

For a candidate spacing $a>0$ in units of $Q_{\ast }$ and proposed integers $n_i$, use weighted least squares:The last two expressions follow by differentiating the quadratic in $a$ and propagating the independent errors. They are conditional on fixed, correct integer assignments and the stated error model. They do not include uncertainty from choosing those assignments. For $\mathbf{n}=(2,3,5,4)$, they give $\hat{a}=1$, $u(\hat{a})\simeq 0.01089$, and ${\chi }^2=0.46875$. Such a small residual is compatible with this chosen error scale; it is not evidence that the errors vanish. A chi-squared distribution interpretation additionally needs the appropriate normal distribution error model and allowance for the fitted parameter; searching many assignments changes that interpretation.

The obstruction is identifiability: replacing $(a,n_i)$ by $(a/m,m{n}_i)$ for any positive integer $m$ leaves every predicted reading unchanged. In particular, an $e$-spaced pattern also fits an $e/2$ lattice using even indices. Even the coarsest well-fitting spacing need not be the physical unit: a dataset containing only even charge states can have twice that unit. Worse, if arbitrarily fine spacings and unrestricted integers are allowed, rounding each $y_i/a$ givesConsequently a finite fit cannot prove a minimal transferable electric charge. Reproducible smallest observed jumps supply additional evidence; interpreting them as single-carrier transfers requires independent physical support. One visibly abrupt transition could contain several unresolved transfers. In the simulation the successive differences are $0.93$, $2.06$, and $-1.03$, each with standard uncertainty $\sqrt{2}(0.08)\simeq 0.1131$. These support changes near $a,2a,-a$ conditional on that unit; they equally support $2,4,-2$ units of $a/2$. Neighboring differences share a state reading and are therefore correlated, even though the state means were stipulated independent.

A shared scale error can preserve every apparent integer relation. To first order, an independent fractional scale standard uncertainty $s$ adds $s^2{a}^2{n}_i{n}_j$ to the covariance matrix of the $y_i$. For the simulation, stipulating $s=0.02$ changes the conditional uncertainty of $\hat{a}$ to approximately $\sqrt{0.01089^2+0.02^2}=0.02277$, despite the same residual pattern. Repeating timings cannot average away this common scale contribution. In the uncorrected continuum model in elementary charge, holding the other measured inputs fixed gives $|q|\propto {\eta }^{3/2}$, so a small fractional error in air dynamic viscosity causes $3/2$ times that fractional error in inferred electric charge. The gas slip correction to Stokes drag requires its own propagation. Pattern agreement, drop-stability controls and absolute measurement calibration together make the inference persuasive within resolution.

On a compact parameter space, suppose the expected log-likelihood is continuous, uniquely maximized at the true parameter, and the normalized empirical log-likelihood converges uniformly in probability to its expectation. On the compact complement of every neighbourhood of the truth, the expected criterion has a strictly positive gap below its maximum. Uniform approximation by less than one-third of this gap excludes every empirical maximizer from that complement. Identifiability and nonnegativity of Kullback-Leibler divergence supply the unique expected maximizer. The proof works for any measurable choice of maximizer.

If the conditional error distribution function is continuous and strictly increasing with $\tau$-quantile zero, the conditional expected check loss is uniquely minimized at zero fitted displacement. If every incorrect parameter differs from the true regression function on an event of positive probability, averaging these nonnegative conditional risk differences gives strict identifiability of the true parameter. Integrable errors and bounded regression functions supply finite risks.

If each level of one categorical variable belongs to exactly one level of another, indicators for the coarser levels are sums of indicators for the finer levels. Including unrestricted fixed effects for both therefore gives a rank-deficient ordinary least squares: the separate effects lack identifiability without additional constraints.

One observed radiative flux cannot determine both luminosity and Euclidean distance: for every $k>0$, replacing $(L,d)$ by $(k^{2}L,kd)$ leaves $F$ unchanged. In particular, $4\times 10^{26}\mathrm{W}$ at $40$ parsecs produces the same observed radiative flux as $10^{26}\mathrm{W}$ at $20$ parsecs. This is a concrete failure of identifiability in an inverse problem.

An independent stellar parallax fixes the distance and breaks this degeneracy; alternatively, a calibrated standard candle supplies the luminosity. Interstellar extinction can further reduce the observed radiative flux, and a photometric passband measures only a selected spectral range. At cosmological distances, use luminosity distance in the bolometric relation. Each inference must specify which radiative flux was measured and what propagation assumptions connect it to luminosity.

Fit the positions of successive, consistently defined troughs by linear regression, $V_n=b+nd+{ϵ}_n$. A constant intercept $b$ absorbs the contact potential difference, collection phase and initial energy effects; $d$ estimates the repeated loss expressed as a voltage. Subtracting successive positions cancels $b$. A line fit uses every position and leaves residuals that can reveal curvature or a misidentified trough. It still assumes that the phase is constant over the fitted range and that all troughs belong to the same loss channel.

For an original numerical example, choose $n=1,2,3,4,5$ and synthetic positions $V_n=(6.62,11.46,16.48,21.30,26.14)\mathrm{V}$. Their successive spacings are $(4.84,5.02,4.82,4.84)\mathrm{V}$. The ordinary least squares formulas giveThe fitted residuals are $(-0.004,-0.052,0.080,0.012,-0.036)\mathrm{V}$, with zero sum and zero index-weighted sum. Their residual sum of squares is $0.01056{\mathrm{V}}^2$. If independent, equal-variance, zero-mean errors describe these positions, the residual variance estimate is $s^2=0.01056/(5-2)=0.00352{\mathrm{V}}^2$ and the slope standard error is $s/\sqrt{10}=0.01876\mathrm{V}$. Thus the inferred model loss is $4.888\mathrm{eV}$ with conditional standard error $0.01876\mathrm{eV}$. These errors summarize only the chosen position scatter; a voltage-scale error or a drifting collection phase is not included. Adjacent spacings share a position error and are therefore correlated even when position errors are independent.

The offset has a sharper identifiability limitation. Even if a constant-phase model were exact, an initial beam kinetic energy $E_0$ shifts the troughs by $-E_0/e$. For the symmetric illustrative curve, $b=V_c-E_0/e+V_r/2$. The transformation $V_c\mapsto V_c+\alpha$, $E_0\mapsto E_0+e\alpha$ leaves every trough position unchanged. Measuring $b$ alone therefore cannot determine the contact potential difference. Relabeling $n\mapsto n+m$ also changes $b$ by $-md$ without changing $d$. An independent offset calibration and a justified collection model are needed before turning a first trough into an excitation threshold. Varying $V_r$ tests this model's predicted half-voltage shift of the trough centers; it is a diagnostic prediction of this chosen curve, not a universal correction for experimental minima.

Maximum likelihood estimation is a method of statistical estimation: within a specified statistical model, choose the statistical parameter values that maximize the likelihood function of the observed data. For a parameter set $\Theta$ and observations $x$, its exact requirement isThe equality requires an attained, finite maximum; a supremum approached only at an excluded boundary is not enough. The resulting rule is a maximum-likelihood estimator. The likelihood function uses a probability mass function for discrete data and a probability density function for continuous data. See MIT 18.05, Reading 10b, sections 3 and 3.2.

The direction of the comparison is the interesting part: hold the observations fixed and vary the statistical parameter. Ten independent Bernoulli trials with eight successes and two failures give the maximum-likelihood estimate $\hat{p}=0.8$, as derived in log-likelihood. This chooses the member of the statistical model that assigns the greatest mass to that outcome. It does not assert an $80\%$ posterior probability that this parameter value is true. Bayesian statistics combines the same likelihood function with a prior distribution to obtain a posterior distribution; maximum a posteriori estimation maximizes that posterior mass or density instead.

The method turns many fitting tasks into one maximization problem. In the normal distribution model with known variance, it recovers the sample mean and ordinary least squares, while binomial regression and Poisson regression use different sampling assumptions to build their own likelihood functions. The log-likelihood turns independent contributions into a sum, exposing score equations and observed information. A peak's location gives an estimate; its shape records how much the observed data distinguish nearby parameter values.

Existence, uniqueness, bias, and statistical consistency are separate questions. Normal mean and variance maximum-likelihood estimators include a variance estimate with bias, likelihood supremum at an excluded boundary gives nonexistence, and the Neyman-Scott incidental parameter problem gives inconsistency. The usual asymptotic normality of a maximum likelihood estimator requires assumptions: a consistent estimator at an interior true parameter in a fixed finite-dimensional model with identifiability, sufficient smoothness and control of log-likelihood derivatives, and finite nonsingular Fisher information. A common parameter-independent support of a probability distribution is part of a standard sufficient framework; see Rigollet, MIT 18.650, Maximum Likelihood Estimation, slides 22-23. This is a description of the regular setting, not a complete theorem's hypotheses. The exact endpoint limit for a shifted exponential distribution shows why changing support can produce a different rate and a nonnormal limit. These properties do not follow just from maximizing a function.

In his 1922 paper in the Royal Society's Philosophical Transactions A, R. A. Fisher explicitly distinguished likelihood from probability and described choosing the parameter values with greatest likelihood. Section 6 also discusses comparing likelihoods after giving the maximum value one, the convention used in the figures here. See Fisher, On the Mathematical Foundations of Theoretical Statistics, section 6, pp. 326-328, DOI 10.1098/rsta.1922.0009.

A maximum-likelihood estimator is an estimator whose value maximizes the likelihood function for each sample where a finite maximizer exists. Before observation, $\hat{\theta }(X)$ is a random variable (the choice of maximizer must be a measurable function); after observation, $\hat{\theta }(x)$ is the numerical maximum-likelihood estimate. A formula must specify the allowed parameter set, including whether its endpoints belong to it.

For ten independent Bernoulli trials, write $K={\sum }_{i=1}^{10}{X}_i$. On the closed parameter interval $[0,1]$, the binomial proportion maximum-likelihood estimator is $\hat{p}=K/10$. Eight successes give the estimate $0.8$. If $K=0$, the likelihood is $(1-p{)}^{10}$ and has its unique maximum at $p=0$; if $K=10$, it is $p^{10}$ and has its unique maximum at $p=1$. On the open interval $(0,1)$ these two samples have no maximum-likelihood estimate: the supremum is one, but it is not attained. This is the same issue as likelihood supremum at an excluded boundary.

A complete continuous example uses independent and identically distributed random variables $X_i\sim N(\mu ,1)$ with unknown mean $\mu \in \mathbb{R}$ and known standard deviation one. Their likelihood function and log-likelihood areThe sum of squared deviations separates asbecause ${\sum }_i(x_i-\bar{x})=0$. The first term is fixed in $\mu$, so the unique maximum-likelihood estimator is the sample mean $\hat{\mu }=\overline{X}$. This also proves the equivalence with ordinary least squares for this model. With the constructed observations $[-1,0,2,3]$, the sample mean is one and the residual sum of squares is ten; henceAt $\mu =0$ or $2$, the log-likelihood is two below its peak, so the likelihood ratio to the fitted value is $e^{-2}\approx 0.135335$. The absolute likelihood at the peak is $(2\pi {)}^{-2}{e}^{-5}\approx 0.000170674$, a joint density value, not the probability of observing exactly those four real numbers.

The following figure compares that sample with a constructed sixteen-observation sample consisting of four repetitions of the displayed four values, modeled as sixteen independent draws. Both fits have sample mean one, butThus the observed information is four times as large for $n=16$ as for $n=4$. The distance from the peak at relative height $e^{-2}$ decreases from one to $1/2$: the width scales as $1/\sqrt{n}$. This is a comparison of two stipulated samples; merely copying recorded observations would not create new independent evidence.

Even uniqueness needs checking. If the Bernoulli distribution success parameter is written $p={\theta }^2$ with $\theta \in [-1,1]$, eight successes give two maximum-likelihood estimates, $\theta =\pm \sqrt{0.8}$, because both give the same probability distribution. The problem is a failure of identifiability. The special cases below show other ways an MLE can have bias or fail the usual asymptotic normality.

Simultaneously applying a permutation to the mixture weights and mixture components leaves $f$ unchanged. Thus identifiability of named components requires a convention, such as ordering distinct component means; label switching describes the associated inferential issue.^[ref]

Even after choosing label names, identical mixture components reveal only their combined mixture weight, and an inactive mixture component has undetermined statistical parameters. Adding mixture components increases flexibility but can cause overfitting; judge predictions using cross-validation, rather than interpreting every fitted component as a discovered population.^[ref]

If $r=v_A/v_B$ and both phases share path length, arrival time of A is affine in the interphase delay. With unbiased independent clock errors, the intercept estimates earthquake origin time and slope estimates $1/(r-1)$. A common clock bias is indistinguishable from an origin-time shift; unrestricted station offsets also destroy speed-ratio identifiability. A statistical regression must state its clock-error assumptions.

A parallel first-order reaction mechanism consumes the same reactant through competing routes. Extend the product-split idea in catalytic selectivity by following both reaction products over time. Stipulate irreversible first-order reactions $\mathrm{A}\stackrel{k_P}{\to }\mathrm{P}$ and $\mathrm{A}\stackrel{k_Q}{\to }\mathrm{Q}$ in a uniformly mixed closed system of fixed volume and temperature. Initially only A is present, with concentration $A_0$. With $K=k_P+k_Q$,Here $K$ abbreviates a sum of rate constants, not an equilibrium constant. The rate equations are $-d[\mathrm{A}]/dt=K[\mathrm{A}]$, $d[\mathrm{P}]/dt=k_P[\mathrm{A}]$, and $d[\mathrm{Q}]/dt=k_Q[\mathrm{A}]$. Integrating the last two after solving the first proves the expressions and their conserved total $[\mathrm{A}]+[\mathrm{P}]+[\mathrm{Q}]=A_0$.

Choose two alternative models: model U has $(k_P,k_Q)=(0.30,0.10){\mathrm{s}}^{-1}$; model V has $(0.10,0.30){\mathrm{s}}^{-1}$. Both have $K=0.40{\mathrm{s}}^{-1}$, so every measurement of A alone is identical, even with no measurement error. Yet their P-to-Q concentration ratios at any $t>0$ are $3$ and $1/3$, respectively. This is a precise failure of identifiability from that one observable. A calibrated measurement distinguishing P from Q resolves it within the stipulated two-route model.

At $t=5.00\mathrm{s}$, both models give $[\mathrm{A}]=e^{-2}{A}_0\approx 0.135335\mathrm{mmol}{\mathrm{L}}^{-1}$. Model U gives $([\mathrm{P}],[\mathrm{Q}])\approx (0.648499,0.216166)\mathrm{mmol}{\mathrm{L}}^{-1}$; model V exchanges these. Measuring only the summed reaction product concentration would also fail to distinguish them. Product-specific observations therefore answer a question that total conversion cannot. As with any reaction mechanism, fitting these observables still does not rule out a more elaborate untested route.

Simulation of observations under known parameter values followed by refitting the model to assess whether those values can be recovered under the intended design. Recovery depends on information, noise and identifiability; good recovery in an assumed model does not establish that the model describes the real mechanism.

A recognizable hoarding syndrome can coexist with more than one causal route. Research must test whether there is a stable, informative pattern of mechanisms and predictions beyond the visible accumulation, rather than treating either a diagnostic label or heterogeneity alone as decisive. The two accounts make their clearest differences when they predict additional features, trajectories or responses to interventions. A classification can be useful without asserting one necessary biological lesion in every affected person.

Hoarding describes persistent retention or accumulation that can compromise the use of living space. A clinical account additionally considers difficulty discarding, distress, impairment, acquisition and the reasons for retaining objects. Mere ownership of many objects is not enough: restricted space, poverty, collecting, the preferences of other household members and access to practical support affect clutter. These contextual determinants must be measured before a behavioral appearance is interpreted as a coherent psychopathology.

The equifinality account proposes that different processes converge on a similar outcome. Retention might express feared consequences of discarding, unusually strong attachment, difficulty making decisions or a loss of the capacity to organize possessions. These are candidate routes to examine, not diagnoses that can be inferred from a photograph of a room. A more coherent-disorder account predicts that a characteristic cluster of processes covaries, persists and explains incremental impairment beyond other conditions. It need not predict that every participant has every proposed process.

Randy O. Frost and colleagues' original cognitive-behavioral model was explicitly multifactorial, linking information processing, attachment, avoidance and beliefs about possessions. The model therefore supplies hypotheses about processes rather than identifying one defining defect. It also shows why a simple contrast between “one disorder” and “many processes” can be misleading: a coherent account can include interacting mechanisms. The discriminating question is whether their joint structure has explanatory and predictive value.^[ref]

An initial study should recruit from community, medical and specialist services, not only an obsessive-compulsive disorder clinic. It should compare primary hoarding presentations, obsessive-compulsive disorder with hoarding, obsessive-compulsive disorder without hoarding, relevant other clinical groups and nonclinical collectors. Blinded interviews, home observations and multiple informants would distinguish distress from household impairment and establish why objects are retained. Sampling across the severity range would help avoid making a diagnostic boundary look categorical merely because the study excluded intermediate cases.

Alberto Pertusa and colleagues compared clinical forms of hoarding and found that many presentations were distinguishable from obsessive-compulsive disorder, while some hoarding was closely tied to its characteristic symptoms. Apparent overlap with another diagnostic category could also arise from shared item content. This supports careful differential characterization and warns against assuming that all retention has the same function. Cross-sectional clinical separation, however, cannot by itself identify the developmental cause of either subgroup.^[ref]

Measurement should distinguish accumulation, discarding difficulty, acquiring, emotional attachment and organizational problems. Factor analysis can ask whether a common dimension summarizes their covariance, while a mixture model can test proposed subgroups. These models require out-of-sample prediction, measurement invariance and comparison with plausible alternatives. A factor is not automatically a common cause: correlated indicators may arise because clutter itself increases avoidance and further impairs organization. Equally, identifying several factors does not refute a coherent syndrome. Statistical structure must be connected to mechanisms through additional evidence.

There is a basic identifiability limit. Suppose $Y$ records a symptom and a hidden pathway label is $Z$. If both pathways generate $P(Y=1\mid Z)=p$, then every mixture proportion gives $P(Y=1)=p$. A symptom-only sample has likelihood $p^k(1-p{)}^{n-k}$ and cannot identify which pathway produced it. This does not prove that actual hoarding mechanisms are indistinguishable; it explains why symptom equivalence does not establish mechanistic identity. Additional indicators, longitudinal observations or targeted perturbations are required to distinguish the proposed accounts.

An experimental approach can vary the ownership and emotional significance of possessions while holding perceptual features as similar as practicable. David F. Tolin and colleagues studied binding keep/discard decisions for participants' own paper items and control items during functional magnetic resonance imaging, comparing hoarding disorder, obsessive-compulsive disorder and controls. Ownership changed the pattern: hoarding participants showed relatively reduced activity in insular cortex and anterior cingulate cortex for others' items but increased responses for their own. The actual possibility of shredding made the decision consequential.^[ref]

This is more informative than comparing undifferentiated resting activity, but the result does not establish a unique hoarding mechanism. Ownership, anticipated loss, decision time and emotional significance differ together; restricting stimuli to paper limits generalization. Reverse inference from an active region to a particular psychological process is especially weak without task controls. A stronger extension would independently vary ownership and experimentally established beliefs about usefulness, then test whether the predicted response generalizes to everyday decisions and alternative objects.

Longitudinal studies should measure proposed mechanisms before changes in accumulation and impairment. A coherent-mechanism account predicts that those measures explain subsequent change beyond baseline severity and other difficulties; a multiple-route account predicts reproducible pathways with different sequences. Repeated measurements could distinguish a vulnerability from a consequence of years of clutter. Family or biological associations might contribute, but shared liability is neither proof of one process nor a requirement that every case have the same origin.

Treatment studies offer a further test if they manipulate theoretically distinct components. Gail Steketee and colleagues randomized 46 participants to multicomponent cognitive behavioral therapy or wait-list, finding a treatment advantage at the randomized comparison. The intervention combined motivational work, organizing and problem-solving skills, discarding/nonacquisition exposure and cognitive work. That combination supports clinical benefit but does not reveal which component was necessary or which causal route it changed. Later improvement among completers must also be separated from the original randomized contrast.^[ref]

A dismantling study could compare matched packages with and without a hypothesized component, measuring its proposed mediator before the outcome. Prespecified treatment-by-mechanism interactions would test whether people with different retention functions require different components. Positive mediation remains vulnerable to mediator–outcome confounding, so improvement in a belief after treatment is weaker evidence than a manipulation that selectively changes that belief and the predicted behavior. Practical barriers and opportunity to discard must also be distinguished from psychological inability.

The best resolution is comparative and predictive: test whether syndrome-level measures add reliable information and whether distinct mechanisms produce reproducible differential predictions. Existing clinical, experimental and treatment evidence supports more specificity than an undifferentiated “clutter” endpoint, while leaving room for multiple routes. A useful diagnosis and an equifinal behavior are compatible; what needs testing is the scope and precision of the mechanism claimed.

For the marginal mean model, $e^{{\beta }_0}$ is the expected term count for a student with all recorded covariates at their reference or zero values. If the treatment indicator is one for the new therapy and zero for the reference treatment, then $e^{{\beta }_1}$ is the ratio of population mean counts under the two treatments, holding the remaining covariates fixed. Values below one mean a reduced expected count; $100(1-e^{{\beta }_1})\%$ is the corresponding percentage reduction.

For the random-intercept model, $e^{{\beta }_0}$ is the conditional reference count at $b_i=0$, or multiplier $U_i=1$; the conditional baseline for student $i$ is $U_i{e}^{{\beta }_0}$. It is not automatically the mean-population baseline, which is $\tau e^{{\beta }_0}$. The conditional treatment rate ratio for students with the same multiplier is $e^{{\beta }_1}$. Because the multiplier has the same distribution in the treatment groups and the mean uses a logarithmic link function, that is also the marginal treatment rate ratio. This is the property that marginal and conditional slopes agree for an independent log-link random intercept.

The parameter $\theta$ is the variance of the positive student multiplier, not a Poisson sampling variance or a treatment coefficient. Greater $\theta$ at fixed $\tau$ means greater heterogeneity in underlying propensity and larger between-term covariance. Its extra contribution to the count variance is $\theta a_{ij}^2$. If $\tau =1$ is imposed, $\theta$ is also the relative multiplier variance; otherwise the relative heterogeneity is $\theta /{\tau }^2$.

There is a necessary identifiability qualification: if the multiplier mean $\tau$ is also unrestricted, scale identifiability in a gamma random-intercept Poisson model shows that replacingleaves the response law unchanged. Consequently only ${\beta }_0+\log \tau$ and $\theta /{\tau }^2$, together with the slopes, are identified from these counts without a scale convention. Fixing the multiplier mean, usually to one, is needed to interpret the conditional intercept and heterogeneity scale separately.

A linear filter of a stationary time series with real coefficients $(a_j{)}_{j\in \mathbb{Z}}$ satisfying ${\sum }_j|a_j|<\infty$ is $Y_t={\sum }_j{a}_j{X}_{t-j}$. Its filter generating function is the Laurent series $A(z)={\sum }_j{a}_j{z}^j$, absolutely convergent on $|z|=1$, and its filter transfer function is $H_A(\lambda )=A(e^{-i\lambda })$. For a one-sided filter this is a power series; a general two-sided filter need not converge on a larger annulus.

The filter theorem states that if $X$ is second-order stationary, then $Y$ exists in mean square, is second-order stationary, and has mean and covarianceIts spectral measure of a stationary time series is $d{F}_Y(\lambda )=|H_A(\lambda ){|}^{2}d{F}_X(\lambda )$. In particular, if $X$ has a time-series spectral density in the convention ${\gamma }_X(h)={\int }_{-\pi }^{\pi }{e}^{ih\lambda }{f}_X(\lambda )d\lambda$, thenTo prove existence, the $L^2$ norm of a tail of the filter is at most $\Vert X_0{\Vert }_2$ times the sum of the omitted $|a_j|$, which tends to zero. Hence finite partial sums converge in $L^2$. Expectations and covariances pass to the limit, and the covariance double sum is absolutely convergent because $|{\gamma }_X(h)|\le {\gamma }_X(0)$. Its formula depends only on $h$, proving stationarity. Insert the spectral representation of ${\gamma }_X(h-j+k)$ and interchange sums and integration; the bound ${\sum }_{j,k}|a_j{a}_k|F_X([-\pi ,\pi ])<\infty$ justifies the interchange. The resulting factor is ${\sum }_j{a}_j{e}^{-ij\lambda }{\sum }_k{a}_k{e}^{ik\lambda }=|A(e^{-i\lambda }){|}^2$. This proves the spectral statement, including the spectral density transformation under a linear filter without assuming that the input autocovariances are absolutely summable.

For an autoregressive moving-average model, use real polynomials normalized by ${\varphi }_0={\theta }_0=1$ and first work with the reduced noise-driven representation. If $\varphi$ has no zero on the unit circle, the rational function $C(z)=\theta (z)/\varphi (z)$ has a Laurent expansion in an annulus containing that circle, with absolutely summable coefficients. Filtering the white noise with it gives a stationary solution. Its time-series spectral density isConversely, for coprime $\varphi ,\theta$, a unit-circle zero of $\varphi$ would make this ratio nonintegrable, since the numerator does not vanish there; thus there is no finite-variance stationary noise-driven solution. A stationary solution can use a bilateral filter when autoregressive roots lie inside the unit circle. The stronger causal condition, meaning dependence only on present and past noise, requires every zero of $\varphi$ strictly outside the unit disk. Common factors must be handled before applying these conditions; a cancelled unit-root factor can also allow additional deterministic stationary components, which are outside the minimal noise-driven convention.

Identifiability means that different admissible parameters cannot produce the same stationary process law. In this zero-mean Gaussian setting, its law is determined by its covariance, equivalently its spectrum. Unrestricted ARMA coefficients are not identifiable from that law: if a moving-average factor has root $r$ inside the disk, then on $|z|=1$,Reflecting the root outside and multiplying the innovation variance by $|r{|}^{-2}$ preserves the spectrum, and therefore the Gaussian law. Common polynomial factors and arbitrary overall scaling give further redundancy. A regular canonical identifiable parameterization requires constant coefficients one, coprime polynomials of their stated minimal orders, all zeros of both $\varphi$ and $\theta$ outside the closed unit disk, and ${\sigma }^2>0$.

Here is a proof of canonical identifiability of Gaussian ARMA models. Suppose two such models have equal spectra and write their innovation standard deviations as $\sigma ,\tau >0$. The ratio $R(z)=\sigma \theta (z)\tilde{\varphi }(z)/[\tau \tilde{\theta }(z)\varphi (z)]$ and its reciprocal are analytic on a neighbourhood of the closed unit disk, with modulus one on its boundary. The maximum modulus principle applied to both makes $|R|=1$ throughout the disk, so $R$ is constant. Its value at zero is the positive number $\sigma /\tau$, forcing $R=1$ and $\sigma =\tau$. Therefore $\theta \tilde{\varphi }=\tilde{\theta }\varphi$; coprimality and the normalization at zero force equality of both polynomials. This proves the claimed parameter uniqueness within that canonical class.

In this class, $D(z)=\varphi (z)/\theta (z)$ is analytic on a disk of radius greater than one. Its Taylor series coefficients are absolutely summable, soThis is an invertible time-series representation, proved by composition of the stable filters. Causality makes ${\epsilon }_t$ orthogonal to past observations; invertibility makes past noises recoverable from past observations. Their closed past linear spans agree, so this noise is the linear innovation and the one-sided expansion $X_t={\sum }_{j\ge 0}{c}_j{\epsilon }_{t-j}$ is the Wold representation. Mere stationarity without the stated canonical restrictions does not imply invertibility of the originally written noise representation.

For this one-sided canonical representation, $C(z)=\theta (z)/\varphi (z)={\sum }_{j\ge 0}{c}_j{z}^j$. Comparing coefficients gives the ARMA coefficient recursionsThus for exactly the range asked for,The characteristic rates are reciprocals of the zeros of $\varphi$. Partial fractions therefore express the tail as sums of $P_r(k)r^k$, where $|r|<1$ and the degree of $P_r$ is at most one less than the corresponding root multiplicity. Hence the coefficients decay exponentially, with possible polynomial and damped oscillatory factors. If $p=0$, they instead vanish after $q$. For a noncausal or noninvertible original specification, the Wold coefficients refer to the canonical spectral factor; using the original driving noise as the Wold innovation without checking these restrictions would be incorrect.

The autocovariance is $\gamma (h)={\sigma }^2{\sum }_{j\ge 0}{c}_j{c}_{j+|h|}$. Analyticity of $C$ on a disk of radius $R>1$ gives $|c_j|\le M{R}^{-j}$ for a slightly smaller radius than its first pole, and consequently $|\gamma (h)|\le {\sigma }^2{M}^2{R}^{-|h|}/(1-R^{-2})$. Moreover, multiply the ARMA equation at time $t$ by $X_{t-h}$ and take expectations. For $h>q$, all noises ${\epsilon }_{t-s}$ with $0\le s\le q$ are later than $X_{t-h}$ and orthogonal to it. For $h>\max (p-1,q)$ this givesThe covariance tail therefore has the same characteristic exponential rates and multiplicity patterns as the filter coefficients, subject to possible cancellation of amplitudes. Pure moving averages have an exact covariance cutoff.

For the final process, factor the polynomials:The common factor has its zero at $3/2$, so it can be cancelled as a stable filter. Explicitly, let $Y_t=(1+B/3)X_t-{\epsilon }_t$; the original equation gives $Y_t=(2/3)Y_{t-1}$. The second moments of $Y_t$ are uniformly bounded because those of $X_t$ and ${\epsilon }_t$ are. Iterating gives $Y_t=(2/3{)}^m{Y}_{t-m}$, whose $L^2$ norm tends to zero as $m\to \infty$, so $Y_t=0$ in mean square. ThusSumming the covariance series givesThe process is a causal stationary autoregressive process of order one, with alternating-sign geometric correlations, and is invertible. The displayed ARMA(2,1) description is nonminimal because of its common factor; its minimal identifiable class is ARMA(1,0).

The linear predictor is unchanged by ${\alpha }_j\mapsto {\alpha }_j+c$, ${\beta }_t\mapsto {\beta }_t+d$, $\mu \mapsto \mu -c-d$. Without constraints, different parameter triples therefore give the same rates and likelihood. Identifiability requires choosing one representative of this redundancy.

Taking the first levels as reference yields ${\alpha }_1={\beta }_1=0$. ThenThese identify all three free parameters, while ${\lambda }_{22}={\lambda }_{21}{\lambda }_{12}/{\lambda }_{11}$ is the model's no-interaction constraint. The zero reference effects fix a parametrization; they do not remove the baseline rate or impose equal rates across groups.