create rasters using pca components of bioclimatic variables
Andre,
The below is slightly modified from the example in
?raster::predict
library(raster)
logo <- brick(system.file("external/rlogo.grd", package="raster"))
sr <- sampleRegular(logo, min(100000, ncell(logo)))
pca <- prcomp(sr)
x <- predict(logo, pca, index=1:3)
plot(x)
Robert
On Thu, Nov 21, 2013 at 7:25 AM, Andre <pintodasilva.a at gmail.com> wrote:
Hi, I performed pca in 4 bioclimatic variables using princomp. ##stack rasters myExpl <- raster::stack (Bio01,Bio02,Bio03,Bio04) ##read as data.frame m <- as.data.frame(myExpl) ##pca (here I used formula to activate na.action) pca <- princomp(formula = ~., data = m, cor = TRUE, na.action=na.exclude) Now I would like to create rasters (with the same size as BIO01) for each component of the PCA to use them as input in biomod2. What is the best way to do this? Many thanks, Andre -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/create-rasters-using-pca-components-of-bioclimatic-variables-tp7585158.html Sent from the R-sig-geo mailing list archive at Nabble.com.
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo