IP Library Granted Patent US 12,724,939
Granted Patent B2
US 12,724,939 · App. 17/810,257 · Granted Sep 1, 2026

Representative simulation results for correlated variables

Inventors: Bradley Curtis Lackey (Redmond, WA); Andrew John McGuinness (Thornhill, GB); Yakoub Hassanov Yakoubov (Reigate, GB)
Assignee: Microsoft Technology Licensing, LLC
G06F30/20G06F17/18G06F2111/10
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,724,939
App. No.
17/810,257
Granted
Sep 1, 2026
Kind
B2
Abstract

A computing system comprises a processor configured to receive, for a plurality of correlated variables, a first predetermined number of simulations from a Monte-Carlo simulation sample, each simulation including a plurality of initial simulation results for the plurality of the variables. A unit interval of a cumulative distribution function (CDF) is segmented into a plurality of bins corresponding to a second predetermined number of strata. An initial discrepancy score is determined based upon a quantity of values in each bin, the first predetermined number, and second predetermined number. At least one of the initial simulation results is removed based upon an initial sum of the initial discrepancy scores. At least one other simulation is added and a plurality of representative simulations is output that represents the CDF across the strata based upon an updated sum of updated discrepancy scores.

Claims (75)

1 . A computing system, comprising:

a processor configured to

receive, for a plurality of correlated variables,

a first predetermined number of simulations from a Monte-Carlo simulation sample, each simulation including a plurality of initial simulation results for the plurality of the variables,

one or more target statistics,

a second predetermined number of strata for each variable and each target statistic, and a cumulative distribution function for each variable and each target statistic;

for each variable and for each of the one or more target statistics,

segment a unit interval of the cumulative distribution function into the second predetermined number of strata and a support of the cumulative distribution function into a plurality of bins such that each bin of the plurality of bins corresponds to one of the strata,

determine an initial discrepancy score based upon a quantity of values in each bin, the first predetermined number of the simulations, and the second predetermined number of the strata for the variable;

determine an initial sum of the initial discrepancy scores;

remove at least one of the plurality of the initial simulations based upon a determination that the initial sum of the initial discrepancy scores is not within an optimization threshold to decrease a sampling error;

add at least one other simulation to a remaining one or more initial simulations based on a policy to minimize the initial discrepancy score;

for each variable and for each of the one or more target statistics, use the quantity of the values in one or more bins corresponding to the at least one of the plurality of the initial simulations and the quantity of the values in one or more bins corresponding to the at least one other simulations to generate an updated discrepancy score;

determine an updated sum of the updated discrepancy scores; and

output a plurality of representative simulations that represent the cumulative distribution function across the strata based upon the updated sum of the updated discrepancy scores.

2 . The computing system of claim 1 , wherein the processor is further configured to accept or reject the at least one other simulation based upon the updated sum of the updated discrepancy scores.

3 . The computing system of claim 1 , wherein the processor is further configured to, for each variable and for each of the one or more target statistics, determine a stratum of a value for each initial simulation result and place the value into one of the plurality of bins based upon the stratum.

4 . The computing system of claim 1 , wherein the at least one other simulation includes a precomputed simulation.

5 . The computing system of claim 1 , wherein the processor is further configured to determine an initial bin-wise discrepancy metric for each bin of the plurality of bins.

6 . The computing system of claim 5 , wherein the initial bin-wise discrepancy metric for a selected bin includes a difference between the quantity of values in the selected bin and the first predetermined number of simulations divided by the second predetermined number of the strata for the variable or target statistic.

7 . The computing system of claim 5 , wherein the initial discrepancy score comprises a maximum bin-wise discrepancy metric or a sum of the initial bin-wise discrepancy metrics for the plurality of bins.

8 . The computing system of claim 5 , wherein the processor is further configured to weight the initial bin-wise discrepancy metric for each bin of the plurality of bins based upon a proximity of a stratum corresponding to the bin to a tail of the cumulative distribution function.

9 . The computing system of claim 5 , wherein the processor is further configured to:

determine an updated bin-wise discrepancy metric for the one or more bins corresponding to the at least one of the plurality of the initial simulations and for the one or more bins corresponding to the at least one other simulation; and

use the updated bin-wise discrepancy metric for the one or more bins corresponding to the at least one of the plurality of the initial simulations and for the one or more bins corresponding to the at least one other simulation, and the initial bin-wise discrepancy metric for each of the remaining one or more initial simulations, to determine the updated discrepancy score.

10 . The computing system of claim 1 , wherein the processor is further configured to, for each variable and target statistic:

decrement the quantity of values in the one or more bins corresponding to the at least one of the plurality of the initial simulations; and

increment the quantity of values in the one or more bins corresponding to the at least one other simulation.

11 . At a computing device, a method comprising:

