24小时热门版块排行榜    

查看: 1014  |  回复: 3
【悬赏金币】回答本帖问题,作者阿毛虫将赠送您 20 个金币

阿毛虫

银虫 (正式写手)

[求助] 想请问一下,多分类逻辑回归,这个代码正确吗 已有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
回复此楼
学习吧
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿毛虫

银虫 (正式写手)

表情为:
学习吧
2楼2022-01-10 09:24:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿毛虫

银虫 (正式写手)

大神,这里是什么意思?Bn=Bo + pinv(X'*W*X)*X'*(Y-Pi);
学习吧
3楼2022-01-10 10:09:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

这里的  pinv(X'*W*X)*X'*(Y-Pi); 是求伪逆矩阵
QQ:2120156492
4楼2022-01-13 16:17:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿毛虫 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 265求调剂 +8 小木虫085600 2026-03-09 11/550 2026-03-10 19:44 by 绿色低碳化学
[考研] 欢迎211本科同学,过A区国家线,A区非偏远一本,交叉学科课题组 +23 lisimayy 2026-03-04 37/1850 2026-03-10 18:10 by a1763909
[考研] 材料工程专硕调剂自荐信,初试323 +6 李白26 2026-03-07 6/300 2026-03-10 16:16 by 无际的草原
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[考研] 307求调剂 +5 超级伊昂大王 2026-03-10 5/250 2026-03-10 14:01 by xxrvvv
[考研] 接受26届调剂生 +21 猪猪猪毛 2026-03-06 21/1050 2026-03-10 13:39 by Jack_333
[考研] 材料调剂,307分 +10 张泳铭1 2026-03-09 10/500 2026-03-10 11:47 by peike
[考研] 中科大材料299求调剂 +10 DAIjiayo 2026-03-05 16/800 2026-03-09 10:45 by 斩魂滴兔子!
[考研] 求调剂,一志愿江南大学,食品科学与工程,总分,320 +3 yyyyyukino 2026-03-07 3/150 2026-03-08 23:07 by 清风月
[考研] 求调剂 +3 拾柒12。 2026-03-08 3/150 2026-03-08 17:33 by bingxueer79
[考研] 347求调剂 +4 浮云满足 2026-03-07 4/200 2026-03-08 16:46 by 星空星月
[考研] 085600材料与化工 292分 求调剂 +6 程晴之 2026-03-06 6/300 2026-03-07 09:22 by 斩魂滴兔子!
[考研] 求调剂 +4 呼呼?~+123456 2026-03-06 4/200 2026-03-06 23:11 by L135790
[考研] 材料277分求调剂 +13 饭饭星球 2026-03-04 14/700 2026-03-06 16:10 by @飒飒飒飒
[考研] 材料与化工304求B区调剂 +4 邱gl 2026-03-06 4/200 2026-03-06 15:51 by 聪明的大松鼠
[考研] 求调剂 +5 danyyyy 2026-03-04 5/250 2026-03-06 08:16 by Leeding1356
[考研] 材料085601一志愿哈工大317 +4 压迫感行 2026-03-04 4/200 2026-03-05 20:08 by 黑衣馒头人
[考研] 求调剂 +3 泡了个椒 2026-03-04 4/200 2026-03-05 14:37 by 泡了个椒
[考研] 301求调剂 +4 李LJR 2026-03-04 4/200 2026-03-05 13:31 by Iveryant
[考研] 320材料一志愿西工大专硕英二数二 有两年光伏工作经验 +5 锤某人 2026-03-04 8/400 2026-03-05 09:24 by oxidpl
信息提示
请填处理意见