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

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

2023-02-01から1ヶ月間の記事一覧

OECD Gender wage pay gap data analysis 4 - Estonia is the outlier from the point of gender wage pay gap view.

UnsplashのRenato Pozziが撮影した写真 www.crosshyou.info This post is following of the above post. I make small data frame that contains only 2010 and 2014 observations. First, I use filter() function to get only 2010 and 2014 data, then I …

読書記録 - 「菅原道真 学者政治家の栄光と没落」 滝川幸司 著 中公新書

菅原道真-学者政治家の栄光と没落 (中公新書 (2559)) 作者:滝川 幸司 中央公論新社 Amazon 菅原道真は太宰府に左遷され、死後、怨霊になったということは知っていいましたが、 その前のことは何も知らなかったです。 この本は、その前のことが書かれていて良…

OECD Gender wage pay gap data analysis 3 - Bootstrapping Confidence Interval and Traditional Method p-value with R

UnsplashのMaria Tejadaが撮影した写真 www.crosshyou.info This post is following of the above post.I will calculate confidence interval. First, I make bootstrapping confidence interval. I use R infer package. I use specify(), generte() and c…

OECD Gender wage gap data analysis 2 - Statistical Inference with Infer package, mean difference, calculation p-value with simulation based method

UnsplashのJames Wainscoatが撮影した写真 www.crosshyou.info This post is following of the above post. In the previous post I see mean of EMPLOYEE is 17.9 and mean of SELFEMPLOYED is 30.6. Let's see whether the difference is statistically si…

OECD Gender wage gap data analysis 1 - Load CSV file data into R

UnsplashのKumiko SHIMIZUが撮影した写真 Hello. In this post, I will analyze Gender wage gap of OECD data with R. First, I download CSV file like belo from OECD web site, https://data.oecd.org/earnwage/gender-wage-gap.htm I use R to analyze …

都道府県別の保護統計調査のデータ分析11 - Rのplmパッケージを使ってパネルデータ分析

UnsplashのKrzysztof Niewolnyが撮影した写真 www.crosshyou.info の続きです。 今回は、Rのplmパッケージを使って、パネルデータで回帰分析をしてみます。 まず、library()関数でplmパッケージを読み込みます。 前回作成した、df3のデータフレームをパネル…

都道府県別の保護統計調査のデータ分析10 - RでPooled Cross-Sectionsのデータ分析

UnsplashのYoksel Zokが撮影した写真 www.crosshyou.info の続きです。 今回からは、Rでパネルデータ分析をしてみます。 まず、パネルデータとして、1997年、2008年、2021年の3年間だけのデータフレームを作ってみます。 filter()関数で1997年、2008年、202…

読書記録 - 「人新世の科学: ニュー・エコロジーが開く地平」 オズワルド・シュミッツ 著 岩波新書

人新世の科学: ニュー・エコロジーがひらく地平 (岩波新書 新赤版 1922) 作者:オズワルド・シュミッツ 岩波書店 Amazon 従来の生態学は自然界と人間界を別々のものとして扱い、調べていたが、新しい生態学、ニュー・エコロジーでは、人間も生態系を構築する…

都道府県別の保護統計調査のデータの分析9 - Rで時系列の回帰分析のSerial Correlationをテストする。

UnsplashのLeo SERRATが撮影した写真 www.crosshyou.info の続きです。 前回はRのdynlmパッケージを使い、時系列データの回帰分析をして、r_one: 1号観察割合をr_four: 4号観察の割合で回帰分析してみました。 今回は、時系列データの回帰分析で考慮しなけ…

都道府県別の保護統計調査のデータの分析8 - R言語のdynlmパッケージで時系列データの回帰分析

UnsplashのTom Podmoreが撮影した写真 www.crosshyou.info 上のブログの続きです。 今回は、保護統計調査のデータを時系列データとして扱って、分析をしてみます。 これまでの回帰分析では、r_one: 1号観察の割合をr_four: 4号観察の割合で回帰分析してい…

東京証券取引所の市場区分再編時の区分選択にガバナンス要因が与える影響

今回の記事は、いつもとは違っています。ある学術誌に投稿したのですが、査読の結果、掲載不可になったものです。掲載不可なので、欠陥や問題点などが多くあるのかと思います。しかし、昨年の3月くらいから、コツコツと書いてきた論文なので、とりあえず、記…

都道府県別の保護統計調査のデータの分析7 - R言語で回帰分析モデルの不均一分散をチェックする。

UnsplashのMaitheli Maitraが撮影した写真 www.crosshyou.info の続きです。 前回は、回帰分析をしました。今回はその続きで、回帰分析したモデルの残差が均一分散かどうかをチェックします。 lmtestパッケージのbptest()関数で簡単にチェックできます。 bpt…

都道府県別の保護統計調査のデータ分析6 - R言語のlm()関数で回帰分析モデルを作成し、stagrgazer()関数でモデルを比較する。

UnsplashのGabriel Garcia Marengoが撮影した写真 www.crosshyou.info の続きです。 今回は回帰分析をしています。被説明変数は、r_one:1号観察の割合、説明変数は、r_two: 2号観察の割合と、r_four: 4号観察の割合にしてみます。 まず、2021年のデータだ…