【自学笔记】基于R语言的copula函数重现期等值线绘制

1.加载相应的包和数据

设置工作路径,并加载相应的边缘分布的pdf(累积分布函数概率值),这期主要讲解重现期等值线,因此具体的边缘分布函数的拟合过程和copula函数的拟合过程就不具体介绍了。

setwd("F:/copula-work")
library(copula) 
library(VineCopula) 
library(MASS)
library(FAdist)
library(fitdistrplus) 
library(copBasic)
data<- read.csv('shuju.csv')
ls<-data[,1]#这里的数据为单变量分布函数F(X)
ld<-data[,2]#这里的数据为单变量分布函数F(Y)

2.绘制联合重现期等值线

根据定义,联合重现期是指P\left ( X>x \, \, \, or\, \, \, Y>y \right ),我们先将其转换为1-P\left ( X<x\, \,\, and\, \,\, Y<y\right ),那么就可以根据copBasic包里的joint.curvesCOP()函数的定义

你可能感兴趣的:(基于R语言的copula函数,r语言)