This function computes Kullback-Leibler divergence (KLD) of quantile spectral estimate.
qkl.divergence(y.qper, qspec, sel.f = NULL, sel.tau = NULL)real number of Kullback-Leibler divergence
matrix or array of quantile spectral estimate from, e.g., qspec.lw()
matrix of array of true quantile spectrum (same dimension as y.qper)
index of selected frequencies for computation (default = NULL: all frequencies)
index of selected quantile levels for computation (default = NULL: all quantile levels)