| 查看: 1133 | 回复: 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 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有46人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有16人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有4人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
阿毛虫
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 549.7
- 散金: 81
- 帖子: 757
- 在线: 117.2小时
- 虫号: 3313820
- 注册: 2014-07-08
- 性别: GG
- 专业: 机械动力学
2楼2022-01-10 09:24:28
阿毛虫
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 549.7
- 散金: 81
- 帖子: 757
- 在线: 117.2小时
- 虫号: 3313820
- 注册: 2014-07-08
- 性别: GG
- 专业: 机械动力学
3楼2022-01-10 10:09:28
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17852.5
- 红花: 53
- 帖子: 2880
- 在线: 606.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
4楼2022-01-13 16:17:28