receiving, for a plurality of correlated variables,

a first predetermined number of simulations from a Monte-Carlo simulation sample, each simulation including a plurality of initial simulation results for the plurality of the variables,

one or more target statistics,

a second predetermined number of strata for each variable and each target statistic, and

a cumulative distribution function for each variable and each target statistic;

for each variable and for each of the one or more target statistics,

segmenting a unit interval of the cumulative distribution function into the second predetermined number of the strata and a support of the cumulative distribution function into a plurality of bins such that each bin of the plurality of bins corresponds to one of the strata,

determining an initial discrepancy score based upon a quantity of values in each bin, the first predetermined number of the simulations, and the second predetermined number of the strata for the variable,

determining an initial sum of the initial discrepancy scores;

removing at least one of the plurality of the initial simulations based upon a determination that the initial sum of the initial discrepancy scores is not within an optimization threshold to decrease a sampling error;

adding at least one other simulation to a remaining one or more initial simulation based on a policy to minimize the initial discrepancy score;

for each variable and for each of the one or more target statistics, using the quantity of the values in one or more bins corresponding to the at least one of the plurality of the initial simulation results and the quantity of the values in one or more bins corresponding to the at least one other simulation result to generate an updated discrepancy score;

determining an updated sum of the updated discrepancy scores; and

outputting a plurality of representative simulations that represent the cumulative distribution function across the strata based upon the updated sum of the updated discrepancy scores.

12 . The method of claim 11 , further comprising accepting or rejecting the at least one other simulation based upon the updated sum of the updated discrepancy scores.

13 . The method of claim 11 , wherein the at least one other simulation includes a precomputed simulation.

14 . The method of claim 11 , wherein determining the initial discrepancy score includes determining an initial bin-wise discrepancy metric for each bin of the plurality of bins.

15 . The method of claim 14 , wherein determining the initial bin-wise discrepancy metric for a selected bin includes determining a difference between the quantity of values in the selected bin and the first predetermined number of the simulations divided by the second predetermined number of the strata for the variable or target statistic.

16 . The method of claim 14 , wherein determining the initial bin-wise discrepancy metric for the selected bin includes determining a maximum bin-wise discrepancy metric or determining a sum of the initial bin-wise discrepancy metrics for the plurality of bins.

17 . The method of claim 14 , wherein determining the initial discrepancy score includes weighting the initial bin-wise discrepancy metric for each bin of the plurality of bins based upon a proximity of a stratum corresponding to the bin to a tail of the cumulative distribution function.

18 . The method of claim 14 , wherein generating the updated discrepancy score includes:

determining an updated bin-wise discrepancy metric for the one or more bins corresponding to the at least one of the plurality of the initial simulation results and for the one or more bins corresponding to the at least one other simulation; and

using the updated bin-wise discrepancy metric for the one or more bins corresponding to the at least one of the plurality of the initial simulation results and for the one or more bins corresponding to the at least one other simulation result, and the initial bin-wise discrepancy metric for each of the remaining one or more initial simulation results, to determine the updated discrepancy score.

19 . The method of claim 11 , wherein generating the updated discrepancy score includes:

decrementing the quantity of values in the one or more bins corresponding to the at least one of the plurality of the initial simulations; and

incrementing the quantity of values in the one or more bins corresponding to the at least one other simulation.

20 . A computing system, comprising:

a processor configured to,

receive, for a plurality of correlated variables,

a first predetermined number of simulations from a Monte-Carlo simulation sample, each simulation including a plurality of initial simulation results for the plurality of the variables;

one or more target statistics,

a second predetermined number of strata for each variable and each target statistic; and

a cumulative distribution function for each variable and each target statistic;

for each variable and for each of the one or more target statistics,

segment a unit interval of the cumulative distribution function into the second predetermined number of strata and a support of the cumulative distribution function into a plurality of bins such that each bin of the plurality of bins corresponds to one of the strata,

count a quantity of values in each bin of the plurality of bins, and

determine an initial discrepancy score based upon a difference between the quantity of values in each bin and the quantity of the initial simulations divided by the second predetermined number of the strata for the variable or target statistic;

determine an initial sum of the initial discrepancy scores;

remove at least one of the plurality of the initial simulation results based upon a determination that the initial sum of the initial discrepancy scores is not within an optimization threshold to decrease a sampling error;

add at least one other simulation to a remaining one or more initial simulations based on a policy to minimize the initial discrepancy score;

for each variable and target statistic,

decrement the quantity of values in one or more bins corresponding to the at least one of the plurality of the initial simulation results,

increment the quantity of values in the one or more bins corresponding to the at least one other simulation result, and

use the quantity of the values in the one or more bins corresponding to the at least one of the plurality of the initial simulation results and the quantity of the values in the one or more bins corresponding to the at least one other simulation result to generate an updated discrepancy score;

