Question

使用R中的Cshapes包，我想创建一个矩阵列表，用于衡量每年是否有两个国家的邻居。

install.packages("cshapes")

运行一年的代码（这里是1990年）可以正常运行：

wmat <- distmatrix(as.Date("1990-1-1"), type="mindist", tolerance=0.5, useGW=FALSE)

这给出了一个具有以下结构的矩阵：

        A   B   C   D
1   A   0   0   210 0
2   B   0   0   637 305
3   C   210 637 0   73
4   D   0   305 73  0

在下一步中，我将两个国家/地区之间0距离的所有组合设置为1，将所有其他组合设置为0，并将对角线设置为0：

wmat[wmat>0]<-5
wmat[wmat==0]<-1
wmat[wmat==5]<-0
diag(wmat)<-0

这给了我以下矩阵：

        A   B   C   D
1   A   0   1   0   1
2   B   1   0   0   0
3   C   0   0   0   0
4   D   1   0   0   0

我很难做的是自动创建1960年至2014年之间的所有年份的矩阵，对每年进行校正，并将结果存储到矩阵列表中，我可以在相应的年份内调用每个矩阵。

非常欢迎任何意见。

Answer 1

你可以尝试

lst <- lapply(1960:2014, function(x) {
           wmat <- distmatrix(as.Date(paste0(x, '-1-1')), 
             type="mindist", tolerance=0.5, useGW=FALSE)
           wmat[wmat>0]<-5
           wmat[wmat==0]<-1
           wmat[wmat==5]<-0
           diag(wmat)<-0
           wmat 

    }
  )

在R中创建距离矩阵列表

1 个答案: