Régression d'ordre 2 avec Apache Commons Math

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:

---

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:

---

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

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:

---

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:

---

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