============================================================================== POWER SIMULATION -- naive_a2_d0_5 [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=+6.962/day, rat SD=9.643, residual SD=12.11, days=6 true stim:day interaction = +6.962 (100% of observed) N/group | per-animal power | LME power ------------------------------------------ 3 | 0.43 | 0.85 <- observed 5 | 0.87 | 0.96 8 | 1.00 | 1.00 12 | 1.00 | 1.00 16 | 1.00 | 1.00 24 | 1.00 | 1.00 true stim:day interaction = +3.481 (50% of observed) N/group | per-animal power | LME power ------------------------------------------ 3 | 0.12 | 0.31 <- observed 5 | 0.33 | 0.42 8 | 0.57 | 0.67 12 | 0.74 | 0.81 16 | 0.90 | 0.89 24 | 0.97 | 0.99 Read the per-animal column as the honest power; LME is optimistic (obs-level DF).