データ分析
www.crosshyou.info の続きです。 今回は、老人福祉費 / 児童福祉費 という比率を計算してみます。 まず、老人福祉費と児童福祉が同じ都道府県の順番で並んでいるか確認します。 == で同じかどうかをテストしました。すべてTRUEなので同じですね。sumでTRUE…
www.crosshyou.info の続きです。 今回は一人当りの児童福祉費(ChildpM)を重回帰分析します。 説明変数は対数をとった面積(logArea)と対数をとった県内総生産(logGDP)です。 まずは、一番複雑なモデルから。 logArea:logGDPの交差項は不要のようです。削除し…
www.crosshyou.info の続きです。 今回は一人当りの老人福祉費を面積(対数をとったもの)と県内総生産(対数をとったもの)の二つの変数で回帰分析してみたいと思います。 まずは、それぞれの変数との散布図を描いてみます。 面積は関係なさそうですが、GDPは関…
www.crosshyou.info の続きです。 今回は棒グラフと信頼区間を表示してみたいと思います。 Statistics: An Introduction Using R 作者: Michael J. Crawley 出版社/メーカー: Wiley 発売日: 2019/12/09 メディア: ペーパーバック この商品を含むブログを見る…
www.crosshyou.info の続きです。 今回の説明変数はIndus, 業種です。summary関数でどういう業種があるか見てみます。 33種類の業種があります。今回は前回までとは違ったアプローチでANOVAをやってみたいと思います。 いつものように参考図書は、 Statistic…
www.crosshyou.info の続きです。 今回は、大企業、中堅企業、中小企業という企業規模の違い、Scaleをexplanatory variableにしてANOVAをしてみます。企業規模の違いで短観の数値に違いはあるでしょうか？ まずは、グラフで様子を確認しましょう。 黒い点が…
www.crosshyou.info の続きです。今回もANOVAをします。今度は、Sector(製造業、非製造業、全産業)を説明変数にして、Value、短観の値を反応変数にします。 またグラフでSectorによって短観の値に違いがあるか見てみましょう。 黒い点が製造業で黒い水平線が…
www.crosshyou.info の続きです。今回もANOVAです。explanatoly variables(説明変数)をType(最近、先行きの２値をとるカテゴリカル変数)にresponse variable(反応変数)はValue(短観の値)です。 まずは、前回と同じくグラフで見てみます。 黒い点が最近、赤い…
今回は先日発表された日銀の短観データを分析してみます。 日銀のウェブサイトからダウンロードしたファイルは、上のようなものでした。 これを こんな感じにデータを組みなおして、read.csv関数でR言語で読み込み、分析してみます。 summary関数でそれぞれ…
www.crosshyou.info の続きです。 児童福祉費と人口は関係ありましたね。児童福祉費と面積はどうでしょうか？ p-valueは0.1507なので有意ではないですね。面積は児童福祉費には関係ないようです。 では、県内総生産を見てみましょう。 p-valueが2.2e-16より…
www.crosshyou.info の続きです。 今回はR言語で回帰分析をします。老人福祉費または児童福祉を人口、面積、県内総生産の3つの変数で回帰分析します。事前の予想は、人口と県内総生産は関係あるけど、面積は関係ないんでは？と思います。 早速やってみます。…
www.crosshyou.info の続きです。 今回は、variance(分散)を調べます。上に画像のある、Michael J. CrawleyのStatistics An Introduction using Rを参考にします。 varianceはvar関数で簡単にわかります。 10年間平均の老人福祉の対数の分散は、 0.07579028…
www.crosshyou.info の続きです。 今回は、老人福祉費と児童福祉費のCentral Dendencyを調べようと思います。 いつものように Michael J. Crawleyの Statistics An introduction using R を参考にします。 Statistics: An Introduction Using R 作者: Michae…
www.crosshyou.info の続きです。前回は都道府県別の10年間の平均値をtapply関数を使って算出しました。 せっかくなので、グラフにしてどのようなデータ分布なのか見てみます。 まずは、いっぺんに三つのグラフを描く関数を作ります。 小さい順に並べたグラ…
今回は、都道府県別の老人福祉費と児童福祉費の分析をしてみようと思います。 データは、政府統計の総合窓口e-stat(www.e-stat.go.jp)から取得しました。 取得したデータは５つです。総人口、総面積、県内総生産額、老人福祉費、児童福祉費です。 ExcelのCSV…
www.crosshyou.info の続きです。 今回は、前回作成したデータフレームの各変数を大きい順、小さい順に並び替えましょう。 まずは、沖縄の平均気温の高い順に並び替えます。 Statistics: An Introduction Using R 作者: Michael J. Crawley 出版社/メーカー:…
www.crosshyou.info の続きです。前回はNAを0に変換しました。 今回は、台風の上陸回数、発生回数をベクトルにします。 まず、現状のデータフレームがどういうふうなのかhead関数で見てみます。 HASSとJOURは1行目が1951年、2行目が1952年と横向きにデータが…
今回は、気象庁のホームページ www.jma.go.jp から、台風の発生数、上陸数のデータを取得して分析しようと思います。 からデータをダウンロードしました。 台風の発生数、上陸数、那覇と東京の月別の平均気温を取得しました。 気温はこんなデータファイルで…
www.crosshyou.info の続きです。 今回は日経平均の騰落を回帰分析したいと思います。ロジスティック回帰分析です。 まずは、日経平均の騰落を表すベクトルを作成します。 まず、vNikkの先頭にNAを追加したベクトル、xを作成しました。 これで、vNikk - x で…
www.crosshyou.info の続きです。というよりは修正です。 前回は、ドル円については、日経平均とラグが無いとして分析しましたが、よく考えたら9月1日時点では、9月末のドル円の値はわかってなくて、8月末の値しかわからないのだから、1か月ずらさないといけ…
www.crosshyou.info の続きです。 最終的な目標は日経平均の値動きをドル円、消費支出、鉱工業生産で説明したい、ということですが、今回はその前段階として利用可能なデータの時点を合わせたいと思います。 どういうことかというと、9月1日の時点で、日経平…
vdata.nikkei.com 今回はこの、日経の経済指標ダッシュボードからデータを保存して分析してみようと思います。 たくさんデータがあるのですが、日経平均(Nikk)、ドル円相場(Yen)、消費支出(Cons)と鉱工業生産指数(IP)の4つを取得しました。最終的には、ドル…
www.crosshyou.info の続きです。 前回は、PC1台当りの生徒数は年度による違いは無いとわかりました。(2013年度と2015年度を比較しました。) 今回は都道府県別のPC1台当りの生徒数を反応変数にして、人口、面積、生徒数の3つの変数で重回帰分析をしてみます…
www.crosshyou.info の続きです。 今回は Statistics: An Introduction Using R 作者: Michael J. Crawley 出版社/メーカー: Wiley 発売日: 2014/11/24 メディア: ペーパーバック この商品を含むブログを見る を参考にして、PC1台当りの生徒数の平均値と95%…
今回は都道府県別の教育用コンピュータ１台当りの生徒数を調べます。 政府統計の総合窓口e-Statからデータを取得します。 www.e-stat.go.jp がe-Statのウェブサイトです。 47都道府県を選択します。 総人口(人), 総面積(ha), 県内総生産(百万円), 教育用コン…
R言語のsummary関数は便利なのですが、標準偏差や変動係数が算出されないことが不満でした。探せばそういう関数あるのかもしれないですが、関数を作る練習として、自作してみました。 # 自動でsummary関数、標準偏差、変動係数を出力する関数 jidou_sammary …
www.crosshyou.info の続きです。 今回は、課税所得伸び率を人口伸び率、面積伸び率、行政職員伸び率の3つの変数で重回帰分析をしてみたいと思います。 まずは、散布図を見てみます。plot関数ですね。 人口伸び率はかなり相関が強そうです。面積伸び率はなさ…
www.crosshyou.info の続きです。今回は、2005年度から2015年度の変化率を調べてみようと思います。 まずは面積の伸びを調べます。 このように、2005年度の面積、2015年度の面積のベクトルを作成します。head関数ではじめの6つを表示しました。五十音順に並…
県庁所在市・政令指定都市・東京都特別区の人口・面積・課税対象所得・一般行政職員人数の分析をします。 まずは、政府統計の総合窓口(www.s-stat.go.jp)からデータを取得します。 このように、特別区、県庁所在市、政令指定都市を指定します。74地域ありま…
www.crosshyou.info の続きです。 前回までの分析で、一人当りの医療費は、一人当り県内総生産がちいさいほど大きくなることがわかりました。今回は、人口密度を説明変数に加えて分析してみます。 まずは、上のように、平均人口密度を表す変数(avgMits)を作…
