首页 > 最新文献

Monte Carlo Methods and Applications最新文献

英文 中文
Monte Carlo method for parabolic equations involving fractional Laplacian 包含分数阶拉普拉斯式的抛物方程的蒙特卡罗方法
IF 0.9 Q3 STATISTICS & PROBABILITY Pub Date : 2022-10-27 DOI: 10.1515/mcma-2022-2129
Caiyu Jiao, Changpin Li
Abstract We apply the Monte Carlo method to solving the Dirichlet problem of linear parabolic equations with fractional Laplacian. This method exploits the idea of weak approximation of related stochastic differential equations driven by the symmetric stable Lévy process with jumps. We utilize the jump-adapted scheme to approximate Lévy process which gives exact exit time to the boundary. When the solution has low regularity, we establish a numerical scheme by removing the small jumps of the Lévy process and then show the convergence order. When the solution has higher regularity, we build up a higher-order numerical scheme by replacing small jumps with a simple process and then display the higher convergence order. Finally, numerical experiments including ten- and one hundred-dimensional cases are presented, which confirm the theoretical estimates and show the numerical efficiency of the proposed schemes for high-dimensional parabolic equations.
摘要应用蒙特卡罗方法求解了具有分数阶拉普拉斯式的线性抛物型方程的Dirichlet问题。该方法利用了由具有跳跃的对称稳定lsamvy过程驱动的相关随机微分方程的弱逼近思想。我们利用跳跃适应方案来近似lsamvy过程,给出了精确的边界退出时间。当解的正则性较低时,通过去掉lsamvy过程的小跳变,建立了数值格式,并给出了收敛阶。当解具有较高的正则性时,我们用一个简单的过程代替小的跳跃,建立一个高阶的数值格式,然后显示更高的收敛阶。最后,给出了十维和一百维情况下的数值实验,验证了理论估计,并证明了所提格式对高维抛物方程的数值效率。
{"title":"Monte Carlo method for parabolic equations involving fractional Laplacian","authors":"Caiyu Jiao, Changpin Li","doi":"10.1515/mcma-2022-2129","DOIUrl":"https://doi.org/10.1515/mcma-2022-2129","url":null,"abstract":"Abstract We apply the Monte Carlo method to solving the Dirichlet problem of linear parabolic equations with fractional Laplacian. This method exploits the idea of weak approximation of related stochastic differential equations driven by the symmetric stable Lévy process with jumps. We utilize the jump-adapted scheme to approximate Lévy process which gives exact exit time to the boundary. When the solution has low regularity, we establish a numerical scheme by removing the small jumps of the Lévy process and then show the convergence order. When the solution has higher regularity, we build up a higher-order numerical scheme by replacing small jumps with a simple process and then display the higher convergence order. Finally, numerical experiments including ten- and one hundred-dimensional cases are presented, which confirm the theoretical estimates and show the numerical efficiency of the proposed schemes for high-dimensional parabolic equations.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"29 1","pages":"33 - 53"},"PeriodicalIF":0.9,"publicationDate":"2022-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44779092","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A random walk algorithm to estimate a lower bound of the star discrepancy 一种估计星差下界的随机游动算法
IF 0.9 Q3 STATISTICS & PROBABILITY Pub Date : 2022-10-21 DOI: 10.1515/mcma-2022-2125
Maryam Alsolami, M. Mascagni
Abstract In many Monte Carlo applications, one can substitute the use of pseudorandom numbers with quasirandom numbers and achieve improved convergence. This is because quasirandom numbers are more uniform that pseudorandom numbers. The most common measure of that uniformity is the star discrepancy. Moreover, the main error bound in quasi-Monte Carlo methods, called the Koksma–Hlawka inequality, has the star discrepancy in the formulation. A difficulty with this bound is that computing the star discrepancy is very costly. The star discrepancy can be computed by evaluating a function called the local discrepancy at a number of points. The supremum of these local discrepancy values is the star discrepancy. If we have a point set in [ 0 , 1 ] s {[0,1]^{s}} with N members, we need to compute the local discrepancy at N s {N^{s}} points. In fact, computing star discrepancy is NP-hard. In this paper, we will consider an approximate algorithm for a lower bound on the star discrepancy based on using a random walk through some of the N s {N^{s}} points. This approximation is much less expensive that computing the star discrepancy, but still accurate enough to provide information on convergence. Our numerical results show that the random walk algorithm has the same convergence rate as the Monte Carlo method, which is O ( N - 1 2 {O(N^{-frac{1}{2}}} ).
摘要在许多蒙特卡罗应用中,可以用拟随机数代替伪随机数的使用,从而达到提高收敛性的目的。这是因为准随机数比伪随机数更均匀。这种均匀性最常见的测量方法是恒星差异。此外,拟蒙特卡罗方法的主要误差界,称为Koksma-Hlawka不等式,在公式中具有星形差异。这个界限的一个困难是,计算恒星差异的成本非常高。星形差异可以通过在若干点上计算一个称为局部差异的函数来计算。这些局部差值的最大值是星形差值。如果我们有一个在[0,1]s {[0,1]^{s}}中有N个成员的点集,我们需要计算N s {N^{s}}点上的局部差异。事实上,计算恒星差异是np困难的。在本文中,我们将考虑一种基于随机遍历一些N s {N^{s}}点的星差下界的近似算法。这种近似方法比计算恒星差异要便宜得多,但仍然足够精确,可以提供关于收敛的信息。我们的数值结果表明,随机漫步算法具有与蒙特卡罗方法相同的收敛速度,即O(N - 1 2 {O(N^{-frac{1}{2}}})。
{"title":"A random walk algorithm to estimate a lower bound of the star discrepancy","authors":"Maryam Alsolami, M. Mascagni","doi":"10.1515/mcma-2022-2125","DOIUrl":"https://doi.org/10.1515/mcma-2022-2125","url":null,"abstract":"Abstract In many Monte Carlo applications, one can substitute the use of pseudorandom numbers with quasirandom numbers and achieve improved convergence. This is because quasirandom numbers are more uniform that pseudorandom numbers. The most common measure of that uniformity is the star discrepancy. Moreover, the main error bound in quasi-Monte Carlo methods, called the Koksma–Hlawka inequality, has the star discrepancy in the formulation. A difficulty with this bound is that computing the star discrepancy is very costly. The star discrepancy can be computed by evaluating a function called the local discrepancy at a number of points. The supremum of these local discrepancy values is the star discrepancy. If we have a point set in [ 0 , 1 ] s {[0,1]^{s}} with N members, we need to compute the local discrepancy at N s {N^{s}} points. In fact, computing star discrepancy is NP-hard. In this paper, we will consider an approximate algorithm for a lower bound on the star discrepancy based on using a random walk through some of the N s {N^{s}} points. This approximation is much less expensive that computing the star discrepancy, but still accurate enough to provide information on convergence. Our numerical results show that the random walk algorithm has the same convergence rate as the Monte Carlo method, which is O ( N - 1 2 {O(N^{-frac{1}{2}}} ).","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"28 1","pages":"341 - 348"},"PeriodicalIF":0.9,"publicationDate":"2022-10-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41818243","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
Superposition of forward and backward motion 前后运动叠加
IF 0.9 Q3 STATISTICS & PROBABILITY Pub Date : 2022-10-19 DOI: 10.1515/mcma-2022-2124
Manfred Harringer
Abstract I consider black body radiation. The wall of the black body exchanges photons with the radiation field in equilibrium, therefore with a common temperature in Planck’s radiation law. The underlying process of radiation consists of creation and annihilation of photons. I want to present an alternate model of motions, where the process of radiation consists of small steps in positive and negative direction, not zero in mean. The detection of radiation consists of storing and restoring of packages of energy. I get an analogue of Planck’s radiation law, where the common temperature emerges from the underlying common model of small steps. The object of the law is not the radiation, but a storage of packages of energy, which belongs to the wall of the black body.
摘要我认为黑体辐射。黑体壁在平衡状态下与辐射场交换光子,因此具有普朗克辐射定律中的共同温度。辐射的基本过程包括光子的产生和湮灭。我想提出一个交替的运动模型,其中辐射过程由正负方向的小步组成,而不是平均值为零。辐射探测包括能量包的储存和恢复。我得到了普朗克辐射定律的类似物,其中共同的温度来自于小台阶的基本共同模型。法律的对象不是辐射,而是能量包的储存,它属于黑体的壁。
{"title":"Superposition of forward and backward motion","authors":"Manfred Harringer","doi":"10.1515/mcma-2022-2124","DOIUrl":"https://doi.org/10.1515/mcma-2022-2124","url":null,"abstract":"Abstract I consider black body radiation. The wall of the black body exchanges photons with the radiation field in equilibrium, therefore with a common temperature in Planck’s radiation law. The underlying process of radiation consists of creation and annihilation of photons. I want to present an alternate model of motions, where the process of radiation consists of small steps in positive and negative direction, not zero in mean. The detection of radiation consists of storing and restoring of packages of energy. I get an analogue of Planck’s radiation law, where the common temperature emerges from the underlying common model of small steps. The object of the law is not the radiation, but a storage of packages of energy, which belongs to the wall of the black body.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"28 1","pages":"329 - 339"},"PeriodicalIF":0.9,"publicationDate":"2022-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45461671","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
Estimation of entropy and extropy based on right censored data: A Bayesian non-parametric approach 基于右截尾数据的熵和熵估计:一种贝叶斯非参数方法
IF 0.9 Q3 STATISTICS & PROBABILITY Pub Date : 2022-09-30 DOI: 10.1515/mcma-2022-2123
L. Al-Labadi, Muhammad Tahir
Abstract Entropy and extropy are central measures in information theory. In this paper, Bayesian non-parametric estimators to entropy and extropy with possibly right censored data are proposed. The approach uses the beta-Stacy process and the difference operator. Examples are presented to illustrate the performance of the estimators.
熵和熵是信息论的核心度量。本文提出了含有可能正确截尾数据的熵和熵的贝叶斯非参数估计。该方法使用了beta-Stacy过程和差分运算符。举例说明了该估计器的性能。
{"title":"Estimation of entropy and extropy based on right censored data: A Bayesian non-parametric approach","authors":"L. Al-Labadi, Muhammad Tahir","doi":"10.1515/mcma-2022-2123","DOIUrl":"https://doi.org/10.1515/mcma-2022-2123","url":null,"abstract":"Abstract Entropy and extropy are central measures in information theory. In this paper, Bayesian non-parametric estimators to entropy and extropy with possibly right censored data are proposed. The approach uses the beta-Stacy process and the difference operator. Examples are presented to illustrate the performance of the estimators.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"28 1","pages":"319 - 328"},"PeriodicalIF":0.9,"publicationDate":"2022-09-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42910811","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
On the practical point of view of option pricing 论期权定价的实用观点
IF 0.9 Q3 STATISTICS & PROBABILITY Pub Date : 2022-09-28 DOI: 10.1515/mcma-2022-2122
N. Halidias
Abstract In this note, we describe a new approach to the option pricing problem by introducing the notion of the safe (and acceptable) price for the writer of an option, in contrast to the fair price used in the Black–Scholes model. Our starting point is that the option pricing problem is closely related with the hedging problem by practical techniques. Recalling that the Black–Scholes model does not give us the price of the option but the initial value of a replicating portfolio, we observe easily that this has a serious disadvantage because it assumes the building of this replicating portfolio continuously in time, and this is a disadvantage of any model that assumes such a construction. Here we study the problem from the practical point of view concerning mainly the over-the-counter market. This approach is not affected by the number of the underlying assets and is particularly useful for incomplete markets. In the usual Black–Scholes or binomial approach or some other approaches, one assumes that one can invest or borrow at the same risk-free rate r > 0 r>0 , which is not true in general. Even if this is the case, one can immediately observe that this risk-free rate is not a universal constant but is different among different people or institutions. So the fair price of an option is not so much fair! Moreover, the two sides are not, in general, equivalent against the risk; therefore, the notion of a fair price has no meaning at all. We also define a variant of the usual binomial model, by estimating safe upward and downward rates u , d u,d , trying to give a cheaper safe or acceptable price for the option.
在这篇文章中,我们描述了一种解决期权定价问题的新方法,通过引入期权出售者的安全(和可接受)价格的概念,与Black-Scholes模型中使用的公平价格形成对比。我们的出发点是,期权定价问题与套期保值问题在实际技术上是密切相关的。回顾Black-Scholes模型并没有给出期权的价格,而是给出了一个可复制投资组合的初始值,我们很容易发现,这个模型有一个严重的缺点,因为它假设这个可复制投资组合在时间上是连续建立的,这是任何假设这种结构的模型的缺点。这里我们主要从实际的角度研究场外交易市场的问题。这种方法不受标的资产数量的影响,对不完全市场特别有用。在通常的布莱克-斯科尔斯方法或二项方法或其他方法中,人们假设人们可以以相同的无风险利率进行投资或借贷,这通常是不正确的。即使是这样,人们也可以立即观察到,这个无风险利率不是一个普遍常数,而是在不同的人或机构之间有所不同。所以期权的公平价格并不是那么公平!此外,一般来说,双方在风险方面并不相等;因此,公平价格的概念根本没有意义。我们还定义了通常的二项模型的一种变体,通过估计安全的上升和下降率u,d, u,d,试图给出一个更便宜的安全或可接受的期权价格。
{"title":"On the practical point of view of option pricing","authors":"N. Halidias","doi":"10.1515/mcma-2022-2122","DOIUrl":"https://doi.org/10.1515/mcma-2022-2122","url":null,"abstract":"Abstract In this note, we describe a new approach to the option pricing problem by introducing the notion of the safe (and acceptable) price for the writer of an option, in contrast to the fair price used in the Black–Scholes model. Our starting point is that the option pricing problem is closely related with the hedging problem by practical techniques. Recalling that the Black–Scholes model does not give us the price of the option but the initial value of a replicating portfolio, we observe easily that this has a serious disadvantage because it assumes the building of this replicating portfolio continuously in time, and this is a disadvantage of any model that assumes such a construction. Here we study the problem from the practical point of view concerning mainly the over-the-counter market. This approach is not affected by the number of the underlying assets and is particularly useful for incomplete markets. In the usual Black–Scholes or binomial approach or some other approaches, one assumes that one can invest or borrow at the same risk-free rate r > 0 r>0 , which is not true in general. Even if this is the case, one can immediately observe that this risk-free rate is not a universal constant but is different among different people or institutions. So the fair price of an option is not so much fair! Moreover, the two sides are not, in general, equivalent against the risk; therefore, the notion of a fair price has no meaning at all. We also define a variant of the usual binomial model, by estimating safe upward and downward rates u , d u,d , trying to give a cheaper safe or acceptable price for the option.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"28 1","pages":"307 - 318"},"PeriodicalIF":0.9,"publicationDate":"2022-09-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46929711","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
Monte Carlo simulation of sensitivity functions for few-view computed tomography of strongly absorbing media 强吸收介质少视点计算机断层扫描灵敏度函数的蒙特卡罗模拟
IF 0.9 Q3 STATISTICS & PROBABILITY Pub Date : 2022-08-25 DOI: 10.1515/mcma-2022-2120
A. Konovalov, V. Vlasov, S. Kolchugin, G. Malyshkin, R. Mukhamadiyev
Abstract The paper describes a sensitivity function calculation method for few-view X-ray computed tomography of strongly absorbing objects. It is based on a probabilistic interpretation of energy transport through the object from a source to a detector. A PRIZMA code package is used to track photons. The code is developed at FSUE “RFNC–VNIITF named after Academ. E. I. Zababakhin” and implements a stochastic Monte Carlo method. The value of the sensitivity function in a discrete cell of the reconstruction region is assumed to be directly proportional to the fraction of photon trajectories which cross the cell from all those recorded by the detector. The method’s efficiency is validated through a numerical experiment on the reconstruction of a section of a spherical heavy-metal phantom with an air cavity and a density difference of 25 Ṫhe proposed method is shown to outperform the method based on projection approximation in case of reconstruction from 9 views.
本文介绍了强吸收物体少视场X射线计算机断层扫描的灵敏度函数计算方法。它基于对物体从源到探测器的能量传输的概率解释。PRIZMA代码包用于跟踪光子。该代码是在以Academy.E命名的FSUE“RFNC–VNIITF”上开发的。 I.Zababakhin”,并实现了一种随机蒙特卡罗方法。假设重建区域的离散单元中的灵敏度函数的值与探测器记录的所有光子轨迹中穿过该单元的光子轨迹的分数成正比。该方法的有效性是通过一个具有空气腔和密度差为25的球形重金属体模截面的重建的数值实验来验证的 Ṫ在从9个视图重建的情况下,所提出的方法优于基于投影近似的方法。
{"title":"Monte Carlo simulation of sensitivity functions for few-view computed tomography of strongly absorbing media","authors":"A. Konovalov, V. Vlasov, S. Kolchugin, G. Malyshkin, R. Mukhamadiyev","doi":"10.1515/mcma-2022-2120","DOIUrl":"https://doi.org/10.1515/mcma-2022-2120","url":null,"abstract":"Abstract The paper describes a sensitivity function calculation method for few-view X-ray computed tomography of strongly absorbing objects. It is based on a probabilistic interpretation of energy transport through the object from a source to a detector. A PRIZMA code package is used to track photons. The code is developed at FSUE “RFNC–VNIITF named after Academ. E. I. Zababakhin” and implements a stochastic Monte Carlo method. The value of the sensitivity function in a discrete cell of the reconstruction region is assumed to be directly proportional to the fraction of photon trajectories which cross the cell from all those recorded by the detector. The method’s efficiency is validated through a numerical experiment on the reconstruction of a section of a spherical heavy-metal phantom with an air cavity and a density difference of 25 Ṫhe proposed method is shown to outperform the method based on projection approximation in case of reconstruction from 9 views.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"28 1","pages":"269 - 278"},"PeriodicalIF":0.9,"publicationDate":"2022-08-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44779220","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Scrambling additive lagged-Fibonacci generators 置乱加性滞后斐波那契生成器
IF 0.9 Q3 STATISTICS & PROBABILITY Pub Date : 2022-08-04 DOI: 10.1515/mcma-2022-2115
Haifa Aldossari, M. Mascagni
Abstract Random numbers are used in a variety of applications including simulation, sampling, and cryptography. Fortunately, there exist many well-established methods of random number generation. An example of a well-known pseudorandom number generator is the lagged-Fibonacci generator (LFG). Marsaglia showed that the lagged-Fibonacci generator using addition failed some of his DIEHARD statistical tests, while it passed all when longer lags were used. This paper presents a scrambler that takes bits from a pseudorandom number generator and outputs (hopefully) improved pseudorandom numbers. The scrambler is based on a modified Feistel function, a method used in the generation of cryptographic random numbers, and multiplication by a chosen multiplier. We show that this scrambler improves the quality of pseudorandom numbers by applying it to the additive LFG with small lags. The scrambler performs well based on its performance with the TestU01 suite of randomness tests. The TestU01 suite of randomness tests is more comprehensive than the DIEHARD tests. In fact, the specific suite of tests we used from TestU01 includes the DIEHARD tests The scrambling of the LFG is so successful that scrambled LFGs with small lags perform as well as unscrambled LFGs with long lags. This comes at the cost of a doubling of execution time, and provides users with generators with small memory footprints that can provide parallel generators like the LFGs in the SPRNG parallel random number generation package.
随机数用于各种应用,包括仿真、采样和密码学。幸运的是,存在许多完善的随机数生成方法。一个众所周知的伪随机数生成器的例子是滞后斐波那契生成器(LFG)。Marsaglia表明,使用加法的滞后斐波那契生成器未能通过他的一些统计测试,而当使用更长的滞后时,它通过了所有测试。本文提出了一种扰频器,它从伪随机数生成器中获取比特并输出(希望)改进的伪随机数。扰频器是基于一个改进的费斯特尔函数,一种用于生成加密随机数的方法,并乘以一个选定的乘数。我们证明了该扰频器通过将其应用于具有小滞后的加性LFG来提高伪随机数的质量。基于其在TestU01随机测试套件中的性能,该扰频器表现良好。test01随机测试套件比DIEHARD测试更全面。实际上,我们从test01中使用的特定测试套件包括DIEHARD测试。对LFG的置乱非常成功,具有小延迟的置乱LFG的性能与具有长延迟的未置乱LFG一样好。这是以双倍的执行时间为代价的,并且为用户提供了内存占用较小的生成器,可以提供类似SPRNG并行随机数生成包中的lfg这样的并行生成器。
{"title":"Scrambling additive lagged-Fibonacci generators","authors":"Haifa Aldossari, M. Mascagni","doi":"10.1515/mcma-2022-2115","DOIUrl":"https://doi.org/10.1515/mcma-2022-2115","url":null,"abstract":"Abstract Random numbers are used in a variety of applications including simulation, sampling, and cryptography. Fortunately, there exist many well-established methods of random number generation. An example of a well-known pseudorandom number generator is the lagged-Fibonacci generator (LFG). Marsaglia showed that the lagged-Fibonacci generator using addition failed some of his DIEHARD statistical tests, while it passed all when longer lags were used. This paper presents a scrambler that takes bits from a pseudorandom number generator and outputs (hopefully) improved pseudorandom numbers. The scrambler is based on a modified Feistel function, a method used in the generation of cryptographic random numbers, and multiplication by a chosen multiplier. We show that this scrambler improves the quality of pseudorandom numbers by applying it to the additive LFG with small lags. The scrambler performs well based on its performance with the TestU01 suite of randomness tests. The TestU01 suite of randomness tests is more comprehensive than the DIEHARD tests. In fact, the specific suite of tests we used from TestU01 includes the DIEHARD tests The scrambling of the LFG is so successful that scrambled LFGs with small lags perform as well as unscrambled LFGs with long lags. This comes at the cost of a doubling of execution time, and provides users with generators with small memory footprints that can provide parallel generators like the LFGs in the SPRNG parallel random number generation package.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"28 1","pages":"199 - 210"},"PeriodicalIF":0.9,"publicationDate":"2022-08-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46604130","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Standard deviation estimation from sums of unequal size samples 不等大小样本和的标准差估计
IF 0.9 Q3 STATISTICS & PROBABILITY Pub Date : 2022-08-04 DOI: 10.1515/mcma-2022-2118
M. Casquilho, J. Buescu
Abstract In numerous industrial and related activities, the sums of the values of, frequently, unequal size samples are systematically recorded, for several purposes such as legal or quality control reasons. For the typical case where the individual values are not or no longer known, we address the point estimation, with confidence intervals, of the standard deviation (and mean) of the individual items, from those sums alone. The estimation may be useful also to corroborate estimates from previous statistical process control. An everyday case of a sum is the total weight of a set of items, such as a load of bags on a truck, which is used illustratively. For the parameters mean and standard deviation of the distribution, assumed Gaussian, we derive point estimates, which lead to weighted statistics, and we derive confidence intervals. For the latter, starting with a tentative reduction to equal size samples, we arrive at a solid conjecture for the mean, and a proposal for the standard deviation. All results are verifiable by direct computation or by simulation in a general and effective way. These computations can be run on public web pages of ours, namely for possible industrial use.
摘要在许多工业和相关活动中,由于法律或质量控制等原因,通常会系统地记录大小不等的样本的值总和。对于单个值未知或不再已知的典型情况,我们仅从这些总和中,用置信区间来处理单个项目的标准偏差(和平均值)的点估计。该估计也可用于证实来自先前统计过程控制的估计。总和的日常情况是一组物品的总重量,例如卡车上的一车袋子,这是示例性使用的。对于分布的参数均值和标准差,假设为高斯,我们导出点估计,这导致加权统计,我们导出置信区间。对于后者,从尝试性地减少到相等大小的样本开始,我们得出了平均值的可靠猜想,以及标准偏差的建议。所有结果都可以通过直接计算或模拟以通用有效的方式进行验证。这些计算可以在我们的公共网页上运行,即用于可能的工业用途。
{"title":"Standard deviation estimation from sums of unequal size samples","authors":"M. Casquilho, J. Buescu","doi":"10.1515/mcma-2022-2118","DOIUrl":"https://doi.org/10.1515/mcma-2022-2118","url":null,"abstract":"Abstract In numerous industrial and related activities, the sums of the values of, frequently, unequal size samples are systematically recorded, for several purposes such as legal or quality control reasons. For the typical case where the individual values are not or no longer known, we address the point estimation, with confidence intervals, of the standard deviation (and mean) of the individual items, from those sums alone. The estimation may be useful also to corroborate estimates from previous statistical process control. An everyday case of a sum is the total weight of a set of items, such as a load of bags on a truck, which is used illustratively. For the parameters mean and standard deviation of the distribution, assumed Gaussian, we derive point estimates, which lead to weighted statistics, and we derive confidence intervals. For the latter, starting with a tentative reduction to equal size samples, we arrive at a solid conjecture for the mean, and a proposal for the standard deviation. All results are verifiable by direct computation or by simulation in a general and effective way. These computations can be run on public web pages of ours, namely for possible industrial use.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"28 1","pages":"235 - 253"},"PeriodicalIF":0.9,"publicationDate":"2022-08-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45218319","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 2
Simulation of transient and spatial structure of the radiative flux produced by multiple recombinations of excitons 激子多次重组产生的辐射通量的瞬态和空间结构的模拟
IF 0.9 Q3 STATISTICS & PROBABILITY Pub Date : 2022-08-04 DOI: 10.1515/mcma-2022-2117
K. Sabelfeld, V. Sapozhnikov
Abstract In this paper, we study the multiple recombination exciton–photon–exciton process governed by a coupled system of the drift-diffusion-recombination equation and the integral radiative transfer equation. We develop a random walk on spheres algorithm for solving this system of equations. The algorithm directly simulates the transient drift-diffusion process of exciton’s motion. Then, at a random time the exciton recombines to a photon that moves in accordance with the radiative transfer equation, which in turn may recombine to an exciton etc. This algorithm is applied to calculate fluxes of excitons and photons as functions of time, and some other characteristics of the process. Calculations have also been carried out to validate the constructed model.
摘要在本文中,我们研究了由漂移-扩散-复合方程和积分辐射传输方程的耦合系统控制的多重复合激子-光子-激子过程。我们开发了一个求解该方程组的随机球上行走算法。该算法直接模拟了激子运动的瞬态漂移扩散过程。然后,在随机时间,激子重组为根据辐射传输方程移动的光子,而光子又可能重组为激子等。该算法用于计算激子和光子的通量作为时间的函数,以及该过程的一些其他特性。还进行了计算,以验证所构建的模型。
{"title":"Simulation of transient and spatial structure of the radiative flux produced by multiple recombinations of excitons","authors":"K. Sabelfeld, V. Sapozhnikov","doi":"10.1515/mcma-2022-2117","DOIUrl":"https://doi.org/10.1515/mcma-2022-2117","url":null,"abstract":"Abstract In this paper, we study the multiple recombination exciton–photon–exciton process governed by a coupled system of the drift-diffusion-recombination equation and the integral radiative transfer equation. We develop a random walk on spheres algorithm for solving this system of equations. The algorithm directly simulates the transient drift-diffusion process of exciton’s motion. Then, at a random time the exciton recombines to a photon that moves in accordance with the radiative transfer equation, which in turn may recombine to an exciton etc. This algorithm is applied to calculate fluxes of excitons and photons as functions of time, and some other characteristics of the process. Calculations have also been carried out to validate the constructed model.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"28 1","pages":"255 - 268"},"PeriodicalIF":0.9,"publicationDate":"2022-08-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47080040","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Approximate bounding of mixing time for multiple-step Gibbs samplers 多阶Gibbs采样器混合时间的近似边界
IF 0.9 Q3 STATISTICS & PROBABILITY Pub Date : 2022-08-04 DOI: 10.1515/mcma-2022-2119
David A. Spade
Abstract Markov chain Monte Carlo (MCMC) methods are important in a variety of statistical applications that require sampling from intractable probability distributions. Among the most common MCMC algorithms is the Gibbs sampler. When an MCMC algorithm is used, it is important to have an idea of how long it takes for the chain to become “close” to its stationary distribution. In many cases, there is high autocorrelation in the output of the chain, so the output needs to be thinned so that an approximate random sample from the desired probability distribution can be obtained by taking a state of the chain every h steps in a process called h-thinning. This manuscript extends the work of [D. A. Spade, Estimating drift and minorization coefficients for Gibbs sampling algorithms, Monte Carlo Methods Appl. 27 2021, 3, 195–209] by presenting a computational approach to obtaining an approximate upper bound on the mixing time of the h-thinned Gibbs sampler.
摘要马尔可夫链蒙特卡罗(MCMC)方法在各种需要从棘手的概率分布中采样的统计应用中是重要的。最常见的MCMC算法是吉布斯采样器。当使用MCMC算法时,重要的是要了解链需要多长时间才能“接近”其平稳分布。在许多情况下,链的输出具有很高的自相关性,因此需要对输出进行细化,以便通过在称为h细化的过程中每h步获取链的状态,可以从所需概率分布中获得近似的随机样本。该手稿扩展了[D.A.Spade,估计吉布斯采样算法的漂移和二阶化系数,蒙特卡罗方法应用27 2021,3195–209]的工作,提出了一种计算方法来获得h稀疏吉布斯采样器混合时间的近似上限。
{"title":"Approximate bounding of mixing time for multiple-step Gibbs samplers","authors":"David A. Spade","doi":"10.1515/mcma-2022-2119","DOIUrl":"https://doi.org/10.1515/mcma-2022-2119","url":null,"abstract":"Abstract Markov chain Monte Carlo (MCMC) methods are important in a variety of statistical applications that require sampling from intractable probability distributions. Among the most common MCMC algorithms is the Gibbs sampler. When an MCMC algorithm is used, it is important to have an idea of how long it takes for the chain to become “close” to its stationary distribution. In many cases, there is high autocorrelation in the output of the chain, so the output needs to be thinned so that an approximate random sample from the desired probability distribution can be obtained by taking a state of the chain every h steps in a process called h-thinning. This manuscript extends the work of [D. A. Spade, Estimating drift and minorization coefficients for Gibbs sampling algorithms, Monte Carlo Methods Appl. 27 2021, 3, 195–209] by presenting a computational approach to obtaining an approximate upper bound on the mixing time of the h-thinned Gibbs sampler.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"28 1","pages":"221 - 233"},"PeriodicalIF":0.9,"publicationDate":"2022-08-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43553227","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 1
期刊
Monte Carlo Methods and Applications
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1