The smallest detectable change

Inverts cpt_power(): searches for the change size at which the detector reaches a target power. The study-design counterpart of a power curve, and the number that belongs in a pre-registration.

Usage

cpt_min_detectable(
  n,
  sigma = 1,
  method = "pelt",
  power = 0.8,
  range = c(0.1, 5),
  n_sim = 100,
  tolerance = 5,
  change_in = "mean",
  location = 0.5,
  noise = "gauss",
  rho = 0,
  df = 3,
  tol = 0.05,
  max_iter = 12,
  seed = NULL,
  ...
)

# S3 method for class 'ggcpt_min_detectable'
print(x, ...)

Arguments

n

Series length.

sigma

Noise standard deviation. Defaults to 1.

method

Detection method. Defaults to "pelt".

power

Target detection probability. Defaults to 0.8.

range

Search range for the change size, in standard deviations. Defaults to c(0.1, 5).

n_sim

Replicates per evaluation. Defaults to 100; the answer is only as precise as this makes it, and the returned object records the Monte Carlo interval at the solution.

tolerance

A detection counts as finding the change when it falls within this many positions of it. Defaults to 5.

change_in

What changes. Defaults to "mean".

location

Changepoint position, as a fraction of n in \((0, 1)\) or an integer position. Defaults to 0.5. A vector runs one scenario per value.

noise

Noise model, passed to cpt_simulate().

rho

AR(1) parameter when noise = "ar1".

df

Degrees of freedom when noise = "t".

tol

Bisection tolerance on the change size. Defaults to 0.05.

max_iter

Maximum bisection steps. Defaults to 12.

seed

Optional seed.

...

Ignored.

x

A ggcpt_min_detectable object.

Value

A list with jump (the smallest change reaching power), achieved_power, mc_se, and the trace of evaluations, with a print() method.

See also

cpt_power().

Examples

# \donttest{
cpt_min_detectable(n = 200, n_sim = 30, max_iter = 4, seed = 1)
#> Smallest detectable change
#>   Target power:      0.8
#>   Change size:       1.02 standard deviations
#>   Achieved power:    0.9 (Monte Carlo SE 0.055)
#> 
#> # A tibble: 6 × 3
#>    jump power  mc_se
#>   <dbl> <dbl>  <dbl>
#> 1 0.1   0     0     
#> 2 0.712 0.5   0.0913
#> 3 1.02  0.9   0.0548
#> 4 1.32  0.967 0.0328
#> 5 2.55  1     0     
#> 6 5     1     0     
# }