==============================================================================
POWER SIMULATION -- naive_boxa_d0_13   [metric: # successes (count)]
==============================================================================
model: behavior ~ stim + day + stim:day + (1|rat)   (behavior = # successes (count); day within-window)
observed groups: stim n=3, control n=8   nrep=120, alpha=0.05
ground truth: stim:day=+1.306/day, rat SD=8.76, residual SD=15.18, days=14

 true stim:day interaction = +1.306 (100% of observed)
  N/group  | per-animal power | LME power
  ------------------------------------------
  3        |  0.12            |  0.32  <- observed
  5        |  0.36            |  0.47
  8        |  0.63            |  0.68  <- observed
  12       |  0.94            |  0.94
  16       |  0.95            |  0.96
  24       |  0.97            |  0.98

 true stim:day interaction = +0.6532 (50% of observed)
  N/group  | per-animal power | LME power
  ------------------------------------------
  3        |  0.07            |  0.17  <- observed
  5        |  0.17            |  0.18
  8        |  0.33            |  0.31  <- observed
  12       |  0.37            |  0.38
  16       |  0.42            |  0.51
  24       |  0.62            |  0.60

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