bon... je me mets à IML....

pour moi, votre problème ne vient pas d'IML...

faisons ensemble un petit test :

data toto;
   do i=1 to 25;
     test=ranuni(-584)/1E12;
     verif=log(test);
     output;
   end;
run;

proc sql noprint ;
   select test,verif into :list separated by ' ',:listl separated by ' ' from toto; 
quit;

proc iml;
   c = {&list};
   b = log(c);
   a = {&listl};
   print b a;
quit;

le calcul des log ne posent pas de problème...

en passant, le passage par une macro variable construite avec PROC SQL est très certainement la chose à pire à faire...

Des modalités de variables numériques, lorsqu'elle sont placées dans une macro variable sont converties au moyen d'un FORMAT BEST8.

Les valeurs que j'envoie donc dans le vecteur a ont donc la forme suivante : 3.91E-13 8.07E-13 9.45E-13 4.17E-13 1.9E-13...

Pour la précision, ce n'est pas gagné... c'est ce que vous montre le contenu du vecteur a (qui lui, contient les log calculés dans l'étape DATA) - mais, franchement parlé, je m'attendais à pire... les différences n'apparaissent qu'à la quatrième décimale...

en conclusion : le calcul de log avec IML, même sur des valeurs très petites ne posent pas de problème.

donc, c'est votre programme...

esa_gc a écrit:

Je ne suis pas un spécialiste, mais il me semble que SAS n'aime pas faire plusieurs choses à la fois

comme les calculatrices quatre opérations qui ne respectent pas les priorités multiplication / addition ?

c'est ça ?

Je me dis que c'est peut être au sein de la routine d'optimisation Newton Raphson

euh... c'est la faute de Newton Raphson ?
les "arrondis" de Newton Raphson ?

allons bon...

quand le sage montre la lune, le fou regarde le doigt ;-)

je demeure convaincu que le problème n'est pas dans le log (le doigt) mais dans votre programme (la lune)

perso, je n'y connais rien à votre routine mais quand je regarde l'aide, j'ai du mal à voir les relations entre vos paramètres et ceux présentés dans l'aide...

1) parmin intervient à un quelconque moment dans ton algorithme d'optimisation? (car on ne voit pas le début de ton programme)

2) dans opt il n'y a pas de virgules mais des espaces, peut être que ça fonctionne avec des virgules, je ne suis pas sûre mais avec des espaces j'en suis sûre :p

bon lui son chiffre est négatif donc normal que son log ne fonctionne pas --' au final ça ne m'avance pas plus ^^

Donc Rick existe vraiment!! Je pensais à un avatar traffiqué pour avoir des cheveux aussi brillant

Je reviens à la charge avec mon problème !

    
 lien       1-lien
0.9970211 0.0029789 
0.7590189 0.2409811 
0.3704836 0.6295164 
0.3580046 0.6419954 
0.3414203 0.6585797 
0.3408354 0.6591646

Mon vecteur lien est trié par ordre décroissant et le 1-lien par ordre croissant. Ce n'est pas 0 ! Le log de 1-lien ne fonctionne pas quand je lance ma routine de maximisation alors que mon programme fonctionne en dehors de celle-ci
J'ai le même problème avec mon programme ou avec celui de Mr alaa-eddine qui a rédiger une macro (bien plus complète que la mienne) sur le sujet qui me concerne http://master-esa.com/viewtopic.php?id=54 ça me redonne espoir en mon programme qui ne doit pas être si mauvais !