============================================================================== 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).