Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/BiCopVuongClarke.r

Based on the Vuong and Clarke tests this function computes a goodness-of-fit score for each bivariate copula family under consideration. For each possible pair of copula families the Vuong and the Clarke tests decides which of the two families fits the given data best and assigns a score—pro or contra a copula family—according to this decision.

1 2 | ```
BiCopVuongClarke(u1, u2, familyset=NA,
correction=FALSE, level=0.05)
``` |

`u1,u2` |
Data vectors of equal length with values in [0,1]. |

`familyset` |
An integer vector of bivariate copula families under consideration, i.e., which are compared in the goodness-of-fit test.
If |

`correction` |
Correction for the number of parameters.
Possible choices: |

`level` |
Numerical; significance level of the tests (default: |

The Vuong as well as the Clarke test compare two models against each other and
based on their null hypothesis, allow for a statistically significant decision among the two models
(see the documentations of `CDVineVuongTest`

and `CDVineClarkeTest`

for descriptions of the two tests).
In the goodness-of-fit test proposed by Belgorodski (2010) this is used for bivariate copula selection.
It compares a model 0 to all other possible models under consideration.
If model 0 is favored over another model, a score of "+1" is assigned and similarly a score of "-1" if the other model is determined to be superior.
No score is assigned, if the respective test cannot discriminate between two models.
Both tests can be corrected for the numbers of parameters used in the copulas.
Either no correction (`correction = FALSE`

), the Akaike correction (`correction = "Akaike"`

)
or the parsimonious Schwarz correction (`correction = "Schwarz"`

) can be used.

The models compared here are bivariate parametric copulas and we would like to determine which family fits the data better than the other families.
E.g., if we would like to test the hypothesis that the bivariate Gaussian copula fits the data best, then we
compare the Gaussian copula against all other copulas under consideration.
In doing so, we investigate the null hypothesis "The Gaussian copula fits the data better than all other copulas under consideration",
which corresponds to *k-1* times the hypothesis "The Gaussian copula *C_j* fits the data better than copula *C_i*" for all *i=1,...,k, i!=j*, where *k* is the number
of bivariate copula families under consideration (length of `familyset`

).
This procedure is done not only for one family but for all families under consideration, i.e., two scores, one based on the Vuong and one based on the Clarke test,
are returned for each bivariate copula family.
If used as a goodness-of-fit procedure, the family with the highest score should be selected.

For more and detailed information about the goodness-of-fit test see Belgorodski (2010).

A matrix with Vuong test scores in the first and Clarke test scores in the second row. Column names correspond to bivariate copula families (see above).

Ulf Schepsmeier, Eike Brechmann, Natalia Belgorodski

Belgorodski, N. (2010) Selecting pair-copula families for regular vines with application to the multivariate analysis of European stock market indices Diploma thesis, Technische Universitaet Muenchen. http://mediatum.ub.tum.de/doc/1079284/1079284.pdf.

Clarke, K. A. (2007). A Simple Distribution-Free Test for Nonnested Model Selection. Political Analysis, 15, 347-363.

Vuong, Q. H. (1989). Ratio tests for model selection and non-nested hypotheses. Econometrica 57 (2), 307-333.

`BiCopGofKendall`

, `CDVineVuongTest`

, `CDVineClarkeTest`

,

`BiCopSelect`

1 2 3 4 5 6 7 8 9 10 11 | ```
# simulate from a t-copula
dat = BiCopSim(500,2,0.7,5)
# apply the test for families 1-10
## Not run:
vcgof = BiCopVuongClarke(dat[,1],dat[,2],familyset=c(1:10))
# display the Vuong test scores
vcgof[1,]
## End(Not run)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.