logistic 回归Matlab代码

function a
alpha = 0.0001;
[m,n] = size(q1x);

max_iters = 500;

X = [ones(size(q1x,1),1), q1x]; % append col of ones for intercept term

theta = zeros(n+1, 1);  % initialize theta
for k = 1:max_iters
  hx = 1./(1+exp(-X*theta));
  theta = theta + alpha * X' * (q1y-hx); 
  ll(k) = sum(q1y .* log(hx) + (1 - q1y) .* log(1 - hx) );
  hold on;
end
plot(q1x(1:50,1),q1x(1:50,2),'o');
hold on;
plot(q1x(51:99,1),q1x(51:99,2),'r*');
hold on;
x1=0:0.01:8;
y=(1-theta(1,1))/theta(3,1)-x1*theta(2,1)/theta(3,1);
plot(x1,y)
q1x.dat,q1y.dat: http://cs229.stanford.edu/
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值