Question

我正在尝试创建一个运行Newton Raphson优化的代码。我正在使用proc iml，但是当我需要评估错误时（e）我需要总结所有的方差并且不知道如何告诉SAS在这种情况下我需要向量的组件的总和而不是矢量。代码如下：

proc iml;  use chap0; read all var{X} into X; 
            read all var{t} into t; 

W=1;           
s= exp(X*w)/(1+ exp(X*w)); print s;
e = (s - t) ** 2; /*here I need the result of the sum for that and not the matrix*/
g=2*(s-t)*s*(1-s);
h=2 * s * (1 - s) * (s * (1 - s) + (s - t) * (1 - 2 * s));

count=0;/*init count number*/


do until (e<1e-8); 

 count=count+1; 
 w=w0-g/h;   /*here I also need the sum of g and h*/
 s= exp(X*w)/(1+ exp(X*w));  
 e = (s - t) ** 2;  
 wo=w;  
end;

谢谢！

Answer 1

您应该使用SSQ函数计算IML中的平方和。

e=ssq(s-t);

以下是其他几种方法。注意第一种方式给出了不同的结果（因为它总结了s和-t的平方），这只是你传递参数的方式的一个例子。

proc iml;
 x = 1:5;
 y = j(1,5,2);
 e1=ssq(x,-y);     *sum of the squares, not actually subtracting note, so it is not the same answer;
 e2=ssq(x-y);      *sum of squares of the differences;
 e3=(x-y)*(x-y)`;  *multiplying a vector by its transpose sums it;
 e4=(x-y)[##];     *summation subscript operator, see note;
 print e1 e2 e3 e4;
quit;

Rick Wicklin关于##运算符有一个post，这非常有用。

proc iml中的平方误差之和

1 个答案: