This paper proposes consistent and asymptotically Gaussian estimators for the drift, the diffusion coefficient and the Hurst exponent of the discretely observed fractional Ornstein-Uhlenbeck process. For the estimation of the drift, the results are obtained only in the case when 1/2 < H < 3/4. This paper also provides ready-to-use software for the R statistical environment based on the YUIMA package.
This paper addresses the problem of estimating drift parameter of the Ornstein - Uhlenbeck type process, driven by the sum of independent standard and fractional Brownian motions. The maximum likelihood estimator is shown to be consistent and asymptotically normal in the large-sample limit, using some recent results on the canonical representation and spectral structure of mixed processes.
In this paper, we will first give the numerical simulation of the sub-fractional Brownian motion through the relation of fractional Brownian motion instead of its representation of random walk. In order to verify the rationality of this simulation, we propose a practical estimator associated with the LSE of the drift parameter of mixed sub-fractional Ornstein-Uhlenbeck process, and illustrate the asymptotical properties according to our method of simulation when the Hurst parameter $H>1/2$.
We calculate the two-point correlation function <x(t2)x(t1)> for a subdiffusive continuous time random walk in a parabolic potential, generalizing well-known results for the single-time statistics to two times. A closed analytical expression is found for initial equilibrium, revealing a clear deviation from a Mittag-Leffler decay.
We study normal diffusive and subdiffusive processes in a harmonic potential (Ornstein-Uhlenbeck process) on a uniformly growing/contracting domain. Our starting point is a recently derived fractional Fokker-Planck equation, which covers both the case of Brownian diffusion and the case of a subdiffusive Continuous-Time Random Walk (CTRW). We find a high sensitivity of the random walk properties to the details of the domain growth rate, which gives rise to a variety of regimes with extremely different behaviors. At the origin of this rich phenomenology is the fact that the walkers still move while they wait to jump, since they are dragged by the deterministic drift arising from the domain growth. Thus, the increasingly long waiting times associated with the ageing of the subdiffusive CTRW imply that, in the time interval between two consecutive jumps, the walkers might travel over much longer distances than in the normal diffusive case. This gives rise to seemingly counterintuitive effects. For example, on a static domain, both Brownian diffusion and subdiffusive CTRWs yield a stationary particle distribution with finite width when a harmonic potential is at play, thus indicating a confinement of the diffusing particle. However, for a sufficiently fast growing/contracting domain, this qualitative behavior breaks down, and differences between the Brownian case and the subdiffusive case are found. In the case of Brownian particles, a sufficiently fast exponential domain growth is needed to break the confinement induced by the harmonic force; in contrast, for subdiffusive particles such a breakdown may already take place for a sufficiently fast power-law domain growth. Our analytic and numerical results for both types of diffusion are fully confirmed by random walk simulations.
Pooled testing (also known as group testing), where diagnostic tests are performed on pooled samples, has broad applications in the surveillance of diseases in animals and humans. An increasingly common use case is molecular xenomonitoring (MX), where surveillance of vector-borne diseases is conducted by capturing and testing large numbers of vectors (e.g. mosquitoes). The R package PoolTestR was developed to meet the needs of increasingly large and complex molecular xenomonitoring surveys but can be applied to analyse any data involving pooled testing. PoolTestR includes simple and flexible tools to estimate prevalence and fit fixed- and mixed-effect generalised linear models for pooled data in frequentist and Bayesian frameworks. Mixed-effect models allow users to account for the hierarchical sampling designs that are often employed in surveys, including MX. We demonstrate the utility of PoolTestR by applying it to a large synthetic dataset that emulates a MX survey with a hierarchical sampling design.
Alexandre Brouste
,Stefano M. Iacus
.
(2011)
.
"Parameter estimation for the discretely observed fractional Ornstein-Uhlenbeck process and the Yuima R package"
.
Alexandre Brouste
هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا