# lm()を使って回帰分析
fit_iris <- lm(Petal.Length ~ Petal.Width_centered * Species, data = mydf)
summary(fit_iris)
# plot()を使って散布図
plot(mydf[,c("Petal.Width_centered","Petal.Length")])