pith. sign in

arxiv: 1907.07317 · v1 · pith:POJWCQZ5new · submitted 2019-07-17 · 🧮 math.OC

Regularized two-stage stochastic variational inequalities for Cournot-Nash equilibrium under uncertainty

Pith reviewed 2026-05-24 20:41 UTC · model grok-4.3

classification 🧮 math.OC
keywords two-stage stochastic variational inequalityCournot-Nash equilibriumsample average approximationregularizationoligopolistic marketproduction planninguncertainty
0
0 comments X

The pith

A regularized sample average approximation solves two-stage stochastic variational inequalities for uncertain Cournot-Nash games.

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

The paper formulates a convex two-stage non-cooperative multi-agent game under uncertainty as a two-stage stochastic variational inequality. It supplies sufficient conditions for the existence of solutions to this SVI and introduces a regularized sample average approximation method to compute them. Convergence of the method is established as the regularization parameter tends to zero and the sample size tends to infinity. The formulation is applied to a two-stage stochastic production and supply planning problem in an oligopolistic market, with numerical tests on crude oil data showing how it captures agents' market shares.

Core claim

Under standard assumptions, the two-stage SVI has solutions, and the regularized sample average approximation method converges to those solutions as the regularization parameter tends to zero and the sample size tends to infinity. The approach is applied to a two-stage stochastic production and supply planning problem with homogeneous commodity in an oligopolistic market.

What carries the argument

The two-stage stochastic variational inequality (SVI) encoding agents' first- and second-stage decisions under uncertainty, solved via the regularized sample average approximation method.

Load-bearing premise

The game is convex and satisfies the standard assumptions required for existence of solutions to the two-stage SVI.

What would settle it

If the solutions produced by the regularized sample average approximation method do not approach a solution of the original two-stage SVI as the regularization parameter tends to zero and the sample size tends to infinity, the convergence claim fails.

Figures

Figures reproduced from arXiv: 1907.07317 by Jie Jiang, Xiaojun Chen, Xiaozhou Wang, Yun Shi.

