3.1 Logistic回归原理,参数推导及Python实现

在这里插入图片描述

import matplotlib.pyplot as plt
import numpy as np
#假设
def logistic(x):
    return 1/(1+np.exp(-x))
#随机梯度下降
def SGD(X,Y,a=0.001):
    theta = np.mat([0, 0], dtype='float64')
    for time in range(3000):
        #用m个样本分别对theta进行更新
        for i in range(30):
            theta += a*( ( Y[i]-logistic(theta*X[i].T) ) * X[i] )
    return theta
#生成30个x
X=np.mat([ [1,x] for x in np.linspace(-5,5,30)])
Y=np.mat(np.append(np.zeros(10),np.ones(20))).T
theta=SGD(X,Y)
# 画出拟合图像
plt.scatter(np.asarray(X[:, 1]), np.asarray(Y))
plt.plot(X[:,1],logistic(theta*X.T).T)
plt.show()
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值