`dss_llogis()` uses the wrong variance. In `R/scores_llogis.R`: `v <- ell^2 * 2 * b / sin(2 * b) - b^2 / sb^2` going by [Wikpedia](https://en.wikipedia.org/wiki/Log-logistic_distribution#Moments) the second part should also be multiplied by ell^2: `v <- ell^2 * (2 * b / sin(2 * b) - b^2 / sb^2)`
dss_llogis()uses the wrong variance. InR/scores_llogis.R:v <- ell^2 * 2 * b / sin(2 * b) - b^2 / sb^2going by Wikpedia the second part should also be multiplied by ell^2:
v <- ell^2 * (2 * b / sin(2 * b) - b^2 / sb^2)