per_change_delta_se <- function()

in utils.R [50:68]

• 15 lines of code
• 1 McCabe index (conditional complexity)

per_change_delta_se <- function(
  m.test, m.control,
  se.test, se.control
) {
  
  mean <- (m.test - m.control) / abs(m.control) -
    se.control ^ 2 * m.test / abs(m.control) ^ 3
  
  var <- (se.test / abs(m.control)) ^ 2 +
    (se.control * abs(m.test)) ^ 2 / abs(m.control) ^ 4
  
  z95 <- stats::qnorm(1 - (1 - 0.95) / 2, 0, 1)
  
  return(data.frame(
    mean     = mean,
    lower95  = mean - z95 * sqrt(var),
    upper95  = mean + z95 * sqrt(var)
  ))
}