关于copula包中的dcopula和dmvdc

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

这是关于绘制一组观测值的等值线命令，为何两部分绘制出来的图像差距那么大？

norm.cop<-ellipCopula(family="normal",param=0.7,dim=2)
norm.mvdc<-mvdc(copula=norm.cop,margins=c("norm","norm"),paramMargins=list(list(mean=0,sd=2),list(mean=1,sd=0.5)))
xx<-rMvdc(10000,norm.mvdc)
zz = HPDregionplot(cbind(xx[,1],xx[,2]),prob = 0.95,xlim=c(-5,5),ylim=c(-5,5),add=FALSE,col="blueviolet",lwd=2)
con.lev1<-zz[[1]]$level
cont.01<-contour(norm.mvdc,dmvdc,xlim=c(-5,5),ylim=c(-5,5),levels= con.lev1,draw = FALSE,add=FALSE)
#/////////////////////////////////////////
norm.cop<-ellipCopula(family="normal",param=0.7,dim=2)
nn<-rCopula(10000,norm.cop)
zz = HPDregionplot(cbind(nn[,1],nn[,2]),prob = 0.95,xlim=c(-5,5),ylim=c(-5,5),add=FALSE,col="blueviolet",lwd=2)
con.lev1<-zz[[1]]$level
cont.01<-contour(norm.cop,dcopula,xlim=c(-5,5),ylim=c(-5,5),levels= con.lev1,draw = FALSE,add=FALSE)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

111111111111111