[求助] 用R 模拟copula函数的随机数的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

x <- mvdc(fgmCopula(c(.2,-.2,-.4,0), dim = 3), c("pareto", "pareto","pareto"), list(list(2, 2), list(2,2),list(2,2)))
x.samp <- rmvdc(x, 100)


上面的模型装了VGAM包和copula包后可以运行,  模型dim是固定的, 但是如果dim是调用一组随机数的话该怎么办?

上面的模拟和我要做的模拟区别在于: 
1. fgmCopula 里面的参数问题, 当 dim 是随机数的时候, c(.2,-.2,-.4,0) 这一部分怎么确定? 括号里面的参数有2^dim-dim-1 个, 如3维时, 参数对应的是{1,2},{1,3},{2,3},{1,2,3}, 现在的想法是只要两两之间的系数不为0, 设置规律: 2^-|i-j|

2.  mvdc() 参数设置问题, 这个和第一个问题差不多.

困扰我很久了, 自己才刚开始学, 实在太菜, 希望有知道的虫友帮帮忙, 不胜感激!

[此贴子已经被作者于2009-4-1 18:55:55编辑过]

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：Copula opula 随机数 函数 模拟 Copula 随机数

555,没人知道吗?

d<-5   #维数
L<-list(2,2)     #子列表
e<-list(L)        #初始的大列表，其中只有一个list（2，2）
 
for(j in 2:d){
e[[j]]<-L
}

#fgmcopula函数中的参数设置
para<-rep(0,2^d-d-1)
k<-1
for(i in 1:(d-1)){
 for(j in (i+1):d){
    para[k]<-0.5^(abs(i-j))
    k<-k+1
  }
}

#X的随机数的产生
x <- mvdc(fgmCopula(para, d),
rep("pareto",d), e)
x.samp <- rmvdc(x, 100)

问题出在给d的赋值上,当d小于等于3是可以做,但是当d>3时就会报错

报错的内容是:

错误于validObject(.Object) :
  invalid class "fgmCopula" object: Bad vector of parameters

而当d=4是,如果将fgmCopula(para, d),写成 fgmCopula(c(0.5,0.5,0.5,0.5,0.5,0.5,0,0,0,0,0),4)
运行可以通过.

上面的问题怎么解决啊????

终于解决了

你好 我想下载一个copula程序包，可是一直找不到，请问要怎么做

limitily 发表于 2009-4-9 10:13
终于解决了

我想下载一个copula程序包，可是一直找不到，请问要怎么做

limitily 发表于 2009-4-1 19:01
d&lt;-5&nbsp;&nbsp; #维数L&lt;-list(2,2)&nbsp;&nbsp;&nbsp;&nbsp; #子列表e&lt;-list(L)&nbsp;&nbsp;&nbsp; ...

我要产生一个gumbel  copula 的随机数对  可是怎么找copula的程序包

mathyanzhen 发表于 2012-12-7 20:51
我要产生一个gumbel  copula 的随机数对  可是怎么找copula的程序包

你好，我要产生一个Clayton copula 的随机数，怎么产生？想请教一下，谢谢

青春之旅 发表于 2015-4-8 17:10
你好，我要产生一个Clayton copula 的随机数，怎么产生？想请教一下，谢谢

查查 VGAM 这个package

青春之旅 发表于 2015-4-8 17:10
你好，我要产生一个Clayton copula 的随机数，怎么产生？想请教一下，谢谢

查查 VGAM 这个package