R中栅格图层的两两相关性

Question

我需要在从WordClim数据库中提取的非洲大陆的19个栅格图层之间完成成对的Pearson相关性分析。我想检查哪些变量层与我的模型更相关/显著。为此，我尝试使用Raster包中的layerStats函数，但在执行后输出不包含数值，所有行和列显示为NA值。 以下是我的脚本。

#Loading raster files from WorldClim database
rastFiles<- list.files(pattern="bil")
a<-stack(rastFiles)

# Adjusting for African Continent
newext<-c(-20, 55, -35, 45)
Africa<-crop(a,newext)
Africa

#Correlation
cor<-layerStats(Africa,'pearson')

- Ricardo Adelino

---

如果您正在使用R，请添加相关标签。 - m69 ''snarky and unwelcoming''

抱歉，我忘记指定编程语言了，我使用的是R。 - Ricardo Adelino

Answer 1

在r中，只需使用以下代码，确保您已经使用na.rm=T来处理跨层的NAs：

library(raster)    
jnk=layerStats(raster_stack, 'pearson', na.rm=T)
corr_matrix=jnk$'pearson correlation coefficient'

- Lucas Fortini

---

你好，你有没有想过为什么我的layerStats函数会返回“Error: Failure during raster IO”错误？在叠加两个栅格图层之前，我必须更改其中一个的范围以匹配另一个。谢谢。 - Sam

“jnk$'mean'”这个是什么意思？我找不到相关信息。谢谢。 - parallax