fitCopula {copula}函数如何的使用？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

fitCopula {copula}的使用介绍，
fitCopula(copula, data, method = c("mpl", "ml", "itau", "irho"),
          start = NULL, lower = NULL, upper = NULL,
          optim.method = "BFGS", optim.control = list(maxit=1000),
          estimate.variance = TRUE, hideWarnings = TRUE)
关于data是这样介绍的：
-matrix of (pseudo-)observations (for "mpl" and "ml" with values necessarily in [0,1]) from the copula to be estimated, where n denotes the sample size and d the dimension. Consider applying the function pobs() first in order to obtain values in [0,1].
请问该解释应该怎么理解？？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享1 收藏3 回帖

关键词：Copula opula ITC TCO fit 如何

以二元正态COPULA为例
norm.cop=normalCopula(0.5,dim=2)#选择COPULA模型，并赋初始值，不是模型参数最终估计值，只是初始迭带点
fit=fitCopula(norm.cop,data=x)#假使数据集为X

X是初始数据通过边缘分布函数映射到(0,1)上的

huyiustc 发表于 2014-5-5 06:27
X是初始数据通过边缘分布函数映射到(0,1)上的

Consider applying the function pobs() first in order to obtain values in [0,1].
这里的pobs() 是一种映射吗？只要保证原数据经过映射取值在（0，1）就可以？而映射可以任意的？

原始数据X，边缘分布函数F（X），映射是u=F（X）

huyiustc 发表于 2014-5-5 09:45
原始数据X，边缘分布函数F（X），映射是u=F（X）

明白了，我还有个问题，如果只是知道两个序列的取值，不知道他们服从什么分布，这里的映射是两个序列的累计分布函数？比如x，y两个序列的累计分布函数ecdf(x),ecdf(y)吗？然后data取值是cbind(ecdf(x),cbind(y))?
还是可以自己设定？我看到一种方法这样写的，z是两序列的原始数据，然后使用的数据是
u <- (apply(x, 2, rank)-0.5) /n，这样写有什么依据吗？

jiaolitao 发表于 2014-5-5 11:27  明白了，我还有个问题，如果只是知道两个序列的取值，不知道他们服从什么分布，这里的映射是两个序列的累 ...

ECDF是用（RanK-1）/n还是RanK/N（忘了，反正是二者之一）作为经验分布，后者是用（RanK-0.5）/n

jiaolitao 发表于 2014-5-5 11:27  明白了，我还有个问题，如果只是知道两个序列的取值，不知道他们服从什么分布，这里的映射是两个序列的累 ...

ECDF是用（RanK-1）/n还是RanK/N（忘了，反正是二者之一）作为经验分布，后者是用（RanK-0.5）/n

错了，ECDF是用（RanK-1）/(n-1)这样会出现0,1而后面的方法不会

huyiustc 发表于 2014-5-5 12:21
错了，ECDF是用（RanK-1）/(n-1)这样会出现0,1而后面的方法不会

明白，非常感谢！