intraBeta {AssetCorr}R Documentation

Parametric Approach of Botha and van Vuuren (2010)- Beta Distribution


The intra asset correlation will be estimated by fitting a beta distribution onto the default rate time series, then calculating the Value-at-Risk (VaR) of this beta distribution and fit it to the theoretical VaR of the Vasicek distribution. The correlation parameter will be backed out numerically. Additionally, bootstrap and jackknife corrections are implemented.


intraBeta(d, n, Quantile=0.999,B = 0, DB=c(0,0), JC = FALSE,
CI_Boot, type="bca", plot=FALSE)



a vector, containing the default time series of the sector.


a vector, containing the number of obligors at the beginning of the period over time.


a number, indicating the desired confidence level of the Value-at-Risk.


an integer, indicating how many bootstrap repetitions should be used for the single bootstrap corrected estimate.


a combined vector, indicating how many bootstrap repetitions should be used for the inner (first entry) and outer loop (second entry) to correct the bias using the double bootstrap.


a logical variable, indicating if the jackknife corrected estimate should be calculated.


a number, indicating the desired confidence interval if the single bootstrap correction is specified. By default, the interval is calculated as the bootstrap corrected and accelerated confidence interval (Bca).


a string, indicating the desired method to calculate the bootstrap confidence intervals. For more details see Studendized confidence intervals are not supported.


a logical variable, indicating whether a plot of the single bootstrap density should be generated.


As stated by Botha and van Vuuren (2010) one can estimate the intra correlation by matching VaR of a parametrized beta distribution onto the VaR of the Vasicek distribution. To do so, the shape parameters (alpha and beta) of the beta distribution are estimated according to Botha and van Vuuren (2010). Afterwards, the VaR_Beta at the confidence level of Quantile will be estimated. In a third step, this VaR_Beta is matched with the theoretical VaR of the Vasicek distribution, given by Vasicek (1991):

VaR_Vasicek= Phi((Phi^-1(PD)+sqrt(rho)*Phi^-1(Quantile))/sqrt(1-rho))

Since Quantile and the corresponding VaR_Beta is known, the intra correlation parameter can be backed out numerically. This estimator is sensitive to the chosen Quantile. Botha and van Vuuren (2010) suggested to use Quantile=0.999, but for validation purposes one may choose different values of Quantile to infer information about the robustness of the correlation estimate.

If DB is specified, the single bootstrap corrected estimate will be calculated by using the bootstrap values of the outer loop (oValues).


The returned value is a list, containing the following components (depending on the selected arguments):


Estimate of the original method


Bootstrap corrected estimate


Double bootstrap corrected estimate


Jackknife corrected estimate


Selected two-sided bootstrap confidence interval


Estimates from the bootstrap resampling


Estimates from the double bootstrap resampling- inner loop


Estimates from the double bootstrap resampling- outer loop


Botha M, van Vuuren G (2010). “Implied asset correlation in retail loan portfolios.” Journal of Risk Management in Financial Institutions, 3(2), 156–173.

Chang J, Hall P (2015). “Double-bootstrap methods that use a single double-bootstrap simulation.” Biometrika, 102(1), 203–214.

Efron B, Tibshirani RJ (1994). An introduction to the bootstrap. CRC press.

Gordy MB (2000). “A comparative anatomy of credit risk models.” Journal of Banking & Finance, 24(1), 119–149.

Vasicek O (1991). “Limiting loan loss probability distribution.” Finance, Economics and Mathematics, 147–148.

See Also

intraFMM, intraJDP1, intraJDP2 intraCMM, intraMLE, intraAMLE intraMode




#Jackknife correction
IntraCorr=intraBeta(d,n, JC=TRUE)

#Bootstrap correction with confidence intervals
IntraCorr=intraBeta(d,n, B=1000, CI_Boot=0.95 )

#Bootstrap correction with confidence intervals and plot
IntraCorr=intraBeta(d,n, B=1000, CI_Boot=0.95, plot=TRUE )

#Double Bootstrap correction with 10 repetitions in the inner loop and 50 in the outer loop
IntraCorr=intraBeta(D1,N1, DB=c(10,50))

[Package AssetCorr version 1.0.4 Index]