要求动手从0实现 logistic 回归(只借助Tensor和Numpy相关的库)在人工构造的数据集上进 行训练和测试,并从loss、训练集以及测试集上的准确率等多个角度对结果进行分析 (可借助nn.BCELoss或nn.BCEWithLogitsLoss作为损失函数,从零实现二元交叉熵为选作)

时间: 2024-04-14 09:26:41 浏览: 22
好的,我将向您展示如何从零开始使用PyTorch和NumPy来实现逻辑回归,并在人工构造的数据集上进行训练和测试。 首先,我们需要导入所需的库: ```python import torch import numpy as np ``` 接下来,我们将创建一个人工构造的数据集。我们假设数据集包含两个特征和一个二元标签(0或1): ```python # 设置随机种子,以便结果可复现 torch.manual_seed(42) np.random.seed(42) # 构造人工数据集 num_samples = 1000 num_features = 2 # 生成特征矩阵 X = np.random.randn(num_samples, num_features) # 生成标签 true_weights = np.array([2, -3]) true_bias = 1.5 logits = np.dot(X, true_weights) + true_bias y = (1 / (1 + np.exp(-logits))) > 0.5 # 将数据转换为Tensor X = torch.from_numpy(X).float() y = torch.from_numpy(y).float() ``` 现在我们已经生成了数据集,接下来是模型的定义和训练过程。 ```python # 定义模型类 class LogisticRegression: def __init__(self, num_features): self.weights = torch.zeros(num_features, 1) self.bias = torch.zeros(1) def forward(self, X): linear = torch.matmul(X, self.weights) + self.bias return torch.sigmoid(linear) def backward(self, X, y, y_pred, learning_rate): gradient_weights = torch.matmul(X.t(), (y_pred - y).reshape(-1, 1)) gradient_bias = torch.sum(y_pred - y) self.weights -= learning_rate * gradient_weights self.bias -= learning_rate * gradient_bias def train(self, X, y, num_epochs, learning_rate): for epoch in range(num_epochs): y_pred = self.forward(X) self.backward(X, y, y_pred, learning_rate) # 创建模型实例并进行训练 model = LogisticRegression(num_features) model.train(X, y, num_epochs=1000, learning_rate=0.03) ``` 现在模型已经训练好了,我们可以使用训练好的模型对训练集和测试集进行预测,并分析结果的准确率。 ```python # 在训练集上进行预测 y_pred_train = model.forward(X) train_accuracy = ((y_pred_train > 0.5) == y).float().mean().item() # 生成测试集数据 X_test = np.random.randn(200, num_features) logits_test = np.dot(X_test, true_weights) + true_bias # 在测试集上进行预测 y_test = (1 / (1 + np.exp(-logits_test))) > 0.5 X_test = torch.from_numpy(X_test).float() y_test = torch.from_numpy(y_test).float() y_pred_test = model.forward(X_test) test_accuracy = ((y_pred_test > 0.5) == y_test).float().mean().item() ``` 最后,我们可以计算损失函数的值并打印出结果。 ```python # 计算损失函数的值 loss_fn = torch.nn.BCELoss() loss_train = loss_fn(y_pred_train, y).item() loss_test = loss_fn(y_pred_test, y_test).item() # 打印结果 print(f"训练集准确率: {train_accuracy}") print(f"测试集准确率: {test_accuracy}") print(f"训练集损失: {loss_train}") print(f"测试集损失: {loss_test}") ``` 通过以上代码,您可以从零开始实现逻辑回归,并在人工构造的数据集上进行训练和测试。您可以根据需要调整模型的超参数,如学习率和训练轮数,以获得更好的结果。

相关推荐

最新推荐

recommend-type

python代码实现逻辑回归logistic原理

主要介绍了python代码实现逻辑回归logistic原理,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
recommend-type

任务三、titanic数据集分类问题

任务说明:分别用下列方法完成titanic数据集分类,给分析每种方法参数设置、方法优缺点分析 logistic回归 决策树 SVM 神经网络
recommend-type

STM32H562实现FreeRTOS内存管理【支持STM32H系列单片机】.zip

STM32H562 FreeRTOS驱动程序,支持STM32H系列单片机。 项目代码可直接运行~
recommend-type

中文翻译Introduction to Linear Algebra, 5th Edition 2.1节

中文翻译Introduction to Linear Algebra, 5th Edition 2.1节 线性代数的核心问题是求解方程组。这些方程都是线性的,即未知数仅与数相乘——我们绝不会 遇见 x 乘以 y。我们的第一个线性方程组较小。接下来你来看看它引申出多远: 两个方程 两个未知数 x − 2y = 1 3x + 2y = 11 (1) 我们一次从一个行开始。第一个方程 x − 2y = 1 得出了 xy 平面的一条直线。由于点 x = 1, y = 0 解 出该方程,因此它在这条直线上。因为 3 − 2 = 1,所以点 x = 3, y = 1 也在这条直线上。若我们选择 x = 101,那我们求出 y = 50。 这条特定直线的斜率是 12,是因为当 x 变化 2 时 y 增加 1。斜率在微积分中很重要,然而这是线 性代数! 图 2.1 将展示第一条直线 x − 2y = 1。此“行图”中的第二条直线来自第二个方程 3x + 2y = 11。你 不能错过两条线的交点 x = 3, y = 1。点 (3, 1) 位于两条线上并且解出两个方程。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

识别MATLAB微分方程求解中的混沌行为:分析非线性方程混沌行为的实用技巧

![matlab求解微分方程](https://img-blog.csdnimg.cn/2021062810300367.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MTQ3OTY3OA==,size_16,color_FFFFFF,t_70) # 1. MATLAB微分方程求解概述 微分方程是描述物理、工程和金融等领域中动态系统的数学模型。MATLAB提供了强大的工具来求解微分方程,包括内置求解器和自定义函数
recommend-type

physon如何做ERP系统

作为一个AI语言模型,我可以给您提供一些基本的建议和思路: 1.明确需求和目标:在开始开发ERP系统之前,需要明确您的业务需求和目标,了解您的业务流程、流程数据、员工角色和权限、财务管理等方面的要求,这样才能更好的设计和开发ERP系统。 2.选择技术框架和开发工具:选择合适的技术框架和开发工具是ERP系统开发的关键。选择一种流行的技术框架和工具可以提高开发效率和质量。 3.设计数据库:ERP系统需要一个功能强大的数据库来存储数据。设计数据库需要考虑数据的完整性、安全性和可扩展性。 4.设计系统架构:系统架构是ERP系统的骨架,需要考虑系统的可扩展性、可维护性和性能。 5.开发和测试:
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

探索MATLAB微分方程求解中的分岔分析:揭示方程动态行为的秘密

![matlab求解微分方程](https://img-blog.csdnimg.cn/2021062810300367.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MTQ3OTY3OA==,size_16,color_FFFFFF,t_70) # 1. MATLAB微分方程求解概述 微分方程在科学、工程和金融等领域有着广泛的应用。MATLAB作为一种强大的数值计算软件,提供了丰富的微分方程求解工具。本章将概述