Nothing
test_that("glmnet and SLOPE return same unpenalized model", {
set.seed(129)
n <- 100
x1 <- rnorm(n)
x2 <- rnorm(n)
prob <- matrix(c(
rep(1, n),
exp(3 + 2 * x1 + x2),
exp(-1 + x1 - 3 * x2)
),
ncol = 3
)
prob <- sweep(prob, 1, apply(prob, 1, sum), "/")
y <- double(n)
for (i in 1:n) {
y[i] <- sample(3, 1, replace = TRUE, prob = prob[i, ])
}
x <- scale(cbind(x1, x2))
y <- factor(y)
library(glmnet)
fit <- glmnet(
x,
y,
family = "multinomial",
lambda = 0,
thresh = 1e-10,
standardize = FALSE
)
g_coef <- as.matrix(do.call(cbind, coef(fit)))
g_coef[, ] <- g_coef[, ] - g_coef[, 3]
g_coef <- g_coef[, 1:2]
ofit <- SLOPE(x, y, family = "multinomial", alpha = 1e-9)
expect_equivalent(g_coef, coef(ofit), tol = 1e-4)
})
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.