Pith. sign in

REVIEW 5 major objections 6 minor 1 cited by

Gamma Distribution for Equilibrium Analysis of Discrete Stochastic Logistic Population Models

T0 review · 5 major / 6 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 For a stochastic logistic model with a gamma-distributed equilibrium, the intrinsic growth rate is pinned to two branches, r+ and r−, lying within the same feasible range as the deterministic logistic map.

desk verdict A clean moment-matching calculation that is sold as a stationary-distribution result; the stationarity claim doesn't survive contact with the map's support, and the feasibility bound is off by an algebra step. read the letter →

arxiv 2411.10167 v1 pith:CU73S4UY submitted 2024-11-15 q-bio.PE

Haiyan Wang This is my paper · ORCID
classification q-bio.PE MSC 92D25
keywords stochasticlogisticmodelgammadistributionstationarymomentmatchingintrinsicgrowthratealternativestablestatesenvironmentalstochasticitydiscrete-timepopulation
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 derives explicit equilibrium conditions for a discrete stochastic logistic model in which the population distribution is assumed to be gamma. It shows that, when the first two moments are stationary, the intrinsic growth rate r must satisfy a quadratic whose two roots, r+ and r−, represent alternative stable states. The result yields a feasible range 1

What carries the argument

The argument's engine is the gamma distribution assumed for $X_t$, combined with the moment-stationarity conditions $E[X_{t+1}]=E[X_t]$ and $\mathrm{Var}[X_{t+1}]=\mathrm{Var}[X_t]$ (equation 11). The gamma's moment formula $E[X^n]=\theta^n\Gamma(k+n)/\Gamma(k)$ turns those two conditions into a closed system in $r$, $k$, and $\theta$. Eliminating $\theta$ produces a quadratic in $r$, whose discriminant gives the noise-variance ceiling $\mathrm{Var}(\epsilon)\le0.5$ and the shape-parameter bound $k\le1/\mathrm{Var}(\epsilon)-2$. The independence of $\epsilon_t$ from $X_t$ is what allows the noise to factor out of the expectations, so only the second moment of $\epsilon_t$ enters.

What would settle it

Simulate the recurrence $X_{t+1}=rX_t(1-X_t)\epsilon_t$ for a chosen $k$ and $\mathrm{Var}(\epsilon)$ with $r$ given by either root in Eq. (39) and $\theta=(r-1)/(r(1+k))$, using a large ensemble and many time steps; if the empirical distribution of $X_t$ does not converge to $\mathrm{Gamma}(k,\theta)$ or the moments drift, the derived equilibrium is not a real stationary distribution.

Watch

Extended reading notes

Core claim

The central claim is that for the stochastic logistic recurrence $X_{t+1}=rX_t(1-X_t)\epsilon_t$ with $E[\epsilon_t]=1$, an equilibrium population that follows a gamma distribution with shape $k$ and scale $\theta$ exists only for certain combinations of $r$, $k$, and the noise variance. Assuming the mean and variance are unchanged from $t$ to $t+1$, the paper derives $\theta=(r-1)/(r(1+k))$ and a quadratic equation for $r$. Its two roots are $r_{\pm}=(2k+4\pm(k+1)\sqrt{1-\mathrm{Var}(\epsilon)(k+2)/(\mathrm{Var}(\epsilon)+1)})/(k+3)$, with real solutions requiring $0\le\mathrm{Var}(\epsilon)\le0.5$ and $0<k\le1/\mathrm{Var}(\epsilon)-2$. The two roots $r_+>r_-$ are interpreted as higher- and lower-density alternative stable states, and all feasible $r$ lie in the deterministic logistic range $1<r<3$. Because $\theta$ cancels from the relation, the growth rate at equilibrium depends on the population's distribution shape and the environmental variance, not on its absolute scale.

Load-bearing premise

The entire derivation rests on treating equality of the first two moments from one time step to the next as the definition of equilibrium, without proving that the gamma distribution itself is invariant under the map; many distributions could share those moments without being stationary.

Editorial extensions

If this is right

  • If the equilibrium is gamma-distributed, the intrinsic growth rate $r$ is fixed by the shape parameter $k$ and the noise variance $\mathrm{Var}(\epsilon)$; the scale $\theta$ does not affect $r$.
  • Environmental noise with $\mathrm{Var}(\epsilon)>0.5$ cannot support this kind of gamma equilibrium, and higher noise shrinks the admissible range of $k$.
  • For fixed $k$ and $\mathrm{Var}(\epsilon)$, the two roots $r_+$ and $r_-$ give two distinct equilibria, implying the model can support alternative stable states at different population densities.
  • The derived feasible range $1<r<3$ matches the deterministic logistic map, so stochastic and deterministic versions share the same stability window.
  • The formula $\theta=(r-1)/(r(1+k))$ links the equilibrium mean $k\theta$ directly to $r$ and $k$, so the population size is fixed once $r$ and $k$ are chosen.

