Rで何かをしたり、読書をするブログ

政府統計の総合窓口のデータや、OECDやUCIやのデータを使って、Rの練習をしています。ときどき、読書記録も載せています。

GDPギャップのデータ分析6 - GDPギャップを潜在成長率と資本投入と労働投入で回帰分析をする。」

 

www.crosshyou.info

 の続きです。

GDPギャップを回帰分析してみましょう。

R言語のlm関数を使います

f:id:cross_hyou:20200530141527j:plain

p-valueは3.126e-07と0.05よりも小さいので有意なモデルです。

InterceptはP値が0.05よりも大きいですね。削除してみます。

f:id:cross_hyou:20200530141911j:plain

update関数で、-1とすれば切片を削除できます。anova関数でmodel1とmodel2を比較しましたが、Pr(>F)が0.1361なので、二つのモデルに有意な違いはありません。

model2をsummary関数で確認してみます。

f:id:cross_hyou:20200530142057j:plain

Multiple R-Squaredが0.3229とmodel1よりもよくなりました。

実際のGDPギャップとモデルの予測値の散布図を見てみましょう。

f:id:cross_hyou:20200530150833j:plain

f:id:cross_hyou:20200530150848j:plain

model2のほうが赤い直線の周りに集まっているように見えます。

model2の残差プロットなどをみてみます。

f:id:cross_hyou:20200530151117j:plain

f:id:cross_hyou:20200530151131j:plain

左上の残差プロットは平均した散らばり具合なのでいい感じだと思います。

GDPギャップは潜在成長率が小さいほど、資本投入と労働投入が大きいほど大きいといえます。

今回は以上です。