BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Xiaolin Luo

Xiaolin Luo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
q-fin.CPq-fin.RMastro-ph.COq-fin.PRgr-qchep-phhep-thmath.NAmath.PRphysics.data-anphysics.geo-phq-fin.ST

Catalog footprint

What is connected

14works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Cosmological constraints from the density gradient weighted correlation function

The mark weighted correlation function (MCF) $W(s,μ)$ is a computationally efficient statistical measure which can probe clustering information beyond that of the conventional 2-point statistics. In this work, we extend the traditional mark weighted statistics by using powers of the density field gradient $|\nabla ρ/ρ|^α$ as the weight, and use the angular dependence of the scale-averaged MCFs to constrain cosmological parameters. The analysis shows that the gradient based weighting scheme is statistically more powerful than the density based weighting scheme, while combining the two schemes together is more powerful than separately using either of them. Utilising the density weighted or the gradient weighted MCFs with $α=0.5,\ 1$, we can strengthen the constraint on $Ω_m$ by factors of 2 or 4, respectively, compared with the standard 2-point correlation function, while simultaneously using the MCFs of the two weighting schemes together can be $1.25$ times more statistically powerful than using the gradient weighting scheme alone. The mark weighted statistics may play an important role in cosmological analysis of future large-scale surveys. Many issues, including the possibility of using other types of weights, the influence of the bias on this statistics, as well as the usage of MCFs in the tomographic Alcock-Paczynski method, are worth further investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Can Non-standard Recombination Resolve the Hubble Tension?

The inconsistent Hubble constant values derived from cosmic microwave background (CMB) observations and from local distance-ladder measurements may suggest new physics beyond the standard $Λ$CDM paradigm. It has been found in earlier works that, at least phenomenologically, non-standard recombination histories can reduce the $\gtrsim 4σ$ Hubble tension to $\sim 2σ$. Following this path, we vary physical and phenomenological parameters in RECFAST, the standard code to compute ionization history of the universe, to explore possible physics beyond standard recombination. We find that the CMB constraint on the Hubble constant is sensitive to the Hydrogen ionization energy and $2s \rightarrow 1s$ two-photon decay rate, both of which are atomic constants, and is insensitive to other details of recombination. Thus, the Hubble tension is very robust against perturbations of recombination history, unless exotic physics modifies the atomic constants during the recombination epoch.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Cosmological constraints from the redshift dependence of the Alcock-Paczynski effect: Fourier space analysis

The tomographic Alcock-Paczynski (AP) method utilizes the redshift evolution of the AP distortion to place constraints on cosmological parameters. It has proved to be a robust method that can separate the AP signature from the redshift space distortion (RSD) effect, and deliver powerful cosmological constraints using the $\lesssim 40h^{-1}\ \rm Mpc$ clustering region. In previous works, the tomographic AP method was performed via the anisotropic 2-point correlation function statistic. In this work we consider the feasibility of conducting the analysis in the Fourier domain and examine the pros and cons of this approach. We use the integrated galaxy power spectrum (PS) as a function of direction, $\hat P_{Δk}(μ)$, to quantify the magnitude of anisotropy in the large-scale structure clustering, and use its redshift variation to do the AP test. The method is tested on the large, high resolution Big-MultiDark Planck (BigMD) simulation at redshifts $z=0-1$, using the underlying true cosmology $Ω_m=0.3071,\ w=-1$. Testing the redshift evolution of $\hat P_{Δk}(μ)$ in the true cosmology and cosmologies deviating from the truth with $δΩ_m=0.1,\ δw=0.3$, we find that the redshift evolution of the AP distortion overwhelms the effects created by the RSD by a factor of $\sim1.7-3.6$. We test the method in the range of $k\in(0.2,1.8)\ h\ \rm Mpc^{-1}$, and find that it works well throughout the entire regime. We tune the halo mass within the range $2\times 10^{13}$ to $10^{14}\ M_{\odot}$, and find that the change of halo bias results in $\lesssim 5 \%$ change in $\hat P_{Δk}(μ)$, which is less significant compared with the cosmological effect. Our work shows that it is feasible to conduct the tomographic AP analysis in the Fourier space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Valuation of Variable Annuities with Guaranteed Minimum Withdrawal and Death Benefits via Stochastic Control Optimization

