교재: Deisenroth·Faisal·Ong, Mathematics for Machine Learning (mml-book.github.io에 공개된 PDF, 2024-01-15 판). 쪽수와 식 번호는 그 PDF 기준이다.
왜 만들었나: 교재 Part II(8~12장)에는 장 끝 연습문제가 없다. 그래서 교재의 절·식·예제를 따라가며 손으로 풀 수 있는 문제를 만들었다. 문제마다 근거가 되는 교재 절과 식 번호를 적었다. 교재에 없는 내용을 묻는 곳은 교재 밖 보충이라고 표시했다.
답과 풀이: 다음 편에 있다. 모든 숫자 답은 numpy·scipy로 검산했다. 개념 해설은 장별 글(8장, 9장, 10장, 11장, 12장)에 있다.
표기: 교재와 같이 공분산의 분모는 NN, 스칼라 함수의 기울기는 행벡터, log⁡\log는 자연로그다. 9장의 데이터 행렬은 N×DN\times D, 10장은 D×ND\times N이다.

문제 종류: [계산] 손으로 숫자를 끝까지 구한다. [유도] 백지에 식을 끝까지 쓴다. [서술] 말로 설명한다.

장문제시간이 없을 때 먼저 풀 문제
8 모델과 데이터8-1 ~ 8-98-3(동전 MLE·MAP·베이즈), 8-5(d-분리), 8-7(베이즈 팩터)
9 선형 회귀9-1 ~ 9-99-1(정규방정식), 9-4(MAP), 9-5(베이즈 회귀)
10 PCA10-1 ~ 10-910-1(네 점 PCA), 10-3(라그랑주 유도), 10-5(N<DN<D)
11 GMM과 EM11-1 ~ 11-811-2(교재 예제 M 단계), 11-3(갱신식 유도), 11-6(EM 한 바퀴)
12 SVM12-1 ~ 12-812-2(하드 마진과 쌍대), 12-3(쌍대 유도), 12-5(XOR 커널)
OX 퀴즈15문항마지막 점검

8장 모델과 데이터

8-1. 경험적 위험과 과적합 [계산]

근거: 8.2.18.2.3절, 식 (8.6)(8.10)

훈련 데이터 (x,y)(x,y)는 (0,0), (1,2), (2,2)(0,0),\ (1,2),\ (2,2), 테스트 데이터는 (3, 2.5), (4, 4.5)(3,\,2.5),\ (4,\,4.5)이다. 예측기 두 개를 비교한다.

  • f1(x)=xf_1(x)=x (직선)
  • f2f_2: 훈련 데이터 세 점을 정확히 지나는 2차식

(a) f2f_2를 구하라. (b) 제곱 손실 ℓ(y,y^)=(y−y^)2\ell(y,\hat y)=(y-\hat y)^2로 두 예측기의 훈련 경험적 위험 Remp(f,Xtrain,ytrain)R_{\text{emp}}(f,X_{\text{train}},y_{\text{train}})을 구하라. (c) 테스트 데이터의 경험적 위험을 구하라. (d) 교재가 정의한 과적합에 비추어 어느 예측기가 과적합인지 설명하라.

8-2. 가우시안 가능도, 최소제곱, 정규화 최소제곱 [유도·계산]

근거: 예제 8.2, 8.3, 8.5, 식 (8.9), (8.12), (8.18)

(a) p(yn∣xn,θ)=N(yn∣xn⊤θ,σ2)p(y_n\mid x_n,\theta)=\mathcal N(y_n\mid x_n^\top\theta,\sigma^2)이고 σ2\sigma^2를 안다고 하자. 음의 로그 가능도 L(θ)\mathcal L(\theta)를 전개하고, L\mathcal L을 최소화하는 θ\theta가 최소제곱 문제 (8.9)의 해와 같음을 보여라. (b) 정규화 최소제곱 (8.12) min⁡θ1N∥y−Xθ∥2+λ∥θ∥2\min_\theta \frac1N\|y-X\theta\|^2+\lambda\|\theta\|^2의 해를 닫힌 꼴로 구하라. (c) 절편 없는 1차원 모델 f(x)=θxf(x)=\theta x와 데이터 x=(1,2,3)x=(1,2,3), y=(2,3,5)y=(2,3,5)에 대해 λ=0\lambda=0일 때와 λ=13\lambda=\frac13일 때 θ\theta를 구하라. (d) 9장의 정규화 최소제곱 (9.32) ∥y−Φθ∥2+λ∥θ∥2\|y-\Phi\theta\|^2+\lambda\|\theta\|^2에 있는 λ\lambda는 (8.12)의 λ\lambda와 어떤 관계인가?

8-3. 동전 던지기: MLE, MAP, 베이즈 추론 [계산]

근거: 식 (8.17), (8.19)(8.23), (8.32)(8.33), 그림 8.10(c), 6.6.1절(켤레 사전분포)

동전을 8번 던져 1,1,0,1,0,1,1,11,1,0,1,0,1,1,1을 얻었다(11=앞면). 모델은 p(x∣μ)=Ber(μ)p(x\mid\mu)=\text{Ber}(\mu)이고 던지기끼리는 독립이다.

(a) 음의 로그 가능도 L(μ)\mathcal L(\mu)를 쓰고 μML\mu_{\text{ML}}을 구하라. (b) 사전분포가 μ∼Beta(2,2)\mu\sim\text{Beta}(2,2)일 때 μMAP\mu_{\text{MAP}}을 구하라. (c) 사후분포 p(μ∣X)p(\mu\mid\mathcal X)를 구하라. (d) 다음 던지기가 앞면일 확률을 식 (8.23)처럼 μ\mu에 대해 적분해서 구하고, μML\mu_{\text{ML}}이나 μMAP\mu_{\text{MAP}}을 그대로 넣은 값과 비교하라.

