码迷,mamicode.com
首页 > 其他好文 > 详细

逻辑(logistic)回归

时间:2016-05-06 16:12:29      阅读:181      评论:0      收藏:0      [点我收藏+]

标签:

技术分享

回忆

----------------------------------------------------------------------------------------------------------------------

训练集合

技术分享


其中技术分享有n个特征,技术分享只有一个,那么假如有如下关系:

技术分享

则:技术分享

-----------------------------------------------------------------------------------------------------------------------

Logistic函数(Sigmoid函数):

技术分享


技术分享则:


技术分享


cost function:


技术分享

                     技术分享

其中:

技术分享

计算技术分享,

技术分享



处理步骤:

技术分享

其中:

技术分享

现在来证明logistic函数和线性回归函数的梯度表达是一样的:

线性回归的梯度:技术分享

logistic的梯度:

技术分享

技术分享

又因为:技术分享

那么:

技术分享

技术分享


其中z就是式中的h,所以和线性回归一样的梯度函数:

                                                     技术分享


实验:

%逻辑回归

%初始数据  
x=[-3;      -2;     -1;     0;      1;      2;     3];  
y=[0.01;    0.1;   0.3;    0.45;   0.8;    0.8;    0.99];  
plot(x,y,'ro');  
hold on  

%拟合
m = length(y);
theta = [0 0];  
a=0.005;  
loss = 1;  
iters = 1;  
eps = 0.0001;  
while loss >eps && iters <100
    loss = 0;
    for i = 1:length(y)
        h = 1./(1+exp(-(theta(1)+theta(2)*x(i))));
        theta(1)=theta(1)+a*(y(i)-h);  
        theta(2)=theta(2)+a*(y(i)-h)*x(i,1);  
        err = theta(1)+theta(2)*x(i,1)-y(i);  
        loss = loss+err*err/m; 
    end
     iters = iters+1;
end
iters  
theta 

%画图对比
for x = -3:0.01:3
     h = 1./(1+exp(-(theta(1)+theta(2)*x)));
     plot(x,h);  
end
hold off



技术分享


逻辑(logistic)回归

标签:

原文地址:http://blog.csdn.net/legotime/article/details/51312393

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!