Discussion :

[lauac]

bonjour à tous ,
j'ai une question concernant la minimisation sous contrainte sous sas iml
j'essaie de faire une minimisation sous contrainte que la variance de certains param (..tous sauf 1) soit égale à 1
j'ai élaboré le code suivant mais ça marche pas.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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 ;

quelqu'un peut m'aider svp ? merci

[Haache]

Bonjour
Ce serait difficile de répondre directement. ça ne marche donc ça affiche quoi?
Il dit qu'il y a erreur dans le code ? Si oui quelle ligne. Essayer de vérifier le journal