合計特殊出生率（１）

奥村晴彦先生が，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）