Сумма квадратичных ошибок в proc iml

Question

Я пытаюсь создать код для запуска оптимизации Ньютона Рафсона. Я использую 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;

У Рика Уиклина есть пост об операторе ##, который весьма полезен.