Reading between the lines

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

  • A natural test of the paper's assumption is to simulate the map with gamma-distributed initial conditions and check whether the distribution at later times remains gamma with the specified parameters; the moment conditions alone do not guarantee invariance of the full distribution.
  • Because only the second moment of $\epsilon_t$ enters the formulas, the same equilibrium relations should hold for any noise distribution with mean 1 and the given variance, such as lognormal or a scaled beta, though the actual stationary law would differ.
  • The two branches $r_+$ and $r_-$ suggest a hysteresis scenario: a population on the high-growth branch might drop to the low-growth branch after a disturbance, and the paper's formulas could be used to test which branch is more resilient to extinction.
  • Extending the moment-matching approach to a beta distribution, as the paper suggests for models with an explicit carrying capacity, would yield analogous conditions and might reveal how the support boundary changes the two-branch structure.
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

5 major / 6 minor

Summary. The paper studies the stochastic logistic recurrence X_{t+1}=rX_t(1-X_t)\epsilon_t with E[\epsilon_t]=1 and nonnegative perturbation. It assumes that, at equilibrium, X_t follows a gamma distribution with shape k and scale \theta (Eq. 4), and imposes equality of the first two moments of X_{t+1} and X_t (Eq. 11). From the mean condition it obtains \theta=(r-1)/(r(1+k)) (Eq. 17), and from the variance condition it derives a quadratic equation for r (Eq. 31). The roots are given in Eq. (39), and the discriminant is used to claim the parameter bounds 0<k\le 1/Var[\epsilon_t]-2 and 0\le Var[\epsilon_t]\le 0.5. The paper interprets the two roots r_+ and r_- as alternative stable states and argues that the feasible range of r is the same as for the deterministic logistic map, namely 1<r<3.

Significance. If the gamma law were shown to be an exact stationary distribution of (2), the explicit r-k-variance relation would be a compact and potentially useful result for stochastic logistic population models, and the analogy with the deterministic interval 1<r<3 would be biologically appealing. The algebraic derivation from Eq. (18) to Eq. (31) is internally consistent, and the mean condition leading to Eq. (17) is straightforward. However, the central claim is not established: the paper verifies only two moment equalities, not invariance of the full gamma distribution, and the gamma's positive support is incompatible with the fact that the map (2) produces negative values whenever X_t>1. No code or numerical simulation is provided; the figures plot Eq. (39) and therefore cannot validate the equilibrium claim independently. The paper is better read as a moment-calibration exercise than as a proof of stationary gamma equilibria.

major comments (5)
  1. [Section 2.1, Eq. (11)] The paper defines equilibrium by equality of the first two moments of X_{t+1} and X_t. For the nonlinear map (2), these equalities are necessary but not sufficient for stationarity of the full distribution. No argument shows that X_{t+1} has the same gamma distribution as X_t, and the paper itself notes immediately after Eq. (10) that X_{t+1} may not be gamma. Consequently, the roots in Eq. (39) are solutions of a moment-calibration equation, not a proof of the existence of a stationary gamma distribution.
  2. [Section 2.2, Eq. (17)] The support of the gamma distribution is incompatible with the map (2). With \theta=(r-1)/(r(1+k)) and r>1, one has \theta<1, so a gamma random variable X_t has positive probability of exceeding 1. On the event X_t>1, the quantity rX_t(1-X_t) is negative, and since \epsilon_t is assumed nonnegative, X_{t+1}<0 with positive probability. A gamma random variable is nonnegative, so X_{t+1} cannot be exactly gamma distributed. This is not a mild approximation issue; the supports are disjoint. The exact stationary gamma claim therefore fails for every r>1.
  3. [Section 2.4, Eqs. (32)-(33)] The feasibility bound is algebraically incorrect. From D=4(k+1)^2[1 - v(k+2)/(v+1)] with v=Var[\epsilon_t], the condition D\ge 0 is equivalent to v(k+1)\le 1, i.e. k\le 1/v - 1, not k\le 1/v - 2. Hence the claimed upper bound Var[\epsilon_t]\le 0.5 does not follow; the discriminant allows values such as v=0.6 with k\le 2/3. Since this bound is used in Section 2.4.2 to restrict the admissible range of k and in Figure 2, the reported parameter range must be corrected.
  4. [Section 2.4, Eq. (39)] The interpretation of r_+ and r_- as "alternative stable states" is unsupported. The two roots are merely the two solutions of a quadratic moment equation; the paper provides no stability analysis, no argument that trajectories converge to either branch, and no demonstration that the two branches correspond to distinct attracting equilibria. Moreover, r is not an externally specified growth parameter in this construction; it is a function of k and v. The statement that the feasible range matches the deterministic condition 1<r<3 is therefore not a dynamical equivalence but a property of Eq. (39).
  5. [Section 3, Conclusion and Discussion] The Conclusion explicitly states that "a more rigorous mathematical framework is needed to prove the relationships between the parameters of stochastic models and their stationary distributions" and suggests that beta distributions might be more suitable because they respect the carrying capacity. These statements are in tension with the Abstract's claims that the paper identifies the feasible range of r at equilibrium and establishes alternative stable states. The paper's own caveats thus acknowledge that the central stationarity result is not proven in the current manuscript.
