generated from LufsX/Experiment-LogReg
This commit is contained in:
+20
-19
@@ -33,28 +33,13 @@ def sigmoid(z):
|
|||||||
return 1 / (1 + np.exp(-z))
|
return 1 / (1 + np.exp(-z))
|
||||||
|
|
||||||
def train_logreg(X, y, learning_rate=0.01, iterations=1000):
|
def train_logreg(X, y, learning_rate=0.01, iterations=1000):
|
||||||
"""
|
"""对率回归训练函数。"""
|
||||||
对率回归训练函数。
|
|
||||||
|
|
||||||
Args:
|
|
||||||
- X (numpy array): 特征数据。
|
|
||||||
- y (numpy array): 标签数据。
|
|
||||||
- learning_rate (float): 学习率,默认值为0.01。
|
|
||||||
- iterations (int): 迭代次数,默认值为1000。
|
|
||||||
|
|
||||||
Returns:
|
|
||||||
- weights (numpy array): 训练得到的权重。
|
|
||||||
"""
|
|
||||||
# 初始化权重
|
|
||||||
weights = np.zeros(X.shape[1])
|
weights = np.zeros(X.shape[1])
|
||||||
|
|
||||||
# 训练模型
|
|
||||||
for _ in range(iterations):
|
for _ in range(iterations):
|
||||||
z = np.dot(X, weights)
|
z = np.dot(X, weights)
|
||||||
predictions = sigmoid(z)
|
predictions = sigmoid(z)
|
||||||
# 计算梯度
|
|
||||||
gradient = np.dot(X.T, (predictions - y)) / y.size
|
gradient = np.dot(X.T, (predictions - y)) / y.size
|
||||||
# 更新权重
|
|
||||||
weights -= learning_rate * gradient
|
weights -= learning_rate * gradient
|
||||||
|
|
||||||
return weights
|
return weights
|
||||||
@@ -89,18 +74,34 @@ def load_data(filepath):
|
|||||||
|
|
||||||
return data_array[:, :-1], data_array[:, -1]
|
return data_array[:, :-1], data_array[:, -1]
|
||||||
|
|
||||||
# 主逻辑
|
def preprocess_data(X, y):
|
||||||
|
"""数据预处理,合并特征23的值,并删除特征3和特征24~28。"""
|
||||||
|
# 合并特征23的值:将2(死亡)和3(安乐死)合并为0(死亡)
|
||||||
|
y = np.where(y == 2, 0, y) # 将值2替换为0
|
||||||
|
y = np.where(y == 3, 0, y) # 将值3替换为0
|
||||||
|
|
||||||
|
# 删除特征3和特征24~28(假设特征索引从0开始)
|
||||||
|
X = np.delete(X, [2] + list(range(23, 28)), axis=1)
|
||||||
|
|
||||||
|
return X, y
|
||||||
|
|
||||||
if __name__ == "__main__":
|
if __name__ == "__main__":
|
||||||
# 加载数据
|
# 加载数据
|
||||||
X, y = load_data("./data/horse-colic.data")
|
X, y = load_data("./data/horse-colic.data")
|
||||||
|
|
||||||
|
# 数据预处理
|
||||||
|
X, y = preprocess_data(X, y)
|
||||||
|
|
||||||
# 训练模型
|
# 训练模型
|
||||||
weights = train_logreg(X, y)
|
weights = train_logreg(X, y)
|
||||||
|
|
||||||
# 加载测试数据并预测
|
# 加载测试数据并进行预处理
|
||||||
X_test, y_test = load_data("./data/horse-colic.test")
|
X_test, y_test = load_data("./data/horse-colic.test")
|
||||||
|
X_test, y_test = preprocess_data(X_test, y_test)
|
||||||
|
|
||||||
|
# 预测
|
||||||
predictions = predict_logreg(X_test, weights)
|
predictions = predict_logreg(X_test, weights)
|
||||||
|
|
||||||
# 打印预测结果和实际结果
|
# 打印预测结果和实际结果
|
||||||
for i in range(min(10, len(y_test))): # 打印前10个结果或全部结果,取决于测试数据量
|
for i in range(min(10, len(y_test))):
|
||||||
print(f"预测结果: {predictions[i]}, 实际结果: {y_test[i]}")
|
print(f"预测结果: {predictions[i]}, 实际结果: {y_test[i]}")
|
||||||
|
|||||||
Reference in New Issue
Block a user