• We study an infinite-horizon discrete-time optimal stopping problem under non-exponential discounting. A new method, which we call the iterative approach, is developed to find subgame perfect Nash equilibria. When the discount function induces decreasing impatience, we establish the existence of an equilibrium through fixed-point iterations. Moreover, we show that there exists a unique optimal equilibrium, which generates larger value than any other equilibrium does at all times. To the best of our knowledge, this is the first time a dominating subgame perfect Nash equilibrium is shown to exist in the literature of time-inconsistency.
  • In this paper, we consider the estimation and inference of precision matrices of a rich class of locally stationary and nonlinear time series assuming that only one realization of the time series is observed. Using a Cholesky decomposition technique, we show that the precision matrices can be directly estimated via a series of least squares linear regressions with smoothly time-varying coefficients. The method of sieves is utilized for the estimation and is shown to be efficient and optimally adaptive in terms of estimation accuracy and computational complexity. We establish an asymptotic theory for a class of ${\cal L}^2$ tests based on the nonparametric sieve estimators. The latter are used for testing whether the precision matrices are diagonal or banded. A high dimensional Gaussian approximation result is established for a wide class of quadratic form of non-stationary and nonlinear processes, which is of interest by itself.
  • Let $\Omega$ to be one of $\X^{N+1},C[0,1],D[0,1]$: a product of Polish spaces, space of continuous functions from $[0,1]$ to a subset of $\mathbb{R}^d$, and space of RCLL (right-continuous with left limits) functions from $[0,1]$ to $\mathbb{R}^d$ respectively. We consider the existence of a probability measure $P$ on $\Omega$ such that $P$ has the given marginals $\alpha$ and $\beta$ and satisfies some other convex transport constraints, which is given by $\Gamma$. The main application we have in mind is the martingale optimal transport problem with when the martingales are assumed to have bounded volatility/quadratic variation. We show that such probability measure exists if and only if the $\alpha$ average of so-called $G$-expectation of bounded uniformly continuous and bounded functions with respect to the measures in $\Gamma$ is less than their $\beta$ average. As a byproduct, we get a necessary and sufficient condition for the Skorokhod embedding for bounded stopping times.
  • Inspired by Strotz's consistent planning strategy, we formulate the infinite horizon mean-variance stopping problem as a subgame perfect Nash equilibrium in order to determine time consistent strategies with no regret. Equilibria among stopping times or randomized stopping times may not exist. This motivates us to consider the notion of liquidation strategies, which lets the stopping right to be divisible. We then argue that the mean-standard deviation variant of this problem makes more sense for this type of strategies in terms of time consistency. It turns out that an equilibrium liquidation strategy always exists. We then analyze whether optimal equilibrium liquidation strategies exist and whether they are unique and observe that neither may hold.
  • A restrictive assumption in change point analysis is "stationarity under the null hypothesis of no change-point", which is crucial for asymptotic theory but not very realistic from a practical point of view. For example, if change point analysis for correlations is performed, it is not necessarily clear that the mean, marginal variance or higher order moments are constant, even if there is no change in the correlation. This paper develops change point analysis for the correlation structures under less restrictive assumptions. In contrast to previous work, our approach does not require that the mean, variance and fourth order joint cumulants are constant under the null hypothesis. Moreover, we also address the problem of detecting relevant change points.
  • Since most of the traded options on individual stocks is of American type it is of interest to generalize the results obtained in semi-static trading to the case when one is allowed to statically trade American options. However, this problem has proved to be elusive so far because of the asymmetric nature of the positions of holding versus shorting such options. Here we provide a unified framework and generalize the fundamental theorem of asset pricing (FTAP) and hedging dualities in arXiv:1502.06681 (to appear in Annals of Applied Probability) to the case where the investor can also short American options. Following arXiv:1502.06681, we assume that the longed American options are divisible. As for the shorted American options, we show that the divisibility plays no role regarding arbitrage property and hedging prices. Then using the method of enlarging probability spaces proposed in arXiv:1604.05517, we convert the shorted American options to European options, and establish the FTAP and sub- and super-hedging dualities in the enlarged space both with and without model uncertainty.
  • A neutron imaging detector based on neutron-sensitive microchannel plates (MCPs) was constructed and tested at beamlines of thermal and cold neutrons. The MCPs are made of a glass mixture containing B-10 and natural Gd, which makes the bulk of the MCP an efficient neutron converter. Contrary to the neutron sensitive scintillator screens normally used in neutron imaging, spatial resolution is not traded off with detection efficiency. While the best neutron imaging scintillators have a detection efficiency around a percent, a detection efficiency of around 50% for thermal neutrons and 70% for cold neutrons has been demonstrated with these MCPs earlier. Our tests show a performance similar to conventional neutron imaging detectors, apart from the orders of magnitude better sensitivity. We demonstrate a spatial resolution better than 150 um. The sensitivity of this detector allows fast tomography and neutron video recording, and will make smaller reactor sites and even portable sources suitable for neutron imaging.
  • In this paper, we introduce a new sparsity-promoting prior, namely, the "normal product" prior, and develop an efficient algorithm for sparse signal recovery under the Bayesian framework. The normal product distribution is the distribution of a product of two normally distributed variables with zero means and possibly different variances. Like other sparsity-encouraging distributions such as the Student's $t$-distribution, the normal product distribution has a sharp peak at origin, which makes it a suitable prior to encourage sparse solutions. A two-stage normal product-based hierarchical model is proposed. We resort to the variational Bayesian (VB) method to perform the inference. Simulations are conducted to illustrate the effectiveness of our proposed algorithm as compared with other state-of-the-art compressed sensing algorithms.
  • We consider the super-hedging price of an American option in a discrete-time market in which stocks are available for dynamic trading and European options are available for static trading. We show that the super-hedging price $\pi$ is given by the supremum over the prices of the American option under randomized models. That is, $\pi=\sup_{(c_i,Q_i)_i}\sum_ic_i\phi^{Q_i}$, where $c_i \in \mathbb{R}_+$ and the martingale measure $Q^i$ are chosen such that $\sum_i c_i=1$ and $\sum_i c_iQ_i$ prices the European options correctly, and $\phi^{Q_i}$ is the price of the American option under the model $Q_i$. Our result generalizes the example given in ArXiv:1604.02274 that the highest model based price can be considered as a randomization over models.
  • On a filtered probability space $(\Omega,\mathcal{F},P,\mathbb{F}=(\mathcal{F}_t)_{t=0,\dotso,T})$, we consider stopper-stopper games $\overline V:=\inf_{\Rho\in\bT^{ii}}\sup_{\tau\in\T}\E[U(\Rho(\tau),\tau)]$ and $\underline V:=\sup_{\Tau\in\bT^i}\inf_{\rho\in\T}\E[U(\Rho(\tau),\tau)]$ in discrete time, where $U(s,t)$ is $\mathcal{F}_{s\vee t}$-measurable instead of $\mathcal{F}_{s\wedge t}$-measurable as is often assumed in the literature, $\T$ is the set of stopping times, and $\bT^i$ and $\bT^{ii}$ are sets of mappings from $\T$ to $\T$ satisfying certain non-anticipativity conditions. We convert the problems into a corresponding Dynkin game, and show that $\overline V=\underline V=V$, where $V$ is the value of the Dynkin game. We also get the optimal $\Rho\in\bT^{ii}$ and $\Tau\in\bT^i$ for $\overline V$ and $\underline V$ respectively.
  • We consider the problem of downlink channel estimation for millimeter wave (mmWave) MIMO-OFDM systems, where both the base station (BS) and the mobile station (MS) employ large antenna arrays for directional precoding/beamforming. Hybrid analog and digital beamforming structures are employed in order to offer a compromise between hardware complexity and system performance. Different from most existing studies that are concerned with narrowband channels, we consider estimation of wideband mmWave channels with frequency selectivity, which is more appropriate for mmWave MIMO-OFDM systems. By exploiting the sparse scattering nature of mmWave channels, we propose a CANDECOMP/PARAFAC (CP) decomposition-based method for channel parameter estimation (including angles of arrival/departure, time delays, and fading coefficients). In our proposed method, the received signal at the BS is expressed as a third-order tensor. We show that the tensor has the form of a low-rank CP decomposition, and the channel parameters can be estimated from the associated factor matrices. Our analysis reveals that the uniqueness of the CP decomposition can be guaranteed even when the size of the tensor is small. Hence the proposed method has the potential to achieve substantial training overhead reduction. We also develop Cramer-Rao bound (CRB) results for channel parameters, and compare our proposed method with a compressed sensing-based method. Simulation results show that the proposed method attains mean square errors that are very close to their associated CRBs, and presents a clear advantage over the compressed sensing-based method in terms of both estimation accuracy and computational complexity.
  • We consider the problem of uplink channel estimation for millimeter wave (mmWave) systems, where the base station (BS) and mobile stations (MSs) are equipped with large antenna arrays to provide sufficient beamforming gain for outdoor wireless communications. Hybrid analog and digital beamforming structures are employed by both the BS and the MS due to hardware constraints. We propose a layered pilot transmission scheme and a CANDECOMP/PARAFAC (CP) decomposition-based method for joint estimation of the channels from multiple users (i.e. MSs) to the BS. The proposed method exploits the sparse scattering nature of the mmWave channel and the intrinsic multi-dimensional structure of the multiway data collected from multiple modes. The uniqueness of the CP decomposition is studied and sufficient conditions for essential uniqueness are obtained. The conditions shed light on the design of the beamforming matrix, the combining matrix and the pilot sequences, and meanwhile provide general guidelines for choosing system parameters. Our analysis reveals that our proposed method can achieve a substantial training overhead reduction by employing the layered pilot transmission scheme. Simulation results show that the proposed method presents a clear advantage over a compressed sensing-based method in terms of both estimation accuracy and computational complexity.
  • We consider a financial market where stocks are available for dynamic trading, and European and American options are available for static trading (semi-static trading strategies). We assume that the American options are infinitely divisible, and can only be bought but not sold. In the first part of the paper, we work within the framework without model ambiguity. We first get the fundamental theorem of asset pricing (FTAP). Using the FTAP, we get the dualities for the hedging prices of European and American options. Based on the hedging dualities, we also get the duality for the utility maximization. In the second part of the paper, we consider the market which admits non-dominated model uncertainty. We first establish the hedging result, and then using the hedging duality we further get the FTAP. Due to the technical difficulty stemming from the non-dominancy of the probability measure set, we use a discretization technique and apply the minimax theorem.
  • We present a comprehensive structural characterization of two different highly pure nuclear graphites that compasses all relevant length scales from nanometers to sub-mm. This has been achieved by combining several experiments and neutron techniques: Small Angle Neutron Scattering (SANS), high-resolution Spin Echo SANS (SESANS) and neutron imaging. In this way it is possible to probe an extraordinary broad range of 6 orders of magnitude in length from microscopic to macroscopic length scales. The results reveal a fractal structure that extends from $ \sim $ 0.6 nm to 0.6 mm and has surface and mass fractal dimensions both very close to 2.5, a value found for percolating clusters and fractured ranked surfaces in 3D.
  • We consider multi-player stopping games in continuous time. Unlike Dynkin games, in our games the payoff of each player is revealed after all the players stop. Moreover, each player can adjust her own stopping strategy by observing other players' behaviors. Assuming the continuity of the payoff functions in time, we show that there always exists an $\eps$-Nash equilibrium in pure stopping strategies for any $\eps>0$.
  • We consider two-player non-zero-sum stopping games in discrete time. Unlike Dynkin games, in our games the payoff of each player is revealed after both players stop. Moreover, each player can adjust her own stopping strategy according to the other player's action. In the first part of the paper, we consider the game where players act simultaneously at each stage. We show that there exists a Nash equilibrium in mixed stopping strategies. In the second part, we assume that one player has to act first at each stage. In this case, we show the existence of a Nash equilibrium in pure stopping strategies.
  • On a filtered probability space $(\Omega ,\mathcal{F}, (\mathcal{F}_t)_{t\in[0,\infty]}, \mathbb{P})$, we consider the two-player non-zero-sum stopping game $u^i := \mathbb{E}[U^i(\rho,\tau)],\ i=1,2$, where the first player choose a stopping strategy $\rho$ to maximize $u^1$ and the second player chose a stopping strategy $\tau$ to maximize $u^2$. Unlike the Dynkin game, here we assume that $U(s,t)$ is $\mathcal{F}_{s\vee t}$-measurable. Assuming the continuity of $U^i$ in $(s,t)$, we show that there exists an $\epsilon$-Nash equilibrium for any $\epsilon>0$.
  • We consider a zero-sum continuous time stopping game in which the pay-off is revealed in the maximum of the two stopping times instead of the minimum, which is the case in Dynkin games.
  • We consider as given a discrete time financial market with a risky asset and options written on that asset and determine both the sub- and super-hedging prices of an American option in the model independent framework of ArXiv:1305.6008. We obtain the duality of results for the sub- and super-hedging prices. For the sub-hedging prices we discuss whether the sup and inf in the dual representation can be exchanged (a counter example shows that this is not true in general). For the super-hedging prices we discuss several alternative definitions and argue why our choice is more reasonable. Then assuming that the path space is compact, we construct a discretization of the path space and demonstrate the convergence of the hedging prices at the optimal rate. The latter result would be useful for numerical computation of the hedging prices. Our results generalize those of ArXiv:1304.3574 to the case when static positions in (finitely many) European options can be used in the hedging portfolio.
  • We consider the optimal stopping problem $v^{(\eps)}:=\sup_{\tau\in\mathcal{T}_{0,T}}\mathbb{E}B_{(\tau-\eps)^+}$ posed by Shiryaev at the International Conference on Advanced Stochastic Optimization Problems organized by the Steklov Institute of Mathematics in September 2012. Here $T>0$ is a fixed time horizon, $(B_t)_{0\leq t\leq T}$ is the Brownian motion, $\eps\in[0,T]$ is a constant, and $\mathcal{T}_{\eps,T}$ is the set of stopping times taking values in $[\eps,T]$. The solution of this problem is characterized by a path dependent reflected backward stochastic differential equations, from which the continuity of $\eps \to v^{(\eps)}$ follows. For large enough $\eps$, we obtain an explicit expression for $v^{(\eps)}$ and for small $\eps$ we have lower and upper bounds. The main result of the paper is the asymptotics of $v^{(\eps)}$ as $\eps\searrow 0$. As a byproduct, we also obtain L\'{e}vy's modulus of continuity result in the $L^1$ sense.
  • An important assumption in the work on testing for structural breaks in time series consists in the fact that the model is formulated such that the stochastic process under the null hypothesis of "no change-point" is stationary. This assumption is crucial to derive (asymptotic) critical values for the corresponding testing procedures using an elegant and powerful mathematical theory, but it might be not very realistic from a practical point of view. This paper develops change point analysis under less restrictive assumptions and deals with the problem of detecting change points in the marginal variance and correlation structures of a non-stationary time series. A CUSUM approach is proposed, which is used to test the "classical" hypothesis of the form $H_0: \theta_1=\theta_2$ vs. $H_1: \theta_1 \not =\theta_2$, where $\theta_1$ and $\theta_2$ denote second order parameters of the process before and after a change point. The asymptotic distribution of the CUSUM test statistic is derived under the null hypothesis. This distribution depends in a complicated way on the dependency structure of the nonlinear non-stationary time series and a bootstrap approach is developed to generate critical values. The results are then extended to test the hypothesis of a {\it non relevant change point}, i.e. $H_0: | \theta_1-\theta_2 | \leq \delta$, which reflects the fact that inference should not be changed, if the difference between the parameters before and after the change-point is small. In contrast to previous work, our approach does neither require the mean to be constant nor - in the case of testing for lag $k$-correlation - that the mean, variance and fourth order joint cumulants are constant under the null hypothesis. In particular, we allow that the variance has a change point at a different location than the auto-covariance.
  • We consider the fundamental theorem of asset pricing (FTAP) and hedging prices of options under non-dominated model uncertainty and portfolio constrains in discrete time. We first show that no arbitrage holds if and only if there exists some family of probability measures such that any admissible portfolio value process is a local super-martingale under these measures. We also get the non-dominated optional decomposition with constraints. From this decomposition, we get duality of the super-hedging prices of European options, as well as the sub- and super-hedging prices of American options. Finally, we get the FTAP and duality of super-hedging prices in a market where stocks are traded dynamically and options are traded statically.
  • We show that the results of ArXiv:1305.6008 on the Fundamental Theorem of Asset Pricing and the super-hedging theorem can be extended to the case in which the options available for static hedging (\emph{hedging options}) are quoted with bid-ask spreads. In this set-up, we need to work with the notion of \emph{robust no-arbitrage} which turns out to be equivalent to no-arbitrage under the additional assumption that hedging options with non-zero spread are \emph{non-redundant}. A key result is the closedness of the set of attainable claims, which requires a new proof in our setting.
  • We investigate the behavior of the Generalized Likelihood Ratio Test (GLRT) (Fan, Zhang and Zhang [Ann. Statist. 29 (2001) 153-193]) for time varying coefficient models where the regressors and errors are non-stationary time series and can be cross correlated. It is found that the GLRT retains the minimax rate of local alternative detection under weak dependence and non-stationarity. However, in general, the Wilks phenomenon as well as the classic residual bootstrap are sensitive to either conditional heteroscedasticity of the errors, non-stationarity or temporal dependence. An averaged test is suggested to alleviate the sensitivity of the test to the choice of bandwidth and is shown to be more powerful than tests based on a single bandwidth. An alternative wild bootstrap method is proposed and shown to be consistent when making inference of time varying coefficient models for non-stationary time series.
  • We investigate the behavior of Fourier transforms for a wide class of nonstationary nonlinear processes. Asymptotic central and noncentral limit theorems are established for a class of nondegenerate and degenerate weighted $V$-statistics through the angle of Fourier analysis. The established theory for $V$-statistics provides a unified treatment for many important time and spectral domain problems in the analysis of nonstationary time series, ranging from nonparametric estimation to the inference of periodograms and spectral densities.