問23.2

Q23.2.R

# -------------
# 第23章 判別分析
# 問23.2 p.213
# -------------
# データ
LD1 <- c(-0.1,-0.1,0.2,0.4,0.9)
LD2 <- c(0.0,0.2,-0.2,0.9,-0.3)

x1 <- c(1.2,1.4,0.3,6.0,4.6)
x2 <- c(1.4,1.4,0.3,0.7,2.1)
x3 <- c(0.8,0.2,1.6,3.6,9.2)

xL <- c(2.5,-1.1,0,-0.8,1.0)

# 射影後の座標
f <- function(x){
  c(LD1 %*% x,LD2 %*% x)
}

# 距離計算
matrix( c(f(x1),f(x2),f(x3),f(xL)), byrow = T, ncol=2 ) |> dist()

> matrix( c(f(x1),f(x2),f(x3),f(xL)), byrow = T, ncol=2 ) |> dist()
         1        2        3
2 5.952520                  
3 5.411248 7.990025         
4 8.052354 2.100405 9.608517
> 

となるので、新しいサンプルから最も近い重心を持つクラスは、x2すなわちクラスFの群となる。