在线性回归模型中使用梯度下降法

上一篇文章记录了梯度下降法模拟实现,下面的代码展示的是在线性回归中使用梯度下降法

import numpy as np
import maltplotlib.pyplot as plt
np.random.seed(666)
x=2*np.random.random(size=100)
y=x*3.+4.+np.random.normal(size=100)
plt.scatter(x,y)
plt.show()

在这里插入图片描述

def J(theta,X_b,y):
	try:
		return np.sum((y-X_b.dot(theta))**2)/len(X_b)
	except:
		return float('inf')
def dJ(theta,X_b,y):
	res=np.empty(len(theta))
	res[0]=np.sum(X_b.dot(theta)-y)
	for i in range(1,len(theta)):
		res[i]=(X_b.dot(theta)-y).dot(X_b[:,i])
	return res*2/len(X_b)
def gradient_descent(X_b,y,initial_theta,eta,n_iters=1e4,epsilon=1e-8)
	theta=initial_theta
	i_iter=0
	while i_iter<n_iters:
		gradient=dJ(theta,X_b,y)
		last_theta=theta
		theta=theta-eta*gradient
		if(abs(J(theta,X_b,y)-J(last_theta,X_b,y))<epsilon):
			break
		i_iter+=1
	return theta

X_b=np.hstack([np.ones((len(x),1)),x.reshape(-1,1)])
initial_theta=np.zeros(X_b.shape[1])
eta=0.01
theta=graident_descent(X_b,y,initial_theta,eta)
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值