全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
8380 9
2009-03-29

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) 这一部分怎么确定? 括号里面的参数有2dim-dim-1 个, 如3维时, 参数对应的是{1,2},{1,3},{2,3},{1,2,3}, 现在的想法是只要两两之间的系数不为0, 设置规律: 2-|i-j|

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

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

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

二维码

扫码加我 拉你入群

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

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

全部回复
2009-3-29 21:54:00
555,没人知道吗?
二维码

扫码加我 拉你入群

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

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

2009-4-1 19:01:00

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)
运行可以通过.

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



二维码

扫码加我 拉你入群

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

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

2009-4-9 10:13:00

终于解决了

二维码

扫码加我 拉你入群

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

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

2012-5-25 11:55:51
你好 我想下载一个copula程序包,可是一直找不到,请问要怎么做
二维码

扫码加我 拉你入群

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

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

2012-12-7 20:50:33
limitily 发表于 2009-4-9 10:13
终于解决了
我想下载一个copula程序包,可是一直找不到,请问要怎么做
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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