In this paper we present a numerical valuation of variable annuities with combined Guaranteed Minimum Withdrawal Benefit (GMWB) and Guaranteed Minimum Death Benefit (GMDB) under optimal policyholder behaviour solved as an optimal stochastic control problem. This product simultaneously deals with financial risk, mortality risk and human behaviour. We assume that market is complete in financial risk and mortality risk is completely diversified by selling enough policies and thus the annuity price can be expressed as appropriate expectation. The computing engine employed to solve the optimal stochastic control problem is based on a robust and efficient Gauss-Hermite quadrature method with cubic spline. We present results for three different types of death benefit and show that, under the optimal policyholder behaviour, adding the premium for the death benefit on top of the GMWB can be problematic for contracts with long maturities if the continuous fee structure is kept, which is ordinarily assumed for a GMWB contract. In fact for some long maturities it can be shown that the fee cannot be charged as any proportion of the account value -- there is no solution to match the initial premium with the fair annuity price. On the other hand, the extra fee due to adding the death benefit can be charged upfront or in periodic instalment of fixed amount, and it is cheaper than buying a separate life insurance.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variable Annuity with GMWB: surrender or not, that is the question

Under the optimal withdrawal strategy of a policyholder, the pricing of variable annuities with Guaranteed Minimum Withdrawal Benefit (GMWB) is an optimal stochastic control problem. The surrender feature available in marketed products allows termination of the contract before maturity, making it also an optimal stopping problem. Although the surrender feature is quite common in variable annuity contracts, there appears to be no published analysis and results for this feature in GMWB under optimal policyholder behaviour - results found in the literature so far are consistent with the absence of such a feature. Also, it is of practical interest to see how the much simpler bang-bang strategy, although not optimal for GMWB, compares with optimal GMWB strategy with surrender option. In this paper we extend our recently developed algorithm (Luo and Shevchenko 2015a) to include surrender option in GMWB and compare prices under different policyholder strategies: optimal, static and bang-bang. Results indicate that following a simple but sub-optimal bang-bang strategy does not lead to significant reduction in the price or equivalently in the fee, in comparison with the optimal strategy. We observed that the extra value added by the surrender option could add very significant value to the GMWB contract. We also performed calculations for static withdrawal with surrender option, which is the same as bang-bang minus the "no-withdrawal" choice. We find that the fee for such contract is only less than 1% smaller when compared to the case of bang-bang strategy, meaning that th "no-withdrawal" option adds little value to the contract.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fast and Simple Method for Pricing Exotic Options using Gauss-Hermite Quadrature on a Cubic Spline Interpolation

