相関係数

同義語:Pearsonの相関係数積率相関係数correlation coefficient

概要

相関係数(correlation coefficient)とは、分散が有限で正の確率変数 $X,Y$ について、共分散を標準偏差の積で割った値 $\rho=\mathrm{Cov}[X,Y]/\sqrt{V[X]V[Y]}$ のことで、$-1\le\rho\le1$ を満たす。$Y$ を 1 次式 $a+bX$ で予測するときの 2 乗誤差の最小値は $V[Y](1-\rho^2)$ であり、$\lvert\rho\rvert=1$ はほとんど確実に 1 次式の関係があることと同値である。共分散は同時分布関数 $H$ と周辺分布関数の積 $FG$ の差の積分に等しく、ここから独立なら無相関であること、正の象限従属な組では無相関と独立が一致することが出る。データの相関係数は、$(X,Y)$ と同じ分布に従う独立な観測を増やすと確率 $1$ で $\rho$ に収束する。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{div}[0]{\mathbin{÷}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 期待値, 分散, Cauchy–Schwarzの不等式

2 つの確率変数 $X,Y$ がどの程度「1 次式の関係」に近いかを $-1$ から $1$ までの数で表したものが 相関係数 である。高校数学ではデータについて定義した(相関係数(高校数学))。この記事では、確率変数の相関係数を一般の確率空間で定め、それが何を測っているかを 2 つの定理で示す。
1 つめは、$Y$ を $X$ の 1 次式 $a+bX$ で予測するときの 2 乗誤差の最小値が $V[Y](1-\rho^2)$ になること(主定理 1)である。$\rho^2$ は「1 次式で説明できる $Y$ の分散の割合」であり、$\lvert\rho\rvert=1$ は誤差 $0$ の 1 次式があることに当たる。2 つめは、共分散が同時分布関数と周辺分布関数の差の積分で書けること(主定理 2)である。ここから、独立なら無相関であること、一方が大きいと他方も大きくなりやすい組(正の象限従属)では相関係数が $0$ 以上で、$0$ なら独立になることが出る。最後に、データの相関係数が確率変数の相関係数に収束することを示す。

定義

この記事を通して、$X,Y$ は確率空間 $(\Omega,\mathcal F,P)$ の上の実数値の確率変数とする。

共分散と相関係数

$E[X^2]<\infty$、$E[Y^2]<\infty$ とする。
$$ \mathrm{Cov}[X,Y]:=E\bigl[(X-E[X])(Y-E[Y])\bigr]=E[XY]-E[X]\,E[Y] $$
を $X$ と $Y$ の 共分散 という。さらに $V[X]>0$、$V[Y]>0$ のとき
$$ \rho(X,Y):=\frac{\mathrm{Cov}[X,Y]}{\sqrt{V[X]}\sqrt{V[Y]}} $$
を $X$ と $Y$ の 相関係数(correlation coefficient)という。Pearson の相関係数、積率相関係数ともいう。$\mathrm{Cov}[X,Y]=0$ のとき、$X$ と $Y$ は 無相関 であるという。

$\lvert XY\rvert\le\frac12(X^2+Y^2)$ から $XY$ は可積分であり、定義の 2 つの式は期待値の線形性で等しい(分散 の定義「共分散」と同じ式である)。共分散と相関係数は $(X,Y)$ の同時分布だけで決まる。
Cauchy–Schwarzの不等式 の命題「期待値の形と相関係数」により、$-1\le\rho(X,Y)\le1$ であり、$\lvert\rho(X,Y)\rvert=1$ となるのは、実数 $c\ne0$、$d$ で、ほとんど確実に $Y=cX+d$ となるものがあるときに限る(このとき $\rho$ は $c$ と同じ符号をもつ)。この記事ではこの命題を引用して使う。

標本相関係数

$n\ge2$ 個の実数の組 $(x_1,y_1),\dots,(x_n,y_n)$ について、$\bar x:=\frac1n\sum_ix_i$、$\bar y:=\frac1n\sum_iy_i$ とおき、$x_i$ がすべて等しくはなく、$y_i$ もすべて等しくはないとき
$$ r:=\frac{\sum_{i=1}^n(x_i-\bar x)(y_i-\bar y)}{\sqrt{\sum_{i=1}^n(x_i-\bar x)^2}\sqrt{\sum_{i=1}^n(y_i-\bar y)^2}} $$
をこのデータの 標本相関係数 という。

各組 $(x_i,y_i)$ に確率 $\frac1n$ を割り当てた分布(経験分布)に従う確率変数の組 $(X,Y)$ をとると、$E[X]=\bar x$、$\mathrm{Cov}[X,Y]=\frac1n\sum_i(x_i-\bar x)(y_i-\bar y)$ などとなるので、$r=\rho(X,Y)$ である。したがって、確率変数について示すことはデータについてもそのまま成り立つ。高校数学の相関係数はこの $r$ である。

1 次変換での不変性

$a,c\ne0$、$b,d$ を実数とすると、$\rho(aX+b,\,cY+d)=\operatorname{sgn}(ac)\,\rho(X,Y)$ である。ここで $\operatorname{sgn}(ac)$ は $ac$ の符号($\pm1$)である。

$aX+b-E[aX+b]=a(X-E[X])$ なので $\mathrm{Cov}[aX+b,cY+d]=ac\,\mathrm{Cov}[X,Y]$、$V[aX+b]=a^2V[X]$、$V[cY+d]=c^2V[Y]$ である。よって $\rho(aX+b,cY+d)=\dfrac{ac}{\lvert a\rvert\lvert c\rvert}\rho(X,Y)$ である。$\square$

たとえば気温を摂氏から華氏に直しても($x\mapsto1.8x+32$)、相関係数は変わらない。一方、$e^X$ のような 1 次でない単調な変換では一般に変わる(下の ex-cor-lognormal)。

主定理 1:1 次式による最良の予測

1 次式による予測の 2 乗誤差

$E[X^2]<\infty$、$E[Y^2]<\infty$、$V[X]>0$ とし、
$$ \beta^*:=\frac{\mathrm{Cov}[X,Y]}{V[X]},\qquad\alpha^*:=E[Y]-\beta^*E[X] $$
とおく。

  1. すべての実数 $a,b$ について $E\bigl[(Y-a-bX)^2\bigr]\ge E\bigl[(Y-\alpha^*-\beta^*X)^2\bigr]$ であり、等号は $(a,b)=(\alpha^*,\beta^*)$ のときに限る。
  2. 誤差 $U:=Y-\alpha^*-\beta^*X$ は $E[U]=0$、$\mathrm{Cov}[U,X]=0$ を満たす。
  3. 最小値は $E[U^2]=V[Y]-\dfrac{\mathrm{Cov}[X,Y]^2}{V[X]}$ である。とくに $V[Y]>0$ なら $E[U^2]=V[Y]\bigl(1-\rho(X,Y)^2\bigr)$ である。

2:$E[U]=E[Y]-\alpha^*-\beta^*E[X]=0$ である。$\mathrm{Cov}[U,X]=\mathrm{Cov}[Y,X]-\beta^*V[X]=0$ である(共分散は各変数について線形で、定数との共分散は $0$)。
1:任意の $a,b$ について $Y-a-bX=U+(\alpha^*-a)+(\beta^*-b)X$ である。$c:=\alpha^*-a+(\beta^*-b)E[X]$、$\delta:=\beta^*-b$ とおくと、$Y-a-bX=U+c+\delta(X-E[X])$ と書ける。2 乗して期待値をとると、2 の $E[U]=0$ と $E[U(X-E[X])]=\mathrm{Cov}[U,X]=0$、$E\bigl[X-E[X]\bigr]=0$ により交差項が消え、
$$ E\bigl[(Y-a-bX)^2\bigr]=E[U^2]+c^2+\delta^2V[X] $$
を得る。$V[X]>0$ なので右辺は $E[U^2]$ 以上で、等号は $c=\delta=0$、すなわち $b=\beta^*$ かつ $a=\alpha^*$ のときに限る。
3:$U=(Y-E[Y])-\beta^*(X-E[X])$ なので
$$ E[U^2]=V[Y]-2\beta^*\mathrm{Cov}[X,Y]+(\beta^*)^2V[X]=V[Y]-\frac{\mathrm{Cov}[X,Y]^2}{V[X]} $$
である。$V[Y]>0$ なら、これは $V[Y]\bigl(1-\frac{\mathrm{Cov}[X,Y]^2}{V[X]V[Y]}\bigr)=V[Y](1-\rho^2)$ である。$\square$

$E[U^2]\ge0$ なので、3 から改めて $\rho^2\le1$ が出る。また $\rho^2=1$ は $E[U^2]=0$、すなわちほとんど確実に $Y=\alpha^*+\beta^*X$ と同値である。$\rho^2$ は、$Y$ の分散 $V[Y]$ のうち $X$ の 1 次式で説明できる部分 $V[\alpha^*+\beta^*X]=\rho^2V[Y]$ の割合であり、決定係数 と呼ばれる。
経験分布に当てはめると、直線 $y=\alpha^*+\beta^*x$ はデータの回帰直線で、残差の 2 乗の和の最小値は $ns_y^2(1-r^2)$ である($s_y^2$ は $y$ の分散)。これは 散布図と回帰直線 の系「回帰直線の性質」の内容であり、説明変数が多い場合への一般化が 回帰分析 の定理「平方和の分解と決定係数」である。

2 つのさいころ

公平なさいころを 2 回独立に投げ、1 回目の目を $X$、2 回目の目を $Z$、和を $Y:=X+Z$ とする。$V[X]=V[Z]=\frac{35}{12}$、$\mathrm{Cov}[X,Z]=0$(分散 の補題「独立な確率変数の積の期待値」)から、$\mathrm{Cov}[X,Y]=V[X]=\frac{35}{12}$、$V[Y]=\frac{35}6$ であり
$$ \rho(X,Y)=\frac{35/12}{\sqrt{35/12}\sqrt{35/6}}=\frac1{\sqrt2}=0.707\ldots $$
である。thm-cor-prediction の最良の予測は $\beta^*=1$、$\alpha^*=7-\frac72=\frac72$ で $Y\approx\frac72+X$、2 乗誤差の最小値は $V[Y](1-\frac12)=\frac{35}{12}=V[Z]$ である。予測できない部分は、2 回目の目の散らばりそのものである。

主定理 2:分布関数による共分散の表示

$(X,Y)$ の 同時分布関数 を $H(s,t):=P(X\le s,\,Y\le t)$、周辺の分布関数を $F(s):=P(X\le s)$、$G(t):=P(Y\le t)$ とする。

共分散の分布関数による表示

$E[X^2]<\infty$、$E[Y^2]<\infty$ とする。関数 $(s,t)\mapsto H(s,t)-F(s)G(t)$ は $\mathbb R^2$ 上で Lebesgue 可積分であり
$$ \mathrm{Cov}[X,Y]=\int_{-\infty}^\infty\int_{-\infty}^\infty\bigl(H(s,t)-F(s)G(t)\bigr)\,ds\,dt $$
である。

この式は Hoeffding の公式と呼ばれる。証明では、2 つの実数の差を指示関数の積分で書く。実数 $x,x'$ について
$$ x-x'=\int_{-\infty}^\infty\bigl(\mathbf 1_{\{x'\le s\}}-\mathbf 1_{\{x\le s\}}\bigr)\,ds $$
である(差の積分表示)。実際、$x'< x$ なら被積分関数は $[x',x)$ で $1$、ほかで $0$ であり、$x'>x$ なら $[x,x')$ で $-1$、ほかで $0$ である。

