Nothing
# test_that(
# desc = "linreg_fit with weights approximately equal to glm()",
# code = {
#
# nrows <- 1000
# ncols <- 20
#
# X <- matrix(data = rnorm(nrows*ncols), nrow = nrows, ncol = ncols)
#
# # X <- cbind(1, X)
#
# colnames(X) <- c(
# # "intercept",
# paste0("x", seq(ncols))
# )
#
# Y <- matrix(rnorm(nrows), ncol = 1)
#
# glm_data <- as.data.frame(cbind(y=as.numeric(Y), X))
#
# # Fit logistic regression using the custom function
#
# W <- sample(1:3, nrow(X), replace=TRUE)
#
# cpp = linreg_fit_exported(X, Y, W, do_scale = TRUE,
# epsilon = 1e-9, iter_max = 20)
#
# R = lm(y ~ ., weights = as.integer(W), data = glm_data)
#
# R_summary <- summary(R)
#
# R_beta_est <- as.numeric(R_summary$coefficients[-1, 'Estimate'])
# R_beta_pvalues <- as.numeric(R_summary$coefficients[-1, 'Pr(>|t|)'])
#
# expect_equal(cpp[,1], R_beta_est, tolerance = 1e-9)
# expect_equal(cpp[,2], R_beta_pvalues, tolerance = 1e-9)
#
#
# }
# )
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.