针对短采样宽带信号近似最大似然(approximated maximum likelihood,AML)方位估计计算量大的问题,将马尔科夫链-蒙特卡罗方法与近似最大似然方位估计相结合,提出一种基于Metropolis-Hastings抽样的近似最大似然方位估计方法(AMLMH)。该...针对短采样宽带信号近似最大似然(approximated maximum likelihood,AML)方位估计计算量大的问题,将马尔科夫链-蒙特卡罗方法与近似最大似然方位估计相结合,提出一种基于Metropolis-Hastings抽样的近似最大似然方位估计方法(AMLMH)。该方法将AML算法的空间谱函数作为信号的概率分布函数,并利用Metropolis-Hastings抽样方法从该概率分布函数中抽样。研究结果表明,AMLMH方法不但保持了原近似最大似然方位估计方法的优良性能,而且减小了计算量。展开更多
This paper mainly talks about a popular approach of volatility of a GARCH-type model in R, while the disturbances are independent and have identical Student-t distribution. It uses the Metropolis-Hastings method to pe...This paper mainly talks about a popular approach of volatility of a GARCH-type model in R, while the disturbances are independent and have identical Student-t distribution. It uses the Metropolis-Hastings method to perform the computations and gives the programs in details in R.展开更多
Metropolis-Hastings algorithms are slowed down by the computation of complex target distributions. To solve this problem, one can use the delayed acceptance Metropolis-Hastings algorithm (MHDA) of Christen and Fox (20...Metropolis-Hastings algorithms are slowed down by the computation of complex target distributions. To solve this problem, one can use the delayed acceptance Metropolis-Hastings algorithm (MHDA) of Christen and Fox (2005). However, the acceptance rate of a proposed value will always be less than in the standard Metropolis-Hastings. We can fix this problem by using the Metropolis-Hastings algorithm with delayed rejection (MHDR) proposed by Tierney and Mira (1999). In this paper, we combine the ideas of MHDA and MHDR to propose a new MH algorithm, named the Metropolis-Hastings algorithm with delayed acceptance and rejection (MHDAR). The new algorithm reduces the computational cost by division of the prior or likelihood functions and increase the acceptance probability by delay rejection of the second stage. We illustrate those accelerating features by a realistic example.展开更多
载荷外推作为载荷谱编制的重要技术手段,当前研究缺乏对于载荷外推总体方法的全面梳理、马尔可夫稳态分布的求解方法适应性不够、缺乏不同非参频次外推方法的比较与选用原则,导致不便生成高精度载荷谱以支撑装备性能设计。围绕坦克在高...载荷外推作为载荷谱编制的重要技术手段,当前研究缺乏对于载荷外推总体方法的全面梳理、马尔可夫稳态分布的求解方法适应性不够、缺乏不同非参频次外推方法的比较与选用原则,导致不便生成高精度载荷谱以支撑装备性能设计。围绕坦克在高机动和极限工况下的载荷谱编制问题,基于某坦克行进间身管位移数据样本,分别使用基于雨流矩阵及核密度估计的非参数外推法、基于马尔可夫链蒙特卡洛(Markov Chain Monte Carlo,MCMC)的信号重构法以及Metropolis-Hastings(简称MH)直接采样法进行了载荷频次外推,并针对MCMC的信号重构法提出了一种改良马尔可夫稳态分布的求解方法。应用所提出的频次-极值相结合的载荷外推总体方法对坦克身管位移进行了频次扩充与极值预测,并结合实车试验结果验证了方法的准确性。研究结果表明:改良的马尔可夫稳态分布求解方法是有效的;在样本长度足够、外推精度要求不甚高的情况下,MH直接采样法可作为一种新的频次外推方法;运用频次-极值相结合的载荷外推总体方法所得结果精度较高;形成的频次外推法选用原则对于载荷谱编制过程中的方法选择具有一定的指导意义。研究工作为装备载荷谱的高质量编制提供了成熟的技术路线和参考。展开更多
The topic of this article is one-sided hypothesis testing for disparity, i.e., the mean of one group is larger than that of another when there is uncertainty as to which group a datum is drawn. For each datum, the unc...The topic of this article is one-sided hypothesis testing for disparity, i.e., the mean of one group is larger than that of another when there is uncertainty as to which group a datum is drawn. For each datum, the uncertainty is captured with a given discrete probability distribution over the groups. Such situations arise, for example, in the use of Bayesian imputation methods to assess race and ethnicity disparities with certain insurance, health, and financial data. A widely used method to implement this assessment is the Bayesian Improved Surname Geocoding (BISG) method which assigns a discrete probability over six race/ethnicity groups to an individual given the individual’s surname and address location. Using a Bayesian framework and Markov Chain Monte Carlo sampling from the joint posterior distribution of the group means, the probability of a disparity hypothesis is estimated. Four methods are developed and compared with an illustrative data set. Three of these methods are implemented in an R-code and one method in WinBUGS. These methods are programed for any number of groups between two and six inclusive. All the codes are provided in the appendices.展开更多
Here the estimating problem of a single sinusoidal signal in the additive symmetricα-stable Gaussian(ASαSG)noise is investigated.The ASαSG noise here is expressed as the additive of a Gaussian noise and a symmetric...Here the estimating problem of a single sinusoidal signal in the additive symmetricα-stable Gaussian(ASαSG)noise is investigated.The ASαSG noise here is expressed as the additive of a Gaussian noise and a symmetricα-stable distributed variable.As the probability density function(PDF)of the ASαSG is complicated,traditional estimators cannot provide optimum estimates.Based on the Metropolis-Hastings(M-H)sampling scheme,a robust frequency estimator is proposed for ASαSG noise.Moreover,to accelerate the convergence rate of the developed algorithm,a new criterion of reconstructing the proposal covar-iance is derived,whose main idea is updating the proposal variance using several previous samples drawn in each iteration.The approximation PDF of the ASαSG noise,which is referred to the weighted sum of a Voigt function and a Gaussian PDF,is also employed to reduce the computational complexity.The computer simulations show that the performance of our method is better than the maximum likelihood and the lp-norm estimators.展开更多
文摘针对短采样宽带信号近似最大似然(approximated maximum likelihood,AML)方位估计计算量大的问题,将马尔科夫链-蒙特卡罗方法与近似最大似然方位估计相结合,提出一种基于Metropolis-Hastings抽样的近似最大似然方位估计方法(AMLMH)。该方法将AML算法的空间谱函数作为信号的概率分布函数,并利用Metropolis-Hastings抽样方法从该概率分布函数中抽样。研究结果表明,AMLMH方法不但保持了原近似最大似然方位估计方法的优良性能,而且减小了计算量。
文摘This paper mainly talks about a popular approach of volatility of a GARCH-type model in R, while the disturbances are independent and have identical Student-t distribution. It uses the Metropolis-Hastings method to perform the computations and gives the programs in details in R.
文摘Metropolis-Hastings algorithms are slowed down by the computation of complex target distributions. To solve this problem, one can use the delayed acceptance Metropolis-Hastings algorithm (MHDA) of Christen and Fox (2005). However, the acceptance rate of a proposed value will always be less than in the standard Metropolis-Hastings. We can fix this problem by using the Metropolis-Hastings algorithm with delayed rejection (MHDR) proposed by Tierney and Mira (1999). In this paper, we combine the ideas of MHDA and MHDR to propose a new MH algorithm, named the Metropolis-Hastings algorithm with delayed acceptance and rejection (MHDAR). The new algorithm reduces the computational cost by division of the prior or likelihood functions and increase the acceptance probability by delay rejection of the second stage. We illustrate those accelerating features by a realistic example.
文摘载荷外推作为载荷谱编制的重要技术手段,当前研究缺乏对于载荷外推总体方法的全面梳理、马尔可夫稳态分布的求解方法适应性不够、缺乏不同非参频次外推方法的比较与选用原则,导致不便生成高精度载荷谱以支撑装备性能设计。围绕坦克在高机动和极限工况下的载荷谱编制问题,基于某坦克行进间身管位移数据样本,分别使用基于雨流矩阵及核密度估计的非参数外推法、基于马尔可夫链蒙特卡洛(Markov Chain Monte Carlo,MCMC)的信号重构法以及Metropolis-Hastings(简称MH)直接采样法进行了载荷频次外推,并针对MCMC的信号重构法提出了一种改良马尔可夫稳态分布的求解方法。应用所提出的频次-极值相结合的载荷外推总体方法对坦克身管位移进行了频次扩充与极值预测,并结合实车试验结果验证了方法的准确性。研究结果表明:改良的马尔可夫稳态分布求解方法是有效的;在样本长度足够、外推精度要求不甚高的情况下,MH直接采样法可作为一种新的频次外推方法;运用频次-极值相结合的载荷外推总体方法所得结果精度较高;形成的频次外推法选用原则对于载荷谱编制过程中的方法选择具有一定的指导意义。研究工作为装备载荷谱的高质量编制提供了成熟的技术路线和参考。
文摘The topic of this article is one-sided hypothesis testing for disparity, i.e., the mean of one group is larger than that of another when there is uncertainty as to which group a datum is drawn. For each datum, the uncertainty is captured with a given discrete probability distribution over the groups. Such situations arise, for example, in the use of Bayesian imputation methods to assess race and ethnicity disparities with certain insurance, health, and financial data. A widely used method to implement this assessment is the Bayesian Improved Surname Geocoding (BISG) method which assigns a discrete probability over six race/ethnicity groups to an individual given the individual’s surname and address location. Using a Bayesian framework and Markov Chain Monte Carlo sampling from the joint posterior distribution of the group means, the probability of a disparity hypothesis is estimated. Four methods are developed and compared with an illustrative data set. Three of these methods are implemented in an R-code and one method in WinBUGS. These methods are programed for any number of groups between two and six inclusive. All the codes are provided in the appendices.
基金supported by National Key R&D Program of China(Grant No.2018YFF01012600)National Natural Science Foundation of China(Grant No.61701021)Fundamental Research Funds for the Central Universities(Grant No.FRF-TP-19-006A3).
文摘Here the estimating problem of a single sinusoidal signal in the additive symmetricα-stable Gaussian(ASαSG)noise is investigated.The ASαSG noise here is expressed as the additive of a Gaussian noise and a symmetricα-stable distributed variable.As the probability density function(PDF)of the ASαSG is complicated,traditional estimators cannot provide optimum estimates.Based on the Metropolis-Hastings(M-H)sampling scheme,a robust frequency estimator is proposed for ASαSG noise.Moreover,to accelerate the convergence rate of the developed algorithm,a new criterion of reconstructing the proposal covar-iance is derived,whose main idea is updating the proposal variance using several previous samples drawn in each iteration.The approximation PDF of the ASαSG noise,which is referred to the weighted sum of a Voigt function and a Gaussian PDF,is also employed to reduce the computational complexity.The computer simulations show that the performance of our method is better than the maximum likelihood and the lp-norm estimators.