wls_qr <- function(X, y, weights) { Phi <- diag(weights) beta <- qr.solve(t(X) %*% Phi %*% X, t(X) %*% Phi %*% y) return(beta) }
Add the following code to your website.
REMOVE THIS Copy to clipboard
For more information on customizing the embed code, read Embedding Snippets.