梯度下降法实现Logistic回归

参考简书中Logistic回归及Python代码实现

Logistic函数的损失函数的偏导数为$\frac{1}{m}\sum_{i=1}^{m}(h_\theta(x_i)-y_i)x_i^j$,所以$\theta$的更新可以写为:$\theta_j=\theta_j-\alpha\frac{1}{m}\sum_{i=1}^{m}(h_\theta(x_i)-y_i)x_i^j$,推导过程中使用了最大似然估计,具体过程见原文

推导过程如下:

1)Logistic函数如下

$$g(z)=\frac{1}{1+e^{-z}}$$

其中

$$z={\theta }_{0}{x }_{0}+{\theta }_{1}{x }_{1}+...+{\theta }_{n}{x }_{n}$$

向量表示为

$$Z={\Theta}^{T}X$$

2)目标函数变形

预测函数为

$$h_\theta(x)=g(\theta^Tx)=\frac{1}{1+e^{-\theta^Tx}}$$

函数$h_\theta(x)$的值有特殊的含义,它表示结果取1的概率,因此,对于输入x分类结果为类别1和类别0的概率分别为:

$$P(y=1|x;\theta)=h_\theta(x)\\ P(y=0|x;\theta)=1-h_\theta(x)$$

将以上两个式子合起来写为
$$P(y|x;\theta)=(h_\theta(x))^y(1-h_\theta(x))^{1-y}$$

3)取似然函数

$$L(\theta)=\prod_{i=1}^mP(y_i|x_i;\theta)=\prod_{i=1}^m(h_\theta(x_i))^{y_i}(1-h_\theta(x_i))^{1-y_i}$$

这是一组观测值的概率,取对数为

$$l(\theta)=\log L(\theta)=\sum_{i=1}^m(y_i\log h_\theta(x_i)+(1-y_i)\log (1-h\theta(x_i)))$$

最大似然估计就是求使$l(\theta)$最大时,$\theta$的值。

 

猜你喜欢

转载自www.cnblogs.com/guesswhy/p/11285753.html