8-4. 인수분해와 그래프 [계산]

근거: 예제 8.7, 8.8, 식 (8.29)~(8.31)

(a) p(a,b,c,d)=p(a) p(b∣a) p(c∣a) p(d∣b,c)p(a,b,c,d)=p(a)\,p(b\mid a)\,p(c\mid a)\,p(d\mid b,c)에 해당하는 방향 그래프의 화살표를 모두 써라. (b) 화살표가 x1→x2x_1\to x_2, x1→x3x_1\to x_3, x2→x4x_2\to x_4, x3→x4x_3\to x_4, x4→x5x_4\to x_5인 그래프의 결합분포를 식 (8.31)로 써라. (c) (b)의 다섯 변수가 모두 0 또는 1만 갖는다면, 결합분포를 표 하나로 저장할 때와 (b)의 인수분해로 저장할 때 필요한 자유 매개변수는 각각 몇 개인가? (교재 밖 보충)

8-5. d-분리 판정 [계산]

근거: 8.5.2절, 예제 8.9, 그림 8.11

그림 8.11의 그래프는 화살표 a→ba\to b, b→cb\to c, a→da\to d, c→dc\to d, d→ed\to e로 되어 있다.

(a) 교재의 네 문장 (8.35) b⊥ ⁣ ⁣ ⁣⊥d∣a,cb\perp\!\!\!\perp d\mid a,c, (8.36) a⊥ ⁣ ⁣ ⁣⊥c∣ba\perp\!\!\!\perp c\mid b, (8.37) b⊥̸ ⁣ ⁣ ⁣⊥d∣cb\not\perp\!\!\!\perp d\mid c, (8.38) a⊥̸ ⁣ ⁣ ⁣⊥c∣b,ea\not\perp\!\!\!\perp c\mid b,e가 맞음을 경로마다 따져 보여라. (b) 다음이 성립하는지 판정하라. (i) a⊥ ⁣ ⁣ ⁣⊥ca\perp\!\!\!\perp c (조건 없음) (ii) b⊥ ⁣ ⁣ ⁣⊥e∣db\perp\!\!\!\perp e\mid d (iii) a⊥ ⁣ ⁣ ⁣⊥e∣b,ca\perp\!\!\!\perp e\mid b,c (iv) b⊥ ⁣ ⁣ ⁣⊥e∣a,cb\perp\!\!\!\perp e\mid a,c

8-6. 잠재변수 모델: 주변화와 사후분포 [계산]

근거: 8.4.3절, 식 (8.24)~(8.28)

잠재변수 z∈{1,2}z\in\{1,2\}의 사전분포가 p(z=1)=0.3p(z=1)=0.3, p(z=2)=0.7p(z=2)=0.7이고, x∣z=1∼N(0,1)x\mid z=1\sim\mathcal N(0,1), x∣z=2∼N(3,1)x\mid z=2\sim\mathcal N(3,1)이다. (N(1∣0,1)≈0.2420\mathcal N(1\mid0,1)\approx0.2420, N(1∣3,1)≈0.0540\mathcal N(1\mid3,1)\approx0.0540)

(a) x=1x=1에서 주변 밀도 p(x=1)p(x=1)을 구하라. 식 (8.25)의 적분이 여기서는 합이 된다. (b) p(z=1∣x=1)p(z=1\mid x=1)을 구하라(식 (8.28)). (c) 사전확률 0.3에 비해 사후확률이 커졌는가, 작아졌는가? 이유는?

8-7. 베이즈 팩터로 모델 비교 [계산]

근거: 8.6.28.6.3절, 식 (8.43)(8.47)

8-3의 동전 데이터(앞면 6, 뒷면 2, 순서까지 주어진 수열)를 두 모델로 설명한다.

  • M1M_1: μ=0.5\mu=0.5로 고정(공정한 동전, 자유 매개변수 없음)
  • M2M_2: μ∼Beta(1,1)\mu\sim\text{Beta}(1,1), 즉 [0,1][0,1] 위의 균등분포

힌트: ∫01μa−1(1−μ)b−1dμ=(a−1)! (b−1)!(a+b−1)!\int_0^1\mu^{a-1}(1-\mu)^{b-1}d\mu=\frac{(a-1)!\,(b-1)!}{(a+b-1)!} (양의 정수 a,ba,b)

(a) p(D∣M1)p(\mathcal D\mid M_1)을 구하라. (b) p(D∣M2)=∫01μ6(1−μ)2dμp(\mathcal D\mid M_2)=\int_0^1\mu^6(1-\mu)^2d\mu를 구하라. (c) 베이즈 팩터 p(D∣M2)p(D∣M1)\frac{p(\mathcal D\mid M_2)}{p(\mathcal D\mid M_1)}를 구하고, 모델 사전확률이 같을 때 어느 모델을 고르는지 말하라. (d) M2M_2의 최대 가능도 max⁡μp(D∣μ)\max_\mu p(\mathcal D\mid\mu)를 p(D∣M1)p(\mathcal D\mid M_1)과 비교하라. (c)와 결론의 세기가 다른 이유를 설명하라. (e) 8번 모두 앞면이었다면 베이즈 팩터는 얼마인가? 모델 사전확률이 p(M1)=0.8p(M_1)=0.8, p(M2)=0.2p(M_2)=0.2라면 사후 오즈 p(M2∣D)p(M1∣D)\frac{p(M_2\mid\mathcal D)}{p(M_1\mid\mathcal D)}는?

8-8. AIC와 BIC [계산]

근거: 8.6.4절, 식 (8.48), (8.49)

데이터 N=100N=100개에 세 모델을 최대가능도로 맞췄다.

모델매개변수 수 MM최대 로그 가능도 log⁡p(x∣θML)\log p(x\mid\theta_{\text{ML}})
AA2−120-120
BB4−116-116
CC10−111-111

