tecktom.com

この章のまとめと試験での出方

2026年10月17日7分で読めますAI生成教科書『資格試験』第14章 相関と回帰

AI 生成・木下が未勉強

統計学の復習 第14章 §14.7

この章の量は、ほとんどすべてが Sxx,Syy,SxyS_{xx},S_{yy},S_{xy} の3つから出た。相関係数 r=Sxy/SxxSyyr=S_{xy}/\sqrt{S_{xx}S_{yy}} は中心化したベクトルのなす角の余弦で、割り方に依存せず −1-1 と 11 の間にある。偏相関係数は第三の変数で説明できる分を除いた残りどうしの相関で、層別解析では偏差積和が層内と層間に分かれることから、全体と層の中で相関の向きが逆になる場合が説明できた。回帰直線は平方完成で求まり、傾き b1=Sxy/Sxxb_1=S_{xy}/S_{xx}、切片 b0=yˉ−b1xˉb_0=\bar y-b_1\bar x である。残差は 11 と xx に直交し、Syy=SR+SeS_{yy}=S_R+S_e、Se=Syy(1−r2)S_e=S_{yy}(1-r^2)、R2=r2R^2=r^2 が出る。誤差に正規分布を仮定すると、b1b_1 は重み付き和として N(β1,σ2/Sxx)\mathrm N(\beta_1,\sigma^2/S_{xx}) に従い、E[Se]=(n−2)σ2E[S_e]=(n-2)\sigma^2 から自由度 n−2n-2 が出て、傾きの tt 検定・無相関の検定・回帰の FF 検定の3つが同じ検定になる。

量式自由度・注意
相関係数Sxy/SxxSyyS_{xy}/\sqrt{S_{xx}S_{yy}}割り方に依存しない
傾き・切片Sxy/SxxS_{xy}/S_{xx}、yˉ−b1xˉ\bar y-b_1\bar x(xˉ,yˉ)(\bar x,\bar y) を通る
残差平方和Syy−Sxy2/Sxx=Syy(1−r2)S_{yy}-S_{xy}^2/S_{xx}=S_{yy}(1-r^2)n−2n-2
傾きの標準誤差Ve/Sxx\sqrt{V_e/S_{xx}}tn−2t_{n-2}
無相関の検定rn−2/1−r2r\sqrt{n-2}/\sqrt{1-r^2}tn−2t_{n-2}(傾きの t0t_0 と同じ)
回帰の FFSR/VeS_R/V_eF1,n−2F_{1,n-2}、F0=t02F_0=t_0^2
ρ\rho の区間z±zα/2/n−3z\pm z_{\alpha/2}/\sqrt{n-3}n−3n-3 はこれだけ
重回帰の Rˉ2\bar R^21−(1−R2)(n−1)/(n−p−1)1-(1-R^2)(n-1)/(n-p-1)残差は n−p−1n-p-1

2級では、5つの和か Sxx,Syy,SxyS_{xx},S_{yy},S_{xy} から相関係数・回帰係数・決定係数を出す問い、統計ソフトの出力(係数・標準誤差・tt 値・pp 値・R2R^2)を読む問い、擬相関・多重共線性・残差プロットの解釈を問う文章の正誤問題が出る。失点は、共分散の割り方の不揃い、自由度 n−2n-2 と n−3n-3 の取り違え、R2R^2 から rr を戻すときの符号に集中する。

1級の統計数理では、回帰は線形モデル Y=Xβ+ε\bm Y=X\bm\beta+\bm\varepsilon として出題される。典型的な小問の連鎖は「最小二乗推定量を導く → 推定量が YiY_i の線形結合であることから平均と分散を求める → 残差平方和の期待値から σ2\sigma^2 の不偏推定量を作る → tt 統計量の分布を示して検定や区間を作る」である。この章の §14.6 の手順はその p=1p=1 の場合で、行列で書いた一般の場合と Se/σ2∼χn−p−12S_e/\sigma^2\sim\chi^2_{n-p-1} の証明は数理統計の解説 第9章 §9.1・§9.2 で扱う。統計応用(理工学)では、回帰の分散分析表が実験計画法の分散分析表と同じ形で現れ、F=t2F=t^2 の関係が1自由度の要因の検定でそのまま使われる。

次の総合問題は、本章の回帰に、第7章の線形結合の分散、第9章の不偏性、第10章の tt 分布をつなぐ。

総合問題 回帰の傾きを、推定量の分布から検定と区間まで組み立てる ★★★

