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

 }

}