Skip to contents

Power calculation for tumor growth

Source: R/tumr_power.R
tumr_power.Rd

For speed, uses a response feature analysis rfeat(); this will likely underestimate power somewhat.

Usage

tumr_power(n, effect_size, N = 1000, ...)

Arguments

n

Sample size per group (integer or vector)

effect_size

As in gendat() (numeric or vector)

N

Number of simulations (default: 1000)

...

Other arguments to gendat()

Value

An array of p-values

Examples

res <- tumr_power(8, 2, 10)
#> 
  |                                                                            
  |                                                                      |   0%
  |                                                                            
  |=======                                                               |  10%
  |                                                                            
  |==============                                                        |  20%
  |                                                                            
  |=====================                                                 |  30%
  |                                                                            
  |============================                                          |  40%
  |                                                                            
  |===================================                                   |  50%
  |                                                                            
  |==========================================                            |  60%
  |                                                                            
  |=================================================                     |  70%
  |                                                                            
  |========================================================              |  80%
  |                                                                            
  |===============================================================       |  90%
  |                                                                            
  |======================================================================| 100%
mean(res$p < 0.05)
#> [1] 1