共分散も右辺も $(X,Y)$ の同時分布だけで決まるので、$(X,Y)$ と同じ分布に従う組 $(X',Y')$ を、$(X,Y)$ と独立になるように積の確率空間の上にとって計算してよい。
段 1。 $E[(X-X')(Y-Y')]=2\,\mathrm{Cov}[X,Y]$ である。実際、展開すると $E[XY]+E[X'Y']-E[XY']-E[X'Y]$ であり、独立性から $E[XY']=E[X]E[Y']=E[X]E[Y]$(積の確率空間の上で Fubiniの定理 を使う)、同様に $E[X'Y]=E[X]E[Y]$、また $E[X'Y']=E[XY]$ である。
段 2。 差の積分表示を $X,X'$ と $Y,Y'$ に使い、積をとると
$$ (X-X')(Y-Y')=\int\!\!\int\bigl(\mathbf 1_{\{X'\le s\}}-\mathbf 1_{\{X\le s\}}\bigr)\bigl(\mathbf 1_{\{Y'\le t\}}-\mathbf 1_{\{Y\le t\}}\bigr)\,ds\,dt $$
である。被積分関数の絶対値の $(s,t)$ についての積分は $\lvert X-X'\rvert\lvert Y-Y'\rvert$ であり、$\lvert X-X'\rvert\lvert Y-Y'\rvert\le\frac12(X-X')^2+\frac12(Y-Y')^2$ の期待値は有限である。よって Fubiniの定理 により、期待値と $(s,t)$ についての積分を入れ替えてよい。
段 3。 各 $(s,t)$ で、被積分関数の期待値は
$$ P(X'\le s,Y'\le t)-P(X'\le s,Y\le t)-P(X\le s,Y'\le t)+P(X\le s,Y\le t)=2H(s,t)-2F(s)G(t) $$
である(2 項目と 3 項目は独立性から $F(s)G(t)$)。段 1・2 と合わせて $2\,\mathrm{Cov}[X,Y]=\int\!\!\int2(H-FG)\,ds\,dt$ を得る。可積分性は、$2\lvert H-FG\rvert$ が被積分関数の絶対値の期待値以下で、その積分が段 2 で有限だったことによる。$\square$

$X$ と $Y$ の関係について何も仮定していないことに注意する。共分散は「同時分布関数が、独立な場合の $F(s)G(t)$ からどちら向きにどれだけずれているか」を全平面で足し合わせたものである。

独立性・正の象限従属と相関

$E[X^2]<\infty$、$E[Y^2]<\infty$ とする。

  1. $X$ と $Y$ が独立なら $\mathrm{Cov}[X,Y]=0$ である。
  2. すべての $s,t$ で $H(s,t)\ge F(s)G(t)$ が成り立つ(正の象限従属)なら $\mathrm{Cov}[X,Y]\ge0$ であり、$\mathrm{Cov}[X,Y]=0$ となるのは $X$ と $Y$ が独立なときに限る。
  3. $g,h$ が単調非減少な関数で $E[g(X)^2]<\infty$、$E[h(X)^2]<\infty$ なら、$\mathrm{Cov}[g(X),h(X)]\ge0$ である。

1:独立なら $H(s,t)=P(X\le s)P(Y\le t)=F(s)G(t)$ なので、thm-cor-hoeffding の被積分関数は $0$ である。
2:$D:=H-FG\ge0$ なので thm-cor-hoeffding から $\mathrm{Cov}\ge0$ である。$\mathrm{Cov}=0$ とすると、$D\ge0$ の積分が $0$ なので、$D=0$ がほとんどいたるところで成り立つ。任意の $(s_0,t_0)$ をとる。正方形 $[s_0,s_0+\frac1k)\times[t_0,t_0+\frac1k)$ は測度が正なので、その中に $D(s_k,t_k)=0$ となる点がある。$H(s_0,t_0)\le H(s_k,t_k)\le H(s_0+\frac1k,t_0+\frac1k)$ であり、事象 $\{X\le s_0+\frac1k,\,Y\le t_0+\frac1k\}$ は $k$ について減少して共通部分が $\{X\le s_0,Y\le t_0\}$ なので、確率空間 の命題「確率の連続性」により $H(s_k,t_k)\to H(s_0,t_0)$ である。同様に $F(s_k)G(t_k)\to F(s_0)G(t_0)$ である。よって $D(s_0,t_0)=0$ がすべての点で成り立つ。2 つの確率変数は、すべての $s,t$ で $H(s,t)=F(s)G(t)$ なら独立である(Dur19 Theorem 2.1.8(p. 47)。この記事では証明しない)。
3:$A:=\{g(X)\le s\}$、$B:=\{h(X)\le t\}$ とおく。$g$ が非減少なので、$g(x)\le s$ となる $x$ の全体 $I$ は「$x\in I$、$x'< x$ なら $x'\in I$」を満たす区間(空集合、$\mathbb R$、$(-\infty,c)$、$(-\infty,c]$ のどれか)であり、$h$ についての $J$ も同様である。このような区間の 2 つは一方が他方に含まれるので、$A=\{X\in I\}$ と $B=\{X\in J\}$ も一方が他方に含まれ、
$$ P(A\cap B)=\min\bigl(P(A),P(B)\bigr)\ge P(A)P(B) $$
である。よって $(g(X),h(X))$ は正の象限従属で、2 から共分散は $0$ 以上である。$\square$

1 の逆は成り立たない。$X$ が $-1,0,1$ を確率 $\frac13$ ずつとり $Y=X^2$ なら無相関だが独立でない(分散 の定義「共分散」の後の例)。2 は、一方の向きにしかずれない組では、無相関と独立が一致することを述べている。3 は Chebyshev の和の不等式の確率版で、たとえば $E[X^6]< \infty$ なら $\mathrm{Cov}[X,X^3]\ge0$ である。

正規分布と順位による相関

$(X,Y)$ が 2 変量の正規分布に従うときは、無相関と独立が同値である(多変量正規分布 の定理「無相関なら独立」)。また、$X=x$ を与えたときの $Y$ の条件付き分布は正規分布で、その平均は thm-cor-prediction の 1 次式 $\alpha^*+\beta^*x$、分散は $V[Y](1-\rho^2)$ である(同記事の定理「条件付き分布」で 2 変量の場合をとったもの)。周辺の分布がそれぞれ正規分布なだけでは足りない(下の ex-cor-normal-margins)。
相関係数は 1 次の関係を測るので、単調だが 1 次でない関係では $1$ にならない。そこで、値を順位に置き換えてから相関係数をとることがある。

Spearman の順位相関と Kendall の順位相関

$X,Y$ の分布関数 $F,G$ が連続であるとする。

  1. $\rho_S(X,Y):=\rho\bigl(F(X),G(Y)\bigr)$ を Spearman の順位相関係数 という。
  2. $(X',Y')$ を $(X,Y)$ と同じ分布に従う独立な組として、$\tau_K(X,Y):=P\bigl((X-X')(Y-Y')>0\bigr)-P\bigl((X-X')(Y-Y')< 0\bigr)$ を Kendall の順位相関係数 という。

$F$ が連続なら $F(X)$ は $[0,1]$ 上の一様分布に従うので(分位点 の定理「分位点変換と確率積分変換」)、$V[F(X)]=\frac1{12}$ であり、$\rho_S=12\,E[F(X)G(Y)]-3$ となる。データでは、$F,G$ を経験分布関数に替えると、同じ値がなければ $F_n(x_i)$ は「$x_i$ の小さい方からの順位を $n$ で割った値」なので、prop-cor-affine により順位どうしの標本相関係数が得られる。これが 順位相関係数(高校数学) の $\rho$ である。$Y=g(X)$($g$ は狭義増加)なら $G(Y)=F(X)$ なので $\rho_S=1$、また $(X-X')(g(X)-g(X'))>0$ が $X\ne X'$ のとき成り立つので $\tau_K=1$ である。

