getHSAF function

High Spectral Anomaly Frequency (HSAF)

Compute High Spectral Anomaly Frequency (HSAF) from the given residual values.

getHSAF(resid, years, for.year = NULL)

Arguments

• resid: A vector of residual values generated with respect to the baseline signal
• years: A vector of years corrosponding to the residual values
• for.year: (optional) Calculate HSAF values only for the given year in this argument. If argument is omitted, HSAF values for all years are calculated.

Returns

Data frame containing two Columns:

year	First column, represents year
HSAF	Second column, represents HSAF values

Examples

# load sample data
data("sycamore")
x = sycamore

# get streamflow object for the sample data
x.streamflow = asStreamflow(x)

# prepare baseline signal 
x.bl = prepareBaseline(x.streamflow)

# HSAF
hsaf = getHSAF(x.bl$resid.sig, x.streamflow$data$year)