Rコード　データ可視化6

今回はmdsrパッケージのSATを使って、散布図をrで書いてみます。

#パッケージのインストール

library(tidyverse)
library(ggplot2)

単純な散布図を書く

g<- ggplot(data = SAT_2010,
      aes(x=expenditure,y=math))+geom_point()
plot(g)

こんなのが書ける

次は、散布図の階層わけしたものを使う

SAT_2010 <- SAT_2010 %>% mutate (SAT_rate=cut(sat_pct,breaks=c(0,30,60,100),labels=c("low","medium","high")))

g<-g %+% SAT_2010

g+aes(color=SAT_rate)

#主効果モデルの重回帰分析を行う

lm_sat<- lm(data = SAT_2010, math~expenditure+ SAT_rate)
summary(lm_sat)

anova(lm_sat)

こんな結果になる

Excelではこんな感じ

#交互作用モデルの重回帰分析

lm_sat<- lm(data= SAT_2010,math~expenditure+ SAT_rate + expenditure*SAT_rate)
summary(lm_sat)
anova(lm_sat)

結果はこんな感じ

Excelではこんな感じ