tmp-tests/faster-IRLS.R
In privefl/bigstatsr: Statistical Tools for Filebacked Big Matrices
N <- 1000
K <- 10
x <- matrix(rnorm(N * K), N)
y <- sample(0:1, N, TRUE)
summary(mod1 <- glm(y ~ x))
x2 <- cbind(1, x); x2[, -2] <- svd(x2[, -2])$u
summary(mod2 <- glm(y ~ x2 + 0))
rbind(mod1$coefficients, mod2$coefficients)
rbind(summary(mod1)$coeff[, 4], summary(mod2)$coeff[, 4])
privefl/bigstatsr documentation built on March 29, 2024, 3:31 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
N <- 1000
K <- 10
x <- matrix(rnorm(N * K), N)
y <- sample(0:1, N, TRUE)
summary(mod1 <- glm(y ~ x))
x2 <- cbind(1, x); x2[, -2] <- svd(x2[, -2])$u
summary(mod2 <- glm(y ~ x2 + 0))
rbind(mod1$coefficients, mod2$coefficients)
rbind(summary(mod1)$coeff[, 4], summary(mod2)$coeff[, 4])
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.