(a) 교재 (8.48)의 AIC log⁡p(x∣θ)−M\log p(x\mid\theta)-M을 계산하고 모델을 골라라. (b) 교재 (8.49)의 BIC log⁡p(x∣θ)−12Mlog⁡N\log p(x\mid\theta)-\frac12M\log N을 계산하고 모델을 골라라(log⁡100≈4.605\log100\approx4.605). (c) 두 기준이 다른 모델을 고른 이유를 설명하라. (d) 다른 책에서는 AIC를 2M−2log⁡p2M-2\log p(작을수록 좋음)로 쓴다. 교재의 정의와 어떤 관계인가?

8-9. 교차 검증과 중첩 교차 검증 [계산·서술]

근거: 8.2.4절, 8.6.1절, 식 (8.39), 8.6.1절 여백 주석(표준오차)

릿지 정규화 세기 λ\lambda 후보 세 개를 5겹 교차 검증(K=5K=5)으로 비교했다. 표는 검증 조각마다의 RMSE다.

λ\lambda조각 1조각 2조각 3조각 4조각 5
0.012.93.43.13.63.0
0.12.62.92.73.02.8
13.03.12.93.23.3

(a) 각 λ\lambda의 교차 검증 추정값 (8.39)와 표준오차 σ/K\sigma/\sqrt K를 구하라. σ\sigma는 다섯 값의 표준편차이고 분모는 KK로 한다. (b) 어떤 λ\lambda를 고르는가? (c) 고른 λ\lambda의 교차 검증 값을 “최종 일반화 성능"으로 보고하면 안 되는 이유를 말하고, 중첩 교차 검증에서 안쪽 고리와 바깥쪽 고리가 각각 무엇을 하는지 설명하라.


9장 선형 회귀

9-1. 직선 맞추기: 정규방정식 [계산]

근거: 9.2.1절, 식 (9.12c), (9.19), (9.22), (9.23)

데이터 (x,y)(x,y): (−1,0), (0,1), (1,1), (2,3)(-1,0),\ (0,1),\ (1,1),\ (2,3). 모델 y=θ0+θ1x+ϵy=\theta_0+\theta_1x+\epsilon, ϵ∼N(0,σ2)\epsilon\sim\mathcal N(0,\sigma^2), 특징 ϕ(x)=[1,x]⊤\phi(x)=[1,x]^\top.

(a) Φ\Phi, Φ⊤Φ\Phi^\top\Phi, Φ⊤y\Phi^\top y를 구하라. (b) θML=(Φ⊤Φ)−1Φ⊤y\theta_{\text{ML}}=(\Phi^\top\Phi)^{-1}\Phi^\top y를 구하라. (c) 잔차 y−ΦθMLy-\Phi\theta_{\text{ML}}이 Φ\Phi의 두 열과 직교함을 확인하라. (d) σML2\sigma^2_{\text{ML}} (9.22)와 RMSE (9.23)을 구하라.

9-2. 노이즈 분산의 최대가능도 [유도]

근거: 식 (9.20)~(9.22)

(a) θ\theta를 고정하고 log⁡p(Y∣X,θ,σ2)\log p(\mathcal Y\mid\mathcal X,\theta,\sigma^2)를 σ2\sigma^2로 미분해 σML2=1N∑n(yn−ϕ⊤(xn)θ)2\sigma^2_{\text{ML}}=\frac1N\sum_n(y_n-\phi^\top(x_n)\theta)^2을 유도하라. (b) 이 값이 최소가 아니라 최대임을 확인하라. (c) 최대 로그 가능도가 −N2log⁡(2πσML2)−N2-\frac N2\log(2\pi\sigma^2_{\text{ML}})-\frac N2임을 보이고, 9-1의 숫자로 값을 구하라.

9-3. 다항 특징과 과적합 [계산]

근거: 예제 9.39.5, 식 (9.16)(9.17), 9.2.2절, 그림 9.6

데이터 (x,y)(x,y): (−1,2), (0,0), (1,1)(-1,2),\ (0,0),\ (1,1).

(a) 2차 다항 특징 ϕ(x)=[1,x,x2]⊤\phi(x)=[1,x,x^2]^\top의 특징 행렬 Φ\Phi를 식 (9.17)처럼 쓰고 θML\theta_{\text{ML}}을 구하라. (b) 1차(직선) 특징 ϕ(x)=[1,x]⊤\phi(x)=[1,x]^\top로 θML\theta_{\text{ML}}과 훈련 RMSE를 구하라. (c) 2차 모델의 훈련 RMSE를 구하라. (d) 새 점 (2, 1.5)(2,\,1.5)에서 두 모델의 예측과 오차를 비교하고, 그림 9.6(차수 MM에 따른 훈련·테스트 RMSE)과 연결해 설명하라.

9-4. MAP 추정과 릿지 회귀 [유도·계산]

근거: 9.2.39.2.4절, 식 (9.28)(9.34)

(a) 가능도 N(y∣Φθ,σ2I)\mathcal N(y\mid\Phi\theta,\sigma^2I), 사전분포 p(θ)=N(0,b2I)p(\theta)=\mathcal N(0,b^2I)일 때 음의 로그 사후분포 (9.28)을 쓰고, 기울기를 0으로 놓아 θMAP\theta_{\text{MAP}} (9.31)을 유도하라. (b) 9-1의 데이터에서 σ2/b2=2\sigma^2/b^2=2일 때 θMAP\theta_{\text{MAP}}을 구하라. (c) ∥θML∥2\|\theta_{\text{ML}}\|^2와 ∥θMAP∥2\|\theta_{\text{MAP}}\|^2를 비교하라. (d) 정규화 최소제곱 (9.32)의 λ\lambda는 σ2\sigma^2, b2b^2와 어떤 관계인가?

9-5. 베이즈 선형 회귀: 사후분포와 예측 [계산]

