liftHurst {liftLRD} | R Documentation |

## Performs (non-decimated) lifting based estimation of the Hurst exponent

### Description

The function exploits the linear relationship in wavelet energy per scale to estimate the long range dependence parameter of a irregular time series.

### Usage

```
liftHurst(x, grid = 1:length(x), model = "FGN", ntraj = 50,
tradonly = FALSE, cutoffs = 0, cut.fine = TRUE, efun = meanmo,
afun = idj, altype = 1, tail = TRUE, normalise = TRUE,
level = 0.05, bc = TRUE, vc = TRUE, jsc = TRUE, BHonly=TRUE,
verbose = FALSE, ...)
```

### Arguments

`x` |
A time series, possibly irregularly spaced. |

`grid` |
The time samples corresponding to the time series |

`model` |
The underlying model the time series is assumed to follow. Possible values are |

`ntraj` |
The number of lifting trajectories (bootstraps) used for the estimation of the Hurst exponent. |

`tradonly` |
A boolean variable, indicating whether the "traditional" minimum integral lifting trajectory should be used for the estimation ( |

`cutoffs` |
A vector indicating if the coarsest scales are to be removed when estimating the linear slope in the energy-scale relationship; for example if |

`cut.fine` |
A boolean variable indicating if the finest (default) or coarsest scales are to be removed when estimating the linear slope in the energy-scale relationship; see |

`efun` |
A function indicating which measure of wavelet energy to use. The default is |

`afun` |
A function indicating which measure of wavelet integral (scale) to use (corresponding to the x-axis of the energy-scale relationship). The default is |

`altype` |
An integer indicating which type of artificial levels to be used, see |

for more information.

`tail` |
A boolean variable indicating whether to amalgamate artificial levels with few coefficients, see |

for more information.

`normalise` |
A boolean variable indicating whether to normalise the detail coefficients by their individual (induced) standard deviations (computed using the diagonal of the lifting matrix). |

`level` |
The significance level for the bootstrap confidence interval of the Hurst exponent estimate. |

`bc` |
A boolean variable indicating whether bias-correction should take place or not, using the approximate chi^2 distribution of the wavelet energies. |

`vc` |
A boolean variable indicating whether a weighted linear regression should be used when estimating the Hurst exponent. If |

`jsc` |
A boolean variable indicating whether the slope of the log-linear relationship between the artificial scales and the log of the integrals should be computed and used to reweight the estimate of the Hurst exponent. |

`BHonly` |
A boolean variable indicating whether only the Hurst estimate should be returned. If |

`verbose` |
A boolean variable indicating whether informative status messages should be printed to the screen. |

`...` |
Any other optional arguments to the function |

### Details

Wavelet lifting is performed on a time series to convert it into a set of wavelet coefficients and corresponding lifting integrals, specific to when the data were "lifted" during the decomposition. The coefficients are then grouped into artificial levels, using the integrals to mimic the support of the wavelets in the classical wavelet setting, and therefore producing a notion of scale. The coefficients in each artificial level are then used to compute values of the wavelet energy for a particular level. The (slope of the) linear relationship between the scales and their energies is then used in computing an estimate of the Hurst exponent for the series. This procedure can be performed for multiple (random) lifting trajectories, each producing a slightly different estimate.

### Value

If `tradonly=TRUE`

, the function returns a matrix of dimension `length(cutoffs)`

x 2. The first column are the slopes of the regression fits for each cutoff, whereas the
second column are the corresponding estimates for the Hurst exponent.

If `tradonly=FALSE`

, the function returns a matrix of dimension `length(cutoffs)`

x 5. The first column are the slopes of the regression fits for each cutoff, where the
average is taken over the `ntraj`

randomly generated lifting trajectories. Similarly, the second column represents the average Hurst exponent for the cutoffs over all
lifting paths. The third column is the standard deviation of the `ntraj`

Hurst estimates through performing non-decimated lifting. The fourth and fifth columns are the lower and upper values of the bootstrap confidence interval of the Hurst exponent estimate.

If `BHonly=FALSE`

, the routine also returns the energies and scales (on a log scale) which are used in the regression to estimate the decay properties of the spectrum (for the last lifting trajectory), as well as the weights used in the regression (if `vc = TRUE`

). if `jsc = TRUE`

, the slope of the integral log-linear relationship is also returned.

### Author(s)

Marina Knight, Matt Nunes

### References

Knight, M. I, Nason, G. P. and Nunes, M. A. (2017) A wavelet lifting approach to long-memory estimation. *Stat. Comput.* **27** (6), 1453–1471. DOI 10.1007/s11222-016-9698-2.

For more details on the weighted linear regression and bias calculations, see e.g.

Veitch, D. and Abry, P. (1999) A Wavelet-Based Joint Estimator of the Parameters of Long-Range Dependence. *IEEE Trans. Info. Theory* **45** (3), 878–897.

### See Also

`artificial.levels`

, `Hfrombeta`

, `fwtnpperm`

### Examples

```
library(fracdiff)
# simulate a long range dependent time series
x<-fracdiff.sim(n = 200, d = 0.3)$series
# perform lifting-based estimation of the Hurst exponent
Hestx<-liftHurst(x, tradonly=TRUE)
# An example with missingness
x1<-fracdiff.sim(n = 500, d = 0.3)$series
gap1<-30:40
gap2<-77:90
gap3<-146:166
timeindex<-setdiff(1:500,c(gap1,gap2,gap3))
Hestx1<-liftHurst(x1[timeindex],grid=timeindex, tradonly=TRUE)
```

*liftLRD*version 1.0-9 Index]