==============================================================================
POWER SIMULATION -- naive_boxa_d6_10   [metric: # successes (count)]
==============================================================================
model: behavior ~ stim + day + stim:day + (1|rat)   (behavior = # successes (count); day within-window)
observed groups: stim n=3, control n=7   nrep=120, alpha=0.05
ground truth: stim:day=-1.619/day, rat SD=11.24, residual SD=10.14, days=5

 true stim:day interaction = -1.619 (100% of observed)
  N/group  | per-animal power | LME power
  ------------------------------------------
  3        |  0.05            |  0.08  <- observed
  5        |  0.06            |  0.07
  8        |  0.17            |  0.21
  12       |  0.26            |  0.27
  16       |  0.28            |  0.30
  24       |  0.41            |  0.45

 true stim:day interaction = -0.8095 (50% of observed)
  N/group  | per-animal power | LME power
  ------------------------------------------
  3        |  0.03            |  0.08  <- observed
  5        |  0.07            |  0.07
  8        |  0.05            |  0.10
  12       |  0.07            |  0.07
  16       |  0.11            |  0.09
  24       |  0.12            |  0.13

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