근거: 9.3절, 식 (9.38), (9.43), (9.57)

매개변수가 하나인 모델 y=θx+ϵy=\theta x+\epsilon, ϵ∼N(0,1)\epsilon\sim\mathcal N(0,1), 사전분포 θ∼N(0,1)\theta\sim\mathcal N(0,1). 데이터 (x,y)(x,y): (1,2), (2,3)(1,2),\ (2,3).

(a) 사후분포 p(θ∣X,Y)=N(mN,SN)p(\theta\mid\mathcal X,\mathcal Y)=\mathcal N(m_N,S_N)을 구하라. (b) θML\theta_{\text{ML}}과 mNm_N을 비교하라. (c) x∗=0x_*=0과 x∗=3x_*=3에서 사후 예측 분포의 평균과 분산을 구하라. (d) x∗=3x_*=3에서 사전 예측 분포 (9.38)과 비교하라. (e) 세 번째 점 (−1,−1)(-1,-1)이 들어왔다. (a)의 사후분포를 새 사전분포로 써서 갱신한 결과와, 세 점을 한꺼번에 써서 구한 사후분포가 같음을 확인하라.

9-6. 주변 가능도로 모델 비교 [계산]

근거: 9.3.5절, 식 (9.60)(9.64), 8.6.28.6.3절

9-5와 같은 데이터 (1,2),(2,3)(1,2),(2,3), σ2=1\sigma^2=1. 상수: log⁡2π≈1.8379\log2\pi\approx1.8379, log⁡6≈1.7918\log6\approx1.7918, log⁡501≈6.2166\log501\approx6.2166.

(a) 사전분포 θ∼N(0,1)\theta\sim\mathcal N(0,1)일 때 p(y∣X)=N(y∣0, s0XX⊤+σ2I)p(y\mid\mathcal X)=\mathcal N(y\mid0,\ s_0XX^\top+\sigma^2I)를 쓰고 로그 값을 구하라. (b) 매개변수가 없는 모델 M0M_0: y=ϵy=\epsilon(즉 θ=0\theta=0 고정)의 log⁡p(y)\log p(y)를 구하라. (c) 사전분포 θ∼N(0,100)\theta\sim\mathcal N(0,100)일 때 log⁡p(y∣X)\log p(y\mid\mathcal X)를 구하라. (d) 세 값을 비교하고 “자동 오컴의 면도날”(8.6.2절), 제프리스-린들리 역설(8.6.3절 Remark)과 연결해 설명하라.

9-7. 최소제곱은 직교사영이다 [계산]

근거: 9.4절, 식 (9.65)~(9.67), 3.8절

X=[1,2,2]⊤X=[1,2,2]^\top, y=[1,1,4]⊤y=[1,1,4]^\top, 모델 y=Xθ+ϵy=X\theta+\epsilon (θ\theta는 스칼라).

(a) θML\theta_{\text{ML}}을 구하라. (b) 사영 XθMLX\theta_{\text{ML}}과 사영 행렬 P=X(X⊤X)−1X⊤P=X(X^\top X)^{-1}X^\top을 구하라. (c) 잔차가 XX와 직교함을 확인하라. (d) P2=PP^2=P, P⊤=PP^\top=P, tr(P)=1\text{tr}(P)=1임을 확인하고 각각의 뜻을 말하라.

9-8. 데이터가 많아지면 베이즈 사후 평균은 MLE로 간다 [유도]

근거: 식 (9.43)

매개변수가 하나인 모델 yn=θxn+ϵny_n=\theta x_n+\epsilon_n, ϵn∼N(0,σ2)\epsilon_n\sim\mathcal N(0,\sigma^2), 사전분포 θ∼N(m0,s0)\theta\sim\mathcal N(m_0,s_0).

(a) 식 (9.43)을 스칼라로 써서 SNS_N, mNm_N을 구하라. (b) mNm_N을 θML\theta_{\text{ML}}과 m0m_0의 가중평균으로 써라. (c) 데이터가 늘어 ∑nxn2→∞\sum_nx_n^2\to\infty이면 SN→0S_N\to0, mN→θMLm_N\to\theta_{\text{ML}}임을 보여라. (d) 반대로 s0→0s_0\to0이면 어떻게 되는가?

9-9. Φ⊤Φ\Phi^\top\Phi의 역행렬이 없을 때 [계산·서술]

근거: 식 (9.12) 옆 여백 주석(rk(X)=D\text{rk}(X)=D 조건), 식 (9.34)

특징 ϕ(x)=[x, 2x]⊤\phi(x)=[x,\ 2x]^\top(둘째 특징이 첫째의 두 배), 데이터 x=(1,2,3)x=(1,2,3), y=(1,2,2)y=(1,2,2).

(a) Φ⊤Φ\Phi^\top\Phi를 구하고 역행렬이 없음을 보여라. (b) 정규방정식의 해가 무수히 많지만 모든 해의 예측값이 같음을 보여라. (c) λ=1\lambda=1일 때 릿지 해 (9.34)를 구하라. (d) λ>0\lambda>0이면 Φ⊤Φ+λI\Phi^\top\Phi+\lambda I가 항상 가역인 이유를 말하라.


10장 PCA

10-1. 네 점 PCA 처음부터 끝까지 [계산]

근거: 10.2절, 식 (10.1), (10.13), (10.15), (10.16), (10.44)

데이터 (3,3), (2,4), (−1,1), (0,0)(3,3),\ (2,4),\ (-1,1),\ (0,0), 주성분 하나(M=1M=1)로 줄인다.

