24小时热门版块排行榜

>论坛更新日志 (3017)
>虫友互识 (233)
>文献求助 (159)
>导师招生 (63)
>休闲灌水 (47)
>考博 (43)
>招聘信息布告栏 (38)
>博后之家 (32)
>硕博家园 (22)
>晶体 (18)
>论文投稿 (18)
>论文道贺祈福 (16)
>公派出国 (13)
>考研 (12)
>教师之家 (11)
>有机资源 (10)

返回列表

【悬赏金币】回答本帖问题，作者阿毛虫将赠送您 20 个金币

阿毛虫

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 90.7
散金: 81
帖子: 748
在线: 112.7小时
虫号: 3313820
注册: 2014-07-08
性别: GG
专业: 机械动力学

[求助] 想请问一下，多分类逻辑回归，这个代码正确吗已有1人参与

%% Training

Betas=MLR(Xtrain,Ytrain);

%% Testing
Ypred_test=getPredMLR(Betas,Xtest);
error_test=getError(Ypred_test,Ytest)
Ypred_train=getPredMLR(Betas,Xtrain);
error_train=getError(Ypred_train,Ytrain)

function [ B ] = MLR( X,Y )
Y_unique=unique(Y);
X=[ones(size(X,1),1),X];
d1=length(Y_unique);
d2=size(X,2);
B=zeros(d1,d2);
% Y_unique=unique(Y);
for i=1:1:d1
Y_blr=getY_blr(Y,Y_unique(i));
B(i,

=getB(B(i,

,X,Y_blr);
end
end

function Bn = getB( X,Y )
Bo=zeros(size(X,2),1);%迭代初始值
for i=1:8
W=getW(X,Bo);
Pi=getPi(Bo,X);
Bn=Bo + pinv(X'*W*X)*X'*(Y-Pi);
Bo=Bn;
end
Bn=Bo;
end

function Pi = getPi(Bo,X)
Pi=[];
S=size(X,1);
% W=zeros(S(1));
for i=1:1:S(1)
P=(exp(Bo'*X(i,

'));
P=P/(1+P);
Pi=[Pi,P];
end
Pi=Pi';
end

function W = getW( X,Bo )
S=size(X);
W=zeros(S(1));
for i=1:1:S(1)
num=(exp(Bo'*X(i,

'));
denom=1+num;
P=num/denom;%sigmond函数
RP=1-P;
% size(W);
W(i,i)=P*RP;
end
end

function [ Y_blr ] = getY_blr(Y,i)
Y_blr=Y;
Y_blr(Y==i)=1;
Y_blr(Y~=i)=0;
end

回复此楼