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编辑过]