(a) 평균을 구하고 데이터를 중심화하라. (b) 공분산 SS를 구하라(분모 NN). (c) SS의 고유값과 단위 고유벡터를 구하라. (d) 좌표 zn=b1⊤(xn−μ)z_n=b_1^\top(x_n-\mu)를 구하고, znz_n의 분산이 λ1\lambda_1과 같음을 확인하라. (e) 복원 x~n=μ+b1zn\tilde x_n=\mu+b_1z_n을 구하고, 평균 제곱 재구성 오차가 λ2\lambda_2와 같음을 확인하라.

10-2. 3차원에서 주성분 개수 고르기 [계산]

근거: 식 (10.24)~(10.25), (10.34), (10.44)

중심화된 데이터의 공분산이 S=[310130001]S=\begin{bmatrix}3&1&0\\1&3&0\\0&0&1\end{bmatrix}이다.

(a) 고유값과 단위 고유벡터를 구하라. (b) M=1M=1, M=2M=2일 때 재구성 오차 JMJ_M과 설명된 분산 비율 VM/VDV_M/V_D를 구하라. (c) 중심화된 점 x=(2,0,1)x=(2,0,1)을 M=1M=1, M=2M=2 주부분공간에 사영하고 각각의 제곱 오차를 구하라.

10-3. 최대 분산 방향은 가장 큰 고유벡터다 [유도]

근거: 식 (10.9)~(10.23)

중심화된 데이터 x1,…,xNx_1,\dots,x_N, S=1N∑nxnxn⊤S=\frac1N\sum_nx_nx_n^\top.

(a) 단위벡터 b1b_1 방향 좌표 z1n=b1⊤xnz_{1n}=b_1^\top x_n의 분산이 V1=b1⊤Sb1V_1=b_1^\top Sb_1임을 보여라. (b) 제약 ∥b1∥=1\|b_1\|=1이 왜 필요한가? (c) 라그랑주 함수로 Sb1=λ1b1Sb_1=\lambda_1b_1을 유도하라. (d) 왜 가장 큰 고유값의 고유벡터를 고르는가? (e) 첫 주성분을 빼낸 데이터 X^=X−b1b1⊤X\hat X=X-b_1b_1^\top X의 공분산 S^\hat S에 대해 S^b1=0\hat Sb_1=0, S^bi=λibi\hat Sb_i=\lambda_ib_i (i≥2i\ge2)임을 보이고, 둘째 주성분이 b2b_2임을 설명하라.

10-4. SVD로 PCA 하기 [계산]

근거: 10.4절, 식 (10.45)~(10.49)

10-1의 중심화된 데이터를 열로 쌓은 D×ND\times N 행렬 X=[21−2−112−1−2]X=\begin{bmatrix}2&1&-2&-1\\1&2&-1&-2\end{bmatrix}.

(a) XX⊤XX^\top을 구하라. (b) XX의 특이값을 구하라. (c) λd=σd2/N\lambda_d=\sigma_d^2/N (10.49)를 확인하라. (d) 9장과 10장의 데이터 행렬 모양이 다른 점을 말하라.

10-5. N<DN<D일 때: 작은 행렬의 고유벡터로 PCA 하기 [계산]

근거: 10.5절, 식 (10.53)~(10.57)

D=4D=4, N=3N=3인 데이터 x1=(1,0,1,0)x_1=(1,0,1,0), x2=(−1,1,0,1)x_2=(-1,1,0,1), x3=(0,−1,−1,−1)x_3=(0,-1,-1,-1).

(a) 평균이 00임을 확인하라. (b) X=[x1,x2,x3]∈R4×3X=[x_1,x_2,x_3]\in\mathbb R^{4\times3}에 대해 X⊤XX^\top X (3×33\times3)를 구하라. (c) 1NX⊤X\frac1NX^\top X의 고유값과 고유벡터 cc를 구하라. (힌트: [1,1,1]⊤[1,1,1]^\top은 고유값 00의 고유벡터다.) (d) XcXc를 정규화해 S=1NXX⊤S=\frac1NXX^\top의 첫 주성분 b1b_1을 구하고 Sb1=λ1b1Sb_1=\lambda_1b_1을 확인하라. (e) 4×44\times4 대신 3×33\times3을 푸는 것이 어떤 이득인지 말하라.

10-6. 표준화 [계산·서술]

근거: 10.6절, 식 (10.58)~(10.61)

두 특징의 평균이 (170, 0)(170,\,0), 표준편차가 (10, 1)(10,\,1), 공분산이 S=[100661]S=\begin{bmatrix}100&6\\6&1\end{bmatrix}이다.

(a) 표준화하지 않은 SS의 고유값을 구하고, 첫 주성분이 거의 (1,0)(1,0)임을 보여라. (b) 표준화한 뒤의 공분산을 구하고 그 고유값과 고유벡터를 구하라. (c) 새 점 x∗=(185, 0.5)x_*=(185,\,0.5)를 (10.58)~(10.61)의 순서대로 표준화 → M=1M=1 사영 → 원래 단위로 되돌리기를 하라.

10-7. 거듭제곱법 [계산]

근거: 10.4절, 식 (10.52)

S=[2112]S=\begin{bmatrix}2&1\\1&2\end{bmatrix}, x0=(1,0)x_0=(1,0).

(a) x1,x2,x3x_1,x_2,x_3을 구하라(매번 정규화). (b) 레일리 몫 xk⊤Sxkx_k^\top Sx_k로 가장 큰 고유값을 추정하라. (c) 수렴 속도가 λ2/λ1\lambda_2/\lambda_1에 달린 이유를 설명하라. (교재 밖 보충)

10-8. 확률적 PCA [계산]

근거: 10.7절, 식 (10.63)~(10.75)

D=2D=2, M=1M=1, B=[2,1]⊤B=[2,1]^\top, μ=0\mu=0, σ2=1\sigma^2=1.

