create a mask of NA values of bioclimatic variables from WorldClim in raster format in R

Question

I have downloaded the 19 bioclimatic variables from Bioclim and I want to create a unique mask for all the 19 variables. The mask should exclude the pixels that have NA value in any of the 19 variables. This procedure is needed for further SDM analysis. I did this with a stack and worked, but it is not working for the Bioclim rasterBrick. In R:

library(raster
bioclim<-raster::getData('worldclim',var='bio', res=10)
s<-sum(bioclim)
MASK[!is.na(s)] <- 1
plot(MASK)

Answer 1

Sticking to your approach, your code worked for me after inserting one line to actually define the "MASK" object.

library(raster)
bioclim<-raster::getData('worldclim',var='bio', res=10)
s<-sum(bioclim)

# create MASK object with NA or other masking value
MASK <- setValues(s[[1]], NA)

MASK[!is.na(s)] <- 1
plot(MASK)

A little simpler, you could use:

MASK2 <- !is.na(s)
plot(MASK2)