大家好，我正在对以下数据进行EOF分析：https://www.dropbox.com/s/b37gzandqv0kdu9/gpcc_spi6_1901_2013.rda?dl=0
这是我创建的一个文件，其中包含某个区域的SPI，我想对它进行EOF分析，这是我到目前为止所做的：
require(utils)
library(colorRamps)
library(sinkr)
Datos <- load("gpcc_spi6_1901_2013.rda")
dim(spi_c)#18,31,1356
#I cut some of the years off
spi1 <- spi_c[ , ,205:1356] #31, 18, 1152 lon ,lat, months
#Then I only take the months from May-Oct for the whole series
reform <- array(spi1, dim = c(31,18,12,96))
verano <- reform[ , ,5:10, ]
#Come back to lon, lat, month format
r_verano <- array(verano, dim=c(31,18,576))
#I know you can only apply EOF to a matrix so I convert the 3D array
#Into a 2D matrix
arreglo <- sapply(1:dim(r_verano)[3], function(x)t(r_verano[,,x])) #558, 576
arreglo <- t(arreglo)
#EOF
neofs <- 20
#I wanted to compare the recursive and non recursive cases
Et <- eof(arreglo, recursive=FALSE, metho="irlba", nu=neofs)
Er <- eof(arreglo, recursive=TRUE, method = "irlba", nu=neofs)
但是现在我不知道如何绘制一个地图中的第一个EOF，然后是另一个地图中的第二个EOF，该地区是墨西哥，所以我想把海岸线放在里面，这样可以看得更清楚
另外，就我对EOF分析的理解而言，时间应该在列上，而行上的天文台却看着帮助它在我看来是另一种方式，这就是为什么我使用{{ 1}}但你能证实吗？感谢。