の続きです。
前回は埼玉県、千葉県、神奈川県は一人当り医療費が少なく、高知県、長崎県、鹿児島県は一人当り医療費が多いことがわかりました。
今回は、一人当り医療費と一人当り県内総生産との関係を調べてみます。
まず、一人当り県内総生産を計算します。
summary関数で平均値を調べると、372万6251円です。
次に、複数年度ある一人当り医療費、一人当り県内総生産の都道府県ごとに平均値を計算します。tapply関数を使います。
plot関数で散布図を見てみます。
右端にポツンと一つはぐれプロットがありますね。一人当り県内総生産が突出して高い都道府県があるのですね。sort関数で並び替えて確認します。
東京都が突出して一人当り県内総生産が大きいですね。
東京都を削除して散布図を描いてみます。
散布図はまんべんなく散らばっている感じですね。なんとなく、一人当り県内総生産が大きいほど、一人当り医療費は小さい感じですね。回帰分析をしてみます。
まずは、東京都を除外した変数を用意します。
length関数でデータの観測数を確認しています。46なので一つ、東京都が除外されていることがわかります。
lm関数で回帰分析をします。
p-valueは0.02855なので0.05よりも小さいですから有意なモデルです。
といってもR-squaredが0.1043ですから、説明力は10%ぐらいしかないです。
avgHGDP2の係数の符号はマイナスですから、やはり一人当り県内総生産が大きいほうが一人当り医療費は少ない、という関係です。
一人当り県内総生産が大きいということは、一人ひとりの生産額が大きい、つまり病気でない、つまり医療費は少ない、ということでしょうかね。
残差のグラフを描いてみます。
最後に散布図を回帰直線を重ねて描いてみます。
今回は以上です。
