2023年6月3日土曜日

合計特殊出生率(1)

奥村晴彦先生が,Rで合計特殊出生率グラフを描いていたのでまね(写経)をしてみた。

厚生労働省の発表では,2022年の合計特殊出生率が1.26となって(前年は1.30),2005年に並んで過去最低を記録した,というのがニュースになったからかな。

写経は成功したが,日本語タイトルが文字化けしていた。ネットで検索すると解決できた。次に,グリッドラインを入れるべく検索したが要領を得ないし,なんだかずれてしまう。GPT-4に相談したところ2回目で正解にたどりつけた。最後に,奥村先生はクリップボードからデータを入力するという高等技術を使っていたので,低レベルのテキストファイル入力に切り替えた。自分で考える力がどんどん喪失されていく・・・

その結果が次のRプログラムと図である。

# X = read.table(pipe("pbpaste"), header=TRUE) (奥村先生のテクニック for macOS)

X <- read.table('/Users/koshi/Desktop/birth.txt', header =TRUE)
par(family = "HiraKakuProN-W3")
plot(X$年, X\$合計特殊出生率, type="o", pch=16, xlab="", ylab="")
t = c(range(X\$合計特殊出生率), X$合計特殊出生率[length(X[,1])])
axis(4, t, t)
title("合計特殊出生率", line=0.5)

# x軸とy軸のメモリの位置を手動で指定
x_ticks <- pretty(X\$年, 11)  # 11個のメモリを生成
y_ticks <- pretty(X\$合計特殊出生率, 9)  # 9個のメモリを生成

# x軸に対してグリッド線を描画
for (i in x_ticks) {
  abline(v = i, col = "lightgray", lty = "dotted")
}

# y軸に対してグリッド線を描画
for (i in y_ticks) {
  abline(h = i, col = "lightgray", lty = "dotted")
}


図:合計特殊出生率の推移(1970-2022)


0 件のコメント: