南瓜书第三章(3) 线性判别分析

本文介绍了线性判别分析(LDA)的基本原理,包括其作为监督降维技术的思路,以及如何通过最大化类间距离和最小化类内距离来提高样本可区分性。文章还涉及了损失函数的推导,使用拉格朗日乘子法求解最优方向w,并讨论了广义特征值和广义瑞利商的概念。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

广义线性模型可以通过“联系函数”完成分类的目的,那线性模型应如何直接完成“分类”呢?

    接下来我们引入最基本的一种线性判别分析(LDA)来解答

 对于西瓜书中部分符号的说明:

1、算法原理

线性判别分析是一个经典的二分类算法。

(1)主要思想:

给定训练样例集,设法将样例投影到一条直线(低维空间)上使得同类样例的投影点尽可能接近、异类样例的投影点尽可能远离:在对新样本进行分类时,将其投影到同样的这条直线上,再根据投影点的位置来确定新样本的类别.图 3.3 给出了一个二维示意图。(因此也可以被看作一种“监督降维技术”)

(2)实例:

       线性判别回归即以一种基于降维的方式将所有的样本映射到一维向量w上,然

<think>好的,我现在需要帮助用户找到《机器学习西瓜第三章的代码实现,并且用户特别提到用Python。首先,我得回忆一下第三章的内容。西瓜第三章主要是关于线性模型,包括线性回归、对数几率回归(也就是逻辑回归)、线性判别分析(LDA)以及多分类学习等。 接下来,我要考虑用户的需求。用户可能已经看过理论部分,现在需要实际的代码来加深理解或者应用。不过西瓜本身并没有提供完整的代码,通常需要参考南瓜或者其他开源项目,比如Datawhale的南瓜项目,他们可能会有配套的代码实现。 然后,我需要检查提供的引用内容。根据用户给出的引用[2],里面提到了第三章的线性模型和类别划分的编码矩阵,还有解码过程。引用[3]提到了第五章的神经网络代码,可能不太相关。引用[1]则涉及聚类算法的代码,属于后面的章节。所以主要的线索在引用[2]里,但具体代码实现可能需要更多信息。 接下来,我需要确定用户是否有其他资源可用。比如,南瓜通常会有配套的代码,或者一些开源社区如GitHub上的项目。可能需要引导用户去这些地方查找。另外,用户可能不知道南瓜的存在,所以需要提及这一点。 然后,考虑用户可能的深层需求。他们可能想通过代码来实践第三章的内容,比如实现线性回归、逻辑回归或者多分类问题。因此,除了直接指出资源的位置,最好也能提供一些简单的示例代码,帮助他们入门。比如,用sklearn实现线性回归或逻辑回归,或者手动编写梯度下降的代码。 另外,用户可能遇到的问题包括数据准备、参数设置、算法理解等,所以示例代码需要简单明了,并附有必要的注释。同时,需要注意代码的正确性,比如在逻辑回归中使用sigmoid函数,处理多分类时的策略如OvR或OvO。 最后,要确保回答结构清晰,分点列出资源链接和示例代码,并提醒用户参考籍的具体页码,如引用中提到的南瓜第56页和西瓜第197页,但第三章可能需要不同的页码。不过用户的问题集中在第三章,可能需要确认正确的引用页码,但根据现有信息,引用[2]提到了第三章,所以可以引用它。 总结步骤: 1. 指出官方代码可能不完整,建议参考南瓜和开源项目。 2. 提供GitHub链接,如Datawhale的南瓜项目。 3. 给出线性回归和逻辑回归的示例代码,使用sklearn和手动实现。 4. 提及多分类学习的实现方法,如OvR和OvO。 5. 添加相关问题,帮助用户深入理解。</think>根据《机器学习》(西瓜第三章“线性模型”的理论内容,以下是代码实现思路及参考资源整理: --- ### 一、第三章核心算法及对应代码实现 第三章主要包含以下模型,代码实现方式如下: 1. **线性回归** - 手动实现梯度下降法: ```python import numpy as np # 生成数据 X = np.array([[1, 2], [2, 3], [3, 4]]) y = np.dot(X, np.array([1, 2])) + 3 + np.random.normal(0, 0.1, 3) # 梯度下降 w = np.random.randn(2) b = 0 lr = 0.01 for _ in range(1000): y_pred = np.dot(X, w) + b grad_w = np.dot(X.T, (y_pred - y)) / len(y) grad_b = np.mean(y_pred - y) w -= lr * grad_w b -= lr * grad_b print(f"权重: {w}, 偏置: {b}") ``` - 使用 `sklearn` 实现: ```python from sklearn.linear_model import LinearRegression model = LinearRegression().fit(X, y) print(f"权重: {model.coef_}, 偏置: {model.intercept_}") ``` 2. **对数几率回归(逻辑回归)** - Sigmoid函数与交叉熵损失实现: ```python def sigmoid(z): return 1 / (1 + np.exp(-z)) # 二分类逻辑回归梯度下降 X = np.array([[1, 2], [2, 3], [3, 4]]) y = np.array([0, 1, 1]) w = np.random.randn(2) b = 0 for _ in range(1000): logits = np.dot(X, w) + b probs = sigmoid(logits) grad_w = np.dot(X.T, (probs - y)) / len(y) grad_b = np.mean(probs - y) w -= lr * grad_w b -= lr * grad_b ``` - 使用 `sklearn` 实现: ```python from sklearn.linear_model import LogisticRegression model = LogisticRegression().fit(X, y) ``` 3. **线性判别分析(LDA)** ```python from sklearn.discriminant_analysis import LinearDiscriminantAnalysis model = LinearDiscriminantAnalysis().fit(X, y) ``` 4. **多分类学习(OvR/OvO策略)** ```python from sklearn.multiclass import OneVsRestClassifier, OneVsOneClassifier from sklearn.svm import SVC ovr_model = OneVsRestClassifier(SVC()).fit(X, y) ovo_model = OneVsOneClassifier(SVC()).fit(X, y) ``` --- ### 二、官方代码资源 1. **南瓜配套代码** - GitHub仓库地址:[https://github.com/datawhalechina/pumpkin-book](https://github.com/datawhalechina/pumpkin-book) - 第三章代码文件通常命名为 `ch03_linear_models.ipynb`,包含线性回归、逻辑回归的推导与实现[^2]。 2. **第三方实现参考** - 线性模型完整实现(含梯度下降可视化):[西瓜线性模型代码示例](https://github.com/Vay-keen/Machine-learning-learning-notes) --- ### 三、关键理论对照 - **编码与解码策略**:多分类问题中,`OneVsRest` 对应“一对其余”(OvR),`OneVsOne` 对应“一对一”(OvO),分类结果通过编码矩阵组合后解码得到最终预测。 - **性能优化**:逻辑回归可通过调整正则化参数(如 `penalty=&#39;l2&#39;`)或使用随机梯度下降(`SGDClassifier`)优化。 ---
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值