r - 面板数据的动态时间规整 (dtw)

Question

对于以下数据框，我使用以下代码计算了成本矩阵。最初，我创建了一个矩阵，其中 [i,j] 条目都是 10。然后我使用 dtw 函数填充下三角部分来计算距离。

如何将“distmx”矩阵的上三角部分设为空白。我的意思是我想删除矩阵“distmx”的上三角部分的 10。

对于这样的数据集（面板数据集），我如何选择查询和参考系列来绘制它？

Lines <- "Ctry  year   Carx 
A 2000   13
A 2001    8
A 2002   12
A 2003   16
A 2004   10
A 2005   13
B 2000   20
B 2001   24
B 2002   32
B 2003   32
B 2004   30
B 2005   36
C 2000   42
C 2001   46
C 2002   45
C 2003   56
C 2004   56
C 2005   58
"
DF <- read.table(text = Lines, header = TRUE)

成本矩阵

library(dtw)

ulke <- unique(DF$Ctry)
distmx <- matrix(rep(10, length(ulke)*length(ulke)), ncol=length(ulke))
colnames(distmx)<-rownames(distmx)<-ulke

for (i in 1:length(ulke)){
  for (j in i:length(ulke)){

    distmx[j,i]<-dtw(DF[DF$Ctry== ulke[i],"Carx"], DF[DF$Ctry== ulke[j],"Carx"],       distance.only=TRUE)$distance

 }

}