在SciPy中创建自定义概率分布以从中绘制随机样本
我想使用蒙特卡罗模拟方法对任意数量的概率分布进行求和。我想要随机抽样连续分布，并将它们与其他连续分布的随机样本相加，最终得到它们组合的概率分布。这些分布本身是经验性的，它们不是函数，而是以P99 = 2.4、P90 = 7.12、P50 = 24.53、P10 = 82.14等形式存在（实际上...
C#蒙特卡罗增量风险计算优化，随机数，并行执行。
我的当前任务是优化一个蒙特卡罗模拟，该模拟计算一组债务人的区域资本充足率数据。 它目前运行速度太慢了，需要在生产中运行并进行每日多次运行。此外，结果数据的粒度将需要逐步提高到桌面甚至可能是书籍级别，我所得到的代码基本上是业务部门在半生产能力下使用的原型。 该应用程序目前是单线程的，因此我需...
如何使用NumPy/SciPy进行简单的高斯混合采样和概率密度函数绘制？
我将三个正态分布相加以得到一个新的分布，如下所示，请问如何在Python中按照这个分布进行采样？ import matplotlib.pyplot as plt import scipy.stats as ss import numpy as np x = np.linspace(0, ...
在Python上无法准确计算圆周率
我是新成员，我将直接着手处理它，因为我整个星期天都在尝试理解它。 我是Python的新手，先前学习了C ++编程的基本中级水平（这是一个10周的大学模块）。 我正在尝试几种迭代技术来计算圆周率，但两种方法都略有不准确，我不确定为什么。 第一种方法是我在大学里学到的-我相信你们中的一些人以...
MCTS UCT与评分系统
我正在尝试通过蒙特卡罗树搜索解决2048的变体。我发现UCT可以在探索/利用之间取得一些权衡。我的唯一问题是，我看到的所有版本都假定分数是胜率。如何使其适应一个分数是最后状态下板子价值的游戏，因此从1-MAX而不是胜利开始的游戏？ 我可以通过除以MAX使用常数c来归一化分数，但这会使游戏...
二叉树存储部分和：名称和现有实现
考虑一个由n个正实数（ai）组成的序列及其部分和序列（si）。给定一个属于(0, sn]的数字x，我们需要找到i使得si−1<x≤si。同时，我们希望能够更改其中一个ai而无需更新所有部分和。通过使用二叉树，将ai作为叶节点值，非叶节点的值为相应子节点值的总和，可以在O(log n)时间内完成...
如何加速蒙特卡罗模拟？
每当我在S-Plus中运行大规模蒙特卡罗模拟时，我总是要等待很长时间，期间我会留起胡须。 在R中运行蒙特卡罗模拟的最佳技巧是什么？有没有好的分布式处理进程的示例？
计算均值置信区间而不存储所有数据点
对于大的n值（查看下面的确定足够大的方法），根据中心极限定理，可以将样本均值分布视为正态分布（高斯分布），但我希望有一种程序可以针对任何n值提供置信区间。方法是使用自由度为n-1的学生T分布。 因此，问题在于，如何在收集或逐个遇到数据点的数据流中计算数据点的平均值的c（例如，c=.95）置信...
PRNG的质量测试
我正在尝试使用伪随机数生成器（如Mersenne Twister和stdlib库中的rand()函数），并希望找到一种好的测试方法，以确定PRNG所产生的随机数据的质量。我已经使用PRNG生成的随机数计算了圆周率的值，并发现rand()和Mersenne Twister非常接近，无法区分（需要...
计算函数可信度的算法 / 蒙特卡罗方法
我正在编写一个程序，试图复制本文开头讨论的算法。 http://www-stat.stanford.edu/~cgates/PERSI/papers/MCMCRev.pdf F是从字符到字符的函数。假设Pl（f）是该函数的“可信度”度量。算法如下： 从初步猜测函数f和一个新函数f*开始 ...