標本相関係数の収束

標本相関係数の一致性

$(X_1,Y_1),(X_2,Y_2),\dots$ が互いに独立で、どれも $(X,Y)$ と同じ分布に従うとし、$E[X^2]<\infty$、$E[Y^2]<\infty$、$V[X]>0$、$V[Y]>0$ とする。最初の $n$ 組の標本相関係数を $r_n$ とすると(定義できない $n$ では $r_n:=0$ とおく)、確率 $1$ で $r_n\to\rho(X,Y)$ である。

$\bar X_n:=\frac1n\sum_{i\le n}X_i$ とし、同様に $\bar Y_n$、$\overline{XY}_n$、$\overline{X^2}_n$、$\overline{Y^2}_n$ を定める。$X_iY_i$ は可積分で、$X_i,Y_i,X_iY_i,X_i^2,Y_i^2$ のそれぞれが独立で同じ分布に従う可積分な確率変数の列なので、大数の強法則(Dur19 Theorem 2.4.1(p. 76)。この記事では証明しない)により、確率 $1$ の事象 $\Omega_0$ の上で 5 つの平均がそれぞれ $E[X]$、$E[Y]$、$E[XY]$、$E[X^2]$、$E[Y^2]$ に収束する。$r_n$ の分子と分母を $n$ で割ると
$$ r_n=\frac{\overline{XY}_n-\bar X_n\bar Y_n}{\sqrt{\overline{X^2}_n-\bar X_n^2}\sqrt{\overline{Y^2}_n-\bar Y_n^2}} $$
である。$\Omega_0$ の上では、分母の 2 つの根号の中身が $V[X]>0$、$V[Y]>0$ に収束するので、十分大きい $n$ で正になり(このとき $r_n$ は定義される)、分子は $\mathrm{Cov}[X,Y]$ に収束する。関数 $(p,q,u,v,w)\mapsto\dfrac{u-pq}{\sqrt{v-p^2}\sqrt{w-q^2}}$ は $v>p^2$、$w>q^2$ の範囲で連続なので、$\Omega_0$ の上で $r_n\to\rho(X,Y)$ である。$\square$