(a) p(x)p(x)의 공분산 BB⊤+σ2IBB^\top+\sigma^2I를 구하라(10.70b). (b) Cov[x,z]\text{Cov}[x,z]와 결합분포 p(x,z)p(x,z)를 써라(10.71)(10.72). (c) x=(3,1)x=(3,1)일 때 사후분포 p(z∣x)=N(m,C)p(z\mid x)=\mathcal N(m,C)를 구하라(10.73)(10.75). (d) σ2→0\sigma^2\to0일 때 사후 평균의 극한을 구하고, 이것이 PCA 좌표와 같아지는 조건을 말하라.

10-9. 재구성 오차는 버린 고유값의 합이다 [유도]

근거: 10.3절, 식 (10.29)~(10.44), (10.25)

(a) 최적 좌표가 zin=bi⊤xnz_{in}=b_i^\top x_n임을 유도하라(10.30)~(10.32). (b) 정규직교 기저 b1,…,bDb_1,\dots,b_D로 xn−x~n=∑j=M+1Dbjbj⊤xnx_n-\tilde x_n=\sum_{j=M+1}^Db_jb_j^\top x_n임을 보여라(10.38). (c) JM=∑j=M+1Dbj⊤SbjJ_M=\sum_{j=M+1}^Db_j^\top Sb_j를 보여라(10.43). (d) JM=∑j=M+1DλjJ_M=\sum_{j=M+1}^D\lambda_j (10.44)와 “전체 분산 = 남긴 분산 + 재구성 오차” (10.25)를 보여라.


11장 GMM과 EM

11-1. 교재 예제의 책임도 [계산]

근거: 예제 11.1, 11.2, 식 (11.6)~(11.8), (11.17), (11.19)

데이터 {−3,−2.5,−1,0,2,4,5}\{-3,-2.5,-1,0,2,4,5\}, 초기 성분 N(x∣−4,1)\mathcal N(x\mid-4,1), N(x∣0,0.2)\mathcal N(x\mid0,0.2), N(x∣8,3)\mathcal N(x\mid8,3), 비율 πk=13\pi_k=\frac13.

x=−1x=-1과 x=2x=2의 책임도 rn1,rn2,rn3r_{n1},r_{n2},r_{n3}을 계산하라.

상수: 12π≈0.3989\frac1{\sqrt{2\pi}}\approx0.3989, 10.4π≈0.8921\frac1{\sqrt{0.4\pi}}\approx0.8921, 16π≈0.2303\frac1{\sqrt{6\pi}}\approx0.2303, e−2.5≈0.08208e^{-2.5}\approx0.08208, e−4.5≈0.01111e^{-4.5}\approx0.01111, e−6≈0.002479e^{-6}\approx0.002479, e−10≈4.54×10−5e^{-10}\approx4.54\times10^{-5}, e−13.5≈1.37×10−6e^{-13.5}\approx1.37\times10^{-6}, e−18≈1.52×10−8e^{-18}\approx1.52\times10^{-8}

11-2. 교재 예제의 M 단계 한 번 [계산]

근거: 예제 11.311.5, 식 (11.20), (11.30), (11.42), (11.27)(11.29), (11.39)(11.41), (11.50)(11.52)

11-1의 데이터에서 교재 (11.19)의 책임도 행렬을 쓴다(행이 데이터, 열이 성분).

xnx_nrn1r_{n1}rn2r_{n2}rn3r_{n3}
−3-31.00.00.0
−2.5-2.51.00.00.0
−1-10.0570.9430.0
000.0010.9990.0
220.00.0660.934
440.00.01.0
550.00.01.0

(a) NkN_k를 구하라. (b) 새 평균 μk\mu_k를 구하라. (c) 새 평균을 써서 새 분산 σk2\sigma_k^2를 구하라. (d) 새 비율 πk\pi_k를 구하라. (e) 교재의 값(평균 −2.7,−0.4,3.7-2.7,-0.4,3.7, 분산 0.14,0.44,1.530.14,0.44,1.53, 비율 0.29,0.29,0.420.29,0.29,0.42)과 비교하라.

11-3. 갱신식 유도: 평균과 분산 [유도]

근거: 정리 11.1, 식 (11.10)~(11.24), 정리 11.2

(a) 로그 가능도 (11.10)을 μk\mu_k로 미분해 μknew=1Nk∑nrnkxn\mu_k^{\text{new}}=\frac1{N_k}\sum_nr_{nk}x_n (11.20)을 유도하라. (b) 이 식이 닫힌 해가 아닌 이유를 말하라. (c) 1차원(Σk=σk2\Sigma_k=\sigma_k^2)에서 σk2\sigma_k^2로 미분해 σk2,new=1Nk∑nrnk(xn−μk)2\sigma_k^{2,\text{new}}=\frac1{N_k}\sum_nr_{nk}(x_n-\mu_k)^2을 유도하라.

11-4. 혼합 비율 갱신: 라그랑주 승수 [유도]

근거: 정리 11.3, 식 (11.42)~(11.49)

(a) 제약 ∑kπk=1\sum_k\pi_k=1을 붙인 라그랑주 함수를 써라. (b) πk\pi_k로 미분해 πk=−Nk/λ\pi_k=-N_k/\lambda를 얻어라. (c) λ=−N\lambda=-N임을 보이고 πk=Nk/N\pi_k=N_k/N을 얻어라.

11-5. 잠재변수 관점과 결정 경계 [계산]

근거: 11.4절, 식 (11.58)~(11.72)

K=2K=2, π=(0.6, 0.4)\pi=(0.6,\,0.4), μ=(0, 2)\mu=(0,\,2), σ12=σ22=1\sigma_1^2=\sigma_2^2=1.

(a) 원-핫 잠재변수 zz의 p(z)p(z)와 p(x∣z)p(x\mid z)에서 p(x)p(x)를 유도하라(11.65)~(11.66). (b) x=1x=1에서 p(z1=1∣x)p(z_1=1\mid x)를 구하라. 결과가 π1\pi_1과 같은 이유는? (c) 두 책임도가 같아지는 xx를 구하라(log⁡1.5≈0.4055\log1.5\approx0.4055).