There is a vast literature on numerical valuation of exotic options using Monte Carlo, binomial and trinomial trees, and finite difference methods. When transition density of the underlying asset or its moments are known in closed form, it can be convenient and more efficient to utilize direct integration methods to calculate the required option price expectations in a backward time-stepping algorithm. This paper presents a simple, robust and efficient algorithm that can be applied for pricing many exotic options by computing the expectations using Gauss-Hermite integration quadrature applied on a cubic spline interpolation. The algorithm is fully explicit but does not suffer the inherent instability of the explicit finite difference counterpart. A `free' bonus of the algorithm is that it already contains the function for fast and accurate interpolation of multiple solutions required by many discretely monitored path dependent options. For illustrations, we present examples of pricing a series of American options with either Bermudan or continuous exercise features, and a series of exotic path-dependent options of target accumulation redemption note (TARN). Results of the new method are compared with Monte Carlo and finite difference methods, including some of the most advanced or best known finite difference algorithms in the literature. The comparison shows that, despite its simplicity, the new method can rival with some of the best finite difference algorithms in accuracy and at the same time it is significantly faster. Virtually the same algorithm can be applied to price other path-dependent financial contracts such as Asian options and variable annuities.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fast Numerical Method for Pricing of Variable Annuities with Guaranteed Minimum Withdrawal Benefit under Optimal Withdrawal Strategy

A variable annuity contract with Guaranteed Minimum Withdrawal Benefit (GMWB) promises to return the entire initial investment through cash withdrawals during the policy life plus the remaining account balance at maturity, regardless of the portfolio performance. Under the optimal withdrawal strategy of a policyholder, the pricing of variable annuities with GMWB becomes an optimal stochastic control problem. So far in the literature these contracts have only been evaluated by solving partial differential equations (PDE) using the finite difference method. The well-known Least-Squares or similar Monte Carlo methods cannot be applied to pricing these contracts because the paths of the underlying wealth process are affected by optimal cash withdrawals (control variables) and thus cannot be simulated forward in time. In this paper we present a very efficient new algorithm for pricing these contracts in the case when transition density of the underlying asset between withdrawal dates or its moments are known. This algorithm relies on computing the expected contract value through a high order Gauss-Hermite quadrature applied on a cubic spline interpolation. Numerical results from the new algorithm for a series of GMWB contract are then presented, in comparison with results using the finite difference method solving corresponding PDE. The comparison demonstrates that the new algorithm produces results in very close agreement with those of the finite difference method, but at the same time it is significantly faster; virtually instant results on a standard desktop PC.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Markov chain Monte Carlo estimation of default and recovery: dependent via the latent systematic factor

It is a well known fact that recovery rates tend to go down when the number of defaults goes up in economic downturns. We demonstrate how the loss given default model with the default and recovery dependent via the latent systematic risk factor can be estimated using Bayesian inference methodology and Markov chain Monte Carlo method. This approach is very convenient for joint estimation of all model parameters and latent systematic factors. Moreover, all relevant uncertainties are easily quantified. Typically available data are annual averages of defaults and recoveries and thus the datasets are small and parameter uncertainty is significant. In this case Bayesian approach is superior to the maximum likelihood method that relies on a large sample limit Gaussian approximation for the parameter uncertainty. As an example, we consider a homogeneous portfolio with one latent factor. However, the approach can be easily extended to deal with non-homogenous portfolios and several latent factors.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bayesian Model Choice of Grouped t-copula

One of the most popular copulas for modeling dependence structures is t-copula. Recently the grouped t-copula was generalized to allow each group to have one member only, so that a priori grouping is not required and the dependence modeling is more flexible. This paper describes a Markov chain Monte Carlo (MCMC) method under the Bayesian inference framework for estimating and choosing t-copula models. Using historical data of foreign exchange (FX) rates as a case study, we found that Bayesian model choice criteria overwhelmingly favor the generalized t-copula. In addition, all the criteria also agree on the second most likely model and these inferences are all consistent with classical likelihood ratio tests. Finally, we demonstrate the impact of model choice on the conditional Value-at-Risk for portfolios of six major FX rates.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dependent default and recovery: MCMC study of downturn LGD credit risk model

There is empirical evidence that recovery rates tend to go down just when the number of defaults goes up in economic downturns. This has to be taken into account in estimation of the capital against credit risk required by Basel II to cover losses during the adverse economic downturns; the so-called "downturn LGD" requirement. This paper presents estimation of the LGD credit risk model with default and recovery dependent via the latent systematic risk factor using Bayesian inference approach and Markov chain Monte Carlo method. This approach allows joint estimation of all model parameters and latent systematic factor, and all relevant uncertainties. Results using Moody's annual default and recovery rates for corporate bonds for the period 1982-2010 show that the impact of parameter uncertainty on economic capital can be very significant and should be assessed by practitioners.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Short Tale of Long Tail Integration

Integration of the form $\int_a^\infty {f(x)w(x)dx} $, where $w(x)$ is either $\sin (ω{\kern 1pt} x)$ or $\cos (ω{\kern 1pt} x)$, is widely encountered in many engineering and scientific applications, such as those involving Fourier or Laplace transforms. Often such integrals are approximated by a numerical integration over a finite domain $(a,\,b)$, leaving a truncation error equal to the tail integration $\int_b^\infty {f(x)w(x)dx} $ in addition to the discretization error. This paper describes a very simple, perhaps the simplest, end-point correction to approximate the tail integration, which significantly reduces the truncation error and thus increases the overall accuracy of the numerical integration, with virtually no extra computational effort. Higher order correction terms and error estimates for the end-point correction formula are also derived. The effectiveness of this one-point correction formula is demonstrated through several examples.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Computing Tails of Compound Distributions Using Direct Numerical Integration

An efficient adaptive direct numerical integration (DNI) algorithm is developed for computing high quantiles and conditional Value at Risk (CVaR) of compound distributions using characteristic functions. A key innovation of the numerical scheme is an effective tail integration approximation that reduces the truncation errors significantly with little extra effort. High precision results of the 0.999 quantile and CVaR were obtained for compound losses with heavy tails and a very wide range of loss frequencies using the DNI, Fast Fourier Transform (FFT) and Monte Carlo (MC) methods. These results, particularly relevant to operational risk modelling, can serve as benchmarks for comparing different numerical methods. We found that the adaptive DNI can achieve high accuracy with relatively coarse grids. It is much faster than MC and competitive with FFT in computing high quantiles and CVaR of compound distributions in the case of moderate to high frequencies and heavy tails.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraining the shape of a gravity anomalous body using reversible jump Markov chain Monte Carlo

Typical geophysical inversion problems are ill-posed, non-linear and non-unique. Sometimes the problem is trans-dimensional, where the number of unknown parameters is one of the unknowns, which makes the inverse problem even more challenging. Detecting the shape of a geophysical object underneath the earth surface from gravity anomaly is one of such complex problems, where the number of geometrical parameters is one of the unknowns. To deal with the difficulties of non-uniqueness, ill-conditioning and nonlinearity, a statistical Bayesian model inference approach is adopted. A reversible jump Markov chain Monte Carlo (RJMCMC) algorithm is proposed to overcome the difficulty of trans-dimensionality. Carefully designed within-model and between-model Markov chain moves are implemented to reduce the rate of generating inadmissible geometries, thus achieving good overall efficiency in the Monte Carlo sampler. Numerical experiments on a 2-D problem show that the proposed algorithm is capable of obtaining satisfactory solutions with quantifiable uncertainty to a challenging trans-dimensional geophysical inverse problem. Solutions from RJMCMC appear to be parsimonious for the given prior, in the sense that among the models satisfactorily represent the true model, models with higher posterior probabilities tend to have fewer number of parameters. The proposed numerical algorithm can be readily adapted to other similar trans-dimensional geophysical inverse applications. Keywords: trans-dimensional geophysical inversion, reversible jump Markov chain Monte Carlo; gravity anomaly.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The t copula with Multiple Parameters of Degrees of Freedom: Bivariate Characteristics and Application to Risk Management

The t copula is often used in risk management as it allows for modelling tail dependence between risks and it is simple to simulate and calibrate. However, the use of a standard t copula is often criticized due to its restriction of having a single parameter for the degrees of freedom (dof) that may limit its capability to model the tail dependence structure in a multivariate case. To overcome this problem, grouped t copula was proposed recently, where risks are grouped a priori in such a way that each group has a standard t copula with its specific dof parameter. In this paper we propose the use of a grouped t copula, where each group consists of one risk factor only, so that a priori grouping is not required. The copula characteristics in the bivariate case are studied. We explain simulation and calibration procedures, including a simulation study on finite sample properties of the maximum likelihood estimators and Kendall's tau approximation. This new copula can be significantly different from the standard t copula in terms of risk measures such as tail dependence, value at risk and expected shortfall. Keywords: grouped t copula, tail dependence, risk management.

0 citations0 reviewsNo code linkedOpen access