# Local methods example

    library(circularEV)
    require(plotly)

### Reading Data

    data(HsSP)
    data(drc)
    timeRange <- 54.5

    idx <- order(drc)
    drc <- drc[idx]
    Data <- HsSP[idx]
    set.seed(1234)
    Data <- Data + runif(length(Data), -0.0001, 0.0001)

    PlotData(Data = Data, drc = drc, thr = NULL, pointSize = 1, cex.axis = 15, cex.lab = 2,
        thrWidth = 2)

    PolarPlotData(Data = Data, drc = drc, thr = NULL, pointSize = 4, fontSize = 14,
        thrWidth = 4, ylim = c(0, max(Data)))

### Threshold selection

Grid values at which the estimation is performed:

    thetaVec <- 1:360

    thrResultMom <- ThrSelection(Data = Data, drc = drc, h = 60, b = 0.35, thetaGrid = thetaVec,
        EVIestimator = "Mom", useKernel = T, concent = 10, bw = 30, numCores = 2)$thr

    PlotData(Data = Data, drc = drc, thr = thrResultMom, pointSize = 1, cex.axis = 15,
        cex.lab = 2, thrWidth = 2)

    PolarPlotData(Data = Data, drc = drc, thr = thrResultMom, pointSize = 4, fontSize = 12,
        thrWidth = 4, ylim = c(0, max(Data)))

### Estimation of EVI, scale and T-year levels

    h <- 60
    useKernel <- TRUE
    concent <- 10
    movThr <- TRUE
    nBoot <- 30
    set.seed(1234)
    output <- LocalEstim(Data = Data, drc = drc, thr = thrResultMom, thetaGrid = thetaVec,
        nBoot = nBoot, EVIestimator = "Mom", h = h, useKernel = useKernel, concent = concent,
        movThr = movThr, TTs = c(100, 10000), timeRange = timeRange)

    RLBoot <- output$RLBoot

    PlotParamEstim(bootEstimates = output$xiBoot, thetaGrid = thetaVec, ylab = bquote(hat(xi)),
        alpha = 0.05, ylim = NULL, cex.axis = 15, cex.lab = 2, thrWidth = 2)

    PlotParamEstim(bootEstimates = output$sigBoot, thetaGrid = thetaVec, ylab = bquote(hat(sigma)),
        alpha = 0.05, ylim = NULL, cex.axis = 15, cex.lab = 2, thrWidth = 2)

    # 100-year level

    PlotRL(RLBootList = RLBoot, thetaGrid = thetaVec, Data = Data, drc = drc, TTs = c(100,
        10000), whichPlot = 1, alpha = 0.05, ylim = NULL, pointSize = 1, cex.axis = 15,
        cex.lab = 2, thrWidth = 2)

    PolarPlotRL(RLBootList = RLBoot, thetaGrid = thetaVec, Data = Data, drc = drc,
        TTs = c(100, 10000), whichPlot = 1, alpha = 0.05, ylim = NULL, pointSize = 4,
        fontSize = 12, lineWidth = 2)

    # 10000-year level

    PlotRL(RLBootList = RLBoot, thetaGrid = thetaVec, Data = Data, drc = drc, TTs = c(100,
        10000), whichPlot = 2, alpha = 0.05, ylim = NULL, pointSize = 1, cex.axis = 15,
        cex.lab = 2, thrWidth = 2)

    PolarPlotRL(RLBootList = RLBoot, thetaGrid = thetaVec, Data = Data, drc = drc,
        TTs = c(100, 10000), whichPlot = 2, alpha = 0.05, ylim = NULL, pointSize = 4,
        fontSize = 12, lineWidth = 2)