b <- brick(system.file("external/rlogo.grd", package="raster"))
layerStats(b, 'pearson')
layerStats(b, 'cov')
# weigh by column number
w <- init(b, v='col')
layerStats(b, 'weighted.cov', w=w)
Run the code above in your browser using DataCamp Workspace