理学療法士がまとめたノート

統計学備忘録(R言語のメモ)

since2016 ときどきTEXのメモ

回帰直線と交互作用

f:id:yoshida931:20210611162531p:plain
パッケージ無しで描く方法

dat <- data.frame(y=c(1,2,3,4), x=c(4,5,8,12), a=c(1,1,0,0))
dat$a <- as.factor(dat$a)

plot(dat$x, dat$y)
fit <- summary(lm(y~x, data=dat))
lines(range(dat$x), fit$coef[1]+ fit$coef[2]*range(dat$x))

pchAB <- ifelse(dat$a == "1", 19, 21)
plot(dat$x, dat$y,  pch=pchAB, cex=1.5 )
legend("bottomright",
    legend=c("a=1", "a=0"),
    pch=c(19, 21),
    lty = c(1, 3)
    )

fit1 <- summary(lm(y ~ x, data=dat[dat$a==1,]))
lines(range(dat$x), fit1$coef[1]+ fit1$coef[2]*range(dat$x))
fit0 <- summary(lm(y ~ x, data=dat[dat$a==0,]))
lines(range(dat$x), fit0$coef[1]+ fit0$coef[2]*range(dat$x), lty=3)