Publications by ayumuR
IV01
1 はじめに Cunningham (2021) の第7章に従って、Card (1993) のデータを用いて、教育年数educが所得lwageに与える影響を推定することを試みる。ここで、教育年数は、内生変数であると考えられ、撹乱項と相関するおそれがある。そのため、操作変数として、自宅がある�...
1641 sym R (5523 sym/23 pcs) 5 img 2 tbl
IV
1 はじめに Cunningham (2021) の第7章に従って、Card (1993) のデータを用いて、教育年数educが所得lwageに与える影響を推定することを試みる。ここで、教育年数は、内生変数であると考えられ、撹乱項と相関するおそれがある。そのため、操作変数として、自宅がある�...
1471 sym R (4929 sym/18 pcs) 4 img 2 tbl
gravity_fixest_es
1 はじめに 本ページでは、政策が２国間の貿易額に与える影響が経時的にどう変化するか分析するため、単純な差の差 (DiD: difference-in-differences) 分析を拡張したevent study分析を実施する。 Anderson and Van Wincoop (2003) が指摘した多角的貿易抵抗指数（物価効果）を制�...
3158 sym R (2252 sym/13 pcs) 2 img 2 tbl
gravity_fixest_twfe
1 はじめに 本ページでは、政策が２国間の貿易額に与える影響を動学的に分析するため、差の差 (DiD: difference-in-differences) 分析、event study分析を実施する。単純なOLSによるDiD分析から始めて、 Anderson and Van Wincoop (2003) が指摘した多角的貿易抵抗指数（物価効果）�...
3432 sym R (4754 sym/22 pcs) 1 img 1 tbl
gravity_fixest_panel
1 はじめに Anderson and Van Wincoop (2003) が指摘した多角的貿易抵抗指数（物価効果）を制御するために、1時点のクロスセクションの貿易データであれば、輸出国のGDPの代わりに輸出国固定効果、輸入国のGDPの代わりに輸入国固定効果を説明変数に加えれば良い。２�...
3876 sym R (3610 sym/20 pcs)
gravity_fixest
1 はじめに 現代的な重力方程式の推定では、多角的貿易抵抗指数（物価効果）を制御するために、輸出国のGDPの代わりに輸出国固定効果、輸入国のGDPの代わりに輸入国固定効果が用いられる。また、マクロショックを制御するために、年固定効果を加えること�...
2994 sym R (3786 sym/15 pcs)
RTips
Rのplot関数でグラフ作成すると、x軸のラベルに不要な小数点がついてしまうことがある。例）2005.0, 2005.5, 2006.0 year <- c(2005,2006,2007) y <- c(1,2,3) plot(year,y) x軸の変数の値を増やすことで、小数点を避けられる場合もあるようだが、以下のようにうまくいかない場合�...
262 sym 3 img
leaflet_basic
1 データ trade.csv 1.1 貿易データ origins: USA, Australia, China, New Zealand, Germanyの5カ国 destinations: Japanのみ import: 上記5カ国から日本への輸出額 longitude.x latitude.x: 輸出国の緯度経度 longitude.y latitude.y: 輸入国（日本）の緯度経度。(137.96, 37.54)。 library(readr) trade <- read_csv("t...
921 sym R (1844 sym/10 pcs) 1 img
Rbasic03
1 棒グラフ library(readxl) asia <- read_excel("asia.xlsx") barplot(asia$ctfp,names.arg = asia$countrycode) barplot(asia$pop,names.arg=asia$countrycode) 2 散布図 人口とGDPの関係 plot(asia$pop,asia$rgdpe) plot(asia$pop,asia$rgdpe,log="xy") plot(asia$pop,asia$rgdpe,log="xy",pch=4) plot(asia$pop,asia$rgdpe,log = "xy",xlab="Populatio...
205 sym R (1621 sym/14 pcs) 8 img
Rbasic01
1 データの読み込み library(readxl) time <- read_excel("time.xlsx") 2 記述統計 # 平均 mean(time$time) ## [1] 62.35294 # 中央値 median(time$time) ## [1] 60 # 最大値 max(time$time) ## [1] 120 # 最小値 min(time$time) ## [1] 25 3 ヒストグラム with hist hist(time$time) 4 データの読み込み library(readxl) time2 <- re...
252 sym R (1196 sym/34 pcs) 12 img