仮定 $E[X^2]<\infty$ は外せない。分散が有限でない分布では、相関係数そのものが定義されない(下の反例の表)。$r_n$ の揺らぎの大きさ($\sqrt n(r_n-\rho)$ の極限分布)は 4 次の積率で決まり、この記事では扱わない。

反例

外す条件反例成り立たなくなること
$V[X]>0$$X=1$(定数)、$Y$ は任意分母が $0$ で、$\rho(X,Y)$ は定義されない(def-cor-rho)
$E[X^2]<\infty$$X$ が標準 Cauchy分布、$Y=X$$E\lvert X\rvert=\infty$ で共分散が定まらず、$\rho$ も定義されない
同時に正規分布に従うex-cor-normal-margins周辺はどちらも正規分布で無相関だが、独立でない
1 次式の関係ex-cor-lognormal($Y=e^X$)単調な関数で決まるのに $\rho=0.762\ldots< 1$($\rho_S=1$)
正の象限従属$X$ が $-1,0,1$ を等確率、$Y=X^2$$H(-1,0)=0< \frac19=F(-1)G(0)$ で正の象限従属でなく、$\rho=0$ なのに独立でない(cor-cor-dependence の 2)
反例:周辺は正規分布だが独立でない無相関の組

$X$ を標準正規分布に従う確率変数、$S$ を $X$ と独立に $\pm1$ を確率 $\frac12$ ずつとる確率変数とし、$Y:=SX$ とおく。$X$ の分布は $0$ について対称なので、$Y$ も標準正規分布に従う。$E[XY]=E[S]E[X^2]=0$ なので無相関である。しかし $\lvert Y\rvert=\lvert X\rvert$ なので、$P(\lvert X\rvert\le1,\lvert Y\rvert\le1)=P(\lvert X\rvert\le1)=0.682\ldots$ は $P(\lvert X\rvert\le1)P(\lvert Y\rvert\le1)=0.466\ldots$ と異なり、独立でない。$(X,Y)$ は 2 変量の正規分布に従わない($X+Y$ は確率 $\frac12$ で $0$ になるので正規分布でない)。満たさない仮定:$(X,Y)$ が同時に正規分布に従うこと。

