myscratchpad

dat <- read.table("https://blog.cnobi.jp/v1/blog/user/89d80905c7038b4121822249e9062fba/1256404631",header=T)

dat1 <- dat

mod1sd <- sd(dat1$mod1)
mod1ab <- dat1$mod1 - mod1sd
mod1bl <- dat1$mod1 - (-1*mod1sd)
mod2sd <- sd(dat1$mod2)
mod2ab <- dat1$mod2 - mod2sd
mod2bl <- dat1$mod2 - (-1*mod2sd)
dat2 <- cbind(dat1, mod1ab, mod1bl, mod2ab, mod2bl)

lmres.whl <- lm(dv ~ idv*mod1*mod2, data = dat2)
summary(lmres.whl)

lmres.1bl.2bl <- lm(dv ~ idv*mod1bl*mod2bl, data = dat2)
lmres.1bl.2ab <- lm(dv ~ idv*mod1bl*mod2ab, data = dat2)
lmres.1ab.2bl <- lm(dv ~ idv*mod1ab*mod2bl, data = dat2)
lmres.1ab.2ab <- lm(dv ~ idv*mod1ab*mod2ab, data = dat2)
summary(lmres.1bl.2bl) #LL
summary(lmres.1bl.2ab) #LH
summary(lmres.1ab.2bl) #HL
summary(lmres.1ab.2ab) #HH

# 標準化値を調べる
dat3 <- dat1
dat3$idv <- scale(dat1$idv)
dat3$mod1 <- scale(dat1$mod1)
dat3$mod2 <- scale(dat1$mod2)
dat3$dv <- scale(dat1$dv)
smod1sd <- sd(dat3$mod1)
smod1ab <- dat3$mod1 - smod1sd
smod1bl <- dat3$mod1 - (-1*smod1sd)
smod2sd <- sd(dat3$mod2)
smod2ab <- dat3$mod2 - smod2sd
smod2bl <- dat3$mod2 - (-1*smod2sd)
dat4 <- cbind(dat3, smod1ab, smod1bl, smod2ab, smod2bl)

slmres.whl <- lm(dv ~ idv*mod1*mod2, data = dat4)
summary(slmres.whl)
slmres.1bl.2bl <- lm(dv ~ idv*smod1bl*smod2bl, data = dat4)
slmres.1bl.2ab <- lm(dv ~ idv*smod1bl*smod2ab, data = dat4)
slmres.1ab.2bl <- lm(dv ~ idv*smod1ab*smod2bl, data = dat4)
slmres.1ab.2ab <- lm(dv ~ idv*smod1ab*smod2ab, data = dat4)
summary(slmres.1bl.2bl) # LL
summary(slmres.1bl.2ab) # LH
summary(slmres.1ab.2bl) # HL
summary(slmres.1ab.2ab) # HH