多元统计分析及R语言建模（第五版）——第6章 判别分析课后习题

第6章 判别分析

文章会用到的数据请在这个网址下下载多元统计分析及R语言建模（第五版）数据
练习题
1）考虑两个数据集x1 = [3 7 2 4 4 7]，x2 = [6 9 5 7 4 8]
（1）计算Fisher线性判别函数

library(MASS)
library(openxlsx)
x1 = c(3,2,4,6,5,4)
x2 = c(7,4,7,9,7,8)
y = c(1,1,1,2,2,2)
(data1 = data.frame(x1,x2,y))
plot(data1[,-3])
text(x1,x2,y,adj = -0.5)

x1_bar = data1[1:3,-3]
x1_bar

x2_bar = data1[4:6,-3]
x2_bar

(x1_mean = apply(x1_bar,2,mean))
(x2_mean = apply(x2_bar,2,mean))

Sp = matrix(c(1,1,1,2),ncol = 2)
Sp_re = solve(Sp)
(2 * cov(x1_bar) +