11-6. EM 한 바퀴: 잘 갈라진 데이터 [계산]

근거: 11.3절, 식 (11.53)~(11.56)

데이터 0,1,5,60,1,5,6, K=2K=2, 초기값 π=(0.5,0.5)\pi=(0.5,0.5), μ=(1,5)\mu=(1,5), σ12=σ22=1\sigma_1^2=\sigma_2^2=1. (e−8≈3.35×10−4e^{-8}\approx3.35\times10^{-4}, e−12≈6.1×10−6e^{-12}\approx6.1\times10^{-6})

(a) E 단계: 네 점의 책임도를 구하라. (b) M 단계: NkN_k, μk\mu_k, σk2\sigma_k^2, πk\pi_k를 구하라. (c) 로그 가능도가 −7.448-7.448에서 −5.676-5.676으로 늘었다. 왜 늘어야 하는가? (d) 책임도가 거의 0 또는 1일 때 이 계산은 K-평균의 어떤 단계와 같아지는가? (교재 밖 보충 포함)

11-7. 최대가능도의 특이점 [계산·서술]

근거: 11.5절(더 읽을거리)

데이터 0,1,5,60,1,5,6, K=2K=2, π=(0.5,0.5)\pi=(0.5,0.5), μ1=0\mu_1=0, σ12=ε\sigma_1^2=\varepsilon, μ2=3\mu_2=3, σ22=10\sigma_2^2=10.

(a) x=0x=0의 로그 밀도 항이 ε→0\varepsilon\to0에서 +∞+\infty로 감을 보여라. (b) 나머지 점의 항은 아래로 유계임을 보이고, 전체 로그 가능도가 +∞+\infty로 감을 결론지어라. (c) 이것이 왜 좋은 모델이 아닌지, 실제로 어떻게 막는지 말하라(교재 밖 보충).

11-8. EM과 GMM 개념 [서술]

근거: 11.2~11.5절

(a) 책임도 rnkr_{nk}가 사후확률 p(znk=1∣xn)p(z_{nk}=1\mid x_n)임을 베이즈 정리로 보여라(11.72). (b) E 단계와 M 단계는 각각 무엇을 고정하고 무엇을 바꾸는가? (c) EM이 전역 최댓값을 보장하지 않는 이유와 실전 대처를 말하라. (d) GMM 최대가능도와 하나의 가우시안 최대가능도의 차이를 “로그 안의 합"으로 설명하라.


12장 SVM

12-1. 초평면까지의 거리와 크기 조정 [계산]

근거: 12.112.2절, 식 (12.2)(12.18)

w=(3,4)w=(3,4), b=−5b=-5.

(a) 점 (3,4)(3,4), (0,0)(0,0), (1, 0.5)(1,\,0.5)의 f(x)=⟨w,x⟩+bf(x)=\langle w,x\rangle+b와 부호 있는 거리 f(x)/∥w∥f(x)/\|w\|를 구하라. (b) (w,b)(w,b)를 (6,8,−10)(6,8,-10)으로 두 배 해도 같은 초평면임을 보이고, yf(x)y f(x) 값은 두 배가 됨을 보여라. (c) 그래서 (12.15)처럼 “가장 가까운 점에서 yn(⟨w,xn⟩+b)=1y_n(\langle w,x_n\rangle+b)=1“로 크기를 고정하면 마진이 1/∥w∥1/\|w\|가 되는 이유를 (12.11)~(12.14)로 설명하라.

12-2. 하드 마진 SVM: 눈에 보이는 답이 최적이 아닐 때 [계산]

근거: 12.2.1~12.2.3절, 식 (12.18), 12.3.2절(볼록 껍질), 식 (12.38), (12.42)

양성 (2,2), (2,3)(2,2),\ (2,3), 음성 (0,0), (0,1)(0,0),\ (0,1).

(a) 수직선 x1=1x_1=1 (w=(1,0)w=(1,0), b=−1b=-1)이 모든 제약 (12.15)를 만족하고 마진이 1임을 확인하라. (b) 두 클래스 볼록 껍질 사이의 가장 가까운 두 점을 찾고, 그로부터 최대 마진 초평면 (w,b)(w,b)를 구하라. (c) (b)의 마진이 (a)보다 큼을 확인하라. (d) 쌍대 변수 α\alpha를 구해 (12.38)로 ww를 다시 만들고, (12.42)로 bb를 구하고, 쌍대 목적값과 원 목적값이 같음을 확인하라.

12-3. 소프트 마진 쌍대 유도 [유도]

근거: 12.3.1절, 식 (12.34)~(12.41)

(a) 원 문제 (12.26a)의 라그랑주 함수를 써라. (b) ww, bb, ξn\xi_n으로 편미분하라. (c) 대입해서 (12.40)을 얻고 쌍대 문제 (12.41)을 써라. (d) 0≤αn≤C0\le\alpha_n\le C가 나오는 이유를 말하라. (e) 하드 마진의 쌍대와 무엇이 다른가?

12-4. 힌지 손실과 여유 변수 [계산]

근거: 12.2.412.2.5절, 식 (12.26), (12.28)(12.33)

w=(1,−1)w=(1,-1), b=0b=0, C=1C=1.

점xnx_nyny_n
1(3,1)(3,1)+1+1
2(1,0)(1,0)+1+1
3(0, 0.5)(0,\,0.5)−1-1
4(2,1)(2,1)−1-1

(a) 각 점의 tn=ynf(xn)t_n=y_nf(x_n), 힌지 손실, 여유 변수 ξn\xi_n을 구하라. (b) 각 점이 마진 밖, 마진 경계 위, 마진 안(올바른 쪽), 잘못 분류 중 어디에 있는지 말하라. (c) 목적함수 (12.31)의 값을 구하라. (d) 0-1 손실과 비교하고, 힌지 손실이 0-1 손실의 상한임을 말하라.

