C# linier regression

Det er vel bare ren matematik ...

Ja det har jeg overvejet.
Jeg kunne jo slå op i min matematikbog
men tænkte jeg kunne slippe lidt nemmere.

Michael

ja okay ...

Det kan jeg godt følge dig i ... skulle bare lige være sikker på hvad du egentlig sprugte om da jeg kom i tvivl ...

Men det er så enten et .NET library eller den matematiske formel som vil virke i .NET ....

efter lidt tid på google ...
Måske ham gutten her har det med i hans Library.
http://community.opensourcedotnet.info/files/Default.aspx?Sort=PostDate

// ouT

Ja måske - ser på det :-)

Jeg giver 200 point, hvis du kan finde ud af om man kan og hvordan.....

Det var ikke i hans Library så ?

Det er pærenemt !

Demo:

using System;
using System.Collections.Generic;

using MathNet.Numerics.LinearAlgebra;

namespace E
{
    public class Program
    {
        private static void Solve(double[] y, out double[] b, double[] x)
        {
            Vector one = Vector.Ones(x.Length);
            Vector x2 = Vector.Create(x);
            Matrix x3 = Matrix.CreateFromColumns(new List<Vector> { one, x2 });
            Vector y2 = Vector.Create(y);
            Matrix y3 = Matrix.CreateFromColumns(new List<Vector> { y2 });
            Matrix b3 = x3.Solve(y3);
            Vector b2 = b3.GetColumnVector(0);
            b = b2.CopyToArray();
        }
        private static readonly Random rng = new Random();
        private static double Eps(double spread)
        {
            return (rng.NextDouble() - 0.5) * spread;   
        }
        public static void Main(string[] args)
        {
            double[] x = { 1, 2, 3, 4, 5 };
            double[] y = new double[x.Length];
            double[] b;
            for(int i = 0; i < y.Length; i++)
            {
                y[i] = 0.456 + 0.123 * x[i] + Eps(0.05);
            }
            Solve(y, out b, x);
            Console.WriteLine("y = " + b[0] + " + " + b[1] + " * x");
            Console.ReadKey();
        }
    }
}

Så du har brugt:
http://mathnet.opensourcedotnet.info/doc/Iridium.ashx

?

Yes.

Der er en linear algrebra del i den.

Og dens Solve bruger ikke simpel b=inv(x'x)x'y men rigtig LU eller QR dekomponering.

Det lyder meget fint arne_V
Nu mangler jeg bare lige og hitte ud af, hvordan jeg ændre point til 200.

Er der nogen der vil have point.
Så må I sende svar:-)
Michael

svar

Skulle buzz ikke have haft point for at finde library'et ?

Ved I hvad, jeg er også lidt flov!!!
Selvfølgelig skulle buzz det!

Hvordan gør jeg nu, hvis buzz skal have point?
Jeg lovede jo 200 for en færdig løsning og du har da, Buzz gjordt en stor del af arbejdet!!!


Michael

Metoden i http://www.eksperten.dk/faq#faq-5-3 kan anvendes.