Figure 5
Figure 5. Figure 5: illustrates the performance of the PHM measured by the number of iterations [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 5.1
Figure 5.1. Figure 5.1: Numerical comparisons among different , J = 10. actual market shares. Based on historical data on crude oil market, we make in-sample back￾tracking test to establish the effectiveness and validity of our model while explaining the market behaviour. Furthermore, the out-of-sample prediction capability of our model is demonstrated when training data is used to specify model parameters. From the results of… view at source ↗
Figure 5.2
Figure 5.2. Figure 5.2: Convergence property of x with increasing ν, J = 10. dramatic changes in oil price. One interesting observation is that the market share behaves rather smoothly even during periods of oil shocks [25]. Majority of the world’s crude oil is sup￾plied by a few large oil exporting countries and they are viewed collectively as a finite number [PITH_FULL_IMAGE:figures/full_fig_p020_5_2.png] view at source ↗
Figure 5.3
Figure 5.3. Figure 5.3: Market shares of different oil-producing countries, 1965-2017. [PITH_FULL_IMAGE:figures/full_fig_p021_5_3.png] view at source ↗
Figure 5.4
Figure 5.4. Figure 5.4: Real, in-sample and out-of-sample market shares results, 2008-2017. [PITH_FULL_IMAGE:figures/full_fig_p023_5_4.png] view at source ↗
read the original abstract

A convex two-stage non-cooperative multi-agent game under uncertainty is formulated as a two-stage stochastic variational inequality (SVI). Under standard assumptions, we provide sufficient conditions for the existence of solutions of the two-stage SVI and propose a regularized sample average approximation method for solving it. We prove the convergence of the method as the regularization parameter tends to zero and the sample size tends to infinity. Moreover, our approach is applied to a two-stage stochastic production and supply planning problem with homogeneous commodity in an oligopolistic market. Numerical results based on historical data in crude oil market are presented to demonstrate the effectiveness of the two-stage SVI in describing the market share of oil producing agents.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 1 minor

Summary. The paper formulates a convex two-stage non-cooperative multi-agent game under uncertainty as a two-stage stochastic variational inequality (SVI). Under standard assumptions, it provides sufficient conditions for the existence of solutions to the two-stage SVI and proposes a regularized sample average approximation (SAA) method for solving it. Convergence of the method is proved as the regularization parameter tends to zero and the sample size tends to infinity. The approach is applied to a two-stage stochastic production and supply planning problem with homogeneous commodity in an oligopolistic market, with numerical results from historical crude oil market data demonstrating effectiveness in describing market shares.

Significance. If the existence conditions and convergence proofs are rigorously established without hidden gaps, the work supplies a theoretically grounded method for computing Cournot-Nash equilibria in uncertain multi-stage settings. The regularized SAA scheme addresses both theoretical solvability and computational tractability, and the oil-market application supplies a concrete, data-driven illustration. No machine-checked proofs or parameter-free derivations are claimed, but the explicit convergence statement under standard assumptions is a clear strength when the derivations are supplied in full.

minor comments (1)
  1. The numerical experiments would be strengthened by including error bars, sensitivity plots with respect to the regularization parameter, and explicit checks that the production-planning instance satisfies the monotonicity/convexity conditions invoked for existence and convergence.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive evaluation of the manuscript and the recommendation for minor revision. The referee's summary accurately describes the formulation of the two-stage SVI, the existence conditions, the regularized SAA convergence analysis, and the crude oil market application.

Circularity Check

0 steps flagged

No significant circularity detected; derivation is self-contained

full rationale

The paper formulates a convex two-stage game as an SVI, states sufficient conditions for existence under standard assumptions, and proves convergence of a regularized SAA scheme as the regularization parameter tends to zero and sample size tends to infinity. These steps are presented as mathematical results (existence via assumptions, convergence via analysis of the regularized approximation), with no equations or claims reducing by construction to fitted parameters, self-definitions, or load-bearing self-citations. The numerical application to the crude oil market is presented as demonstration rather than a central derivation. No patterns of self-definitional claims, fitted inputs renamed as predictions, or ansatzes smuggled via prior self-work appear in the provided claims. This is the expected outcome for a standard existence-and-convergence analysis in stochastic variational inequalities.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract, the central claims rest on unspecified standard assumptions for SVI existence and on the introduction of a regularization parameter that is driven to zero.

free parameters (1)
  • regularization parameter
    The method introduces a regularization parameter that is taken to zero in the convergence statement.
axioms (1)
  • domain assumption Standard assumptions guaranteeing existence of solutions to the two-stage SVI
    Invoked to establish existence and to support the convergence proof.

pith-pipeline@v0.9.0 · 5647 in / 1175 out tokens · 19837 ms · 2026-05-24T20:41:17.116022+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages

  1. [1]

    Bornstein, P

    G. Bornstein, P. Krusell and S. Rebelo, Lags, costs, and shocks: An equilibrium model of the oil industry, National Bureau of Economic Research, 2017

  2. [2]

    Burke, X

    J.V. Burke, X. Chen and H. Sun, Subdifferentiation and smoothing of nonsmooth integral func- tionals, Department of applied mathematics, The Hong Kong Polytechnic University, 2017

  3. [3]

    Chen, Smoothing methods for nonsmooth, nonconvex minimization, Math

    X. Chen, Smoothing methods for nonsmooth, nonconvex minimization, Math. Program., 134 (2012), 71-99

  4. [4]

    Chen and M

    X. Chen and M. Fukushima, A smoothing method for a mathematical program with P-matrix linear complementarity constraints, Comput. Optim. Appl., 27 (2004), 223-246

  5. [5]

    Chen and M

    X. Chen and M. Fukushima, Expected residual minimization method for stochastic linear comple- mentarity problems, Math. Oper. Res., 30 (2005), 1022-1038

  6. [6]

    Chen, T.K

    X. Chen, T.K. Pong and R.J.B. Wets, Two-stage stochastic variational inequalities: an ERM- solution procedure, Math. Program., 165 (2017), 71-111

  7. [7]

    X. Chen, H. Sun and R.J.B. Wets, Regularized mathematical programs with stochastic equilibrium constraints: estimating structural demand models, SIAM J. Optim., 25 (2015), 53-75

  8. [8]

    X. Chen, A. Shapiro and H. Sun, Convergence analysis of sample average approximation of two-state stochastic generalized equations, SIAM J. Optim., 29 (2019), 135-161

  9. [9]

    X. Chen, H. Sun and H. Xu, Discrete approximation of two-stage stochastic and distributionally robust linear complementarity problems, Math. Program., (2018), to appear. Two-stage stochastic variational inequalities for Nash equilibrium 25

  10. [10]

    Chen, R.J.B

    X. Chen, R.J.B. Wets and Y. Zhang, Stochastic variational inequalities: residual minimization smoothing sample average approximations, SIAM J. Optim., 22 (2012), 649-673

  11. [11]

    Chen and S

    X. Chen and S. Xiang, Newton iterations in implicit time-stepping scheme for differential linear complementarity systems, Math. Program., 138 (2013), 579-606

  12. [12]

    Chen and Y

    X. Chen and Y. Ye, On homotopy-smoothing methods for box-constrained variational inequalities, SIAM J. Control. Optim., 37 (1999), 589-616

  13. [13]

    Chen and Y

    X. Chen and Y. Ye, On smoothing methods for the P0 matrix linear complementarity problem, SIAM J. Optim., 11 (2000), 341-363

  14. [14]

    Facchinei and J.S

    F. Facchinei and J.S. Pang, Finite-dimensional variational inequalities and complementarity prob- lems, Springer Science & Business Media, 2007

  15. [15]

    Gabriel and J.J

    S.A. Gabriel and J.J. More, Smoothing of mixed complementarity problems, in Complementarity and Variational Problems: State of the Art, M.C. Ferris and J.S. Pang, eds., SIAM, (1997), 105-116

  16. [16]

    Genc, S.S

    T.S. Genc, S.S. Reynolds and S. Sen, Dynamic oligopolistic games under uncertainty: A stochastic programming approach, J. Econ. Dyn. Control., 31 (2007), 55-80

  17. [17]

    Hamilton, This is what happened to the oil price-macroeconomy relationship, J

    J.D. Hamilton, This is what happened to the oil price-macroeconomy relationship, J. Monet. Econ., 38 (1996), 215-220

  18. [18]

    Hamilton, What is an oil shock?, J

    J.D. Hamilton, What is an oil shock?, J. Econom., 113:2 (2003), 363-398

  19. [19]

    Hamilton, Historical oil shocks, National Bureau of Economic Research, 2011

    J.D. Hamilton, Historical oil shocks, National Bureau of Economic Research, 2011

  20. [20]

    Harker and J.S

    P.T. Harker and J.S. Pang, Finite-dimensional variational inequality and nonlinear complementar- ity problems: a survey of theory, algorithms and applications, Math. Program., 48 (1990), 161-200

  21. [21]

    Hobbs and J.S

    B.F. Hobbs and J.S. Pang, Nash-Cournot equilibria in electric power markets with piecewise linear demand functions and joint constraints, Oper. Res., 55 (2007), 113-127

  22. [22]

    Hotelling, The economics of exhaustible resources, J

    H. Hotelling, The economics of exhaustible resources, J. Polit. Econ., 39 (1931), 137-175

  23. [23]

    Huang, M

    J. Huang, M. Leng and M. Parlar, Demand functions in decision modeling: A comprehensive survey and research directions, Decis. Sci., 44 (2013), 557-609

  24. [24]

    Jofr´ e, R.T

    A. Jofr´ e, R.T. Rockafellar and R.J.B. Wets, Variational inequalities and economic equilibrium, Math. Oper. Res., 32 (2007), 32-50

  25. [25]

    Kilian, Not all oil price shocks are alike: Disentangling demand and supply shocks in the crude oil market, Am

    L. Kilian, Not all oil price shocks are alike: Disentangling demand and supply shocks in the crude oil market, Am. Econ. Rev., 99 (2009), 1053-1069

  26. [26]

    Kilian and D.P

    L. Kilian and D.P. Murphy, The role of inventories and speculative trading in the global market for crude oil, J. Appl. Econom., 29 (2014), 454-478

  27. [27]

    Murphy, H.D

    F.H. Murphy, H.D. Sherali and A.L. Soyster, A mathematical programming approach for deter- mining oligopolistic market equilibrium, Math. Program., 24 (1982), 92-105

  28. [28]

    J.S. Pang, S. Sen and U.V. Shanbhag, Two-stage non-cooperative games with risk-averse players, Math. Program., 165 (2017), 235-290

  29. [29]

    Pang, C.L

    J.S. Pang, C.L. Su, Y.C. Lee, A constructive approach to estimating pure characteristics demand models with pricing, Oper. Res., 63 (2015), 639-659

  30. [30]

    Philpott, M.C

    A. Philpott, M.C. Ferris and R.J.B. Wets, Equilibrium, uncertainty and risk in hydro-thermal electricity systems, Math. Program., 157 (2016), 483-513

  31. [31]

    Outrata, M.C

    J. Outrata, M.C. Ferris, M. ˇCervinka and M. Outrata, On Cournot-Nash-Walras equilibria and their computation, Set-Valued Var. Anal., 24 (2016), 387-402

  32. [32]

    Richard, J.S

    C.W. Richard, J.S. Pang and S.E. Richard, The linear complementarity problem, Academic Press, 1992

  33. [33]

    Rockafellar and J

    R.T. Rockafellar and J. Sun, Solving monotone stochastic variational inequalities and complemen- tarity problems by progressive hedging, Math. Program., 174 (2019), 453-471

  34. [34]

    Rockafellar and R.J.B

    R.T. Rockafellar and R.J.B. Wets, Stochastic convex programming: Kuhn-Tucker conditions, J. Math. Econ., 2 (1975), 349-370

  35. [35]

    Rockafellar and R.J.B

    R.T. Rockafellar and R.J.B. Wets, Stochastic convex programming: Basic duality, Pac. J. Math., 62 (1976), 173-195

  36. [36]

    Rockafellar and R.J.B

    R.T. Rockafellar and R.J.B. Wets, Stochastic convex programming: Singular multipliers and ex- 26 J. JIANG, Y. SHI, X. WANG AND X. CHEN tended duality singular multipliers and duality, Pac. J. Math., 62 (1976), 173-195

  37. [37]

    Rockafellar and R.J.B

    R.T. Rockafellar and R.J.B. Wets, Stochastic convex programming: Relatively complete recourse and induced feasibility, SIAM. J. Control. Optim., 14 (1976), 574-589

  38. [38]

    Rockafellar and R.J.B

    R.T. Rockafellar and R.J.B. Wets, Scenarios and policy aggregation in optimization under uncer- tainty, Math. Oper. Res., 16 (1991), 119-147

  39. [39]

    Rockafellar and R.J.B

    R.T. Rockafellar and R.J.B. Wets, Variational analysis, Springer Science & Business Media, 2009

  40. [40]

    Rockafellar and R.J.B

    R.T. Rockafellar and R.J.B. Wets, Stochastic variational inequalities: single-stage to multistage, Math. Program., 165, (2017), 331-360

  41. [41]

    Rosen, Existence and uniqueness of equilibrium points for concave n-person games, Econo- metrica, (1965), 520-534

    B.J. Rosen, Existence and uniqueness of equilibrium points for concave n-person games, Econo- metrica, (1965), 520-534

  42. [42]

    Salant, Exhaustible resources and industrial structure: A Nash-Cournot approach to the world oil market, J

    S.W. Salant, Exhaustible resources and industrial structure: A Nash-Cournot approach to the world oil market, J. Polit. Econ., 84 (1976), 1079-1093

  43. [43]

    Shapiro, D

    A. Shapiro, D. Dentcheva and A. Ruszczy´ nski,Lectures on stochastic programming: modeling and theory, SIAM, 2009

  44. [44]

    Shapiro, Monte Carlo sampling methods, Handbooks in operations research and management science, 10 (2003), 353-425

    A. Shapiro, Monte Carlo sampling methods, Handbooks in operations research and management science, 10 (2003), 353-425

  45. [45]

    Smith, World oil: market or mayhem?, J

    J.L. Smith, World oil: market or mayhem?, J. Econ. Perspect., 23 (2009), 145-164

  46. [46]

    Sun, C.L

    H. Sun, C.L. Su and X. Chen, SAA-regularized methods for multiproduct price optimization under the pure characteristics demand model, Math. Program., 165 (2017), 361-389

  47. [47]

    H. Xu, J.S. Pang, F. Ord´ onez and M. Dessouky, Complementarity models for traffic equilibrium with ridesharing, Transport. Res. B-meth., 81 (2015), 161-182

  48. [48]

    J. Yao, I. Adler and S.S. Oren, Modeling and computing two-settlement oligopolistic equilibrium in a congested electricity network, Oper. Res., 56 (2008), 34-47

  49. [49]

    Juvenal and I

    L. Juvenal and I. Petrella, Speculation in the oil market, J. Appl. Econom., 30 (2015), 621-649. A. Appendix Proof of Theorem 3.1. Proof. By direct computation, we have that ( Mϵ(ξ) (yϵ λϵ ) +q(x,ξ ) ) j = { γ(ξ)(yϵ)j +γ(ξ)Tϵ + (λϵ)j−pj(ξ), j = 1,...,J ; xj−J− (yϵ)j−J +ϵ(λϵ)j−J, j =J + 1,..., 2J. Then, we can rewrite problem (3.2) as below: { 0≤ (yϵ)j⊥γ(ξ...