Discussion :

[mydomdom]

Bonjour,

Je cherche à faire une régression polynomiale d'ordre 2 avec la bibliothèque Apache.Commons.Math.
A titre d'exemple, le code suivant (adapté de http://commons.apache.org/proper/com...e/fitting.html) devrait logiquement aboutir à cette équation:
y = 39.79 x² - 497.66 x + 997.45

Code :

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

1
2
3
4
5
6
7
8
double[] y = { 540.0, 160.0, -140.0, -360.0, -480.0, -560.0, -540.0, -440.0, -260.0, 0.0, 340.0};              
        final WeightedObservedPoints obs = new WeightedObservedPoints();
        for (double figure:y){
            obs.add(1.0, figure);// 1.0 = le poids de chaque donnee. La doc precise de mettre 1.0 si on ne dispose pas de l'info
        }
        final PolynomialCurveFitter fitter = PolynomialCurveFitter.create(2);
        final double[] coeff = fitter.fit(obs.toList());
        System.out.println("coef="+Arrays.toString(coeff));

Or, il n'en est rien puisqu'il me donne les coefficients suivants:

Code :

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

coef=[-53.73522460839947, -52.22329678670934, -52.22329678670934]

J'apprécierais toute aide me permettant d'utiliser correctement cette belle bibliothèque ...
Merci d'avance
Dominique

[mydomdom]

J'ai la réponse. Il ne faut pas oublier de mettre les x ... Gloups...
Le code devient si on part des x avec 1 et qu'on incrémente de 1 à chaque fois:

Code :

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

1
2
3
4
5
6
7
8
9
10
double[] y = { 540.0, 160.0, -140.0, -360.0, -480.0, -560.0, -540.0, -440.0, -260.0, 0.0, 340.0};              
        final WeightedObservedPoints obs = new WeightedObservedPoints();
        double m = 1.0;
        	for (double fig:y4){
        		obs.add(m, fig);
        		m++;
        	}
        final PolynomialCurveFitter fitter = PolynomialCurveFitter.create(2);
        final double[] coeff = fitter.fit(obs.toList());
        System.out.println("coef="+Arrays.toString(coeff));

La sortie devient:

Code :

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

1
2
coef=[997.4545454545454, -497.6643356643356, 39.790209790209786]
//constante, coef 1er degre, coef second degre

soit: y = 39.79 x² - 497.66 x + 997.45
C'est tout bon .
Si ça peut aider un jour quelqu'un...
Dominique