Pith. sign in

REVIEW 4 major objections 5 minor 25 references

Equilibrium Analysis of Discrete Stochastic Population Models with Gamma Distribution

T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read At equilibrium, the intrinsic growth rate $r$ in stochastic logistic and Ricker models is fixed by the gamma shape parameter $k$ and the perturbation variance, splitting into two branches $r_+$ and $r_-$.

desk verdict Moment-consistency formulas for gamma-shaped noise in two discrete population models, but the "equilibrium" and "alternative stable states" framing outruns the evidence. read the letter →

arxiv 2411.15859 v1 pith:YPLQMREJ submitted 2024-11-24 q-bio.PE

Haiyan Wang This is my paper · ORCID
classification q-bio.PE MSC 39A5092D25
keywords gammadistributionstochasticlogisticequationRickerstationaryintrinsicgrowthratemomentmatchingalternativestablestatespopulationdynamics
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper aims to establish that at equilibrium, $r$ in two canonical discrete stochastic population models is not a free parameter but is tied to the shape of the population-size distribution and the strength of environmental noise. Assuming the equilibrium population size $X_t$ is gamma-distributed, and matching only the mean and variance at a single time step, the paper derives explicit formulas for $r$: a closed-form quadratic-root formula for the modified stochastic logistic equation $X_{t+1}=rX_t(1-X_t^n)\varepsilon_t$ and a transcendental equation for the stochastic Ricker equation $X_{t+1}=X_t e^{r(1-X_t)}\varepsilon_t$. Both formulas yield two growth-rate branches, $r_+$ and $r_-$, interpreted as alternative stable states. A further conclusion is that $r$ is independent of the gamma scale parameter $\theta$, so growth depends on population structure, not absolute abundance. A reader should care because these relations turn a purely statistical description of population fluctuations into testable constraints on a core ecological parameter.

What carries the argument

The load-bearing object is the gamma distribution $\mathrm{Gamma}(k,\theta)$ — the two-parameter positive-density law with mean $k\theta$ and variance $k\theta^2$ — used as the assumed stationary distribution of population size. The mechanism is moment consistency: the paper imposes that one stochastic update leaves both the mean and variance unchanged, $E[X_{t+1}]=E[X_t]$ and $\mathrm{Var}(X_{t+1})=\mathrm{Var}(X_t)$. All calculations pass through the identity $E[X_t^n e^{-sX_t}] = \frac{\Gamma(k+n)}{\Gamma(k)} \frac{\theta^n}{(1+s\theta)^{k+n}}$, proved in the appendix, which gives the higher moments needed to write both conditions. For the logistic model, eliminating $\theta$ between the mean and variance equations produces a quadratic in $r$; for the Ricker model, it produces a transcendental equation in $r$ and $k$ alone.

What would settle it

Simulate equations (1) and (2) numerically from a non-gamma initial distribution with independent mean-one multiplicative noise, run to stationarity, and test whether the empirical stationary distribution is gamma; if the fitted shape $k$ and noise variance satisfy (34) or (53) with the simulated mean per-capita growth rate, the relations hold, and a systematic mismatch for many parameter settings would falsify them.

Watch

Extended reading notes

Core claim

The central claim is that the equilibrium conditions $E[X_{t+1}]=E[X_t]$ and $\mathrm{Var}(X_{t+1})=\mathrm{Var}(X_t)$, together with the gamma assumption $X_t\sim\mathrm{Gamma}(k,\theta)$, force the intrinsic growth rate $r$ in the modified logistic model to satisfy a quadratic relation whose roots are $r_\pm = \frac{A+n \pm \sqrt{n^2+AB+2nB}}{A+2n}$ with $A=k(1-Q)+1$ and $B=\frac{\mathrm{Var}(\varepsilon_t)}{\mathrm{Var}(\varepsilon_t)+1}(k+1)$, and in the Ricker model to satisfy the transcendental relation $2 e^{r/(k+1)} - \bigl((1+\mathrm{Var}(\varepsilon_t)) e^{2r}\bigr)^{1/(k+2)} = 1$. In both models, positive noise variance admits two positive growth rates, whereas zero variance recovers a single branch as a limiting case. Because $\theta$ is eliminated between the mean and variance equations, $r$ depends only on $k$ and $\mathrm{Var}(\varepsilon_t)$, not on the scale of the population distribution.

Load-bearing premise

