On convergence of the inexact Rayleigh quotient iteration with MINRES

For the Hermitian inexact Rayleigh quotient iteration (RQI), we present a new general theory, independent of iterative solvers for shifted inner linear systems. The theory shows that the method converges at least quadratically under a new condition, called the uniform positiveness condition, that may allow the residual norm _ξk≥1${\xi }_k\ge 1$ of the inner linear system at outer iteration k+1$k+1$ and can be considerably weaker than the condition _ξk≤ξ<1${\xi }_k\le \xi <1$ with ξ$\xi$ a constant not near one commonly used in the literature. We consider the convergence of the inexact RQI with the unpreconditioned and tuned preconditioned MINRES methods for the linear systems. Some attractive properties are derived for the residuals obtained by MINRES. Based on them and the new general theory, we make a refined analysis and establish a number of new convergence results. Let ‖_rk‖$\Vert r_k\Vert$ be the residual norm of approximating eigenpair at outer iteration k$k$. Then all the available cubic and quadratic convergence results require _ξk=O(‖_rk‖)${\xi }_k=O\left(\Vert r_k\Vert \right)$ and _ξk≤ξ${\xi }_k\le \xi$ with a fixed ξ$\xi$ not near one, respectively. Fundamentally different from these, we prove that the inexact RQI with MINRES generally converges cubically, quadratically and linearly provided that _ξk≤ξ${\xi }_k\le \xi$ with a constant ξ<1$\xi <1$ not near one, _ξk=1−O(‖_rk‖)${\xi }_k=1-O\left(\Vert r_k\Vert \right)$ and _ξk=1−O(‖_rk‖²)${\xi }_k=1-O\left({\Vert r_k\Vert }^2\right)$, respectively. The new convergence conditions are much more relaxed than ever before. The theory can be used to design practical stopping criteria to implement the method more effectively. Numerical experiments confirm our results.     

Estimation for stochastic differential equations with a small diffusion coefficient

We consider a multidimensional diffusion X$X$ with drift coefficient b(_Xt,α)$b(X_t,\alpha )$ and diffusion coefficient εa(_Xt,β)$\epsilon a(X_t,\beta )$ where α$\alpha$ and β$\beta$ are two unknown parameters, while ε$\epsilon$ is known. For a high frequency sample of observations of the diffusion at the time points k/n$k/n$, k=1,…,n$k=1,\dots ,n$, we propose a class of contrast functions and thus obtain estimators of (α,β)$(\alpha ,\beta )$. The estimators are shown to be consistent and asymptotically normal when n→∞$n\to \infty$ and ε→0$\epsilon \to 0$ in such a way that ε⁻¹n^−ρ${\epsilon }^{-1}{n}^{-\rho }$ remains bounded for some ρ>0$\rho >0$. The main focus is on the construction of explicit contrast functions, but it is noted that the theory covers quadratic martingale estimating functions as a special case. In a simulation study we consider the finite sample behaviour and the applicability to a financial model of an estimator obtained from a simple explicit contrast function.     

Parametric inference for discretely observed multidimensional diffusions with small diffusion coefficient

We consider a multidimensional diffusion X$X$ with drift coefficient b(α,_Xt)$b(\alpha ,X_t)$ and diffusion coefficient ?σ(β,_Xt)$?\sigma (\beta ,X_t)$. The diffusion sample path is discretely observed at times _tk=kΔ$t_k=k\Delta$ for k=1…n$k=1\dots n$ on a fixed interval [0,T]$[0,T]$. We study minimum contrast estimators derived from the Gaussian process approximating X$X$ for small ?$?$. We obtain consistent and asymptotically normal estimators of α$\alpha$ for fixed Δ$\Delta$ and ?→0$?\to 0$ and of (α,β)$(\alpha ,\beta )$ for Δ→0$\Delta \to 0$ and ?→0$?\to 0$ without any condition linking ?$?$ and Δ$\Delta$. We compare the estimators obtained with various methods and for various magnitudes of Δ$\Delta$ and ?$?$ based on simulation studies. Finally, we investigate the interest of using such methods in an epidemiological framework.     

Frobenius groups and retract rationality

Let k$k$ be any field, G$G$ be a finite group acting on the rational function field k(_xg:g∈G)$k(x_g:g\in G)$ by h⋅_xg=_xhg$h\cdot x_g=x_{hg}$ for any h,g∈G$h,g\in G$. Define k(G)=k(_xg:g∈G)^G$k\left(G\right)=k{(x_g:g\in G)}^G$. Noether’s problem asks whether k(G)$k\left(G\right)$ is rational (= purely transcendental) over k$k$. A weaker notion, retract rationality introduced by Saltman, is also very useful for the study of Noether’s problem. We prove that, if G$G$ is a Frobenius group with abelian Frobenius kernel, then k(G)$k\left(G\right)$ is retract k$k$-rational for any field k$k$ satisfying some mild conditions. As an application, we show that, for any algebraic number field k$k$, for any Frobenius group G$G$ with Frobenius complement isomorphic to S_L2(_F5)$S{L}_2\left({\mathbb{F}}_5\right)$, there is a Galois extension field K$K$ over k$k$ whose Galois group is isomorphic to G$G$, i.e. the inverse Galois problem is valid for the pair (G,k)$(G,k)$. The same result is true for any non-solvable Frobenius group if k(_ζ8)$k\left({\zeta }_8\right)$ is a cyclic extension of k$k$.     

Regularity in the obstacle problem for parabolic non-divergence operators of Hörmander type

