View source: R/NonparTrawlEstimation.R
c4est | R Documentation |
This function estimates the fourth cumulant of the trawl process.
c4est(data, Delta)
data |
The data set used to estimate the fourth cumulant |
Delta |
The width Delta of the observation grid |
According to Sauri and Veraart (2022), estimator based on X_0, X_{Δ_n}, …, X_{(n-1)Δ_n} is given by
\hat c_4(L')=RQ_n/\hat a(0),
where
RQ_n:=\frac{1}{√{2 nΔ_{n}}} ∑_{k=0}^{n-2}(X_{(k+1)Δ_n}-X_{kΔ_n})^4,
and
\hat a(0)=\frac{1}{2Δ_{n}n} ∑_{k=0}^{n-2}(X_{(k+1)Δ_n}-X_{kΔ_n})^{2}.
The function returns the estimated fourth cumulant of the Levy seed: \hat c_4(L').
##Simulate a trawl process ##Determine the sampling grid my_n <- 1000 my_delta <- 0.1 my_t <- my_n*my_delta ###Choose the model parameter #Exponential trawl function: my_lambda <- 2 #Poisson marginal distribution trawl my_v <- 1 #Set the seed set.seed(123) #Simulate the trawl process Poi_data<-ambit::sim_weighted_trawl(my_n, my_delta, "Exp", my_lambda, "Poi", my_v)$path #Estimate the fourth cumulant of the trawl process c4est(Poi_data, my_delta)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.