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

 true stim:day interaction = +10.2 (100% of observed)
  N/group  | per-animal power | LME power
  ------------------------------------------
  3        |  0.55            |  0.91  <- observed
  5        |  0.93            |  0.98
  8        |  1.00            |  1.00
  12       |  1.00            |  1.00
  16       |  1.00            |  1.00
  24       |  1.00            |  1.00

 true stim:day interaction = +5.1 (50% of observed)
  N/group  | per-animal power | LME power
  ------------------------------------------
  3        |  0.11            |  0.33  <- observed
  5        |  0.42            |  0.57
  8        |  0.67            |  0.76
  12       |  0.89            |  0.90
  16       |  0.97            |  0.97
  24       |  1.00            |  1.00

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