교재: 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장 )에 있다.표기: 교재와 같이 공분산의 분모는 N N N , 스칼라 함수의 기울기는 행벡터, log \log log 는 자연로그다. 9장의 데이터 행렬은 N × D N\times D N × D , 10장은 D × N D\times N D × N 이다.
문제 종류: [계산] 손으로 숫자를 끝까지 구한다. [유도] 백지에 식을 끝까지 쓴다. [서술] 말로 설명한다.
장 문제 시간이 없을 때 먼저 풀 문제 8 모델과 데이터 8-1 ~ 8-9 8-3(동전 MLE·MAP·베이즈), 8-5(d-분리), 8-7(베이즈 팩터) 9 선형 회귀 9-1 ~ 9-9 9-1(정규방정식), 9-4(MAP), 9-5(베이즈 회귀) 10 PCA 10-1 ~ 10-9 10-1(네 점 PCA), 10-3(라그랑주 유도), 10-5(N < D N<D N < D ) 11 GMM과 EM 11-1 ~ 11-8 11-2(교재 예제 M 단계), 11-3(갱신식 유도), 11-6(EM 한 바퀴) 12 SVM 12-1 ~ 12-8 12-2(하드 마진과 쌍대), 12-3(쌍대 유도), 12-5(XOR 커널) OX 퀴즈 15문항 마지막 점검
8장 모델과 데이터# 8-1. 경험적 위험과 과적합 [계산]# 근거: 8.2.18.2.3절, 식 (8.6)(8.10)
훈련 데이터 ( x , y ) (x,y) ( x , y ) 는 ( 0 , 0 ) , ( 1 , 2 ) , ( 2 , 2 ) (0,0),\ (1,2),\ (2,2) ( 0 , 0 ) , ( 1 , 2 ) , ( 2 , 2 ) , 테스트 데이터는 ( 3 , 2.5 ) , ( 4 , 4.5 ) (3,\,2.5),\ (4,\,4.5) ( 3 , 2.5 ) , ( 4 , 4.5 ) 이다. 예측기 두 개를 비교한다.
f 1 ( x ) = x f_1(x)=x f 1 ( x ) = x (직선)f 2 f_2 f 2 : 훈련 데이터 세 점을 정확히 지나는 2차식(a) f 2 f_2 f 2 를 구하라.
(b) 제곱 손실 ℓ ( y , y ^ ) = ( y − y ^ ) 2 \ell(y,\hat y)=(y-\hat y)^2 ℓ ( y , y ^ ) = ( y − y ^ ) 2 로 두 예측기의 훈련 경험적 위험 R emp ( f , X train , y train ) R_{\text{emp}}(f,X_{\text{train}},y_{\text{train}}) R emp ( f , X train , y train ) 을 구하라.
(c) 테스트 데이터의 경험적 위험을 구하라.
(d) 교재가 정의한 과적합에 비추어 어느 예측기가 과적합인지 설명하라.
8-2. 가우시안 가능도, 최소제곱, 정규화 최소제곱 [유도·계산]# 근거: 예제 8.2, 8.3, 8.5, 식 (8.9), (8.12), (8.18)
(a) p ( y n ∣ x n , θ ) = N ( y n ∣ x n ⊤ θ , σ 2 ) p(y_n\mid x_n,\theta)=\mathcal N(y_n\mid x_n^\top\theta,\sigma^2) p ( y n ∣ x n , θ ) = N ( y n ∣ x n ⊤ θ , σ 2 ) 이고 σ 2 \sigma^2 σ 2 를 안다고 하자. 음의 로그 가능도 L ( θ ) \mathcal L(\theta) L ( θ ) 를 전개하고, L \mathcal L L 을 최소화하는 θ \theta θ 가 최소제곱 문제 (8.9)의 해와 같음을 보여라.
(b) 정규화 최소제곱 (8.12) min θ 1 N ∥ y − X θ ∥ 2 + λ ∥ θ ∥ 2 \min_\theta \frac1N\|y-X\theta\|^2+\lambda\|\theta\|^2 min θ N 1 ∥ y − Xθ ∥ 2 + λ ∥ θ ∥ 2 의 해를 닫힌 꼴로 구하라.
(c) 절편 없는 1차원 모델 f ( x ) = θ x f(x)=\theta x f ( x ) = θ x 와 데이터 x = ( 1 , 2 , 3 ) x=(1,2,3) x = ( 1 , 2 , 3 ) , y = ( 2 , 3 , 5 ) y=(2,3,5) y = ( 2 , 3 , 5 ) 에 대해 λ = 0 \lambda=0 λ = 0 일 때와 λ = 1 3 \lambda=\frac13 λ = 3 1 일 때 θ \theta θ 를 구하라.
(d) 9장의 정규화 최소제곱 (9.32) ∥ y − Φ θ ∥ 2 + λ ∥ θ ∥ 2 \|y-\Phi\theta\|^2+\lambda\|\theta\|^2 ∥ y − Φ θ ∥ 2 + λ ∥ θ ∥ 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 , 1 1,1,0,1,0,1,1,1 1 , 1 , 0 , 1 , 0 , 1 , 1 , 1 을 얻었다(1 1 1 =앞면). 모델은 p ( x ∣ μ ) = Ber ( μ ) p(x\mid\mu)=\text{Ber}(\mu) p ( x ∣ μ ) = Ber ( μ ) 이고 던지기끼리는 독립이다.
(a) 음의 로그 가능도 L ( μ ) \mathcal L(\mu) L ( μ ) 를 쓰고 μ ML \mu_{\text{ML}} μ ML 을 구하라.
(b) 사전분포가 μ ∼ Beta ( 2 , 2 ) \mu\sim\text{Beta}(2,2) μ ∼ Beta ( 2 , 2 ) 일 때 μ MAP \mu_{\text{MAP}} μ MAP 을 구하라.
(c) 사후분포 p ( μ ∣ X ) p(\mu\mid\mathcal X) p ( μ ∣ X ) 를 구하라.
(d) 다음 던지기가 앞면일 확률을 식 (8.23)처럼 μ \mu μ 에 대해 적분해서 구하고, μ ML \mu_{\text{ML}} μ ML 이나 μ MAP \mu_{\text{MAP}} μ 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) p ( a , b , c , d ) = p ( a ) p ( b ∣ a ) p ( c ∣ a ) p ( d ∣ b , c ) 에 해당하는 방향 그래프의 화살표를 모두 써라.
(b) 화살표가 x 1 → x 2 x_1\to x_2 x 1 → x 2 , x 1 → x 3 x_1\to x_3 x 1 → x 3 , x 2 → x 4 x_2\to x_4 x 2 → x 4 , x 3 → x 4 x_3\to x_4 x 3 → x 4 , x 4 → x 5 x_4\to x_5 x 4 → x 5 인 그래프의 결합분포를 식 (8.31)로 써라.
(c) (b)의 다섯 변수가 모두 0 또는 1만 갖는다면, 결합분포를 표 하나로 저장할 때와 (b)의 인수분해로 저장할 때 필요한 자유 매개변수는 각각 몇 개인가? (교재 밖 보충 )
8-5. d-분리 판정 [계산]# 근거: 8.5.2절, 예제 8.9, 그림 8.11
그림 8.11의 그래프는 화살표 a → b a\to b a → b , b → c b\to c b → c , a → d a\to d a → d , c → d c\to d c → d , d → e d\to e d → e 로 되어 있다.
(a) 교재의 네 문장 (8.35) b ⊥ ⊥ d ∣ a , c b\perp\!\!\!\perp d\mid a,c b ⊥ ⊥ d ∣ a , c , (8.36) a ⊥ ⊥ c ∣ b a\perp\!\!\!\perp c\mid b a ⊥ ⊥ c ∣ b , (8.37) b ⊥̸ ⊥ d ∣ c b\not\perp\!\!\!\perp d\mid c b ⊥ ⊥ d ∣ c , (8.38) a ⊥̸ ⊥ c ∣ b , e a\not\perp\!\!\!\perp c\mid b,e a ⊥ ⊥ c ∣ b , e 가 맞음을 경로마다 따져 보여라.
(b) 다음이 성립하는지 판정하라.
(i) a ⊥ ⊥ c a\perp\!\!\!\perp c a ⊥ ⊥ c (조건 없음) (ii) b ⊥ ⊥ e ∣ d b\perp\!\!\!\perp e\mid d b ⊥ ⊥ e ∣ d (iii) a ⊥ ⊥ e ∣ b , c a\perp\!\!\!\perp e\mid b,c a ⊥ ⊥ e ∣ b , c (iv) b ⊥ ⊥ e ∣ a , c b\perp\!\!\!\perp e\mid a,c b ⊥ ⊥ e ∣ a , c
8-6. 잠재변수 모델: 주변화와 사후분포 [계산]# 근거: 8.4.3절, 식 (8.24)~(8.28)
잠재변수 z ∈ { 1 , 2 } z\in\{1,2\} z ∈ { 1 , 2 } 의 사전분포가 p ( z = 1 ) = 0.3 p(z=1)=0.3 p ( z = 1 ) = 0.3 , p ( z = 2 ) = 0.7 p(z=2)=0.7 p ( z = 2 ) = 0.7 이고, x ∣ z = 1 ∼ N ( 0 , 1 ) x\mid z=1\sim\mathcal N(0,1) x ∣ z = 1 ∼ N ( 0 , 1 ) , x ∣ z = 2 ∼ N ( 3 , 1 ) x\mid z=2\sim\mathcal N(3,1) x ∣ z = 2 ∼ N ( 3 , 1 ) 이다. (N ( 1 ∣ 0 , 1 ) ≈ 0.2420 \mathcal N(1\mid0,1)\approx0.2420 N ( 1 ∣ 0 , 1 ) ≈ 0.2420 , N ( 1 ∣ 3 , 1 ) ≈ 0.0540 \mathcal N(1\mid3,1)\approx0.0540 N ( 1 ∣ 3 , 1 ) ≈ 0.0540 )
(a) x = 1 x=1 x = 1 에서 주변 밀도 p ( x = 1 ) p(x=1) p ( x = 1 ) 을 구하라. 식 (8.25)의 적분이 여기서는 합이 된다.
(b) p ( z = 1 ∣ x = 1 ) p(z=1\mid x=1) p ( z = 1 ∣ x = 1 ) 을 구하라(식 (8.28)).
(c) 사전확률 0.3에 비해 사후확률이 커졌는가, 작아졌는가? 이유는?
8-7. 베이즈 팩터로 모델 비교 [계산]# 근거: 8.6.28.6.3절, 식 (8.43)(8.47)
8-3의 동전 데이터(앞면 6, 뒷면 2, 순서까지 주어진 수열)를 두 모델로 설명한다.
M 1 M_1 M 1 : μ = 0.5 \mu=0.5 μ = 0.5 로 고정(공정한 동전, 자유 매개변수 없음)M 2 M_2 M 2 : μ ∼ Beta ( 1 , 1 ) \mu\sim\text{Beta}(1,1) μ ∼ Beta ( 1 , 1 ) , 즉 [ 0 , 1 ] [0,1] [ 0 , 1 ] 위의 균등분포힌트: ∫ 0 1 μ a − 1 ( 1 − μ ) b − 1 d μ = ( 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)!} ∫ 0 1 μ a − 1 ( 1 − μ ) b − 1 d μ = ( a + b − 1 )! ( a − 1 )! ( b − 1 )! (양의 정수 a , b a,b a , b )
(a) p ( D ∣ M 1 ) p(\mathcal D\mid M_1) p ( D ∣ M 1 ) 을 구하라.
(b) p ( D ∣ M 2 ) = ∫ 0 1 μ 6 ( 1 − μ ) 2 d μ p(\mathcal D\mid M_2)=\int_0^1\mu^6(1-\mu)^2d\mu p ( D ∣ M 2 ) = ∫ 0 1 μ 6 ( 1 − μ ) 2 d μ 를 구하라.
(c) 베이즈 팩터 p ( D ∣ M 2 ) p ( D ∣ M 1 ) \frac{p(\mathcal D\mid M_2)}{p(\mathcal D\mid M_1)} p ( D ∣ M 1 ) p ( D ∣ M 2 ) 를 구하고, 모델 사전확률이 같을 때 어느 모델을 고르는지 말하라.
(d) M 2 M_2 M 2 의 최대 가능도 max μ p ( D ∣ μ ) \max_\mu p(\mathcal D\mid\mu) max μ p ( D ∣ μ ) 를 p ( D ∣ M 1 ) p(\mathcal D\mid M_1) p ( D ∣ M 1 ) 과 비교하라. (c)와 결론의 세기가 다른 이유를 설명하라.
(e) 8번 모두 앞면이었다면 베이즈 팩터는 얼마인가? 모델 사전확률이 p ( M 1 ) = 0.8 p(M_1)=0.8 p ( M 1 ) = 0.8 , p ( M 2 ) = 0.2 p(M_2)=0.2 p ( M 2 ) = 0.2 라면 사후 오즈 p ( M 2 ∣ D ) p ( M 1 ∣ D ) \frac{p(M_2\mid\mathcal D)}{p(M_1\mid\mathcal D)} p ( M 1 ∣ D ) p ( M 2 ∣ D ) 는?
8-8. AIC와 BIC [계산]# 근거: 8.6.4절, 식 (8.48), (8.49)
데이터 N = 100 N=100 N = 100 개에 세 모델을 최대가능도로 맞췄다.
모델 매개변수 수 M M M 최대 로그 가능도 log p ( x ∣ θ ML ) \log p(x\mid\theta_{\text{ML}}) log p ( x ∣ θ ML ) A A A 2 − 120 -120 − 120 B B B 4 − 116 -116 − 116 C C C 10 − 111 -111 − 111
(a) 교재 (8.48)의 AIC log p ( x ∣ θ ) − M \log p(x\mid\theta)-M log p ( x ∣ θ ) − M 을 계산하고 모델을 골라라.
(b) 교재 (8.49)의 BIC log p ( x ∣ θ ) − 1 2 M log N \log p(x\mid\theta)-\frac12M\log N log p ( x ∣ θ ) − 2 1 M log N 을 계산하고 모델을 골라라(log 100 ≈ 4.605 \log100\approx4.605 log 100 ≈ 4.605 ).
(c) 두 기준이 다른 모델을 고른 이유를 설명하라.
(d) 다른 책에서는 AIC를 2 M − 2 log p 2M-2\log p 2 M − 2 log p (작을수록 좋음)로 쓴다. 교재의 정의와 어떤 관계인가?
8-9. 교차 검증과 중첩 교차 검증 [계산·서술]# 근거: 8.2.4절, 8.6.1절, 식 (8.39), 8.6.1절 여백 주석(표준오차)
릿지 정규화 세기 λ \lambda λ 후보 세 개를 5겹 교차 검증(K = 5 K=5 K = 5 )으로 비교했다. 표는 검증 조각마다의 RMSE다.
λ \lambda λ 조각 1 조각 2 조각 3 조각 4 조각 5 0.01 2.9 3.4 3.1 3.6 3.0 0.1 2.6 2.9 2.7 3.0 2.8 1 3.0 3.1 2.9 3.2 3.3
(a) 각 λ \lambda λ 의 교차 검증 추정값 (8.39)와 표준오차 σ / K \sigma/\sqrt K σ / K 를 구하라. σ \sigma σ 는 다섯 값의 표준편차이고 분모는 K K K 로 한다.
(b) 어떤 λ \lambda λ 를 고르는가?
(c) 고른 λ \lambda λ 의 교차 검증 값을 “최종 일반화 성능"으로 보고하면 안 되는 이유를 말하고, 중첩 교차 검증에서 안쪽 고리와 바깥쪽 고리가 각각 무엇을 하는지 설명하라.
9장 선형 회귀# 9-1. 직선 맞추기: 정규방정식 [계산]# 근거: 9.2.1절, 식 (9.12c), (9.19), (9.22), (9.23)
데이터 ( x , y ) (x,y) ( x , y ) : ( − 1 , 0 ) , ( 0 , 1 ) , ( 1 , 1 ) , ( 2 , 3 ) (-1,0),\ (0,1),\ (1,1),\ (2,3) ( − 1 , 0 ) , ( 0 , 1 ) , ( 1 , 1 ) , ( 2 , 3 ) . 모델 y = θ 0 + θ 1 x + ϵ y=\theta_0+\theta_1x+\epsilon y = θ 0 + θ 1 x + ϵ , ϵ ∼ N ( 0 , σ 2 ) \epsilon\sim\mathcal N(0,\sigma^2) ϵ ∼ N ( 0 , σ 2 ) , 특징 ϕ ( x ) = [ 1 , x ] ⊤ \phi(x)=[1,x]^\top ϕ ( x ) = [ 1 , x ] ⊤ .
(a) Φ \Phi Φ , Φ ⊤ Φ \Phi^\top\Phi Φ ⊤ Φ , Φ ⊤ y \Phi^\top y Φ ⊤ y 를 구하라.
(b) θ ML = ( Φ ⊤ Φ ) − 1 Φ ⊤ y \theta_{\text{ML}}=(\Phi^\top\Phi)^{-1}\Phi^\top y θ ML = ( Φ ⊤ Φ ) − 1 Φ ⊤ y 를 구하라.
(c) 잔차 y − Φ θ ML y-\Phi\theta_{\text{ML}} y − Φ θ ML 이 Φ \Phi Φ 의 두 열과 직교함을 확인하라.
(d) σ ML 2 \sigma^2_{\text{ML}} σ ML 2 (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) log p ( Y ∣ X , θ , σ 2 ) 를 σ 2 \sigma^2 σ 2 로 미분해 σ ML 2 = 1 N ∑ n ( y n − ϕ ⊤ ( x n ) θ ) 2 \sigma^2_{\text{ML}}=\frac1N\sum_n(y_n-\phi^\top(x_n)\theta)^2 σ ML 2 = N 1 ∑ n ( y n − ϕ ⊤ ( x n ) θ ) 2 을 유도하라.
(b) 이 값이 최소가 아니라 최대임을 확인하라.
(c) 최대 로그 가능도가 − N 2 log ( 2 π σ ML 2 ) − N 2 -\frac N2\log(2\pi\sigma^2_{\text{ML}})-\frac N2 − 2 N log ( 2 π σ ML 2 ) − 2 N 임을 보이고, 9-1의 숫자로 값을 구하라.
9-3. 다항 특징과 과적합 [계산]# 근거: 예제 9.39.5, 식 (9.16)(9.17), 9.2.2절, 그림 9.6
데이터 ( x , y ) (x,y) ( x , y ) : ( − 1 , 2 ) , ( 0 , 0 ) , ( 1 , 1 ) (-1,2),\ (0,0),\ (1,1) ( − 1 , 2 ) , ( 0 , 0 ) , ( 1 , 1 ) .
(a) 2차 다항 특징 ϕ ( x ) = [ 1 , x , x 2 ] ⊤ \phi(x)=[1,x,x^2]^\top ϕ ( x ) = [ 1 , x , x 2 ] ⊤ 의 특징 행렬 Φ \Phi Φ 를 식 (9.17)처럼 쓰고 θ ML \theta_{\text{ML}} θ ML 을 구하라.
(b) 1차(직선) 특징 ϕ ( x ) = [ 1 , x ] ⊤ \phi(x)=[1,x]^\top ϕ ( x ) = [ 1 , x ] ⊤ 로 θ ML \theta_{\text{ML}} θ ML 과 훈련 RMSE를 구하라.
(c) 2차 모델의 훈련 RMSE를 구하라.
(d) 새 점 ( 2 , 1.5 ) (2,\,1.5) ( 2 , 1.5 ) 에서 두 모델의 예측과 오차를 비교하고, 그림 9.6(차수 M M M 에 따른 훈련·테스트 RMSE)과 연결해 설명하라.
9-4. MAP 추정과 릿지 회귀 [유도·계산]# 근거: 9.2.39.2.4절, 식 (9.28)(9.34)
(a) 가능도 N ( y ∣ Φ θ , σ 2 I ) \mathcal N(y\mid\Phi\theta,\sigma^2I) N ( y ∣ Φ θ , σ 2 I ) , 사전분포 p ( θ ) = N ( 0 , b 2 I ) p(\theta)=\mathcal N(0,b^2I) p ( θ ) = N ( 0 , b 2 I ) 일 때 음의 로그 사후분포 (9.28)을 쓰고, 기울기를 0으로 놓아 θ MAP \theta_{\text{MAP}} θ MAP (9.31)을 유도하라.
(b) 9-1의 데이터에서 σ 2 / b 2 = 2 \sigma^2/b^2=2 σ 2 / b 2 = 2 일 때 θ MAP \theta_{\text{MAP}} θ MAP 을 구하라.
(c) ∥ θ ML ∥ 2 \|\theta_{\text{ML}}\|^2 ∥ θ ML ∥ 2 와 ∥ θ MAP ∥ 2 \|\theta_{\text{MAP}}\|^2 ∥ θ MAP ∥ 2 를 비교하라.
(d) 정규화 최소제곱 (9.32)의 λ \lambda λ 는 σ 2 \sigma^2 σ 2 , b 2 b^2 b 2 와 어떤 관계인가?
9-5. 베이즈 선형 회귀: 사후분포와 예측 [계산]# 근거: 9.3절, 식 (9.38), (9.43), (9.57)
매개변수가 하나인 모델 y = θ x + ϵ y=\theta x+\epsilon y = θ x + ϵ , ϵ ∼ N ( 0 , 1 ) \epsilon\sim\mathcal N(0,1) ϵ ∼ N ( 0 , 1 ) , 사전분포 θ ∼ N ( 0 , 1 ) \theta\sim\mathcal N(0,1) θ ∼ N ( 0 , 1 ) . 데이터 ( x , y ) (x,y) ( x , y ) : ( 1 , 2 ) , ( 2 , 3 ) (1,2),\ (2,3) ( 1 , 2 ) , ( 2 , 3 ) .
(a) 사후분포 p ( θ ∣ X , Y ) = N ( m N , S N ) p(\theta\mid\mathcal X,\mathcal Y)=\mathcal N(m_N,S_N) p ( θ ∣ X , Y ) = N ( m N , S N ) 을 구하라.
(b) θ ML \theta_{\text{ML}} θ ML 과 m N m_N m N 을 비교하라.
(c) x ∗ = 0 x_*=0 x ∗ = 0 과 x ∗ = 3 x_*=3 x ∗ = 3 에서 사후 예측 분포의 평균과 분산을 구하라.
(d) x ∗ = 3 x_*=3 x ∗ = 3 에서 사전 예측 분포 (9.38)과 비교하라.
(e) 세 번째 점 ( − 1 , − 1 ) (-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) ( 1 , 2 ) , ( 2 , 3 ) , σ 2 = 1 \sigma^2=1 σ 2 = 1 . 상수: log 2 π ≈ 1.8379 \log2\pi\approx1.8379 log 2 π ≈ 1.8379 , log 6 ≈ 1.7918 \log6\approx1.7918 log 6 ≈ 1.7918 , log 501 ≈ 6.2166 \log501\approx6.2166 log 501 ≈ 6.2166 .
(a) 사전분포 θ ∼ N ( 0 , 1 ) \theta\sim\mathcal N(0,1) θ ∼ N ( 0 , 1 ) 일 때 p ( y ∣ X ) = N ( y ∣ 0 , s 0 X X ⊤ + σ 2 I ) p(y\mid\mathcal X)=\mathcal N(y\mid0,\ s_0XX^\top+\sigma^2I) p ( y ∣ X ) = N ( y ∣ 0 , s 0 X X ⊤ + σ 2 I ) 를 쓰고 로그 값을 구하라.
(b) 매개변수가 없는 모델 M 0 M_0 M 0 : y = ϵ y=\epsilon y = ϵ (즉 θ = 0 \theta=0 θ = 0 고정)의 log p ( y ) \log p(y) log p ( y ) 를 구하라.
(c) 사전분포 θ ∼ N ( 0 , 100 ) \theta\sim\mathcal N(0,100) θ ∼ N ( 0 , 100 ) 일 때 log p ( y ∣ X ) \log p(y\mid\mathcal X) log p ( y ∣ 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 X = [ 1 , 2 , 2 ] ⊤ , y = [ 1 , 1 , 4 ] ⊤ y=[1,1,4]^\top y = [ 1 , 1 , 4 ] ⊤ , 모델 y = X θ + ϵ y=X\theta+\epsilon y = Xθ + ϵ (θ \theta θ 는 스칼라).
(a) θ ML \theta_{\text{ML}} θ ML 을 구하라.
(b) 사영 X θ ML X\theta_{\text{ML}} X θ ML 과 사영 행렬 P = X ( X ⊤ X ) − 1 X ⊤ P=X(X^\top X)^{-1}X^\top P = X ( X ⊤ X ) − 1 X ⊤ 을 구하라.
(c) 잔차가 X X X 와 직교함을 확인하라.
(d) P 2 = P P^2=P P 2 = P , P ⊤ = P P^\top=P P ⊤ = P , tr ( P ) = 1 \text{tr}(P)=1 tr ( P ) = 1 임을 확인하고 각각의 뜻을 말하라.
9-8. 데이터가 많아지면 베이즈 사후 평균은 MLE로 간다 [유도]# 근거: 식 (9.43)
매개변수가 하나인 모델 y n = θ x n + ϵ n y_n=\theta x_n+\epsilon_n y n = θ x n + ϵ n , ϵ n ∼ N ( 0 , σ 2 ) \epsilon_n\sim\mathcal N(0,\sigma^2) ϵ n ∼ N ( 0 , σ 2 ) , 사전분포 θ ∼ N ( m 0 , s 0 ) \theta\sim\mathcal N(m_0,s_0) θ ∼ N ( m 0 , s 0 ) .
(a) 식 (9.43)을 스칼라로 써서 S N S_N S N , m N m_N m N 을 구하라.
(b) m N m_N m N 을 θ ML \theta_{\text{ML}} θ ML 과 m 0 m_0 m 0 의 가중평균으로 써라.
(c) 데이터가 늘어 ∑ n x n 2 → ∞ \sum_nx_n^2\to\infty ∑ n x n 2 → ∞ 이면 S N → 0 S_N\to0 S N → 0 , m N → θ ML m_N\to\theta_{\text{ML}} m N → θ ML 임을 보여라.
(d) 반대로 s 0 → 0 s_0\to0 s 0 → 0 이면 어떻게 되는가?
9-9. Φ ⊤ Φ \Phi^\top\Phi Φ ⊤ Φ 의 역행렬이 없을 때 [계산·서술]# 근거: 식 (9.12) 옆 여백 주석(rk ( X ) = D \text{rk}(X)=D rk ( X ) = D 조건), 식 (9.34)
특징 ϕ ( x ) = [ x , 2 x ] ⊤ \phi(x)=[x,\ 2x]^\top ϕ ( x ) = [ x , 2 x ] ⊤ (둘째 특징이 첫째의 두 배), 데이터 x = ( 1 , 2 , 3 ) x=(1,2,3) x = ( 1 , 2 , 3 ) , y = ( 1 , 2 , 2 ) y=(1,2,2) y = ( 1 , 2 , 2 ) .
(a) Φ ⊤ Φ \Phi^\top\Phi Φ ⊤ Φ 를 구하고 역행렬이 없음을 보여라.
(b) 정규방정식의 해가 무수히 많지만 모든 해의 예측값이 같음을 보여라.
(c) λ = 1 \lambda=1 λ = 1 일 때 릿지 해 (9.34)를 구하라.
(d) λ > 0 \lambda>0 λ > 0 이면 Φ ⊤ Φ + λ I \Phi^\top\Phi+\lambda I Φ ⊤ Φ + λ 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) ( 3 , 3 ) , ( 2 , 4 ) , ( − 1 , 1 ) , ( 0 , 0 ) , 주성분 하나(M = 1 M=1 M = 1 )로 줄인다.
(a) 평균을 구하고 데이터를 중심화하라.
(b) 공분산 S S S 를 구하라(분모 N N N ).
(c) S S S 의 고유값과 단위 고유벡터를 구하라.
(d) 좌표 z n = b 1 ⊤ ( x n − μ ) z_n=b_1^\top(x_n-\mu) z n = b 1 ⊤ ( x n − μ ) 를 구하고, z n z_n z n 의 분산이 λ 1 \lambda_1 λ 1 과 같음을 확인하라.
(e) 복원 x ~ n = μ + b 1 z n \tilde x_n=\mu+b_1z_n x ~ n = μ + b 1 z n 을 구하고, 평균 제곱 재구성 오차가 λ 2 \lambda_2 λ 2 와 같음을 확인하라.
10-2. 3차원에서 주성분 개수 고르기 [계산]# 근거: 식 (10.24)~(10.25), (10.34), (10.44)
중심화된 데이터의 공분산이 S = [ 3 1 0 1 3 0 0 0 1 ] S=\begin{bmatrix}3&1&0\\1&3&0\\0&0&1\end{bmatrix} S = 3 1 0 1 3 0 0 0 1 이다.
(a) 고유값과 단위 고유벡터를 구하라.
(b) M = 1 M=1 M = 1 , M = 2 M=2 M = 2 일 때 재구성 오차 J M J_M J M 과 설명된 분산 비율 V M / V D V_M/V_D V M / V D 를 구하라.
(c) 중심화된 점 x = ( 2 , 0 , 1 ) x=(2,0,1) x = ( 2 , 0 , 1 ) 을 M = 1 M=1 M = 1 , M = 2 M=2 M = 2 주부분공간에 사영하고 각각의 제곱 오차를 구하라.
10-3. 최대 분산 방향은 가장 큰 고유벡터다 [유도]# 근거: 식 (10.9)~(10.23)
중심화된 데이터 x 1 , … , x N x_1,\dots,x_N x 1 , … , x N , S = 1 N ∑ n x n x n ⊤ S=\frac1N\sum_nx_nx_n^\top S = N 1 ∑ n x n x n ⊤ .
(a) 단위벡터 b 1 b_1 b 1 방향 좌표 z 1 n = b 1 ⊤ x n z_{1n}=b_1^\top x_n z 1 n = b 1 ⊤ x n 의 분산이 V 1 = b 1 ⊤ S b 1 V_1=b_1^\top Sb_1 V 1 = b 1 ⊤ S b 1 임을 보여라.
(b) 제약 ∥ b 1 ∥ = 1 \|b_1\|=1 ∥ b 1 ∥ = 1 이 왜 필요한가?
(c) 라그랑주 함수로 S b 1 = λ 1 b 1 Sb_1=\lambda_1b_1 S b 1 = λ 1 b 1 을 유도하라.
(d) 왜 가장 큰 고유값의 고유벡터를 고르는가?
(e) 첫 주성분을 빼낸 데이터 X ^ = X − b 1 b 1 ⊤ X \hat X=X-b_1b_1^\top X X ^ = X − b 1 b 1 ⊤ X 의 공분산 S ^ \hat S S ^ 에 대해 S ^ b 1 = 0 \hat Sb_1=0 S ^ b 1 = 0 , S ^ b i = λ i b i \hat Sb_i=\lambda_ib_i S ^ b i = λ i b i (i ≥ 2 i\ge2 i ≥ 2 )임을 보이고, 둘째 주성분이 b 2 b_2 b 2 임을 설명하라.
10-4. SVD로 PCA 하기 [계산]# 근거: 10.4절, 식 (10.45)~(10.49)
10-1의 중심화된 데이터를 열로 쌓은 D × N D\times N D × N 행렬 X = [ 2 1 − 2 − 1 1 2 − 1 − 2 ] X=\begin{bmatrix}2&1&-2&-1\\1&2&-1&-2\end{bmatrix} X = [ 2 1 1 2 − 2 − 1 − 1 − 2 ] .
(a) X X ⊤ XX^\top X X ⊤ 을 구하라.
(b) X X X 의 특이값을 구하라.
(c) λ d = σ d 2 / N \lambda_d=\sigma_d^2/N λ d = σ d 2 / N (10.49)를 확인하라.
(d) 9장과 10장의 데이터 행렬 모양이 다른 점을 말하라.
10-5. N < D N<D N < D 일 때: 작은 행렬의 고유벡터로 PCA 하기 [계산]# 근거: 10.5절, 식 (10.53)~(10.57)
D = 4 D=4 D = 4 , N = 3 N=3 N = 3 인 데이터 x 1 = ( 1 , 0 , 1 , 0 ) x_1=(1,0,1,0) x 1 = ( 1 , 0 , 1 , 0 ) , x 2 = ( − 1 , 1 , 0 , 1 ) x_2=(-1,1,0,1) x 2 = ( − 1 , 1 , 0 , 1 ) , x 3 = ( 0 , − 1 , − 1 , − 1 ) x_3=(0,-1,-1,-1) x 3 = ( 0 , − 1 , − 1 , − 1 ) .
(a) 평균이 0 0 0 임을 확인하라.
(b) X = [ x 1 , x 2 , x 3 ] ∈ R 4 × 3 X=[x_1,x_2,x_3]\in\mathbb R^{4\times3} X = [ x 1 , x 2 , x 3 ] ∈ R 4 × 3 에 대해 X ⊤ X X^\top X X ⊤ X (3 × 3 3\times3 3 × 3 )를 구하라.
(c) 1 N X ⊤ X \frac1NX^\top X N 1 X ⊤ X 의 고유값과 고유벡터 c c c 를 구하라. (힌트: [ 1 , 1 , 1 ] ⊤ [1,1,1]^\top [ 1 , 1 , 1 ] ⊤ 은 고유값 0 0 0 의 고유벡터다.)
(d) X c Xc X c 를 정규화해 S = 1 N X X ⊤ S=\frac1NXX^\top S = N 1 X X ⊤ 의 첫 주성분 b 1 b_1 b 1 을 구하고 S b 1 = λ 1 b 1 Sb_1=\lambda_1b_1 S b 1 = λ 1 b 1 을 확인하라.
(e) 4 × 4 4\times4 4 × 4 대신 3 × 3 3\times3 3 × 3 을 푸는 것이 어떤 이득인지 말하라.
10-6. 표준화 [계산·서술]# 근거: 10.6절, 식 (10.58)~(10.61)
두 특징의 평균이 ( 170 , 0 ) (170,\,0) ( 170 , 0 ) , 표준편차가 ( 10 , 1 ) (10,\,1) ( 10 , 1 ) , 공분산이 S = [ 100 6 6 1 ] S=\begin{bmatrix}100&6\\6&1\end{bmatrix} S = [ 100 6 6 1 ] 이다.
(a) 표준화하지 않은 S S S 의 고유값을 구하고, 첫 주성분이 거의 ( 1 , 0 ) (1,0) ( 1 , 0 ) 임을 보여라.
(b) 표준화한 뒤의 공분산을 구하고 그 고유값과 고유벡터를 구하라.
(c) 새 점 x ∗ = ( 185 , 0.5 ) x_*=(185,\,0.5) x ∗ = ( 185 , 0.5 ) 를 (10.58)~(10.61)의 순서대로 표준화 → M = 1 M=1 M = 1 사영 → 원래 단위로 되돌리기를 하라.
10-7. 거듭제곱법 [계산]# 근거: 10.4절, 식 (10.52)
S = [ 2 1 1 2 ] S=\begin{bmatrix}2&1\\1&2\end{bmatrix} S = [ 2 1 1 2 ] , x 0 = ( 1 , 0 ) x_0=(1,0) x 0 = ( 1 , 0 ) .
(a) x 1 , x 2 , x 3 x_1,x_2,x_3 x 1 , x 2 , x 3 을 구하라(매번 정규화).
(b) 레일리 몫 x k ⊤ S x k x_k^\top Sx_k x k ⊤ S x k 로 가장 큰 고유값을 추정하라.
(c) 수렴 속도가 λ 2 / λ 1 \lambda_2/\lambda_1 λ 2 / λ 1 에 달린 이유를 설명하라. (교재 밖 보충 )
10-8. 확률적 PCA [계산]# 근거: 10.7절, 식 (10.63)~(10.75)
D = 2 D=2 D = 2 , M = 1 M=1 M = 1 , B = [ 2 , 1 ] ⊤ B=[2,1]^\top B = [ 2 , 1 ] ⊤ , μ = 0 \mu=0 μ = 0 , σ 2 = 1 \sigma^2=1 σ 2 = 1 .
(a) p ( x ) p(x) p ( x ) 의 공분산 B B ⊤ + σ 2 I BB^\top+\sigma^2I B B ⊤ + σ 2 I 를 구하라(10.70b).
(b) Cov [ x , z ] \text{Cov}[x,z] Cov [ x , z ] 와 결합분포 p ( x , z ) p(x,z) p ( x , z ) 를 써라(10.71)(10.72).
(c) x = ( 3 , 1 ) x=(3,1) x = ( 3 , 1 ) 일 때 사후분포 p ( z ∣ x ) = N ( m , C ) p(z\mid x)=\mathcal N(m,C) p ( z ∣ x ) = N ( m , C ) 를 구하라(10.73)(10.75).
(d) σ 2 → 0 \sigma^2\to0 σ 2 → 0 일 때 사후 평균의 극한을 구하고, 이것이 PCA 좌표와 같아지는 조건을 말하라.
10-9. 재구성 오차는 버린 고유값의 합이다 [유도]# 근거: 10.3절, 식 (10.29)~(10.44), (10.25)
(a) 최적 좌표가 z i n = b i ⊤ x n z_{in}=b_i^\top x_n z in = b i ⊤ x n 임을 유도하라(10.30)~(10.32).
(b) 정규직교 기저 b 1 , … , b D b_1,\dots,b_D b 1 , … , b D 로 x n − x ~ n = ∑ j = M + 1 D b j b j ⊤ x n x_n-\tilde x_n=\sum_{j=M+1}^Db_jb_j^\top x_n x n − x ~ n = ∑ j = M + 1 D b j b j ⊤ x n 임을 보여라(10.38).
(c) J M = ∑ j = M + 1 D b j ⊤ S b j J_M=\sum_{j=M+1}^Db_j^\top Sb_j J M = ∑ j = M + 1 D b j ⊤ S b j 를 보여라(10.43).
(d) J M = ∑ j = M + 1 D λ j J_M=\sum_{j=M+1}^D\lambda_j J M = ∑ j = M + 1 D λ 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\} { − 3 , − 2.5 , − 1 , 0 , 2 , 4 , 5 } , 초기 성분 N ( x ∣ − 4 , 1 ) \mathcal N(x\mid-4,1) N ( x ∣ − 4 , 1 ) , N ( x ∣ 0 , 0.2 ) \mathcal N(x\mid0,0.2) N ( x ∣ 0 , 0.2 ) , N ( x ∣ 8 , 3 ) \mathcal N(x\mid8,3) N ( x ∣ 8 , 3 ) , 비율 π k = 1 3 \pi_k=\frac13 π k = 3 1 .
x = − 1 x=-1 x = − 1 과 x = 2 x=2 x = 2 의 책임도 r n 1 , r n 2 , r n 3 r_{n1},r_{n2},r_{n3} r n 1 , r n 2 , r n 3 을 계산하라.
상수: 1 2 π ≈ 0.3989 \frac1{\sqrt{2\pi}}\approx0.3989 2 π 1 ≈ 0.3989 , 1 0.4 π ≈ 0.8921 \frac1{\sqrt{0.4\pi}}\approx0.8921 0.4 π 1 ≈ 0.8921 , 1 6 π ≈ 0.2303 \frac1{\sqrt{6\pi}}\approx0.2303 6 π 1 ≈ 0.2303 , e − 2.5 ≈ 0.08208 e^{-2.5}\approx0.08208 e − 2.5 ≈ 0.08208 , e − 4.5 ≈ 0.01111 e^{-4.5}\approx0.01111 e − 4.5 ≈ 0.01111 , e − 6 ≈ 0.002479 e^{-6}\approx0.002479 e − 6 ≈ 0.002479 , e − 10 ≈ 4.54 × 10 − 5 e^{-10}\approx4.54\times10^{-5} e − 10 ≈ 4.54 × 1 0 − 5 , e − 13.5 ≈ 1.37 × 10 − 6 e^{-13.5}\approx1.37\times10^{-6} e − 13.5 ≈ 1.37 × 1 0 − 6 , e − 18 ≈ 1.52 × 10 − 8 e^{-18}\approx1.52\times10^{-8} e − 18 ≈ 1.52 × 1 0 − 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)의 책임도 행렬을 쓴다(행이 데이터, 열이 성분).
x n x_n x n r n 1 r_{n1} r n 1 r n 2 r_{n2} r n 2 r n 3 r_{n3} r n 3 − 3 -3 − 3 1.0 0.0 0.0 − 2.5 -2.5 − 2.5 1.0 0.0 0.0 − 1 -1 − 1 0.057 0.943 0.0 0 0 0 0.001 0.999 0.0 2 2 2 0.0 0.066 0.934 4 4 4 0.0 0.0 1.0 5 5 5 0.0 0.0 1.0
(a) N k N_k N k 를 구하라.
(b) 새 평균 μ k \mu_k μ k 를 구하라.
(c) 새 평균을 써서 새 분산 σ k 2 \sigma_k^2 σ k 2 를 구하라.
(d) 새 비율 π k \pi_k π k 를 구하라.
(e) 교재의 값(평균 − 2.7 , − 0.4 , 3.7 -2.7,-0.4,3.7 − 2.7 , − 0.4 , 3.7 , 분산 0.14 , 0.44 , 1.53 0.14,0.44,1.53 0.14 , 0.44 , 1.53 , 비율 0.29 , 0.29 , 0.42 0.29,0.29,0.42 0.29 , 0.29 , 0.42 )과 비교하라.
11-3. 갱신식 유도: 평균과 분산 [유도]# 근거: 정리 11.1, 식 (11.10)~(11.24), 정리 11.2
(a) 로그 가능도 (11.10)을 μ k \mu_k μ k 로 미분해 μ k new = 1 N k ∑ n r n k x n \mu_k^{\text{new}}=\frac1{N_k}\sum_nr_{nk}x_n μ k new = N k 1 ∑ n r nk x n (11.20)을 유도하라.
(b) 이 식이 닫힌 해가 아닌 이유를 말하라.
(c) 1차원(Σ k = σ k 2 \Sigma_k=\sigma_k^2 Σ k = σ k 2 )에서 σ k 2 \sigma_k^2 σ k 2 로 미분해 σ k 2 , new = 1 N k ∑ n r n k ( x n − μ k ) 2 \sigma_k^{2,\text{new}}=\frac1{N_k}\sum_nr_{nk}(x_n-\mu_k)^2 σ k 2 , new = N k 1 ∑ n r nk ( x n − μ k ) 2 을 유도하라.
11-4. 혼합 비율 갱신: 라그랑주 승수 [유도]# 근거: 정리 11.3, 식 (11.42)~(11.49)
(a) 제약 ∑ k π k = 1 \sum_k\pi_k=1 ∑ k π k = 1 을 붙인 라그랑주 함수를 써라.
(b) π k \pi_k π k 로 미분해 π k = − N k / λ \pi_k=-N_k/\lambda π k = − N k / λ 를 얻어라.
(c) λ = − N \lambda=-N λ = − N 임을 보이고 π k = N k / N \pi_k=N_k/N π k = N k / N 을 얻어라.
11-5. 잠재변수 관점과 결정 경계 [계산]# 근거: 11.4절, 식 (11.58)~(11.72)
K = 2 K=2 K = 2 , π = ( 0.6 , 0.4 ) \pi=(0.6,\,0.4) π = ( 0.6 , 0.4 ) , μ = ( 0 , 2 ) \mu=(0,\,2) μ = ( 0 , 2 ) , σ 1 2 = σ 2 2 = 1 \sigma_1^2=\sigma_2^2=1 σ 1 2 = σ 2 2 = 1 .
(a) 원-핫 잠재변수 z z z 의 p ( z ) p(z) p ( z ) 와 p ( x ∣ z ) p(x\mid z) p ( x ∣ z ) 에서 p ( x ) p(x) p ( x ) 를 유도하라(11.65)~(11.66).
(b) x = 1 x=1 x = 1 에서 p ( z 1 = 1 ∣ x ) p(z_1=1\mid x) p ( z 1 = 1 ∣ x ) 를 구하라. 결과가 π 1 \pi_1 π 1 과 같은 이유는?
(c) 두 책임도가 같아지는 x x x 를 구하라(log 1.5 ≈ 0.4055 \log1.5\approx0.4055 log 1.5 ≈ 0.4055 ).
11-6. EM 한 바퀴: 잘 갈라진 데이터 [계산]# 근거: 11.3절, 식 (11.53)~(11.56)
데이터 0 , 1 , 5 , 6 0,1,5,6 0 , 1 , 5 , 6 , K = 2 K=2 K = 2 , 초기값 π = ( 0.5 , 0.5 ) \pi=(0.5,0.5) π = ( 0.5 , 0.5 ) , μ = ( 1 , 5 ) \mu=(1,5) μ = ( 1 , 5 ) , σ 1 2 = σ 2 2 = 1 \sigma_1^2=\sigma_2^2=1 σ 1 2 = σ 2 2 = 1 . (e − 8 ≈ 3.35 × 10 − 4 e^{-8}\approx3.35\times10^{-4} e − 8 ≈ 3.35 × 1 0 − 4 , e − 12 ≈ 6.1 × 10 − 6 e^{-12}\approx6.1\times10^{-6} e − 12 ≈ 6.1 × 1 0 − 6 )
(a) E 단계: 네 점의 책임도를 구하라.
(b) M 단계: N k N_k N k , μ k \mu_k μ k , σ k 2 \sigma_k^2 σ k 2 , π k \pi_k π k 를 구하라.
(c) 로그 가능도가 − 7.448 -7.448 − 7.448 에서 − 5.676 -5.676 − 5.676 으로 늘었다. 왜 늘어야 하는가?
(d) 책임도가 거의 0 또는 1일 때 이 계산은 K-평균의 어떤 단계와 같아지는가? (교재 밖 보충 포함)
11-7. 최대가능도의 특이점 [계산·서술]# 근거: 11.5절(더 읽을거리)
데이터 0 , 1 , 5 , 6 0,1,5,6 0 , 1 , 5 , 6 , K = 2 K=2 K = 2 , π = ( 0.5 , 0.5 ) \pi=(0.5,0.5) π = ( 0.5 , 0.5 ) , μ 1 = 0 \mu_1=0 μ 1 = 0 , σ 1 2 = ε \sigma_1^2=\varepsilon σ 1 2 = ε , μ 2 = 3 \mu_2=3 μ 2 = 3 , σ 2 2 = 10 \sigma_2^2=10 σ 2 2 = 10 .
(a) x = 0 x=0 x = 0 의 로그 밀도 항이 ε → 0 \varepsilon\to0 ε → 0 에서 + ∞ +\infty + ∞ 로 감을 보여라.
(b) 나머지 점의 항은 아래로 유계임을 보이고, 전체 로그 가능도가 + ∞ +\infty + ∞ 로 감을 결론지어라.
(c) 이것이 왜 좋은 모델이 아닌지, 실제로 어떻게 막는지 말하라(교재 밖 보충 ).
11-8. EM과 GMM 개념 [서술]# 근거: 11.2~11.5절
(a) 책임도 r n k r_{nk} r nk 가 사후확률 p ( z n k = 1 ∣ x n ) p(z_{nk}=1\mid x_n) p ( z nk = 1 ∣ 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) w = ( 3 , 4 ) , b = − 5 b=-5 b = − 5 .
(a) 점 ( 3 , 4 ) (3,4) ( 3 , 4 ) , ( 0 , 0 ) (0,0) ( 0 , 0 ) , ( 1 , 0.5 ) (1,\,0.5) ( 1 , 0.5 ) 의 f ( x ) = ⟨ w , x ⟩ + b f(x)=\langle w,x\rangle+b f ( x ) = ⟨ w , x ⟩ + b 와 부호 있는 거리 f ( x ) / ∥ w ∥ f(x)/\|w\| f ( x ) /∥ w ∥ 를 구하라.
(b) ( w , b ) (w,b) ( w , b ) 를 ( 6 , 8 , − 10 ) (6,8,-10) ( 6 , 8 , − 10 ) 으로 두 배 해도 같은 초평면임을 보이고, y f ( x ) y f(x) y f ( x ) 값은 두 배가 됨을 보여라.
(c) 그래서 (12.15)처럼 “가장 가까운 점에서 y n ( ⟨ w , x n ⟩ + b ) = 1 y_n(\langle w,x_n\rangle+b)=1 y n (⟨ w , x n ⟩ + b ) = 1 “로 크기를 고정하면 마진이 1 / ∥ w ∥ 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) ( 2 , 2 ) , ( 2 , 3 ) , 음성 ( 0 , 0 ) , ( 0 , 1 ) (0,0),\ (0,1) ( 0 , 0 ) , ( 0 , 1 ) .
(a) 수직선 x 1 = 1 x_1=1 x 1 = 1 (w = ( 1 , 0 ) w=(1,0) w = ( 1 , 0 ) , b = − 1 b=-1 b = − 1 )이 모든 제약 (12.15)를 만족하고 마진이 1임을 확인하라.
(b) 두 클래스 볼록 껍질 사이의 가장 가까운 두 점을 찾고, 그로부터 최대 마진 초평면 ( w , b ) (w,b) ( w , b ) 를 구하라.
(c) (b)의 마진이 (a)보다 큼을 확인하라.
(d) 쌍대 변수 α \alpha α 를 구해 (12.38)로 w w w 를 다시 만들고, (12.42)로 b b b 를 구하고, 쌍대 목적값과 원 목적값이 같음을 확인하라.
12-3. 소프트 마진 쌍대 유도 [유도]# 근거: 12.3.1절, 식 (12.34)~(12.41)
(a) 원 문제 (12.26a)의 라그랑주 함수를 써라.
(b) w w w , b b b , ξ n \xi_n ξ n 으로 편미분하라.
(c) 대입해서 (12.40)을 얻고 쌍대 문제 (12.41)을 써라.
(d) 0 ≤ α n ≤ C 0\le\alpha_n\le C 0 ≤ α n ≤ C 가 나오는 이유를 말하라.
(e) 하드 마진의 쌍대와 무엇이 다른가?
12-4. 힌지 손실과 여유 변수 [계산]# 근거: 12.2.412.2.5절, 식 (12.26), (12.28)(12.33)
w = ( 1 , − 1 ) w=(1,-1) w = ( 1 , − 1 ) , b = 0 b=0 b = 0 , C = 1 C=1 C = 1 .
점 x n x_n x n y n y_n y n 1 ( 3 , 1 ) (3,1) ( 3 , 1 ) + 1 +1 + 1 2 ( 1 , 0 ) (1,0) ( 1 , 0 ) + 1 +1 + 1 3 ( 0 , 0.5 ) (0,\,0.5) ( 0 , 0.5 ) − 1 -1 − 1 4 ( 2 , 1 ) (2,1) ( 2 , 1 ) − 1 -1 − 1
(a) 각 점의 t n = y n f ( x n ) t_n=y_nf(x_n) t n = y n f ( x n ) , 힌지 손실, 여유 변수 ξ n \xi_n ξ 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 ) (1,-1),\ (-1,1) ( 1 , − 1 ) , ( − 1 , 1 ) .
(a) 원래 2차원에서 직선으로 분리할 수 없음을 볼록 껍질로 보여라.
(b) 다항 커널 k ( x , x ′ ) = ( x ⊤ x ′ ) 2 k(x,x')=(x^\top x')^2 k ( x , x ′ ) = ( x ⊤ x ′ ) 2 의 특징 사상이 ϕ ( x ) = ( x 1 2 , 2 x 1 x 2 , x 2 2 ) \phi(x)=(x_1^2,\ \sqrt2x_1x_2,\ x_2^2) ϕ ( x ) = ( x 1 2 , 2 x 1 x 2 , x 2 2 ) 임을 보이고, x = ( 1 , 2 ) x=(1,2) x = ( 1 , 2 ) , x ′ = ( 3 , − 1 ) x'=(3,-1) x ′ = ( 3 , − 1 ) 로 숫자를 확인하라.
(c) 네 점의 그람 행렬 K K K 를 구하고 양의 준정부호인지 확인하라(12.53).
(d) 특징 공간에서 하드 마진 SVM의 w w w , b b b , 마진을 구하고 결정 함수를 원래 좌표로 써라.
(e) 쌍대 변수 α \alpha α 를 구하고, 커널만으로 결정 함수를 다시 써서 (d)와 같음을 확인하라.
(f) k ( x , x ′ ) = − ∥ x − x ′ ∥ 2 k(x,x')=-\|x-x'\|^2 k ( x , x ′ ) = − ∥ x − x ′ ∥ 2 이 커널이 될 수 없음을 점 0 0 0 과 1 1 1 (1차원)로 보여라.
12-6. C C C 와 정규화 계수 λ \lambda λ 의 관계 [유도·서술]# 근거: 12.2.5절, 식 (12.31), 8.2.3절, 식 (8.6), (8.12)
(a) (12.31)을 C N CN CN 으로 나누어 “평균 손실 + λ × \lambda\times λ × 정규화항” 꼴로 바꾸고 λ \lambda λ 와 C C C 의 관계를 구하라.
(b) C C C 가 클 때와 작을 때 마진과 훈련 오류가 어떻게 되는지 말하라.
(c) C → ∞ C\to\infty C → ∞ 이면 어떤 문제가 되는가?
12-7. 커널 SVM으로 예측하기 [계산]# 근거: 식 (12.38), 12.4절
RBF 커널 k ( x , x ′ ) = exp ( − 1 2 ∥ x − x ′ ∥ 2 ) k(x,x')=\exp\!\big(-\frac12\|x-x'\|^2\big) k ( x , x ′ ) = exp ( − 2 1 ∥ x − x ′ ∥ 2 ) . 서포트 벡터 x 1 = ( 0 , 0 ) x_1=(0,0) x 1 = ( 0 , 0 ) , y 1 = + 1 y_1=+1 y 1 = + 1 , α 1 = 1 \alpha_1=1 α 1 = 1 과 x 2 = ( 1 , 1 ) x_2=(1,1) x 2 = ( 1 , 1 ) , y 2 = − 1 y_2=-1 y 2 = − 1 , α 2 = 1 \alpha_2=1 α 2 = 1 , b = 0 b=0 b = 0 이 주어졌다고 하자. 결정 함수는 f ( x ) = ∑ n α n y n k ( x n , x ) + b f(x)=\sum_n\alpha_ny_nk(x_n,x)+b f ( x ) = ∑ n α n y n k ( x n , x ) + b 다.
(a) x ∗ = ( 0 , 1 ) x_*=(0,1) x ∗ = ( 0 , 1 ) , ( 0 , 0.5 ) (0,\,0.5) ( 0 , 0.5 ) , ( 2 , 2 ) (2,2) ( 2 , 2 ) 에서 f f f 와 예측 레이블을 구하라. (e − 0.5 ≈ 0.6065 e^{-0.5}\approx0.6065 e − 0.5 ≈ 0.6065 , e − 0.125 ≈ 0.8825 e^{-0.125}\approx0.8825 e − 0.125 ≈ 0.8825 , e − 0.625 ≈ 0.5353 e^{-0.625}\approx0.5353 e − 0.625 ≈ 0.5353 , e − 1 ≈ 0.3679 e^{-1}\approx0.3679 e − 1 ≈ 0.3679 , e − 4 ≈ 0.0183 e^{-4}\approx0.0183 e − 4 ≈ 0.0183 )
(b) 결정 경계를 구하라.
(c) 왜 w w w 를 직접 계산하지 않는가?
12-8. SVM 개념 [서술]# 근거: 12.2~12.4절
(a) 서포트 벡터란 무엇인가? α n = 0 \alpha_n=0 α n = 0 인 점을 지워도 해가 그대로인 이유는?
(b) 0 < α n < C 0<\alpha_n<C 0 < α n < C 인 점과 α n = C \alpha_n=C α n = C 인 점은 기하학적으로 어떻게 다른가?
(c) 원 문제 대신 쌍대 문제를 푸는 이유 두 가지를 말하라.
(d) 소프트 마진 원 문제 (12.26a)가 볼록 최적화 문제인 이유는?
OX 퀴즈# 8~12장 공통. 맞으면 O, 틀리면 X. 틀린 것은 고쳐 말해 본다.
MML의 PCA에서 데이터 공분산은 1 N − 1 \frac1{N-1} N − 1 1 로 나눈다. 교재 (8.48)의 AIC는 값이 작을수록 좋은 모델이다. MAP 추정은 매개변수의 사후분포 전체를 준다. 사전분포 N ( 0 , b 2 I ) \mathcal N(0,b^2I) N ( 0 , b 2 I ) 의 MAP 추정은 λ = σ 2 / b 2 \lambda=\sigma^2/b^2 λ = σ 2 / b 2 인 정규화 최소제곱 (9.32)과 같다. 베이즈 선형 회귀의 사후 예측 분산은 노이즈 분산 σ 2 \sigma^2 σ 2 보다 작아질 수 있다. a → d ← c a\to d\leftarrow c a → d ← c 에서 d d d 를 관측하면 a a a 와 c c c 가 독립이 된다.PCA의 첫 주성분은 공분산의 가장 작은 고유값의 고유벡터다. PCA의 평균 제곱 재구성 오차는 버린 고유값의 합이다. GMM의 책임도 r n k r_{nk} r nk 를 k k k 에 대해 더하면 1이다. EM 알고리즘은 항상 로그 가능도의 전역 최댓값에 수렴한다. 가장 가까운 점에서 y n ( ⟨ w , x n ⟩ + b ) = 1 y_n(\langle w,x_n\rangle+b)=1 y n (⟨ w , x n ⟩ + b ) = 1 로 크기를 고정하면 하드 마진 SVM의 마진은 1 / ∥ w ∥ 1/\|w\| 1/∥ w ∥ 다. 소프트 마진 SVM에서 α n = 0 \alpha_n=0 α n = 0 인 점은 w w w 에 기여하지 않는다. 유효한 커널의 그람 행렬은 양의 준정부호다. 9장의 데이터 행렬은 D × N D\times N D × N , 10장은 N × D N\times D N × D 다. 주변 가능도(증거)는 매개변수가 많은 모델을 항상 선호한다.