minor comments (6)
  1. [Abstract and Section 1] The phrase "analysis the impact" should be "analyze the impact," and the repeated phrase "sufficient small" in Sections 2.4.1 and 2.4.2 should be "sufficiently small."
  2. [Section 2.1, Eqs. (9)-(10)] Equation (9) is followed by an empty numbered equation (10); the numbering should be cleaned up so that the display is not followed by a blank line.
  3. [Figure 2] Figure 2 is described as containing four panels, but the panels are not labeled and the axis labels and parameter values used in the plots are not specified, which makes the figure difficult to interpret.
  4. [References] Reference [9] is a Wikipedia article; a standard textbook or peer-reviewed reference for the gamma distribution would be more appropriate in a journal article.
  5. [Throughout] The notation V ar[\epsilon_t] appears with a space and inconsistent formatting; the right-hand side of Eq. (26) would be clearer with explicit parentheses around the product of the variance fraction, (k+1)^2, and 1/(r-1).
  6. [References [19]-[20]] The reference list contains apparent typos: "Dirersity" should be "Diversity" and "Stochustic Abundunce" should be "Stochastic Abundance."

Circularity Check

1 steps flagged · score 2.0 of 10

Moment-matching algebra is self-contained; only a minor self-referential validation of the gamma ansatz and an overstated equilibrium interpretation.

  1. other [Section 2.1, Eq. (11) and following sentence]
    "Therefore, at equilibrium, we assume that, for a specific time t, E[Xt+1] = E[Xt], V ar[Xt+1] = V ar[Xt] (11) ... These mathematical relations further confirm the gamma distribution’s flexibility in accommodating skewed distributions and its ability to model variability around an equilibrium make it particularly suitable for representing the outcomes of stochastic processes in ecology."

    The gamma law is posited at Eq. (4), and 'equilibrium' is defined by the moment conditions (11). Equations (17) and (39) are algebraic consequences of that posited form plus those moment conditions, so they cannot independently confirm that the gamma distribution is suitable: the output of the assumption is being used as evidence for the assumption. The later description of the two roots r+ and r− as 'alternative stable states' also imports stability language without an invariance or stability analysis, so that conclusion is no stronger than the input moment conditions.

full rationale

The central derivation is not circular in the fitting or self-citation sense. Starting from the explicit gamma assumption (Eq. 4), the independence and unit-mean assumptions on ε, and the two moment equalities (Eq. 11), the paper obtains θ=(r−1)/(r(1+k)) (Eq. 17), a quadratic in r (Eq. 31), and the roots r± (Eq. 39) by direct algebra. These results are internally derived, not fitted to a subset of data, and no load-bearing claim rests on a self-citation. The comparison with the deterministic range 1<r<3 is an external consistency statement, not an input. However, the paper overstates the status of these formulas: equality of the first two moments is necessary, not sufficient, for a stationary distribution, and the paper does not prove that a gamma law is invariant under (2); Section 3 concedes that a rigorous framework for stationary distributions is future work. There is one rhetorical circularity: the text says the derived relations 'further confirm' the gamma distribution's suitability, but those relations were generated from the gamma assumption, so they cannot confirm it. This is minor and not load-bearing for the algebraic results. Separately, the feasibility claim Var(ε)≤0.5 rests on a discriminant simplification error: D≥0 gives k≤1/Var(ε)−1, not k≤1/Var(ε)−2, so that cap is a correctness issue rather than a circularity.

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

No data are fitted; k and Var(ε) are symbolic inputs. The core of the paper is two imposed assumptions: the gamma law and moment-based stationarity. Both are stated but not proven. The paper introduces no new particles, forces, or entities.

