创建一个矩阵列表，以度量R中增长率的差异

Question

假设我有以下列表：

$`1990`
  cntry year growth
1     1 1990      1
2     2 1990      2
3     3 1990     14

$`1991`
  cntry year growth
4     1 1991      6
5     2 1991      5
6     3 1991      7

创建列表的代码：

cntry<-c(1,2,3,1,2,3)
year<-c(1990,1990,1990,1991,1991,1991)
growth<-c(1,2,14,6,5,7)

a<-data.frame(cntry,year,growth)
a.split<-split(a, as.factor(a$year))

我想要创建的是一个矩阵列表，用来衡量两国之间增长率的差异。

1. 在各自的矩阵中，只有当列中国家的增长率高于行中国家的增长率时，才应计算增长率的差异。
2. 如果列中的增长率等于或小于行中国家的增长率，则矩阵应显示值为0。
3. 如果组合中至少有一个国家具有NA值，则矩阵应在各自的字段中显示0。

结果如下：

$`1990`
  1 2  3
1 0 1 13
2 0 0 12
3 0 0  0

$`1991`
  1 2 3
1 0 0 1
2 1 0 2
3 0 0 0

任何帮助都是非常欢迎的！

Answer 1

使用outer

> lapply(a.split, function(x) t(outer(x$growth, x$growth, function(u,v) ifelse(u>v, u-v, 0))))
#$`1990`
#     [,1] [,2] [,3]
#[1,]    0    1   13
#[2,]    0    0   12
#[3,]    0    0    0

#$`1991`
#     [,1] [,2] [,3]
#[1,]    0    0    1
#[2,]    1    0    2
#[3,]    0    0    0