ccov.np(data, score, bwscale = 1.1, thetagrid = seq(-3, 3, len = 200), progress = TRUE, scale_score = TRUE)
bwscale
times $N^{-1/5}$.
score
should be z standardized in advance of the calculation of
conditional covariancesZhang, J., & Stout, W. (1999). Conditional covariance structure of generalized compensatory multidimensional items, Psychometrika, 64, 129-152.