Description Usage Arguments Details Value Author(s) References See Also Examples
Robust estimator of the Hurst parameter of a fractional Brownian possibly contaminated by additive outliers and/or an additive noise.
1 2 3 |
fbm |
data |
nma |
name of the filter used for filtering the data. See |
M1 |
Minimum value of the dilatation factor. Default is |
M2 |
Maximum value of the dilatation factor. Default is |
method |
Type of the discrete variations method. |
par |
Parameters depending on |
llplot |
If true a plot of \log(U_n^{a^m}) against \log(m) for m=M_1,...,M_2 is produced. |
An estimate of the Hurst exponent parameter is provided without estimating the scaling coefficient C and σ (parameter related to an additive noise). The standard method ST is based on filtering the data with dilated versions of the initial filter (whose name is nma
). Other methods are improvements. Methods TM and Q are based on trimmed-means and sample quantiles respectively. Methods B0 and B1 exploit the fact that the contamination is a Brownian motion or a Gaussian white noise. Other methods are combinations of the two last classes. See Achard and Coeurjolly (2009) for more details.
Returns the Hurst parameter estimate
J.-F. Coeurjolly
S. Achard and J.-F. Coeurjolly (2009). Discrete variations of the fractional Brownian in the presence of outliers and an additive noise. Submitted
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | n<-10000;H<-.8
## no
z<-perturbFBM(n,H,type="no",plot=FALSE)
dvFBM(z,method="ST")
dvFBM(z,method="TM",par=list(beta1=.1,beta2=.1))
dvFBM(z,method="B0-Q",par=list(vecp=.5,vecc=1))
dvFBM(z,method="B1-ST")
## AO
z<-perturbFBM(n,H,type="AO",SNR=-20,plot=FALSE)
dvFBM(z,nma="d4",method="ST")
dvFBM(z,nma="d4",method="TM",par=list(beta1=.1,beta2=.1))
## B0
z<-perturbFBM(n,H,type="B0",SNR=0,plot=FALSE)
dvFBM(z,M2=10,method="ST")
dvFBM(z,M2=10,method="B0-ST")
## B1
z<-perturbFBM(n,H,type="B1",SNR=0,plot=FALSE)
dvFBM(z,method="ST")
dvFBM(z,method="B1-ST")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.