==============================================================================
POWER SIMULATION -- boxa_b2_d0_5   [metric: # successes (count)]
==============================================================================
model: behavior ~ stim + day + stim:day + (1|rat)   (behavior = # successes (count); day within-window)
observed groups: stim n=4, control n=3   nrep=120, alpha=0.05
ground truth: stim:day=+4.707/day, rat SD=5.872, residual SD=11.3, days=6

 true stim:day interaction = +4.707 (100% of observed)
  N/group  | per-animal power | LME power
  ------------------------------------------
  3        |  0.27            |  0.59  <- observed
  5        |  0.63            |  0.72
  8        |  0.91            |  0.93
  12       |  1.00            |  1.00
  16       |  0.99            |  0.99
  24       |  1.00            |  1.00

 true stim:day interaction = +2.354 (50% of observed)
  N/group  | per-animal power | LME power
  ------------------------------------------
  3        |  0.10            |  0.20  <- observed
  5        |  0.20            |  0.24
  8        |  0.31            |  0.33
  12       |  0.52            |  0.54
  16       |  0.62            |  0.62
  24       |  0.85            |  0.87

Read the per-animal column as the honest power; LME is optimistic (obs-level DF).
