tests/testthat/test-error_length_mu_sd.R

#context("Error lenght mu sigma")
test_that("error difference lenght of mu and sigma", {
  local_edition(3)
  mu_T <- c(37162.0, 35702.0)
  mu_R1 <- c(37705.0, 35930.0, 125.1)
  mu_R2 <- c(37702.8, 35862.4, 126.9)
  mu_list <- list(mu_T,mu_R1,mu_R2)

  # Save variance covariance matrix in a list
  sigma_T <- c(11113.62172, 9132.75342, 16.89586)
  sigma_R1 <- c(12332.41615, 9398.42182, 17.88151)
  sigma_R2 <- c(12113.72, 9098.42182, 17.1586)

  sigma_list <- list(sigma_T,sigma_R1,sigma_R2)

  # Same treatment allocation rate
  TAR = c(1,1,1) # we assume same allocation rate in both arms

  expect_error(SimTOST::sampleSize( mu_list = mu_list,
                                      sigma_list = sigma_list,
                                      varcov_list = NA,
                                      power = 0.9,
                                      dtype = "parallel",
                                      ctype = "ROM",
                                      vareq = T,
                                      lognorm = T,
                                      k=1,
                                      list_comparator =list(c("T","R1"),c("T","R2")),
                                      arm_names=c("T","R_1","R2"),
                                      ncores=1),"In each arm, 'mu', 'sigma', and 'y_name' must have the same length.")
})

Try the SimTOST package in your browser

Any scripts or data that you put into this service are public.

SimTOST documentation built on April 3, 2025, 9:05 p.m.