ot_statistique.cpp

Aller à la documentation de ce fichier.

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       ot_statistique.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:54:00 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
#pragma hdrstop
 
#include "ot_statistique.h"
#pragma package(smart_init)
 
 
double OT_STATISTIQUE::moyenne (double * x, int n)
{
    int i;
    double moyenne_x = 0;
    for (i=0;i<n;i++)
        moyenne_x += x[i];
    moyenne_x /= n;
    return moyenne_x;
}
 
double OT_STATISTIQUE::ecart_type (double * x, int n)
{
    int i;
    double moyenne_x = moyenne(x,n);
 
    double moyenne_x2 = 0;
    for (i=0;i<n;i++)
        moyenne_x2 += x[i]*x[i];
    moyenne_x2 /= n;
 
    return sqrt(moyenne_x2-moyenne_x*moyenne_x);
}
 
double OT_STATISTIQUE::covariance (double * x, double * y, int n)
{
    int i;
 
    double moyenne_xy = 0;
    for (i=0;i<n;i++)
        moyenne_xy += x[i]*y[i];
    moyenne_xy /= n;
 
    double moyenne_x = 0;
    for (i=0;i<n;i++)
        moyenne_x += x[i];
    moyenne_x /= n;
 
    double moyenne_y = 0;
    for (i=0;i<n;i++)
        moyenne_y += y[i];
    moyenne_y /= n;
 
    return moyenne_xy - moyenne_x * moyenne_y;
}
 
double OT_STATISTIQUE::variance (double * x, int n)
{
    int i;
 
    double moyenne_carres = 0;
    for (i=0;i<n;i++)
        moyenne_carres += x[i]*x[i];
    moyenne_carres /= n;
 
    double moyenne = 0;
    for (i=0;i<n;i++)
        moyenne += x[i];
    moyenne /= n;
 
    return moyenne_carres - moyenne*moyenne;
}
 
void OT_STATISTIQUE::regression_lineaire(double *x, double *y, int n, double *pente, double *valY0)
{
    int i;
 
    // coef directeur de la droite
    double a = variance(y,n) / covariance(x,y,n);
 
    // b = valeur de y (0) par moindre carr�s
    double b = 0;
    for (i=0; i<n; i++)
        b += y[i] - a*x[i];
    b /= n;
 
    *pente = a;
    *valY0 = b;
}
 
double OT_STATISTIQUE::max(double *x, int n)
{
    if (n==0)
        return 0;
 
    double max_x = x[0];
    for (int i=0;i<n;i++)
        if (x[i] > max_x)
            max_x = x[i];
 
    return max_x;
}
 
double OT_STATISTIQUE::min(double *x, int n)
{
    if (n==0)
        return 0;
 
    double min_x = x[0];
    for (int i=0;i<n;i++)
        if (x[i] < min_x)
            min_x = x[i];
 
    return min_x;
}