vine {pencopulaCond} | R Documentation |
"Estimating Non-Simplified Vine Copulas Using Penalized Splines"
Description
Estimating Non-Simplified Vine Copulas Using Penalized Splines
Usage
vine(data,d=2,d2=2,D=4,D3=6,lambda=c(100,50),type="Rvine",order.Dvine=FALSE,m=2,
cores=NULL,q=1,mod.cond=TRUE,max.iter=51,fix.lambda=FALSE,RVM=NULL,cal.cond=FALSE,
id=NULL,test.ind=FALSE,test.cond=2,lambda.search=FALSE,lam1.vec=NULL,lam2.vec=NULL)
Arguments
data |
'data' contains the data. 'data' has to be a matrix or a data.frame with two columns. |
d |
refers to the hierachy level of the marginal hierarchical B-spline for copulas in the first tree of the vine, default is d=2. |
d2 |
refers to the hierachy level of the marginal hierarchical B-spline for copulas in the second tree and in the following trees of the vine, default is d2=2. |
D |
referes to the maximum hierachy level for copulas in the first tree of the vine, default is D=4. If D<d, it follows D<-d. |
D3 |
referes to the maximum hierachy level for copulas in the second tree and in the following trees of the vine, default is D3=6. |
lambda |
Starting values for lambda, first start values for copulas in the first tree, second start value for copulas in the second tree and in the following trees of the vine, default is lambda=c(100,50). |
type |
Default is type="Rvine", fitting a regular vine copula. An alternative is type="Dvine", fitting a D-vine copula. |
order.Dvine |
Only relevant for type="Dvine". Indicating if the first level of the Dvine is ordered, default order.Dvine=TRUE. |
m |
Indicating the order of differences to be penalised. Default is "m=2". |
cores |
Default=NULL, the number of cpu cores used for parallel computing can be specified. |
q |
Degree of B-splines. Default is q=1. |
mod.cond |
Default=TRUE. If mod.cond=FALSE each pair-copula in the vine is estimated as simplified copula. The argument test.cond varies the test for the simplyfing assumption, which is imported from the R-package pacotest. |
max.iter |
maximum number of iteration, the default is max.iter=51. |
fix.lambda |
Default=FALSE, using the algorithm in the paper for estimating the optimal penalty parameter. If fix.lambda=TRUE, lambda is constant throughout the estimation. |
RVM |
Default=NULL. If RVM is a RVine-Matrix, this matrix determines the structure of the vine. |
cal.cond |
Default=FALSE. If cal.cond=TRUE each copula in the second tree and in the following trees of the vine is estimated as conditional copula. |
id |
Optional, one set id to any value. Especially important for simulations, starting with several starting values for lambda. |
test.ind |
Default=FALSE. If test.ind=TRUE, the fitted log-likelihood of each pair-copula is evaluated. If ("log.like"/"n"<0.001), where "n" is the sample size, the program set the corresponding pair copula as independence copula. We do not use this in our simulations or applications in the article. |
test.cond |
If test.cond=2, testType='ECORR' is chosen for the test of the simplyfing assumption as proposed in the article. There is an additional second test available in the R-package pactotest. testType="VI" is chosen with test.cond=1. |
lambda.search |
TRUE/FALSE, indicating if a search about several starting values for lambda should be performed. If search is selected, the starting value 'lambda' does not work anymore. |
lam1.vec |
Vector of candidate values for penalty parameter lambda for copulas in the first tree of the vine |
.
lam2.vec |
Vector of candidate values for penalty parameter lambda for copulas in the second tree and in the following trees of the vine |
.
Details
The calculation of the vine is done stepwise. The specifications in 'vine' are done for every paircopula in the vine with the identical specification. There is no option to change parameters for some pair-copulas.
Value
Returning a list containing
vine |
The estimated vine copula, an object of class 'pencopulaCond' |
log.like |
the estimated log-likelihood |
log.like.vec |
A vector with the estimated log.like.vec of each pair-copula |
AIC |
AIC value |
AIC.vec |
A vector with the estimated AIC of each pair-copula |
cAIC |
corrected AIC value |
cAIC.vec |
A vector with the estimated cAIC of each pair-copula |
d |
Used d |
d2 |
Used d2 |
D |
Used D |
D3 |
Used D3 |
order |
the used order of the first level (reported only for D-vines) |
S |
Sequence seq(1:(dim(data)[2])) |
N |
Number of observations, that is dim(data)[1] |
base |
Used basis function |
q |
Used degree of the B-spline basis |
no.cond.dens |
Estimated number of condtional copulas |
pca |
Indicating the used number of pca |
D.struc |
Used D.struc |
type |
Selected type of the vine copula |
VineMatrix |
VineMatrix, reported for type="Rvine" |
Author(s)
Christian Schellhase <cschellhase@wiwi.uni-bielefeld.de>
References
Flexible Copula Density Estimation with Penalized Hierarchical B-Splines, Kauermann G., Schellhase C. and Ruppert, D. (2013), Scandinavian Journal of Statistics 40(4), 685-705.
Estimating Non-Simplified Vine Copulas Using Penalized Splines, Schellhase, C. and Spanhel, F. (2017), Statistics and Computing.