反例:単調だが 1 次でない関係

$X$ を標準正規分布に従う確率変数とし、$Y:=e^X$ とする。標準正規分布の密度を $\varphi$ とすると $e^x\varphi(x)=e^{1/2}\varphi(x-1)$(平方完成)なので、$E[e^X]=e^{1/2}$、$E[Xe^X]=e^{1/2}\int x\varphi(x-1)\,dx=e^{1/2}$ であり、同様に $E[e^{2X}]=e^2$ である。したがって $\mathrm{Cov}[X,Y]=E[Xe^X]-E[X]E[e^X]=e^{1/2}$、$V[Y]=E[e^{2X}]-E[e^X]^2=e^2-e$ であり
$$ \rho(X,Y)=\frac{e^{1/2}}{\sqrt{e^2-e}}=\frac1{\sqrt{e-1}}=0.762\ldots $$
である。$Y$ は $X$ の狭義増加関数なので $\rho_S(X,Y)=\tau_K(X,Y)=1$ だが、1 次式ではないので $\lvert\rho\rvert< 1$ である。$X$ を $e^X$ に替えるだけで相関係数が変わることは、prop-cor-affine が 1 次変換に限ることを示している。満たさない性質:$Y$ が $X$ の 1 次式であること。

関連項目

参考文献

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する