r - `freeParam` 拟合正常 copula 时出错 [R]

Question

我正在使用 R 版本 3.3.2 和包 copula 版本 0.999-15 来评估正常 copula 与我的数据的拟合。我的数据和代码是:

数据:https://www.dropbox.com/s/tdg8bfzmy4nd1dd/jumps.dat?dl=0

library(copula)
data <- read.csv(file="jumps.dat", head=F, sep="")

cop_model <- ellipCopula("normal", dim = 2)
m <- pobs(as.matrix(data))
fitCopula(cop_model, m, method = 'mpl')

运行代码后，我收到以下错误:

Error in `freeParam<-`(`*tmp*`, value = estimate) : the length of 'value' is not equal to the number of free parameters
Calls: fitCopula ... fitCopula.ml -> fitCopStart -> fitCopula.icor ->      freeParam<-
Execution halted

我不知道这里发生了什么。 Clayton 和 Gumbel 的合身非常好。在网上搜索类似的错误，我一无所获。阅读文档 ( https://www.rdocumentation.org/packages/copula/versions/0.999-15/topics/fitCopula ?) 以了解 ellipCopula 的一些特殊性，我找到了 posDef 的特定选项，但它根本没有返回任何解决方案。

Answer 1

老问题，但我发现了这个，所以会分享我的解决方案。

尝试运行以下命令，这是一个最小的工作示例:

library(copula)
print("-----------")
mycop <- ellipCopula("normal", dim=4)
data <- matrix(runif(400), nrow=4)
fitCopula(mycop, t(data))
print("-----------")

对我来说，如果我打开 R 并逐行输入，这会很好地工作，但如果我使用 Rscript 作为脚本运行，则会失败。解决方案是您还需要 library(methods)。

出于某种原因，这适用于 copula v0.999-v14，但被 v0.999-v16 破坏了。唉。