オリジナルある装置の設定値 xx と出力 yy について n=10n=10 組のデータを取り、xˉ=5\bar x=5、yˉ=12\bar y=12、Sxx=40S_{xx}=40、Syy=50S_{yy}=50、Sxy=36S_{xy}=36 を得た。モデルは Yi=β0+β1xi+εiY_i=\beta_0+\beta_1x_i+\varepsilon_i で、εi\varepsilon_i は独立に N(0,σ2)\mathrm N(0,\sigma^2) に従うとする。

  1. (1) 回帰直線を求めよ。(使う道具:本章 定理 14.6)
  2. (2) 残差平方和を2通りの形で求め、決定係数と相関係数を求めよ。(使う道具:本章 式 (14.16)・定義 14.9)
  3. (3) b1=∑wiYib_1=\sum w_iY_i と書けることを使い、E[b1]=β1E[b_1]=\beta_1、V[b1]=σ2/SxxV[b_1]=\sigma^2/S_{xx} を示せ。(使う道具:第7章 線形結合の分散、本章 式 (14.24))
  4. (4) H0:β1=0H_0:\beta_1=0 を有意水準 5%5\% の両側で検定し、β1\beta_1 の 95%95\% 信頼区間を求めよ(t8(0.025)=2.306t_8(0.025)=2.306)。t0t_0 を相関係数からも計算して確かめよ。(使う道具:第9章 不偏推定量、第10章 tt 分布、本章 定理 14.12)
  5. (5) x0=8x_0=8 での yy の平均 β0+8β1\beta_0+8\beta_1 の 95%95\% 信頼区間を求めよ。(使う道具:本章 式 (14.31))
解答

(1) b1=36/40=0.9b_1=36/40=0.9、b0=12−0.9×5=7.5b_0=12-0.9\times5=7.5 で、回帰直線は y^=7.5+0.9x\hat y=7.5+0.9x である。

(2) Se=Syy−Sxy2/Sxx=50−1296/40=50−32.4=17.6S_e=S_{yy}-S_{xy}^2/S_{xx}=50-1296/40=50-32.4=17.6、また Se=Syy−b1Sxy=50−0.9×36=17.6S_e=S_{yy}-b_1S_{xy}=50-0.9\times36=17.6 である。決定係数は R2=SR/Syy=32.4/50=0.648R^2=S_R/S_{yy}=32.4/50=0.648、相関係数は傾きが正なので r=+0.648=0.805r=+\sqrt{0.648}=0.805 である(36/40×50=36/44.72=0.80536/\sqrt{40\times50}=36/44.72=0.805 と一致)。

(3) wi=(xi−xˉ)/Sxxw_i=(x_i-\bar x)/S_{xx} は ∑wi=0\sum w_i=0、∑wixi=1\sum w_ix_i=1、∑wi2=1/Sxx\sum w_i^2=1/S_{xx} を満たす。期待値の線形性から

E[b1]=∑wi(β0+β1xi)=β0⋅0+β1⋅1=β1E[b_1]=\sum w_i(\beta_0+\beta_1x_i)=\beta_0\cdot0+\beta_1\cdot1=\beta_1

である。YiY_i は独立なので、第7章の線形結合の分散の公式で共分散の項がすべて 00 になり、

V[b1]=∑wi2V[Yi]=σ2∑wi2=σ2SxxV[b_1]=\sum w_i^2V[Y_i]=\sigma^2\sum w_i^2=\frac{\sigma^2}{S_{xx}}

である。

(4) 定理 14.11 から Ve=17.6/8=2.2V_e=17.6/8=2.2 が σ2\sigma^2 の不偏推定値で、SE(b1)=2.2/40=0.055=0.2345\mathrm{SE}(b_1)=\sqrt{2.2/40}=\sqrt{0.055}=0.2345 である。t0=0.9/0.2345=3.838t_0=0.9/0.2345=3.838 で、自由度 88 の tt 分布の両側 5%5\% 点 2.3062.306 を超えるので、H0H_0 を棄却する。相関係数からは t0=0.805×8/1−0.648=2.2768/0.5933=3.838t_0=0.805\times\sqrt8/\sqrt{1-0.648}=2.2768/0.5933=3.838 で一致する。β1\beta_1 の区間は 0.9±2.306×0.2345=0.9±0.5410.9\pm2.306\times0.2345=0.9\pm0.541 で [ 0.359, 1.441 ][\,0.359,\ 1.441\,] であり、00 を含まないことが棄却と対応している。

(5) 予測値は y^0=7.5+0.9×8=14.7\hat y_0=7.5+0.9\times8=14.7 である。式 (14.31) の σ2\sigma^2 を VeV_e で置き換えると、標準誤差は

2.2(110+(8−5)240)=2.2×0.325=0.715=0.8456\sqrt{2.2\left(\frac1{10}+\frac{(8-5)^2}{40}\right)}=\sqrt{2.2\times0.325}=\sqrt{0.715}=0.8456

で、半幅は 2.306×0.8456=1.9502.306\times0.8456=1.950 である。区間は [ 12.75, 16.65 ][\,12.75,\ 16.65\,] である。Y^0−(β0+8β1)\hat Y_0-(\beta_0+8\beta_1) は正規分布に従い、VeV_e と独立(VeV_e は Yˉ\bar Y とも b1b_1 とも独立)なので、標準誤差で割った量は t8t_8 に従う。区間の作り方は (4) と同じである。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『資格試験』の記事一覧へ