今回は、上の画像、UCI Machine Learning Repositoryの中にあった、Real Estate valuation data setのデータを分析しようと思います。 URLは、 https://archive.ics.uci.edu/ml/datasets/Real+estate+valuation+data+set です。 Citation requestがYeh, I. C…
www.crosshyou.info の続きです。 今回は前日比を規制有りの空売り比率と規制無しの空売り比率で回帰分析をしてみたいと思います。R言語のlm関数を使います。 規制無しの空売り比率と前日比の相関係数が、0.3965843ですので、この相関係数をどれだけ上回る回…
www.crosshyou.info の続きです。 今回は、空売り比率と前日比の相関関係を調べてみます。 まずは、R言語のplot関数で散布図を描いてみます。 規制有りと前日比は相関なさそうで、規制有りと前日比は少し相関ありそうです。 cor.test関数で相関係数を調べて…
www.crosshyou.infoの続きです。 今回は、規制有り、規制無し、前日比の各データの分布が正規分布かどうかを調べます。 skewとkurtosisを計算してこの値が有意に0と違っていれば正規分布では無いとわかります。 skewは分子は、 m3 = sum((y - ybar)^3) / n …
www.crosshyou.info の続きです。 今回は各データのVariance(分散)を調べてみます。 まずは、規制有りのVarianceと規制無しのVarianceを調べてみます。 var関数でVarianceはわかります。規制有りのVarianceは33.4、規制無しのVarianceは4.87です。かなり違い…
www.crosshyou.infoの続きです。 今回は、Michael J. Crawley著 Statistics An introduction using Rの第3章、Central Tendencyを参考にして分析の練習をしてみます。 まずは、hist関数でヒストグラムを描きます。 規制有りは左の裾野が広く、規制無しと前日…
今回は東証の業種別空売りデータを分析しようと思います。 空売り集計 | 日本取引所グループ このサイトからPDFファイルをダウンロードしました。 こういうデータです。 それと、 リアルタイム株価指数値一覧 | 日本取引所グループ このサイトの業種別指数の…
www.crosshyou.info の続きです。 今回は、R言語のglm関数で退院済フラグをresponse variableに、性別と年代をexplanatory variablesにしてロジスティクス回帰分析をしてみます。 まず、退院済フラグ、性別と年代だけのデータフレームを作ります。 head関数…
www.crosshyou.info の続きです。 今回は、性別と退院済フラグのクロス表分析をします。 性別の調査中と不明を削除しましょう。 table関数でクロス表を作ります。 男性のほうが退院済フラグが1の比率が多いようです。 prop.table関数で見てみましょう。 男性…
www.crosshyou.info の続きです。 今回は年代と退院済フラグでクロス表を作成します。 年代の不明を削除して、100歳以上を90代と一緒にして、90代以上にしましょう。 クロス表を作成します。 prop.table関数で比率をみてみます。 年代が上がるほど、退院済フ…
www.crosshyou.info の続きです。 今回は、年代と性別でクロス表分析をします。 性別の調査中と不明、年代の不明を削除しましょう。 まず、年代の不明を削除しました。 次に、性別の調査中と不明を削除します。 性別の100歳以上を90代と合計して、90代以上と…
www.crosshyou.infoの続きです。 今回は曜日 x 退院済フラグでクロス表を作ります。 これは、不要なデータを削除する必要はないですね。 prop.table関数で比率を見ましょう。 金曜日が退院済フラグが0、退院していない人が多いです。 barplot関数でグラフに…
www.crosshyou.infoの続きです。 今回は、曜日と性別でクロス表分析をします。 まずは普通にtable関数で集計します。 調査中と不明は削除しましょう。 集計用に作成した曜日と性別で集計します。 prop.table関数で男女の比率を見てみます。 prop.table関数で…
www.crosshyou.info の続きです。 前回は、1変数について陽性者数を集計しました。 今回は、2変数で集計してみます。 曜日 x 年代、曜日 x 性別、曜日 x 退院済、年代 x 性別、年代 x 退院済、性別 x 退院済、の6つです。 10歳未満と10代を統合、90代と100歳…
www.crosshyou.infoの続きです。 今回はR言語のtable関数を使って集計をしてみます。 まずは公表_日付で集計しました。 数が多いですから、plot関数でグラフにします。 テレビでよく見るようなグラフになりました。 曜日で集計してみます。 曜日の順番がバラ…
知人から「日本学術会議ニュース・メール」というのが転送されてきました。 各都道府県データリストのURLが記載されていましたので、アクセスしてみました。 Links to Raw Case Data by Each Prefecture ｜日本学術会議 このように、東京都の陽性者数のデー…
www.crosshyou.info の続きです。 GDPギャップを回帰分析してみましょう。 R言語のlm関数を使います p-valueは3.126e-07と0.05よりも小さいので有意なモデルです。 InterceptはP値が0.05よりも大きいですね。削除してみます。 update関数で、-1とすれば切片…
www.crosshyou.infoの続きです。 各変数間の散布図マトリックスをみてみます。 R言語のpairs関数で散布図マトリックスを描きます。散らばっているというよりは、ぐにゃぐにゃしている感じですね。 cor関数で相関係数マトリックスをみてみます。 GDPギャップ…
www.crosshyou.infoの続きです。 前回は四半期ごとの違いはありませんでした。 今回は、1980年台、1990年台、2000年台、2010年台と10年ごとに区切った年台ごとに差があるかどうかを調べます。 まず、年台を表す変数を作成します。 こうして作成したdf$Nendai…
www.crosshyou.infoの続きです。今回は第1四半期、第2四半期、第3四半期、第4四半期で差があるのかどうかを調べようと思います。 R言語のaov関数とsummary関数を使います。 GDPギャップからみていきます。 Pr(>F)が0.997ですから四半期によってGDPギャップに…
www.crosshyou.info の続きです。 今回はそれぞれのデータのグラフを描いてみます。plot関数を使います。 GDPギャップのグラフです。マイナスの時が多いようです。 GDPの潜在成長率のグラフです。最近は1%のちょっと下あたりです。 TFPのグラフです。長期的…
今回から、GDPギャップのデータを分析しようと思います。 月例経済報告 - 内閣府 のサイトからExcelファイルを取得しました。 このようなファイルでした。 このようなCSVファイルにしました。これをR言語のread.csv関数で読み込みます。 read.csv関数でデー…
www.crosshyou.info 今回はchgperInc, 1人当り課税所得の変化幅を他の変数で回帰分析してみようと思います。 前回の散布図マトリックスから、直線的な相関関係というよりは複雑な関係がありそうなので、交差項や2乗項を入れて回帰分析してみます。 step関数…
www.crosshyou.info の続きです。 今回は前回作成した変化幅の変数をヒストグラムや箱ひげ図でグラフにします。 まず、作成した5つの変化幅の変数を一つのデータフレームにまとめます。R言語のdata.frame関数を使います。 summary関数で各変数の平均値などを…
www.crosshyou.info の続きです。 今回は2015年と昔の年を比較して、どのように変化してきたか調べてみます。 まずはYearを確認します。 1985年度が一番古いです。1985年度と2015年度でどれだけ変化があったかを調べます。 まずは人口密度から調べます。 人…
www.crosshyou.infoの続きです。 今回は、R言語のglm関数でロジスティクス回帰分析をして、東京23区かそうじゃないかを見分けようと思います。 predict関数で予測をして、table関数で実際の値と比較してみます。 予測で1と予測したのが23区です。1と1のとこ…
www.crosshyou.info の続きです。 今回はlogperInc(1人当り課税所得を対数変換した値)をlogMitsu(人口密度を対数変換した値), Ratio1(第1次産業就業者数の割合), Ratio2(第2次産業就業者数の割合), Tko(東京23区かどうか)の4つの変数で回帰分析をしてみます…
www.crosshyou.info の続きです。 前回の分析で東京都の23区は他とはかなり違いがあることがわかりましたので、まず、東京都23区かどうかという変数を作ろうと思います。R言語のgrepl関数を使いました。 東京都の中には八王子市も含まれていますので、八王子…
www.crosshyou.info の続きです。 今回は前回作成した、人口密度、1人当り課税所得、就業者の割合のランキングをみてみようと思います。最新の年度でみてみようと思いますので、まず、table関数で年度を調べます。 2015年度が最新です。2015年度のデータでラ…
今回からは、全国主要都市の第1次、第2次、第3次産業の就業者数のデータを分析しようと思います。 データは、政府統計の総合窓口e-statから取得しました。 www.e-stat.go.jp選択した都市は、 特別区、県庁所在市、政令指定都市、中核市の合計102都市です。 …
