df <- read.csv("data09.csv", header=T, fileEncoding="utf-8")

res.1 <- lm(T2_SDS ~ scale(T1_SDS), data=df)

res.2 <- lm(T2_SDS ~ scale(T1_SDS) + scale(Rumination) + scale(Stress), data=df)

res.3 <- lm(T2_SDS ~ scale(T1_SDS) + scale(Rumination) * scale(Stress), data=df)

anova(res.1, res.2, res.3)