“…According to Sklar's theorem where F represents the joint cdf of ( x l,t , x j,t ), and C the copula function (Zimmer ). For each state‐pair combination, we proceed to implementing Racine's () version of Fermanian and Scaillet's () nonparametric estimator of F , that is, where , and the integral corresponds to a standard kernel density estimator that utilizes an adaptive nearest‐neighbor bandwidth with a second‐order Gaussian kernel. In terms of empirical implementation, the interesting feature in this type of modeling is that the bandwidth changes with each sample realization in the dataset when estimating the cdfs at any point of the multivariate (copula) distribution and it demonstrates superior numerical performance when estimating extreme tails (Li and Racine ).…”