查看: 247  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 给考研学子的忠告 +4 HiRabit 2022-05-14 8/400 2022-05-18 13:47 by 欢欢鱼yhh
[考博] 申请中山大学博士有多难 +9 cogito7 2022-05-17 9/450 2022-05-18 13:47 by DODOLUCK
[硕博家园] 替羞涩的西安妹妹找靠谱老公(人在西安的山西女子) +5 飞雪飘云 2022-05-17 7/350 2022-05-18 12:16 by wx24002381
[公派出国] 有去美国的小伙伴吗?你们拿了几个国家的offer ,一起交流 +3 lunl115 2022-05-13 6/300 2022-05-18 11:54 by jiachui
[高分子] 有机硅树脂 +4 金金金迷茫 2022-05-14 9/450 2022-05-18 11:16 by 金金金迷茫
[基金申请] E04还没有送审 +6 xmccnmmmm 2022-05-12 7/350 2022-05-18 09:47 by 木命克火
[论文投稿] CEJ投稿交流 +5 eric1130 2022-05-17 6/300 2022-05-18 01:17 by mxnjfu
[硕博家园] 该怎么选 +7 1131113644 2022-05-14 9/450 2022-05-18 00:55 by 为了科研丶
[论文投稿] 本科生毕设论文投稿 +3 Melody_Li 2022-05-17 3/150 2022-05-17 19:38 by georose
[硕博家园] 黑眼圈大怎么办??? +7 倩女离魂 2022-05-11 7/350 2022-05-17 11:51 by 杶白
[基金申请] 今年的医学口青基涵评资助指标好像只有1/3不到,太卷了! +6 535708218 2022-05-16 7/350 2022-05-17 05:47 by liuqingv
[考研] 寻找研友 +12 Seaforest 2022-05-13 20/1000 2022-05-16 16:36 by 常州大学666
[硕博家园] 经费少 是做生信还是自学ai还是写综述 +7 将心 2022-05-11 9/450 2022-05-16 14:48 by university
[论文投稿] 想请教一下各位实时影响因子一般准确吗 +3 拉架隆多 2022-05-15 8/400 2022-05-15 23:01 by 拉架隆多
[基金申请] 还有打招呼不让支持某某人的 +42 babu2015 2022-05-13 48/2400 2022-05-15 18:19 by lee86208
[高分子] 不饱和树脂固化问题 +3 金金金迷茫 2022-05-11 7/350 2022-05-14 18:49 by 金金金迷茫
[基金申请] 基金参与者收不到当年评审本子吗? +14 小谦谦2008 2022-05-11 14/700 2022-05-14 09:39 by wang1325
[论文投稿] 投MDPI期刊被编辑拒掉,是不是意味着文章太次了 +19 对不齐呀 2022-05-11 20/1000 2022-05-14 08:21 by 重大粉
[基金申请] 请问数理学部力学口的青年项目有几位评审专家? +6 DS_Audio 2022-05-13 7/350 2022-05-13 20:59 by 黑色甜心123
[基金申请] 收到“i科研人”的电子邮件,推荐一批科研课题。 +4 sunshineSJ 2022-05-11 8/400 2022-05-12 09:17 by sunshineSJ
信息提示
请填处理意见