Calculo de muestra

# Supuestos
alpha <- 0.05
power <- 0.80
HR <- 0.60
p_exp <- 0.30        # prevalencia de isquemia

# Z scores
Z_alpha <- qnorm(1 - alpha/2)
Z_beta  <- qnorm(power)

# Eventos requeridos
E <- (Z_alpha + Z_beta)^2 /
  ((log(HR))^2 * p_exp * (1 - p_exp))
E

## [1] 143.233

# Pasar a tamaño muestral
inc_event <- 0.55
inc_comp  <- 0.03

N <- E / inc_event
N_adj <- N / (1 - inc_comp)

# Ajuste conservador por covariables
N_final <- N_adj * 1.20
N_final

## [1] 322.1735