tests/testthat/test_stdmod_selected_print_summary.R

skip("Test in interactive sessions")
library(testthat)
library(stdmod)

set.seed(64252)
n <- 1000
x1 <- rnorm(n)
x2 <- rnorm(n) / 1000
x3 <- rnorm(n)
y <- .2 * x1 + .3 * x2 + 0 * x3 + rnorm(n, 0, 1)

dat <- data.frame(x1, x2, x3 = x3 * 100, y)

lm_out <- lm(y ~ x1 + x2 + x3, dat)
lm_std <- std_selected(lm_out, to_standardize = ~ y + x1 + x2)
lm_std_boot <- std_selected_boot(lm_out, to_standardize = ~ y + x1 + x2)
lm_out2 <- std_selected(lm_out, to_standardize = ~ y + x1)
lm_std2 <- std_selected(lm_out, to_standardize = ~ y + x1)
lm_std_boot2 <- std_selected_boot(lm_out2, to_standardize = ~ y + x1)

summary(lm_std)
print(summary(lm_std), est_digits = 4, pvalue_less_than = 1e-3)
print(summary(lm_std), est_digits = 2, pvalue_less_than = 1e-3)
print(summary(lm_std), est_digits = 2, pvalue_less_than = 1e-4, t_digits = 2)

summary(lm_std_boot)
print(summary(lm_std_boot), est_digits = 4, pvalue_less_than = 1e-3)
print(summary(lm_std_boot), est_digits = 2, pvalue_less_than = 1e-2, t_digits = 1)
print(summary(lm_std_boot), est_digits = 2, pvalue_less_than = 1e-2, t_digits = 5)

print(summary(lm_std_boot), default_style = TRUE)

summary(lm_std2)
print(summary(lm_std2), est_digits = 4, pvalue_less_than = 1e-3)
print(summary(lm_std2), est_digits = 2, pvalue_less_than = 1e-3)
print(summary(lm_std2), est_digits = 2, pvalue_less_than = 1e-4, t_digits = 2)

summary(lm_std_boot2)
print(summary(lm_std_boot2), est_digits = 4, pvalue_less_than = 1e-3)
print(summary(lm_std_boot2), est_digits = 2, pvalue_less_than = 1e-2, t_digits = 1)
print(summary(lm_std_boot2), est_digits = 2, pvalue_less_than = 1e-2, t_digits = 5)

print(summary(lm_std_boot2), default_style = TRUE)

Try the stdmod package in your browser

Any scripts or data that you put into this service are public.

stdmod documentation built on Sept. 30, 2024, 9:42 a.m.