The assumption that the equilibrium population size is exactly gamma-distributed and that stationarity is fully captured by matching the first two moments across one time step; the paper does not prove convergence to a gamma distribution, so if the true stationary law is not gamma, the derived $r(k, \mathrm{Var}\,\varepsilon)$ relations need not describe actual equilibria.

Editorial extensions

If this is right

  • At equilibrium, $r$ is fixed by the gamma shape parameter $k$ and the noise variance, so observed growth rates can be checked against distributional shape.
  • For positive noise variance, two growth rates $r_+$ and $r_-$ coexist, representing alternative stable states with different resilience and extinction risk.
  • Because $r$ does not depend on $\theta$, populations of very different average sizes can share the same equilibrium growth rate if their shape parameters match.
  • In the $n=1$ logistic case, the formula reduces to the earlier result with feasible ranges $\mathrm{Var}(\varepsilon_t) \le 0.5$ and $k \le 1/\mathrm{Var}(\varepsilon_t) - 2$, tying feasible shapes to noise strength.
  • For the Ricker model, the upper branch can exceed the deterministic chaotic threshold near $2.69$ at small $k$, suggesting that noise can stabilize dynamics that would be chaotic without perturbations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same two-moment matching procedure could be rerun with other stationary families, such as log-normal or negative binomial laws, generating parallel $r(k, \mathrm{Var}\,\varepsilon)$ constraints that could be compared across species to distinguish which distributional assumption best fits a given time series.
  • The two-branch structure implies a testable dynamical prediction: a large transient disturbance could move a population from the high-growth branch $r_+$ to the low-growth branch $r_-$ without any change in gamma shape or noise variance, producing apparent regime shifts with no underlying parameter change.
  • Because the derivation matches exactly two moments while the gamma family has two parameters plus $r$, the independence of $\theta$ may be an artifact of the matching scheme; matching a third moment or the full distribution would reveal whether the $r(k, \mathrm{Var}\,\varepsilon)$ relation survives stronger consistency requirements.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper postulates that at equilibrium the population size in the stochastic logistic map (1) and the stochastic Ricker map (2) follows a gamma distribution with shape parameter k and scale parameter θ, and it replaces stationarity by equality of the first two moments across one time step. Under that ansatz, the paper derives explicit relations for the intrinsic growth rate r: a quadratic formula, Eq. (34), for the logistic map and a transcendental equation, Eq. (53), for the Ricker map. Each relation gives two branches r± and eliminates θ. The paper interprets the two branches as alternative stable states and interprets the absence of θ as evidence that r is governed by internal population dynamics rather than by population scale.

Significance. If the claims were established, the paper would provide a compact moment-consistency criterion linking the intrinsic growth rate to the gamma shape parameter and the noise variance, with a direct ecological interpretation. The algebraic derivations are mostly sound: the elimination of θ is explicit, the Ricker calculation uses the Laplace-transform identity (54) elegantly, and the parameter-free status of the final relations is a strength. However, the central claim is contingent on an unproved gamma-stationarity ansatz, and the paper's own Section 4 acknowledges the need for a rigorous foundation. The contribution is therefore best assessed as a two-moment matching exercise rather than an equilibrium or stability analysis of the original stochastic processes.

