1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
| START F_td (P) GLOBAL ( y1, y );
%macro essai ;
%do i=1 %to 29 ;
F&i. = (y[&i.]- probnorm((y1-sqrt(1*p[30])*p[&i.])/sqrt(1-p[30])))##2 ;
%end ;
F=F1+F2+F3+F4+F5+F6+F7+F8+F9+F10+F11+F12+F13+F14+F15+F16+F17+F18+F19+F20+F21+F22+F23+F24+F25+F26+F27+F28+F29 ;
%mend ;
%essai ;
RETURN (F);
FINISH F_td ;
START c_td (P);
c= 1-(1/29)*(p[1]*p[1]+p[2]*p[2]+p[3]*p[3]+p[4]*p[4]+p[5]*p[5]+p[6]*p[6]+p[7]*p[7]+p[8]*p[8]+p[9]*p[9]+
p[10]*p[10]+p[11]*p[11]+p[12]*p[12]+p[13]*p[13]+p[14]*p[14]+p[15]*p[15]+p[16]*p[16]+p[17]*p[17]+p[18]*p[18]+p[19]*p[19]+
p[20]*p[20]+p[21]*p[21]+p[22]*p[22]+p[23]*p[23]+p[24]*p[24]+p[25]*p[25]+p[26]*p[26]+p[27]*p[27]+p[28]*p[28]+p[29]*p[29]) ;
RETURN (c);
FINISH c_td ;
P = {0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.1};
optn= j(1,30,.);OPTN[2]=3;
CALL NLPnms (RC,"F_td", P,OPTN) nlc="c_td" ;
quit ; |
Partager