Function to get the frequency of the climate value, which will be used to
provide fx
correction for WA-PLS and TWA-PLS.
fx(x, bin, show_plot = FALSE)
Numeric vector with the modern climate values.
Binwidth to get the frequency of the modern climate values.
Boolean flag to show a plot of fx ~ x
.
Numeric vector with the frequency of the modern climate values.
cv.w
, cv.pr.w
, and
sse.sample
if (FALSE) {
# Load modern pollen data
modern_pollen <- read.csv("/path/to/modern_pollen.csv")
# Get the frequency of each climate variable fx
fx_Tmin <- fxTWAPLS::fx(modern_pollen$Tmin, bin = 0.02, show_plot = TRUE)
fx_gdd <- fxTWAPLS::fx(modern_pollen$gdd, bin = 20, show_plot = TRUE)
fx_alpha <- fxTWAPLS::fx(modern_pollen$alpha, bin = 0.002, show_plot = TRUE)
}