tests/testthat/profile_chisq.R

In miheerdew/bmdupdate: This package implements the pvalue calculation for the ChiSq test

source("helper.R")
library(bmdupdate)
library(rbenchmark)

sim <- generate_data(1000, 1000, rho, 500, beta, s2)
b <- new(BackendChisq, sim$X, sim$Y)

A <- 1:400
testX <- TRUE

show(benchmark(
          Cpar = b$trace_multi(A, testX, TRUE),
          Cseq = b$trace_multi(A, testX, FALSE),
          RseqX = {
            c <- precomputations(sim$X[,A], sim$Y)
            t(sapply(1:sim$dy, function(i) traceXR(i,c)))
        },
        RseqM = {
          c <- precomputations(sim$X[,A], sim$Y)
          t(sapply(1:sim$dy, function(i) traceMR(i,c)))
        },
        replications=1))