www.crosshyou.infowww.crosshyou.info政府統計の総合窓口のデータや、OECDのデータを使って、Rの練習をしています。ときどき、読書記録も載せています。R言語のsample関数とfor関数でモンテカルロシミュレーションデータ分析www.crosshyou.infoの続きです。今回はR言語のlm関数で回帰分析をしてみます。variableはtravel:旅行・行楽の年間行動者率にしてみます。Y軸をforeignにして散布図を描きます。右肩上がりの散布図ですね。cor.test関数で相関があるかどうかの検定をしてみます。p値は1.047e-10と0.05よりも小さいです。正の相関があると言えます。この相関係数の信頼区間をモンテカルロで求めてみましょう。横道にそれますがやってみます。180個あるdfのデータの中からランダムにn個(今回は100個)の行をとって、そのn個の行でtravelとforeignの相関係数を計算します。for関数でN回(今回は10万回)繰り返して、soukanという箱に入れました。hist関数でヒストグラムを描いてみましょう。正規分布っぽいヒストグラムですね。quantile関数で下位0.025と上位0.975の値を求めます。0.5364974の間にある、ということです。0.5562115とよく似た値ですね。回帰分析するつもりがモンテカルロになってしましました。次回はwww.crosshyou.infoです。はじめから読むには、www.crosshyou.infoです。