major comments (4)
  1. [§2.1, Eqs. (4)–(5)] The gamma-stationarity ansatz is assumed, not derived, and stationarity is enforced only through equality of the first two moments. A true stationary distribution must have all moments invariant, so Eqs. (34) and (53) are at most necessary conditions for first-two-moment invariance, not characterizations of stationary distributions. The third moment provides a direct falsification test. For the logistic map with n=1, E[X_{t+1}^3] = r^3 E[ε^3] (E[X_t^3] - 3E[X_t^4] + 3E[X_t^5] - E[X_t^6]). Taking k=1 and Var(ε)=0.01, the r+ branch from Eq. (34) gives r≈1.99 and θ≈0.249, and the bracketed expression is negative (≈ -0.011), so E[X_{t+1}^3] < 0 while E[X_t^3] > 0. This is impossible for a nonnegative random variable X_{t+1}. The same insufficiency applies to Eq. (53): the Ricker derivation never imposes third- or higher-moment stationarity. The abstract and Section 4 overstate the results. The paper should either prove gamma stationarity, or at least moment-closedness of the maps, or explicitly reframe the results as two-moment matching conditions.
  2. [§2.4 and §3.2] The two roots r± are algebraic solutions for the parameter r in a model with fixed dynamics; they are not alternative stable states of a dynamical system. In Eqs. (1) and (2), r is a constant parameter, so r+ and r− correspond to different models with different growth rates, not to two coexisting attractors of one system. Statements such as "alternative stable states of the intrinsic growth rate" and the associated resilience and extinction discussion should be replaced by wording about two admissible parameter branches of the moment-consistency relation.
  3. [Eqs. (35)–(36)] In the n=1 reduction, substituting B = Var(ε)/(Var(ε)+1)(k+1) into Eq. (35) gives sqrt(1 - Var(ε)(k+3)/(Var(ε)+1)), not sqrt(1 - Var(ε)(k+2)/(Var(ε)+1)) as printed in Eq. (36). The printed formula is algebraically inconsistent with the preceding line. This also affects the parameter bounds quoted in Section 2.4, since the discriminant condition changes unless the bound is derived from the correct expression.
  4. [§2.4 and §4] The claimed independence of r from θ is an artifact of eliminating θ via the mean condition Eq. (16), which fixes θ in terms of r, k, and n. Since θ is not a free parameter under the equilibrium conditions, the interpretive statement that r is "primarily governed by internal dynamics rather than population scale" overreads a parameterization choice. The manuscript should state that the moment-consistency relation does not involve θ because θ is slaved to the other parameters by mean invariance.
minor comments (5)
  1. [Abstract] The abstract contains grammatical errors, for example "at equilibrium examines with the gamma distribution"; the sentence should be rewritten.
  2. [Figures 1–2] Figure 1 contains a placeholder caption, "(a) Caption 4", and neither figure states which panel corresponds to which n and Var(ε); all panels should be fully labeled and captioned.
  3. [§2.4] The paper uses the companion preprint [26] as the source of the feasibility bounds, but since Eq. (36) is derived in the present paper, the bounds should be rederived directly from Eq. (36) rather than cited from [26].
  4. [References] References [20]–[22] are duplicates of [12]–[14], and the titles in [24]–[25] contain typos; the bibliography should be cleaned up.
  5. [§2.2, Eq. (10)] The symbol n is used both for the nonlinearity exponent in Eqs. (1)–(2) and as a generic moment index in Eq. (10), which is potentially confusing; a different index, such as m, should be used for moments.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the r–k–Var(ε) relations are derived algebraically from explicitly stated moment assumptions, and the sole self-citation is not load-bearing.

full rationale

The derivation chain is conditional but not circular. The paper explicitly postulates X_t ~ Gamma(k, θ) at equilibrium (Eq. 4) and replaces stationarity with first-two-moment matching (Eq. 5). All subsequent equations, including (16), (34), and (53), are algebraic consequences of these stated assumptions; no parameter is fitted to data and then relabeled as a prediction. Eq. (34) solves a quadratic for r that follows from the mean and variance conditions, and the claimed independence of r from θ arises from eliminating θ via the mean condition, not from assuming the conclusion. The only self-citation, [26], is used to check the n=1 special case ('This is exactly the formula in [26]') and to supply parameter ranges, but the same formula is derived independently within this paper, so the citation is not load-bearing. The manuscript itself flags its main limitation in Section 2.1: 'Because of the nonlinearity and random perturbation in (1), Xt+1 may not necessarily follow the gamma distribution,' and in the Conclusion: 'establishing a robust mathematical foundation to rigorously prove the relationships between model parameters and their stationary distributions is crucial.' These are honest caveats about sufficiency and rigorous stationarity, not evidence that the derivation reduces to its inputs. The interpretive claim that r+ and r− are 'alternative stable states' is an overreach beyond the algebra, but overinterpretation is not circularity. Therefore the paper receives a circularity score of 0.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The paper introduces no new particles, forces, dimensions, or other invented entities. It imports the gamma distribution as an assumed stationary law, imposes first-two-moment matching, and varies standard model parameters k, theta, Var(epsilon), and n, none of which are estimated from data.

free parameters (3)
  • Gamma shape parameter k
    Assumed shape of the equilibrium gamma distribution; r is expressed as a function of k and the feasible range of k is explored numerically in Figures 1 and 2.
  • Gamma scale parameter theta
    Assumed scale of the equilibrium gamma distribution; it cancels out of the central r relations, so it does not enter the final formulas.
  • Perturbation variance Var(epsilon)
    Variance of the multiplicative noise, treated as an input parameter; the paper explores values from 0 to about 0.5 and shows it constrains the existence of real r branches.
