偏最小二乘logistic回归

xcsyb

偏最小二乘logistic回归

哪位有偏最小二乘logistic回归的程序？
想用它来做分类
可是第一步求logistic的系数就不知道怎么求了
郁闷

现在只编写了一般的线性偏最小二乘程序

哪位有的提供一下或给个思路
谢谢！！！

风花雪月

1. % PLS 单因变量偏最小二乘回归
   
2. clear
   
3. % 原始数据：自变量 X, 因变量 Y
   
4. load XYpls
   
5. % 标准化 X--E0 , Y--F0
   
6. [E0,meanX,stdX]=prestd(X') ; E0=E0'
   
7. [F0,meanY,stdY]=prestd(Y') ; F0=F0'
   
8. 
   
9. E=E0; F=F0;
   
10. beta=0 ; Fg=0;
    
11. for ii=1:3
    
12. 
    
13. EF=[E,F] ; [n,p]=size(E)
    
14. EFcov=cov(EF) ; EFcov=EFcov(1:end-1,:) %协方差矩阵
    
15. w=EFcov(:,p+1)/sqrt(sum(EFcov(:,p+1).^2))
    
16. t=E*w
    
17. q=E0'*t/(t'*t)
    
18. r=F0'*t/(t'*t)
    
19. ww=pinv(E0)*t % w*
    
20. 
    
21. beta=r*ww+beta % 回归系数β
    
22. Fg=t*r+Fg % 预测值
    
23. 
    
24. E=E-t*q' ; F=F-t*r % 残差矩阵
    
25. end
    
26. 
    
27. Yg=poststd(Fg,meanY,stdY) % 预测值（还原值）
    
28. 
    
29. % 模型检验
    
30. % 复测定系数 R^2
    
31. SSE=sum((Y-Yg).^2)
    
32. SST=sum((Y-ones(n,1)*meanY).^2)
    
33. R2 = 1-SSE/SST % 复测定系数 R^2
    
34. R2tz = 1-(1-R2)*(n-1)/(n-p-1) % 调整复测定系数 R^2tz
    
35. % F 检验
    
36. F = ((SST-SSE)/p)/(SSE/(n-p-1))

复制代码

来自：研学

aika1985

您好，这段程序中对初始的数据是不是有什么要求，我做的输入是800*4的矩阵，输出是800*5的矩阵，运行后不行，t=E*w这步不对呢，w是一个8*1的矩阵，E是4*800的矩阵，对程序的算法也不十分明白，能给解释一下不？:loveliness:
敬待回信