============================================================================== POWER SIMULATION -- boxa_b2_d0_10 [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=+1.604/day, rat SD=6.047, residual SD=12.42, days=11 true stim:day interaction = +1.604 (100% of observed) N/group | per-animal power | LME power ------------------------------------------ 3 | 0.17 | 0.33 <- observed 5 | 0.46 | 0.58 8 | 0.67 | 0.70 12 | 0.90 | 0.93 16 | 0.95 | 0.97 24 | 1.00 | 1.00 true stim:day interaction = +0.802 (50% of observed) N/group | per-animal power | LME power ------------------------------------------ 3 | 0.11 | 0.11 <- observed 5 | 0.09 | 0.14 8 | 0.23 | 0.25 12 | 0.33 | 0.35 16 | 0.51 | 0.53 24 | 0.67 | 0.66 Read the per-animal column as the honest power; LME is optimistic (obs-level DF).