In this paper we continue the study initiated in [15] concerning the obstacle problem for a class of parabolic non-divergence operators structured on a set of vector fields X={_X1,…,_Xq}$X=\{X_1,\dots ,X_q\}$ in Rⁿ${\mathbf{R}}^n$ with C^∞$C^{\infty }$-coefficients satisfying Hörmander?s finite rank condition, i.e., the rank of Lie[_X1,…,_Xq]$\mathrm{Lie}[X_1,\dots ,X_q]$ equals n   at every point in Rⁿ${\mathbf{R}}^n$. In [15] we proved, under appropriate assumptions on the operator and the obstacle, the existence and uniqueness of strong solutions to a general obstacle problem. The main result of this paper is that we establish further regularity, in the interior as well as at the initial state, of strong solutions. Compared to [15] we in this paper assume, in addition, that there exists a homogeneous Lie group G=(Rⁿ,°,_δλ)$\mathbf{G}=({\mathbf{R}}^n,°,{\delta }_{\lambda })$ such that _X1,…,_Xq$X_1,\dots ,X_q$ are left translation invariant on G and such that _X1,…,_Xq$X_1,\dots ,X_q$ are _δλ${\delta }_{\lambda }$-homogeneous of degree one.     

On the removal lemma for linear systems over Abelian groups

In this paper we present an extension of the removal lemma to integer linear systems over abelian groups. We prove that, if the k$k$-determinantal of an integer (k×m)$(k\times m)$ matrix A$A$ is coprime with the order n$n$ of a group G$G$ and the number of solutions of the system Ax=b$Ax=b$ with _x1∈_X1,…,_xm∈_Xm$x_1\in X_1,\dots ,x_m\in X_m$ is o(n^m−k)$o\left(n^{m-k}\right)$, then we can eliminate o(n)$o\left(n\right)$ elements in each set to remove all these solutions.     

Invariance principles for generalized domains of semistable attraction

Let X,_X1,_X2,…$X,X_1,X_2,\dots$ be independent and identically distributed R^d${\mathbb{R}}^d$-valued random vectors and assume X$X$ belongs to the generalized domain of attraction of some operator semistable law without normal component. Then without changing its distribution, one can redefine the sequence on a new probability space such that the properly affine normalized partial sums converge in probability and consequently even in L^p$L^p$ (for some p>0$p>0$) to the corresponding operator semistable Lévy motion.     

On weighted Hilbert spaces and integration of functions of infinitely many variables

We study aspects of the analytic foundations of integration and closely related problems for functions of infinitely many variables _x1,_x2,…∈D$x_1,x_2,\dots \in D$. The setting is based on a reproducing kernel k$k$ for functions on D$D$, a family of non-negative weights _γu${\gamma }_u$, where u$u$ varies over all finite subsets of N$\mathbb{N}$, and a probability measure ρ$\rho$ on D$D$. We consider the weighted superposition K=_∑uγuku$K={\sum }_u{\gamma }_u{k}_u$ of finite tensor products _ku$k_u$ of k$k$. Under mild assumptions we show that K$K$ is a reproducing kernel on a properly chosen domain in the sequence space D^N$D^{\mathbb{N}}$, and that the reproducing kernel Hilbert space H(K)$H\left(K\right)$ is the orthogonal sum of the spaces H(_γuku)$H\left({\gamma }_u{k}_u\right)$. Integration on H(K)$H\left(K\right)$ can be defined in two ways, via a canonical representer or with respect to the product measure ρ^N${\rho }^{\mathbb{N}}$ on D^N$D^{\mathbb{N}}$. We relate both approaches and provide sufficient conditions for the two approaches to coincide.     

On functional limits of short- and long-memory linear processes with GARCH(1,1) noises

This paper considers the short- and long-memory linear processes with GARCH (1,1) noises. The functional limit distributions of the partial sum and the sample autocovariances are derived when the tail index α$\alpha$ is in (0,2)$(0,2)$, equal to 2, and in (2,∞)$(2,\infty )$, respectively. The partial sum weakly converges to a functional of α$\alpha$-stable process when α<2$\alpha <2$ and converges to a functional of Brownian motion when α≥2$\alpha \ge 2$. When the process is of short-memory and α<4$\alpha <4$, the autocovariances converge to functionals of α/2$\alpha /2$-stable processes; and if α≥4$\alpha \ge 4$, they converge to functionals of Brownian motions. In contrast, when the process is of long-memory, depending on α$\alpha$ and β$\beta$ (the parameter that characterizes the long-memory), the autocovariances converge to either (i) functionals of α/2$\alpha /2$-stable processes; (ii) Rosenblatt processes (indexed by β$\beta$, 1/2<β<3/4$1/2<\beta <3/4$); or (iii) functionals of Brownian motions. The rates of convergence in these limits depend on both the tail index α$\alpha$ and whether or not the linear process is short- or long-memory. Our weak convergence is established on the space of càdlàg functions on [0,1]$[0,1]$ with either (i) the _J1$J_1$ or the _M1$M_1$ topology (Skorokhod, 1956); or (ii) the weaker form S$S$ topology (Jakubowski, 1997). Some statistical applications are also discussed.     

An extremal problem for polynomials

We give a solution to an extremal problem for polynomials, which asks for complex numbers α₀,…,α_n${\alpha }_0,\dots ,{\alpha }_n$ of unit magnitude that minimise the largest supremum norm on the unit circle for all polynomials of degree n whose k  -th coefficient is either α_k${\alpha }_k$ or −α_k$-{\alpha }_k$.