本人采用R语言计算gumbel copula的Kendall 秩相关系数τ
data<-read.csv("C:/Users/pc/Desktop/1.csv",header=F)
gumbel.cop<- gumbelCopula ()
cg<-iTau(gumbel.cop,cor(cbind(data[1],data[2],data[3],data[4]), method=”kendall”))
计算结果如下:
V1 V2 V3 V4
V1 Inf 1.000882 1.014005 1.000000
V2 1.000882 Inf 1.708333 1.000000
V3 1.014005 1.708333 Inf 1.006209
V4 1.000000 1.000000 1.006209 Inf
但是Kendall 秩相关系数τ应该是在[-1,1]之间。
求哪位大神指导一下,谢谢