データ分析
www.crosshyou.info の続きです。 今回はウグイスとスズメについて調べてみます。 まずは、ウグイス、スズメだけの作業用のデータフレームを作成します。 ウグイスは124、スズメは121の標本があります。 それでは、都道府県別の標本数を見てみます。table関…
www.crosshyou.info の続きです。 今回はメジロとヒヨドリでは、標本の年、月、都道府県に差があるか調べてみます。 まずは、メジロとヒヨドリだけの作業用のデータフレームを作ります。 Nameには、アオアシシギなど不要のファクタ水準があるので、削除しま…
www.crosshyou.info の続きです。 今回は、Moku(目), Zoku(属), Shu(種)について調べます。 まず、Mokuは何種類あるか、標本数の多いMokuは何か調べてみましょう。 levels関数とlevels関数でMokuの種類数を調べます。 Mokuは214種類ですね。どのMokuが一番多…
国立科学博物館の標本データベースからCSVファイルをダウンロードできることを発見しましたので、データをダウンロードして遊んでみようと思います。 詳細検索(鳥類) | 標本・資料統合データベース から日本の標本を検索してみました。 こんなCSVファイルが…
www.crosshyou.info の続きです。 今回はbarplot関数の練習をしてみようと思います。 習うより慣れろの精神でいろいろやってみます。 まずはデータを確認します。 １列目は文字(ファクター)で、２、３、４列目が数値データですね。 まずは、barplot(df)とし…
安倍晋三首相の在任日数が憲政史上最長になりました。 そこで今回は、各報道機関の世論調査、内閣支持率がどうなっているかを調べてみました。 こちらは、NHKの世論調査です。 http://www.nhk.or.jp/senkyo/shijiritsu/ 2019年11月8日から10日の調査で、1253…
www.crosshyou.info の続きです。 tapply関数でメーカー別の販売台数を見てみましょう。 トヨタが204,628台でダントツですね。2位がホンダの63,507台、3位が日産の62,943台です。 メーカーごとの平均順位はどうでしょうか？ トヨタが15.78で1位ですね。2位が…
www.crosshyou.info の続きです。 今回は販売台数をグラフにしてみましょう。 まずは、小さい順グラフ、ヒストグラム、箱ひげ図の３つのグラフを一度に作成する関数を作ります。 gpという名前の関数を作りました。 それでは販売台数をグラフにしてみます。 …
www.crosshyou.info の続きです。 今回は標準偏差、標準誤差、変動係数などを計算する関数を作成してみます。 まずは、前回保存したCSVファイルを読み込んでsummary関数で様子を見てみます。 こうなりました。summary関数だと数値データは最小値、第1分位、…
今回は、R言語のXMLパッケージのreadHTMLTable関数でWebの表形式のデータを直接読込んでみます。 参考図書は、Rクックブックです。 Rクックブック 作者: Paul Teetor,大橋真也,木下哲也 出版社/メーカー: オライリージャパン 発売日: 2011/12/22 メディア: …
www.crosshyou.info の続きです。 今回は、県内総生産当りの凶悪犯認知件数を人口と可住地面積で回帰分析しようと思います。 まずは、県内総生産当りの凶悪犯認知件数を算出します。 avgGDPは百万円単位なので、百万をかけて1円当りの件数にしています。富山…
www.crosshyou.info の続きです。 前回作成した、可住地面積10万ha当りの凶悪犯認知件数を人口と県内総生産で回帰分析してみます。 まずはplot関数で散布図を描いてみます。 右側の対数変換した後の散布図のほうが分析するにはよさそうですね。 県内総生産と…
www.crosshyou.info の続きです。 今回は可住地面積当りの凶悪犯認知件数を計算してみます。 10万ha当りの件数を算出しました。秋田県は8.4件、大阪府は780.5件です。100倍近い差がありますね。 グラフで分布をみてみます。 大きな値に外れ値がいっぱいある…
www.crosshyou.info の続きです。 今回は、ファイ係数を計算する関数を作ってみます。 実は、 www.crosshyou.info のときにも関数を作ってみたのですが、何故かエラーになったんです。 そのとき作った関数が これです。 そして、この関数を実行したら、 と「…
www.crosshyou.info の続きです。 今回は対数値のavgArea, avgGDPで回帰分析してみます。 log(avgArea):log(avgGDP)の項目はいらないですね。削除します。 model2を採用します。 I(log(avgGDP)^2)はいらないですね。 model3を採用します。 I(log(avgArea)^2)…
www.crosshyou.info の続きです。 今回は、人口10万人当りの凶悪犯認知件数を算出して、これを可住地面積と県内総生産額で回帰分析してみます。 大阪府は人口10万人当りの件数でも一番多いですね。 一番少ないのは山形県です。 グラフにしてみます。 大阪府…
今回は都道府県別の凶悪犯認知件数を分析します。 e-stat, 政府統計の総合窓口からデータを取得しました。 説明変数として、総人口、可住地面積、県内総生産額も一緒に取得しました。 このようなExcelファイルを作成しました。これをread.csv関数でR言語に読…
www.crosshyou.info の続きです。 今回は勝利数を引き分け数で回帰分析してみます。 引き分け数が多いチームほど勝利数が多ければ、引き分けが多いのは粘り強さの証になりますし、逆に引き分けが多いチームほど勝利数が少なければ、引き分けが多いのは勝ちき…
www.crosshyou.info の続きです。 今回は 2 x 2 のクロス表で関連性の強さを見る指標、ファイ係数を計算しようと思います。 bellcurve.jpこのサイトに このように計算式が載っていました。これを使って計算します。 まずは、Aクラス・Bクラス x セリーグ・パ…
www.crosshyou.info の続きです。 今回はクロス表を作って、カイ二乗検定をしてみます。 勝利数を Aクラス・Bクラス x セリーグ・パリーグ、 Aクラス・Bクラス x 関東・非関東 セリーグ・パリーグ x 関東・非関東 この３つのクロス表で分類してカイ二乗検定…
www.crosshyou.info の続きです。 前回は、セリーグ、パリーグで勝利数の平均値に有意な違いは無いことがわかりました。今回は、関東の球団かそうじゃないか、AクラスかBクラスかで勝利数に違いがあるかを調べます。 var.test関数でvarianceが同じかどうか確…
今回はプロ野球の勝敗データを分析してみようと思います。 npb.jpこのサイトにあった の画像にある勝敗データを使おうと思います。 まず、Excelにデータを転記しました。 関東にある球団か、そうでないかという変数、Kantouと各リーグ上位をAクラス、下位をB…
www.crosshyou.info 今回は前回の続きで、こんどは人口、面積、県内総生産の対数値、logPop, logArea, logGDPで回帰分析をしてみたいと思います。前回は北海道の面積が他と比べて圧倒的に大きかったので、その影響が大きく出ていました。今回はどうでしょう…
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によって短観の値に違いがあるか見てみましょう。 黒い点が製造業で黒い水平線が…
