tests/testthat/test-d.spls.norm.R

In dual.spls: Dual Sparse Partial Least Squares Regression

test_that("d.spls.norm works", {

#Data1
vec=1:100

#Positive value
expect_gt(d.spls.norm1(vec),0)
expect_gt(d.spls.norm2(vec),0)

#equal norm1 and norm2
expect_equal(d.spls.norm1(vec),sum(abs(vec)))
expect_equal(d.spls.norm2(vec),sqrt(sum(vec^2)))

#Data2
n <- 100
p <- 50
nondes <- 20
sigmaondes <- 0.5
data=d.spls.simulate(n=n,p=p,nondes=nondes,sigmaondes=sigmaondes)
X <- data$X
y <- data$y

#Positive value
expect_gt(d.spls.norm1(y),0)
expect_gt(d.spls.norm2(y),0)

#equal norm1 and norm2
expect_equal(d.spls.norm1(y),sum(abs(y)))
expect_equal(d.spls.norm2(y),sqrt(sum(y^2)))})