assumptions (4)
  • ad hoc to paper At equilibrium, X_t follows a gamma distribution with shape k and scale θ (Eq. 4).
    This is imposed, not derived; the nonlinear map does not preserve the gamma family, and the paper provides no invariance proof.
  • domain assumption At equilibrium, it is sufficient to require E[X_{t+1}]=E[X_t] and Var[X_{t+1}]=Var[X_t] (Eq. 11).
    These are necessary conditions for stationarity but not sufficient for a nonlinear map; no proof is given that X_{t+1} has the same distribution as X_t.
  • domain assumption The perturbation ε_t has mean 1, is nonnegative, and is independent of X_t.
    This is the model definition in Eq. (2); the independence is used to factor expectations in Eq. (12).
  • ad hoc to paper Population values X_t can be treated as positive real numbers with gamma moments, despite the logistic map's (0,1) domain.
    Gamma support extends beyond 1, so the moment formulas assume positivity but not the upper bound required by the map; the paper does not restrict the state space.

how reviews work

0 comments
Cite this review

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

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

Stochastic models play an essential role in accounting for the variability and unpredictability seen in real-world. This paper focuses on the application of the gamma distribution to analysis of the stationary distributions of populations governed by the discrete stochastic logistic equation at equilibrium. It is well known that the population dynamics of deterministic logistic models are dependent on the range of intrinsic growth rate. In this paper, we identify the same feasible range of the intrinsic growth rate for the stochastic model at equilibrium and establish explicit mathematical relation among the parameters of the gamma distribution and the stochastic models. We analyze the biological implications of these relationships, with particular emphasis on how the shape and scale parameters of the gamma distribution reflect population dynamics at equilibrium. These mathematical relations describe the impact of the variance of the stochastic perturbation on the intrinsic growth rate, and, in particular, reveal that there are two branches of the intrinsic growth rates representing alternative stable states at equilibrium.

Figures

Figures reproduced from arXiv: 2411.10167 by the authors.

Figure 1
Figure 1. Plot of r in terms of k It is well known that the population dynamics of the deterministic logistic model (1), such as stable growth, periodic oscillations, chaos, or extinction, are dependent on the value of r. If 1 < r < 3, the population grows and eventually reaches a nonzero steady state[13, 14, 15]. 10 [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Plots of r in term of k The condition 0 < k ≤ 1 Var[ϵt] − 2 suggests that the feasible range of k values is restricted by the environmental variance Var(ϵt). Smaller values of Var(ϵt) extend the range of k, which implies that populations with larger k (more symmetrical and regular population distribution) are better supported in stable environments with low variability. As we can see from [PITH_FULL_IMAGE:figures/f… view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Equilibrium Analysis of Discrete Stochastic Population Models with Gamma Distribution

    q-bio.PE 2024-11 conditional novelty 5.0 of 10

    Assuming a gamma stationary distribution, stochastic logistic and Ricker models yield explicit moment-consistency formulas with two intrinsic-growth-rate branches.

Reference graph

Works this paper leans on

21 extracted references · 7 canonical work pages · cited by 1 Pith paper

  1. [1]

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

    2010

  2. [2]

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

    1991

  3. [3]

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

    work page 1974

  4. [4]

    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

  5. [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

  6. [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

  7. [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

    1981

  8. [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 the stochastic logistic model of population growth, Ecological Modelling, 161 (2003) 239–248

    2003

Show all 21 references
  1. [9]

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

    2024

  2. [10]

    Kot, Elements of Mathematical Ecology

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

    2001

  3. [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/s12080-015-0288-2

    2016 doi

  4. [12]

    Density Dependence in Rangeland Grasshoppers (Orthoptera: Acrididae)

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

    1993

  5. [13]

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

    1976

  6. [14]

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

    2018

  7. [15]

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

    2000

  8. [16]

    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

  9. [17]

    Random varia- tion in model parameters: A comprehensive review of stochastic 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 stochastic logistic growth equation , Ecological Modelling, Volume 484, 2023, 110475

    2023

  10. [18]

    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

  11. [19]

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

    1975

  12. [20]

    Engen, (1978)

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

    1978

  13. [21]

    and Huang, Shuo and Jiang, Jifa and Wang, Hao (2021)

    Schreiber, Sebastian J. and Huang, Shuo and Jiang, Jifa and Wang, Hao (2021). Ex- tinction and Quasi-Stationarity for Discrete-Time, Endemic SIS and SIR Models , SIAM Journal on Applied Mathematics, 81(5) 2195-2217. 15

    2021

Pith tools

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