WAPLS.predict.w | R Documentation |
WA-PLS predict function
WAPLS.predict.w(WAPLSoutput, fossil_taxa)
WAPLSoutput |
The output of the |
fossil_taxa |
Fossil taxa abundance data to reconstruct past climates, each row represents a site to be reconstructed, each column represents a taxon. |
A list of the reconstruction results. Each element in the list is described below:
fit
The fitted values using each number of components.
nPLS
The total number of components extracted.
WAPLS.w
## Not run: # Load modern pollen data modern_pollen <- read.csv("/path/to/modern_pollen.csv") # Extract taxa taxaColMin <- which(colnames(modern_pollen) == "taxa0") taxaColMax <- which(colnames(modern_pollen) == "taxaN") taxa <- modern_pollen[, taxaColMin:taxaColMax] # Load reconstruction data Holocene <- read.csv("/path/to/Holocene.csv") taxaColMin <- which(colnames(Holocene) == "taxa0") taxaColMax <- which(colnames(Holocene) == "taxaN") core <- Holocene[, taxaColMin:taxaColMax] ## Train fit_Tmin <- fxTWAPLS::WAPLS.w(taxa, modern_pollen$Tmin, nPLS = 5) fit_f_Tmin <- fxTWAPLS::WAPLS.w( taxa, modern_pollen$Tmin, nPLS = 5, usefx = TRUE, fx_method = "bin", bin = 0.02 ) fit_Tmin2 <- fxTWAPLS::WAPLS.w2(taxa, modern_pollen$Tmin, nPLS = 5) fit_f_Tmin2 <- fxTWAPLS::WAPLS.w2( taxa, modern_pollen$Tmin, nPLS = 5, usefx = TRUE, fx_method = "bin", bin = 0.02 ) ## Predict fossil_Tmin <- fxTWAPLS::WAPLS.predict.w(fit_Tmin, core) fossil_f_Tmin <- fxTWAPLS::WAPLS.predict.w(fit_f_Tmin, core) fossil_Tmin2 <- fxTWAPLS::WAPLS.predict.w(fit_Tmin2, core) fossil_f_Tmin2 <- fxTWAPLS::WAPLS.predict.w(fit_f_Tmin2, core) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.