Accueil | Cast3M

   1 : $$$$ AGRE     NOTICE  FD218221  25/03/12    21:15:01     12191          
   2 :                                              DATE     25/03/12
   3 : 
   4 :   Operateur AGRE                       Voir aussi : @STAT SOMM DIME
   5 :       --------------  
   6 :       RES = AGRE MOT1 LIST1 (P) ('ROBU') ;
   7 :                                                                          
   8 :     
   9 :       Objet :                                                            
  10 :       ______
  11 :                                                                          
  12 :       Cet operateur calcule une fonction d'agregation d'une liste de valeurs.
  13 : 
  14 :       Commentaires :
  15 :       _____________
  16 :                                                                          
  17 :       LIST1  : Liste de reels (type LISTREEL).
  18 : 
  19 :       MOT1   : Mot clef indiquant la fonction a appliquer
  20 :                a choisir dans la liste suivante :
  21 : 
  22 :       Statistiques, normes, moyennes
  23 :         'SOMM' : Somme
  24 :                  = Sum[i=1;n](x_i)
  25 :         'PROD' : Produit
  26 :                  = Pro[i=1;n](x_i)
  27 :         'MOYE' : Moyenne arithmetique (moment d'ordre 1)
  28 :                  = SOMM/n
  29 :         'MOHA' : Moyenne harmonique
  30 :                  = n/(Sum[i=1;n](1/x_i))
  31 :         'MOGE' : Moyenne geometrique
  32 :                  = PROD**(1/n)
  33 :         'VARI' : Variance (moment centre d'ordre 2)
  34 :                  = (1/n) * Sum[i=1;n]((x_i-MOYE)**2)
  35 :         'ECTY' : Ecart type
  36 :                  = VARI**0.5
  37 :         'ASYM' : Coefficient d'asymetrie (moment centre reduit d'ordre 3)
  38 :                  = (1/n) * Sum[i=1;n](((x_i-MOYE)/ECTY)**3)
  39 :         'KURT' : Kurtosis (moment centre reduit d'ordre 4)
  40 :                  = (1/n) * Sum[i=1;n](((x_i-MOYE)/ECTY)**4)
  41 :         'MEDI' : Mediane
  42 : 
  43 :       Normes/moyennes/maximums generalises (de parametre P)
  44 :       (le mot clef 'ROBU' est disponible pour ces mots clefs)
  45 :         'PMOM' : Moment d'ordre P
  46 :                  = Sum[i=1;n](x_i**P)
  47 :         'PMOY' : Moyenne generalisee d'ordre P
  48 :                  = (PMOM/n)**(1/P)
  49 :         'PNOR' : Norme generalisee d'ordre P
  50 :                  = (Sum[i=1;n](|x_i|**P))**(1/P)
  51 :         'LEHM' : Fonction de Lehmer d'ordre P
  52 :                  = (Sum[i=1;n](x_i**P) / (Sum[i=1;n](x_i**(P-1))))
  53 :         'KSL'  : Fonction de Kreisselmeir Steinhauser inferieure d'ordre P (ou MellowMax)
  54 :                  = (1/P) * ln((1/n) * Sum[i=1;n](exp(P*x_i)))
  55 :         'KSU'  : Fonction de Kreisselmeir Steinhauser superieure d'ordre P (ou LogSumExp)
  56 :                  = (1/P) * ln(Sum[i=1;n](exp(P*x_i)))
  57 :         'BOLT' : Fonction de Boltzmann d'ordre P
  58 :                  = (Sum[i=1;n](x_i*exp(P*x_i))) / (Sum[i=1;n](exp(P*x_i)))
  59 : 
  60 :       P      : Parametre d'exponentiation pour certaines fonctions (type FLOTTANT).
  61 : 
  62 :       'ROBU' : Mot clef, optionnel, pour utiliser une version "robuste" des fonctions parametrees.
  63 :                Ceci est utile lorsque les valeurs de LIST1, ou de P, sont grandes et conduisent
  64 :                a des debordements (overflow) lors de l'application des fonctions puissance ou exponentielle.
  65 :                On procede en normalisant les valeurs x avec la norme infinie VINF ou bien le maximum VMAX :
  66 :                - P moyenne robuste :  PMOY(x,p) = VINF * PMOY(x/VINF,p)
  67 :                - P norme   robuste :  PNOR(x,p) = VINF * PNOR(x/VINF,p)
  68 :                - P LEHM    robuste :  LEHM(x,p) = VINF * LEHM(x/VINF,p)
  69 :                - KSL       robuste :  KSL(x,p)  = VMAX + KSL(x-VMAX,p)
  70 :                - KSU       robuste :  KSU(x,p)  = VMAX + KSU(x-VMAX,p)
  71 :                - BOLT      robuste :  BOLT(x,p) = VMAX + BOLT(x-VMAX,p)
  72 : 
  73 :       RES    : Resultat (type FLOTTANT).
  74 : 
  75 :       Remarques :
  76 :       __________
  77 : 
  78 :       1) La fonction PMOY (P moyenne) est prise egale a la moyenne geometrique MOGE si P=0.
  79 : 
  80 :       2) Certaines fonctions d'agregation necessittent de calculer de tres grands nombres,
  81 :          une prudence s'impose donc sur le choix de l'exposant P lorsque les valeurs de LIST1
  82 :          sont elevees afin d'eviter les debordements de nombre flottant.
  83 :