determine an updated sum of the updated discrepancy scores; and

output a plurality of representative simulations that represent the cumulative distribution function across the strata based upon the updated sum of the updated discrepancy scores.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jul 15, 2022
From: LACKEY, BRADLEY CURTIS; MCGUINNESS, ANDREW JOHN; YAKOUBOV, YAKOUB HASSANOV
To: MICROSOFT TECHNOLOGY LICENSING, LLC
Reel/Frame 060520/0629 →
Continuity (1)
Related Publication 20240005057A1 · Jan 4, 2024
References Cited (32)
US 9805158B1 · Liu · 2017 [cited by applicant]
US 10628539B2 · Cox · 2020 [cited by applicant]
US 10831951B1 · Kossow · 2020 [cited by applicant]
US 10867249B1 · Bosnjakovic · 2020 [cited by applicant]
US 20080126024A1 · Gottsman · 2008 [cited by examiner]
US 20110054869A1 · Li · 2011 [cited by applicant]
US 20120173456A1 · Hirl · 2012 [cited by examiner]
US 20160283629A1 · Weckx · 2016 [cited by examiner]
US 20170109323A9 · Chen et al. · 2017 [cited by applicant]
US 20190377774A1 · Joshi · 2019 [cited by applicant]
US 20210049242A1 · Halle · 2021 [cited by examiner]
US 20210319158A1 · Bhattacharyya · 2021 [cited by applicant]
US 20240005056A1 · McGuinness · 2024 [cited by applicant]
US 20240005058A1 · McGuinness · 2024 [cited by applicant]
EP 2565830A1 · 2013 [cited by applicant]
WO 2017134045A1 · 2017 [cited by applicant]
Packham, et al., “Latin Hypercube Sampling with Dependence and Applications in Finance”, In Journal of SSRN, Oct. 1, 2008, 29 Pages. [cited by applicant]
“International Search Report and Written Opinion Issued in PCT Application No. PCT/US23/022624”, Mailed Date: Aug. 29, 2023, 10 Pages. [cited by applicant]
Scheingraber, et al., “Uncertainty Quantification for Seismic Risk Assessment Using Latin Hypercube Sampling and Quasi Monte Carlo Simulation”, In Proceedings of 16th European Conference on Earthquake Engineering, Jun. … [cited by applicant]
Gui, et al., “Fitting the Erlang Mixture Model to Data via a GEM-CMM algorithm”, In Journal of Computational and Applied Mathematics, vol. 343, Dec. 1, 2018, 26 Pages. [cited by applicant]
Arbenz, et al., “Copula based hierarchical risk aggregation through sample reordering”, In Journal of Insurance: Mathematics and Economics, vol. 51, Issue 1, Mar. 19, 2012, pp. 122-133. [cited by applicant]
Bon, et al., “Accelerating sequential Monte Carlo with surrogate likelihoods”, In Repository of arXiv:2009.03699v2, Jul. 20, 2021, 835 Pages. [cited by applicant]
Booth, et al., “Probabilistic analysis of a coal mine roadway including correlation control between model input parameters”, In Journal of Computers and Geotechnics, vol. 74, Feb. 4, 2016, pp. 151-162. [cited by applicant]
Cabri, Alberto, “Quantum inspired approach for early classification of time series”, In PhD Theses Series DIBRIS-TH-2019-21, Dec. 2019, 97 Pages. [cited by applicant]
Communication pursuant to Article 94(3) EPC received for EP Application No. 23741103.8, mailed on Nov. 26, 2025, 9 Pages. [cited by applicant]
Communication pursuant to Article 94(3) EPC received for EP Application No. 23741441.2, mailed on Nov. 26, 2025, 8 Pages. [cited by applicant]
International Search Report and Written Opinion received for PCT Application No. PCT/GB2023/051726, mailed on Sep. 21, 2023, 14 pages. [cited by applicant]
International Search Report and Written Opinion received for PCT Application No. PCT/GB2023/051736, mailed on Sep. 22, 2023, 15 pages. [cited by applicant]
Lee, et al., “Bayesian variable selection under the proportional hazards mixed-effects model”, In Journal of Computational Statistics and Data Analysis, vol. 75, Jul. 2014, pp. 53-65. [cited by applicant]
Non-Final Office Action mailed on Aug. 20, 2025, in U.S. Appl. No. 17/810,276, 12 pages. [cited by applicant]
Non-Final Office Action mailed on Oct. 1, 2025, in U.S. Appl. No. 17/810,251, 19 pages. [cited by applicant]
Notice of Allowance mailed on May 11, 2026, in U.S. Appl. No. 17/810,251, 09 pages. [cited by applicant]