12-5. 커널로 XOR 풀기 [계산]

근거: 12.4절, 식 (12.52)~(12.53), 12.3.2절

양성 (1,1), (−1,−1)(1,1),\ (-1,-1), 음성 (1,−1), (−1,1)(1,-1),\ (-1,1).

(a) 원래 2차원에서 직선으로 분리할 수 없음을 볼록 껍질로 보여라. (b) 다항 커널 k(x,x′)=(x⊤x′)2k(x,x')=(x^\top x')^2의 특징 사상이 ϕ(x)=(x12, 2x1x2, x22)\phi(x)=(x_1^2,\ \sqrt2x_1x_2,\ x_2^2)임을 보이고, x=(1,2)x=(1,2), x′=(3,−1)x'=(3,-1)로 숫자를 확인하라. (c) 네 점의 그람 행렬 KK를 구하고 양의 준정부호인지 확인하라(12.53). (d) 특징 공간에서 하드 마진 SVM의 ww, bb, 마진을 구하고 결정 함수를 원래 좌표로 써라. (e) 쌍대 변수 α\alpha를 구하고, 커널만으로 결정 함수를 다시 써서 (d)와 같음을 확인하라. (f) k(x,x′)=−∥x−x′∥2k(x,x')=-\|x-x'\|^2이 커널이 될 수 없음을 점 00과 11(1차원)로 보여라.

12-6. CC와 정규화 계수 λ\lambda의 관계 [유도·서술]

근거: 12.2.5절, 식 (12.31), 8.2.3절, 식 (8.6), (8.12)

(a) (12.31)을 CNCN으로 나누어 “평균 손실 + λ×\lambda\times정규화항” 꼴로 바꾸고 λ\lambda와 CC의 관계를 구하라. (b) CC가 클 때와 작을 때 마진과 훈련 오류가 어떻게 되는지 말하라. (c) C→∞C\to\infty이면 어떤 문제가 되는가?

12-7. 커널 SVM으로 예측하기 [계산]

근거: 식 (12.38), 12.4절

RBF 커널 k(x,x′)=exp⁡ ⁣(−12∥x−x′∥2)k(x,x')=\exp\!\big(-\frac12\|x-x'\|^2\big). 서포트 벡터 x1=(0,0)x_1=(0,0), y1=+1y_1=+1, α1=1\alpha_1=1과 x2=(1,1)x_2=(1,1), y2=−1y_2=-1, α2=1\alpha_2=1, b=0b=0이 주어졌다고 하자. 결정 함수는 f(x)=∑nαnynk(xn,x)+bf(x)=\sum_n\alpha_ny_nk(x_n,x)+b다.

(a) x∗=(0,1)x_*=(0,1), (0, 0.5)(0,\,0.5), (2,2)(2,2)에서 ff와 예측 레이블을 구하라. (e−0.5≈0.6065e^{-0.5}\approx0.6065, e−0.125≈0.8825e^{-0.125}\approx0.8825, e−0.625≈0.5353e^{-0.625}\approx0.5353, e−1≈0.3679e^{-1}\approx0.3679, e−4≈0.0183e^{-4}\approx0.0183) (b) 결정 경계를 구하라. (c) 왜 ww를 직접 계산하지 않는가?

12-8. SVM 개념 [서술]

근거: 12.2~12.4절

(a) 서포트 벡터란 무엇인가? αn=0\alpha_n=0인 점을 지워도 해가 그대로인 이유는? (b) 0<αn<C0<\alpha_n<C인 점과 αn=C\alpha_n=C인 점은 기하학적으로 어떻게 다른가? (c) 원 문제 대신 쌍대 문제를 푸는 이유 두 가지를 말하라. (d) 소프트 마진 원 문제 (12.26a)가 볼록 최적화 문제인 이유는?


OX 퀴즈

8~12장 공통. 맞으면 O, 틀리면 X. 틀린 것은 고쳐 말해 본다.

  1. MML의 PCA에서 데이터 공분산은 1N−1\frac1{N-1}로 나눈다.
  2. 교재 (8.48)의 AIC는 값이 작을수록 좋은 모델이다.
  3. MAP 추정은 매개변수의 사후분포 전체를 준다.
  4. 사전분포 N(0,b2I)\mathcal N(0,b^2I)의 MAP 추정은 λ=σ2/b2\lambda=\sigma^2/b^2인 정규화 최소제곱 (9.32)과 같다.
  5. 베이즈 선형 회귀의 사후 예측 분산은 노이즈 분산 σ2\sigma^2보다 작아질 수 있다.
  6. a→d←ca\to d\leftarrow c에서 dd를 관측하면 aa와 cc가 독립이 된다.
  7. PCA의 첫 주성분은 공분산의 가장 작은 고유값의 고유벡터다.
  8. PCA의 평균 제곱 재구성 오차는 버린 고유값의 합이다.
  9. GMM의 책임도 rnkr_{nk}를 kk에 대해 더하면 1이다.
  10. EM 알고리즘은 항상 로그 가능도의 전역 최댓값에 수렴한다.
  11. 가장 가까운 점에서 yn(⟨w,xn⟩+b)=1y_n(\langle w,x_n\rangle+b)=1로 크기를 고정하면 하드 마진 SVM의 마진은 1/∥w∥1/\|w\|다.
  12. 소프트 마진 SVM에서 αn=0\alpha_n=0인 점은 ww에 기여하지 않는다.
  13. 유효한 커널의 그람 행렬은 양의 준정부호다.
  14. 9장의 데이터 행렬은 D×ND\times N, 10장은 N×DN\times D다.
  15. 주변 가능도(증거)는 매개변수가 많은 모델을 항상 선호한다.