This is an internal helper function to compute the working correlation matrix based on standardized residuals.
.update_correlation(residuals, nobs, corstr, tau, cluster_indices)
The updated working correlation matrix `R`.
A numeric vector of residuals (y - mu).
A numeric vector of observations per subject.
A string specifying the correlation structure.
The quantile being estimated.
A numeric vector of cumulative subject observation counts.