通过循环列来拟合双变量copula

时间:2014-10-23 16:43:46

标签: r

我试图将copula模型拟合为双变量数据,u2是一个向量,u1的每一列是一个包含2个变量的数据框。这是我的代码

library(copula)
normal.cop <- normalCopula(c(0.6,0.3,0.6),dim=3,dispstr="un")
n=100
set.seed(2)
Dataset=rcopula(normal.cop, n)
u <- apply(Dataset, 2, rank) / (n + 1)
u1=u[,2:3]
u2=u[,1]

for(m in 1:ncol(u1)) 
{ 
  LIST[[m]]<-fitCopula(normal.cop1, cbind(u2,m), method="mpl") 
  LIST[[m]]=coef(m)
} 

不幸的是,这不起作用,我收到了这个错误:

  

特征(mat)中的错误:&#39; x&#39;

中的无限或缺失值

另外:警告信息:

  

在cor(x,method =&#34; kendall&#34;)中:标准差为零。

我感谢您的帮助

0 个答案:

没有答案