有谁用SQL存储过程写过Logisitic回归的算法吗

shiyiming

如题，希望能多多请教，谢谢。

shiyiming

应该很难，反而在EXCEL中或者有更多相关的函数可以用。
以前在图书馆看到过在SQL SEVER上做建模
兄弟也做logistic regression? 多多交流

shiyiming

这里采用一篇文章中的移栽后不同天数的生长数据。
根据k=(y2^2(y1+y3)-2y1y2y3)/(y2^2-y1y3)公式等距计算k值，调用
data m;                                                                                                                                 
input x y@@;                                                                                                                           
yp=log((4.59-y)/y);                                                                                                                     
cards;                                                                                                                                 
30      0.06                                                                                                                           
45      0.71                                                                                                                           
60      3.31                                                                                                                           
75      4.47                                                                                                                           
90      4.58                                                                                                                           
;                                                                                                                                       
proc reg;                                                                                                                              
model yp=x;                                                                                                                             
run;                                                                                                                                    
可计算出a、b值，再调用下面语句进行曲线模拟：
data n;                                                                                                                                 
set m;                                                                                                                                 
yp=4.59/(1+14083.51*exp(-0.17479*x));                                                                                                   
proc gplot data=n;                                                                                                                     
symbol1 i=join v=dot;                                                                                                                  
symbol2 i=join v=none co=red;                                                                                                           
plot y*x=1 yp*x=2/overlay;                                                                                                              
run;
从拟合图可以看出拟合效果非常好！