This R package implements the generalized additive modeling framework for bivariate copulas introduced by Vatter and Chavez-Demoulin (2015) and its extension to Pair-Copula Constructions (or Vine Copulas) by Vatter and Nagler (2017). It includes functions for parameter estimation, model selection, simulation, and visualization. The package is still under development. Please see the API documentation for a detailed description of all functions.

You can install:

the stable release on CRAN:

`install.packages("gamCopula")`

the latest development version:

`devtools::install_github("tvatter/gamCopula")`

Below, we list most functions and features you should know about. As usual in copula models, data are assumed to be serially independent and lie in the unit hypercube.

`gamBiCop`

: Creates a GAM bivariate copula by specifying a family and model, namely a`gamObject`

as return by the`gam`

function from the`mgcv`

package. Returns an object of class`gamBiCop`

. The class has the following methods:`print`

,`summary`

: a brief or comprehensive overview of the bivariate copula, respectively.`plot`

: plot method based on`plot.gam`

from`mgcv`

.`logLik`

,`AIC`

,`BIC`

,`nobs`

: usual fit statistics.`EDF`

: Equivalent degrees of freedom for the components of the model.

`gamBiCopSimulate`

: Simulates from a bivariate GAM copula.`gamBiCopFit`

: Estimates parameters of a bivariate copula with a prespecified family. Returns an object of class`gamBiCop`

.`gamBiCopSelect`

: Estimates the parameters of a bivariate copula for a set of families and selects the best fitting model (using either AIC or BIC). Returns an object of class`gamBiCop`

.`gamBiCopPredict`

,`gamBiCopPDF`

,`gamBiCopCDF`

: Predict and PDF/CDF methods for the GAM copula model.

`gamVine`

: Creates a GAM vine copula model by specifying a tree structure and list of`gamBicop`

objects corresponding to each edge. Returns an object of class`gamVine`

. The class has the following methods:`print`

,`summary`

: a brief or comprehensive overview of the bivariate copula, respectively.`plot`

: plots based on`plot.gamBiCop`

.

`gamVineSimulate`

: Simulates from a GAM vine copula model.`gamVineSeqFit`

: Estimates the parameters of a GAM vine copula model with prespecified structure and families.`gamVineCopSelect`

: Estimates the parameters and selects the best family for a GAM vine copula model with prespecified structure matrix.`gamVineStructureSelect`

: Fits a GAM vine copula model assuming no prior knowledge. It selects the R-vine structure using Dissmann et al. (2013)’s method, estimates parameters for various families, and selects the best family for each pair.`RVM2GVC`

: converts an`RVineMatrix`

object from the`VineCopula`

package into a`gamVine`

In this package several bivariate copula families are included for bivariate and multivariate analysis using vine copulas. It provides functionality of elliptical (Gaussian and Student-t) as well as Archimedean (Clayton, Gumbel, Frank) copulas to cover a large range of dependence patterns. For the Clayton and Gumbel copula families, rotated versions are included to cover negative dependence as well.

A copula family: 1 Gaussian, 2 Student t, 5 Frank, 301 Double Clayton type I (standard and rotated 90 degrees), 302 Double Clayton type II (standard and rotated 270 degrees), 303 Double Clayton type III (survival and rotated 90 degrees), 304 Double Clayton type IV (survival and rotated 270 degrees), 401 Double Gumbel type I (standard and rotated 90 degrees), 402 Double Gumbel type II (standard and rotated 270 degrees), 403 Double Gumbel type III (survival and rotated 90 degrees), 404 Double Gumbel type IV (survival and rotated 270 degrees).

The following table shows the parameter ranges of bivariate copula families with parameters `par`

and `par2`

and internal coding `family`

:

Copula family | `family` |
`par` |
`par2` |
---|---|---|---|

Gaussian | `1` |
`(-1, 1)` |
- |

Student t | `2` |
`(-1, 1)` |
`(2,Inf)` |

Double Clayton type I (standard and 90 degrees) | `301` |
`(-Inf, Inf)` |
- |

Double Clayton type II (standard and 270 degrees) | `302` |
`(-Inf, Inf)` |
- |

Double Clayton type III (survival and 90 degrees) | `303` |
`(-Inf, Inf)` |
- |

Double Clayton type IV (survival and 270 degrees) | `304` |
`(-Inf, Inf)` |
- |

Double Gumbel type I (standard and 90 degrees) | `401` |
`(-Inf, Inf)` |
- |

Double Gumbel type II (standard and 270 degrees) | `402` |
`(-Inf, Inf)` |
- |

Double Gumbel type III (survival and 90 degrees) | `403` |
`(-Inf, Inf)` |
- |

Double Gumbel type IV (survival and 270 degrees) | `404` |
`(-Inf, Inf)` |
- |

Frank | `5` |
`R \ {0}` |
- |

Vatter, T., Nagler, T. (2017)

*Generalized Additive Models for Pair-Copula Constructions*.

Preprint available at arXiv:1608.01593.

Vatter, T., Chavez-Demoulin, V. (2015).

*Generalized additive models for conditional dependence structures*.

Journal of Multivariate Analysis, 141: 147-167, http://dx.doi.org/10.1016/j.jmva.2015.07.003.