我试图将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;)中:标准差为零。
我感谢您的帮助