# 第1章 モデル # 第2節 単回帰分析 # 福元健太郎 # データの読み込み # データは森健太郎氏作成(2010年度受講生、本人了承済み) データ <- read.csv("abortion.csv") # オブジェクトの内容を表示 データ # データの登録 attach(データ) # 散布図(plot)の作成 plot(大統領D,中絶率,main="リベラルな州ほど中絶の割合が高い",xlab="民主党得票率(%)",ylab="千人当たり中絶件数") # 千人当たり中絶件数を民主党得票率で説明する単回帰分析 # lmはLinear Model(線形モデル)の略 回帰分析結果 <- lm(中絶率 ~ 大統領D) # 係数や有意水準の要約(summary) summary(回帰分析結果) # 回帰直線 # y=a+bxをabline(a,b)で描くのでablineか 回帰係数 <- 回帰分析結果$coef 回帰係数 回帰係数[1] 回帰係数[2] abline(回帰係数[1],回帰係数[2]) ###################### # 平均値の比較との関係 # 平均値の比較 (共和党知事州の平均中絶率 <- mean(中絶率[知事==1])) (民主党知事州の平均中絶率 <- mean(中絶率[知事==0])) 共和党知事州の平均中絶率 - 民主党知事州の平均中絶率 # 単回帰分析 回帰分析結果2 <- lm(中絶率 ~ 知事) summary(回帰分析結果2) # 両者の違い (共和党知事州の平均中絶率 - 民主党知事州の平均中絶率) - 回帰分析結果2$coef[2]