test <- lmer(ln_auc ~ zNSE_mean + waketimeday1 + zsoc_negaff + zp_soc_negaff + black + female + zincome + zpared + zstress_sum
            + (1 | id), data = na.omit(MADIC_hlm))

fixef(test)["zNSE_mean"] <- .2

powerSim(test, nsim = 1000)