test_that("omega inv sqrt monovar", {
expect_equal(omega_inv_sqrt(c(1,5,6,7),0.5) %>% unname,
c(sqrt(0.75),4.5,3.5,4))
})
test_that("error praislm", {
expect_error(praislm(ts(matrix(1:10,ncol=1),start=2011),ts(matrix(1:10,ncol=1),start=2010),FALSE,FALSE,numeric(),NULL),
"same windows and frequencies")
expect_error(praislm(1:120,ts(1:10,frequency=1,start=2010),
TRUE,numeric(),NULL),"ts object")
expect_error(praislm(ts(1:120,frequency=12,start=2010),1:10,
TRUE,numeric(),NULL),"ts object")
expect_error(praislm(ts(matrix(1,1,1),frequency=1,start=2010),
ts(1,frequency=1,start=2010),
FALSE, TRUE, numeric(), NULL),"less than two observations")
expect_error(praislm(ts(structure(c(1, 1, pi, 17 * pi),
dim = c(2L, 2L),
dimnames = list(NULL, c("constant", "hfserie"))),
frequency=1,start=2010),
ts(c(99,-4),frequency=1,start=2010),
TRUE, FALSE, numeric(), NULL),"rank becomes imperfect")
})
test_that("max iter",{
expect_warning(praislm_impl(matrix(c(1.324398899494450,-0.953033704627492,0.153551486265972),ncol=1),
c(26.080220590193644,-20.148635100560220,2.606520694629608),TRUE),
"Maximum iterations")
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.