Nothing
library(fdasrvf)
testthat::test_that("Verify pair_align_functions_bayes() works as intended", {
set.seed(1)
data('simu_data')
out <- pair_align_functions_bayes(f1 = simu_data$f[1:100,1], f2 = simu_data$f[1:100,2],
timet = simu_data$time[1:100], iter = 5000, times = 5, tau = 2, powera = 1,
showplot = FALSE, extrainfo = TRUE)
# verify functions match approximately
# testthat::expect_equal(sum(simu_data$f[1:100,1] - out$f2_a), -10.6867105)
testthat::expect_equal(sum(simu_data$f[1:100,1] - out$f2_q), -10.58674295)
# verify warning issued for times = 2
testthat::expect_warning(pair_align_functions_bayes(f1 = simu_data$f[,1],
f2 = simu_data$f[,2], timet = simu_data$time, iter = 1000, times = 2,
tau = 2, powera = 1, showplot = FALSE, extrainfo = FALSE))
# verify acceptance rate
# testthat::expect_equal(mean(out$pct_accept), 0.0168)
})
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.