Calcul de régression linéaire en C#

Version imprimable

28/12/2010, 10h00
olibara

Calcul de régression linéaire en C#

Bonjour

Je cherche une source CSharp permetant de faire le calcul de régression linéaire d'un ensemple de points (le but et de recuperer a et b pour evaluer la tendance d'évolution dans le temps de mesures quotidiennes)

J'avoue que mes notions de statistique et l'application des formules se sont un peu estompées et qu'a priori j'ai un peu peur d'y plonger tout seul

Donc je cherche un petit bout de code ou une recette pratique qui m'aideront plus que des formules abstraites et indigeste ou qu'une librairie hermétique !

Une idée ?
Un conseil ?

Merci de votre aide

N.B. Du code dans un autre language est aussi bon, j'aurais plus facile a traduire qu'une formule abstraite :roll:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
void Main()
{
	var points = new List<Point>()
						{
							new Point{X=1,Y=8},
							new Point{X=3,Y=1},
							new Point{X=2,Y=4}
						};
 
	var a = points.Variance(p=>p.X,p=>p.Y)/points.Variance(p=>p.X,p=>p.X);
	var b = points.Average(p=>p.Y)-a*points.Average(p=>p.X);
}
 
public static class Extension
{
	public static double Variance<T>(this IEnumerable<T> list, Func<T,double> selectA, Func<T,double> selectB)
	{
		return list.Average(p=>selectA(p)*selectB(p))-(list.Average(p=>selectA(p))*list.Average(p=>selectB(p)));
	}
}
 
public struct Point
{
	public double X {get;set;}
	public double Y {get;set;}
}

Vérifié sur Excel :)

http://www.thesmileysign.com/temp/54...3f9e226afa.gif

28/12/2010, 12h22
olibara

Merci PitMaverick78

Je teste
J'ai déja renommé Point en PointD pour lever l'ambiguité sur Point
J'ai juste un peu de mal a capter dans ton exemple comment la liste points peut recevoir la methode Variance

Et le compilateur a du mal aussi

Citation:

'System.Collections.Generic.List<AdmStat.Common.C_Regression.PointD>' does not contain a definition for 'Variance' and no extension method 'Variance' accepting a first argument of type 'System.Collections.Generic.List<AdmStat.Common.C_Regression.PointD>' could be found (are you missing a using directive or an assembly reference?)

Salut PitMaverick78

Ce ne serais pas plutot ceci ?

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
    static public void Mainu()
    {
      var points = new List<PointD>()
      {
        new PointD{X=1,Y=8},
        new PointD{X=3,Y=1},
        new PointD{X=2,Y=4}
      };
 
      var a = Variance(points, p => p.X, p => p.Y) / Variance(points, p => p.X, p => p.X);
      var b = points.Average(p => p.Y) - a * points.Average(p => p.X);
    }
 
    public static double Variance<T>(this IEnumerable<T> list, Func<T, double> selectA, Func<T, double> selectB)
    {
      return list.Average(p => selectA(p) * selectB(p)) - (list.Average(p => selectA(p)) * list.Average(p => selectB(p)));
    }
 
    public struct PointD
    {
      public double X { get; set; }
      public double Y { get; set; }
    }

28/12/2010, 13h00
Nathanael Marchand

Oui alors enfait Variance est une methode d'extension qui doit être mise dans une classe statique en effet. J'ai un peu résumé le code.