全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
4888 4
2014-10-15

在R中用fitCopula()估计二元正态copula的参数,

cop.nor<- normalCopula(param=c(1), dim = 2, dispstr = "ex")

para.cop<-fitCopula(cop.nor,cbind(x,y),method="mpl",start=NULL,lower=NULL,upper=NULL,optim.control = list(NULL), optim.method = "BFGS", estimate.variance = TRUE)

得到以下结果:


fitCopula() estimation based on 'maximum pseudo-likelihood'

and a sample of size 501.

       Estimate Std. Error z value Pr(>|z|)   

rho.1 0.9917392  0.0007053    1406   <2e-16 ***

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

The maximized loglikelihood is  977

Optimization converged

Number of loglikelihood evaluations:

function gradient

      38        4




现在想把其中的参数提取出来,

最开始用$,结果被告知“para.cop$param : $ operator not defined for this S4 class”,因此换成了@,如下:


cop.nor<- normalCopula(param= para.cop@param, dim = 2, dispstr = "ex")


结果被告知:错误于stopifnot((pdim <- length(param)) >= 1, is.numeric(param)) :

"param"槽名不存在于"fitCopula"类别对象中,

现向各位大神请教:怎样才能顺利提取参数?









二维码

扫码加我 拉你入群

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

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

全部回复
2014-10-15 22:18:52
试试coef(para.cop)
二维码

扫码加我 拉你入群

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

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

2014-10-16 10:33:55
掘墓者 发表于 2014-10-15 22:18
试试coef(para.cop)
跪谢大神啊~~~
二维码

扫码加我 拉你入群

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

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

2014-10-16 18:28:29
写得真好                                                                        
                                       
                                                     
                                                     
                                             
                                                                 
                                                                                 
二维码

扫码加我 拉你入群

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

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

2015-10-10 21:06:53
请问,Number of loglikelihood evaluations 和 function gradient 是什么意思?
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群