Fitting Linear Quantile Regression Mixed Models with
Relationship Matrix
Fit a quantile regression mixed model involved Relationship Matrix using a sparse implementation of the
Frisch-Newton interior-point algorithm as described in
Portnoy and Koenker (1977, Statistical Science) < https://www.jstor.org/stable/2246216>.