Learn R Programming

Sim.DiffProc (version 2.5)

PEOUG: Parametric Estimation of Hull-White/Vasicek (HWV) Gaussian Diffusion Models(Exact likelihood inference)

Description

Parametric estimation of Hull-White/Vasicek Model.

Usage

PEOUG(X, delta, starts = list(r= 1, theta= 1, sigma= 1), leve = 0.95)

Arguments

X
a numeric vector of the observed time-series values.
delta
the fraction of the sampling period between successive observations.
starts
named list. Initial values for optimizer.
leve
the confidence level required.

Value

  • coefCoefficients extracted from the model.
  • AICA numeric value with the corresponding AIC.
  • vcovA matrix of the estimated covariances between the parameter estimates in the linear or non-linear predictor of the model.
  • confintA matrix (or vector) with columns giving lower and upper confidence limits for each parameter. These will be labelled as (1-level)/2 and 1 - (1-level)/2.

Details

the Vasicek or Ornstein-Uhlenbeck model solves the stochastic differential equation : $$dX(t) = r * (theta - X(t)) * dt + sigma * dW(t)$$ It is ergodic for r > 0. We have also shown its exact conditional and stationary densities. In particular, the conditional density p(t,.|x) is the density of a Gaussian law with mean = theta+(x0-theta)*exp(-r*t) and variance = (sigma^2/(2*r))*(1-exp(-2*r*t)). R has the [dqpr]norm functions to evaluate the density, the quantiles, and the cumulative distribution or generate pseudo random numbers from the normal distribution.

See Also

PEABM Parametric Estimation of Arithmetic Brownian Motion, PEOUexp Explicit Estimators of Ornstein-Uhlenbeck Model, PEOU Parametric Estimation of Ornstein-Uhlenbeck Model, PEBS Parametric Estimation of model Black-Scholes.

Examples

Run this code
## example 1
## t0 = 0 ,T = 10
 data(DATA1)
 res <- PEOUG(DATA1,delta=0.01,starts=list(r=2,theta=0,sigma=1))
 res
 HWVF(N=1000,M=10,t0=0,T=10,x0=40,r=0.9979465,theta=16.49602,sigma=3.963486)
 points(seq(0,10,length=length(DATA1)),DATA1,type="l",lwd=3,col="blue")

Run the code above in your browser using DataLab