assumptions (4)
  • domain assumption X_t follows a Gamma(k, theta) distribution at equilibrium
    Assumed in Section 2.1, Eq. (4). This is the load-bearing modeling choice that makes the moment formulas possible.
  • domain assumption At equilibrium E[X_{t+1}] = E[X_t] and Var(X_{t+1}) = Var(X_t)
    Eq. (5) in Section 2.1. Used to close the moment system, but equality of the first two moments is weaker than full stationarity of the distribution.
  • domain assumption The multiplicative noise epsilon_t is independent of X_t with E[epsilon_t] = 1
    Stated in Section 1 and used to factor E[epsilon_t] out of moment expressions in the logistic and Ricker derivations.
  • standard math Gamma integral identity E[X^n e^{-sX}] = Gamma(k+n)/Gamma(k) * theta^n / (1 + s*theta)^{k+n}
    Proved in the Appendix, Eq. (54), and used to evaluate moments in the Ricker model.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Equilibrium Analysis of Discrete Stochastic Population Models with Gamma Distribution." pith.science (2026). https://pith.science/paper/YPLQMREJ

@misc{pith2026241115859,
  author       = {Pith},
  title        = {Pith review of: Equilibrium Analysis of Discrete Stochastic Population Models with Gamma Distribution},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YPLQMREJ}},
  note         = {Machine review of arXiv:2411.15859}
}
read the original abstract

This paper analyzes the stationary distributions of populations governed by the discrete stochastic logistic and Ricker difference equations at equilibrium examines with the gamma distribution. We identify mathematical relationships between the intrinsic growth rate in the stochastic equations and the parameters of the gamma distribution with a small stochastic perturbation. We present the biological significance of these relationships, emphasizing how the stochastic perturbation and shape parameter of the gamma distribution influence population dynamics at equilibrium. Furthermore, we identify two branches of the intrinsic growth rate, representing alternative stable states corresponding to higher and lower growth rates. This duality provides deeper insights into population stability and resilience under stochastic conditions.

Figures

Figures reproduced from arXiv: 2411.15859 by the authors.

Figure 1
Figure 1. Plots of r in term of k for the modified logistic equation it is easy to see that a population with a larger gamma shape parameter has a more sym￾metrical and regular population distribution. From [PITH_FULL_IMAGE:figures/full_fig_p011_1.png] view at source ↗
Figure 2
Figure 2. Plots of r in term of k for Ricker equation For the deterministic Ricker model given by (2), population dynamics, including stable growth, periodic oscillations, chaos, or extinction, are closely tied to the value of r [17, 10]. For 0 < r < 2.69, the population typically grows and stabilizes at a nonzero steady state or exhibits damped oscillations or populations alternating between large and small sizes in successi… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

25 extracted references · 18 canonical work pages

  1. [26]

    Gamma Distribution for Equilibrium Analysis of Discrete Stochastic Logistic Population Models

    Haiyan Wang, (2024). Gamma Distribution for Equilibrium Analysis of Discrete Stochastic Logistic Population Models . arXiv:2411.10167

    work page Pith review arXiv 2024

  2. [1]

    Allen, L. J. S. (2010). An Introduction to Stochastic Processes with Applications to Biol- ogy. Chapman and Hall/CRC

    work page 2010

  3. [3]

    R. M. May, (1974) Stability and Complexity in Model Ecosystems , Princeton U.P., Prince- ton

    work page 1974

  4. [5]

    C. S. Peters, M. Mangel, R. F. Costantino, Stationary Distribution of Population Size in Tribolium, Bulletin of Mathematical Biology, 51(5), 625–638, 1989

    1989

  5. [6]

    Dennis and R

    B. Dennis and R. F. Costantino, Analysis of Steady-State Populations with the Gamma Abundance Model and its Application to Tribolium , Ecology, 69(4), 1200–1213, 1988

    1988

  6. [7]

    R. F. Costantino and R. A. Desharnais, Gamma Distributions of Adult Numbers for Tribolium Populations in the Region of their Steady-States , Journal of Animal Ecology, 50, 667–681, 1981. 19

    work page 1981

  7. [8]

    Matis, Thomas R

    James H. Matis, Thomas R. Kiffe b, Eric Renshawc, Janet Hassan, A simple saddle- point approximation for the equilibrium distribution of th e stochastic logistic model of population growth, Ecological Modelling, 161 (2003) 239–248

    work page 2003

  8. [9]

    Wikipedia contributors, ”Gamma distribution,” Wikipedia, The Free Encyclopedia, 2024

    2024

Show all 25 references
  1. [10]

    Kot, Elements of Mathematical Ecology

    M. Kot, Elements of Mathematical Ecology . Cambridge University Press, 2001

    2001

  2. [11]

    Dennis, B., Assas, L., Elaydi, S. et al. Allee effects and resilience in stochastic popula- tions. Theor Ecol, 9, 323–335 (2016). https://doi.org/10.1007/s1208 0-015-0288-2

    2016 doi

  3. [12]

    Braverman, A

    E. Braverman, A. Rodkina (2013). Difference equations of Ricker and logistic types under bounded stochastic perturbations with positive mean , Computers and Mathematics with Applications, Volume 66, Issue 11, 2013, Pages 2281-2294

    2013

  4. [16]

    Density Dependence in Rangelan d Grasshoppers (Orthoptera: Acrididae)

    Kemp, William P., and Brian Dennis. “Density Dependence in Rangelan d Grasshoppers (Orthoptera: Acrididae).” Oecologia 96, no. 1 (199 3): 1–8. http://www.jstor.org/stable/4220493

  5. [17]

    May, R. M. (1976). Simple mathematical models with very complicated dynamics . Nature, 261(5560), 459-467

    1976

  6. [18]

    Strogatz, S. H. (2018). Nonlinear Dynamics and Chaos: With Applications to Physics , Biology, Chemistry, and Engineering . CRC Press

    2018

  7. [19]

    Hilborn, R. C. (2000). Chaos and Nonlinear Dynamics: An Introduction for Scientis ts and Engineers . Oxford University Press

    2000

  8. [20]

    Braverman, A

    E. Braverman, A. Rodkina (2013). Difference equations of Ricker and logistic types under bounded stochastic perturbations with positive mean , Computers and Mathematics with Applications, Volume 66, Issue 11, 2013, Pages 2281-2294. 20

    2013

  9. [21]

    Random varia- tion in model parameters: A comprehensive review of stochas tic logistic growth equation , Ecological Modelling, Volume 484, 2023, 110475

    Md Aktar, Ul Karim, Vikram Aithal, Amiya Ranjan Bhowmick, (2023 ). Random varia- tion in model parameters: A comprehensive review of stochas tic logistic growth equation , Ecological Modelling, Volume 484, 2023, 110475

    2023

  10. [22]

    Cheke, Qianqian Zhang, San yi Tang, (2024)

    Dingding Yan, Mengqi He, Robert A. Cheke, Qianqian Zhang, San yi Tang, (2024). A stochastic hormesis Ricker model and its application to mul tiple fields . Chaos, Solitons & Fractals, Volume 185, 2024, 115131

    2024

  11. [23]

    Theoretical Population Biology, Volume 63, Issue 2, 2003, Pages 159-168,

    Ingemar N ˚ asell,(2003).Moment closure and the stochastic logistic model . Theoretical Population Biology, Volume 63, Issue 2, 2003, Pages 159-168,

    2003

  12. [24]

    E. C. Pielou (1975). Ecological Dirersity, Wiley, New York, 1975

    1975

  13. [25]

    Engen, (1978)

    S. Engen, (1978). Stochustic Abundunce Models , Chapman and Hall, London, 1978

    1978

  14. [27]

    Renshaw, E. (1991). Modelling Biological Populations in Space and Time . Cambridge University Press

    1991

  15. [28]

    R. M. May, (1974) Stability and Complexity in Model Ecosystems , Princeton U.P., Princeton

    1974

  16. [29]

    Dennis and G

    B. Dennis and G. P. Patil, The Gamma Distribution and Weighted Multimodal Gamma Distributions as Models of Population Abundance , Mathematical Biosciences, 68, 187–212, 1984

    1984

  17. [30]

    C. S. Peters, M. Mangel, R. F. Costantino, Stationary Distribution of Population Size in Tribolium, Bulletin of Mathematical Biology, 51(5), 625–638, 1989. 21

    1989

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.