Nothing
library("SISIR")
context("Testing ridge regression `ridgeSIR`")
test_that("Normalization in ridgeSIR works as expected.", {
set.seed(1140)
tsteps <- seq(0, 1, length = 200)
nsim <- 100
simulate_bm <- function() return(c(0, cumsum(rnorm(length(tsteps)-1, sd=1))))
x <- t(replicate(nsim, simulate_bm()))
beta <- cbind(sin(tsteps*3*pi/2), sin(tsteps*5*pi/2))
y <- log(abs(x %*% beta[ ,1])) + sqrt(abs(x %*% beta[ ,2])) +
rnorm(nsim, sd = 0.1)
# perform ridge regression
res_ridge <- ridgeSIR(x, y, H = 10, d = 5, mu2 = 10^8)
# test if normalization is OK
norm_EDR <- res_ridge$utils$norm_EDR
expect_equal(diag(res_ridge$parameters$d),
crossprod(res_ridge$EDR, norm_EDR) %*% res_ridge$EDR,
tolerance = 10^(-10))
})
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.