Linearni regrese – C / C++ – Fórum – Programujte.com

Tohle je dotázek na pokročilé, podařilo se mi najít kód na linearni regresi: http://www.oocities.org/…p/linreg.htm . Vzorec na linearni regresi je y=a+b*x. Já ale potřebuji vytvořit třídu na mnohonasobnou linearni regresi, kde mam treba 50 hodnot x. Tzn: y=a+b1*x1+b2*x2+b3*x3+..+b50*x50

Díky za jakýkoliv help, nedávám tomu ale moc, že se někdo ozve :)

První nástřel: z x a b bych vytvořil strukturu, ty naskládal do vektoru a ten projel cyklem

hu (= hlucheucho)

Založil jsem git: https://github.com/Netix/Multiple-Linear-Regression

Zatím obsahuje třídu LinearRegression sloužící jako předloha.

#3 fnenks
Stacilo dat link na zdroj a nemusel si ni zakladat  :

http://www.oocities.org/…/linregh.htm

mňo, dával jsem link v prvním příspěvku..

#5 fnenks
A funguje to?
 

Funguje, ale problém je že to je jenom y=a+b*x. Já potřebuji vytvořit třídu na mnohonásobnou lineárni regresi.

Mno, takže jsem přemýšlel nad implementací, v konstruktoru bude pouze jeden parametr a to počet x, dále funkce pro přidání dat bude mít dva parametry, první parametr bude vektor, který musí obsahovat stejný počet hodnot, jak bylo uvedeno v konstruktoru, jinak to hodí výjmku, druhý parametr bude y. Dále pak další dvě funkce: getA(), a getB(), přičemž getB bude vrace vektor. V tomhle mam jasno. Problém je s celým algoritmem, který sám nezvládnu.

https://github.com/…10fac9e8ca7c

pořád nevyřešené :(

nechápu nad čím bádáš, mě to připomíná polynom na výpočet teploty z napětí u termočlánků.

hu

#include <vector>

typedef struct {         //struktura predstavuje bn, xn
           float b;
           float x;
        } clen;

class Regrese
{
  private:
    float a;
    std::vector<clen> cleny;  //zde jsou ulozeny cleny b0, x0 ... bn, xn

  public:
    //potrebujes metody:
    int RekniPocetClenu();
    void Nastav_a(float a);
    int PridejClen(float b, float x);   //vrati pocet clenu po vlozeni
    float Regrese_y();  //vrati regresi
}

a implementace metody float Regrese_y():

float Regrese::Regrese_y()
{
  float y;
  std::vector<clen>::iterator it;  //iterator funguje jako pointer na prvky v kontejneru

  y = a;
  for(it = cleny.begin(); it != cleny.end(); it++)
  {
    y = y + it->b * it->x;  //zde se pristupuje k polozkam ve strukture
  }
  return y;
}

Přidat prvek do vektoru a určit, kolik jich je by nemělo být problém

hu

tohle neni linearni regrese

ja vkladam data x1, x2, x3, x4, x5, x6, x7, y

posilam ukazku v excelu: https://dl.dropboxusercontent.com/u/97350425/linreg.xlsx

první tabulka je naplněna daty, druha tabulka z dat vytvori hodnoty do vzorce, 3ti tabulka aplikuje vzorec a dopocita Y

ja potrebuji ziskat ty hodnoty v druhe tabulce z dat prvni tabulky, ale ne v excelu, ale v c++

#15 fnenks
Mas problem s implementaciou alebo s tym ze nevies ako sa pocita mnohonasobna lin. regresia?

podle vzorce co jsi sem dal, je y výsledek, ten nevkládáš. Někde musíš vzít ty koeficienty a, b a někam je uložit. Stejně tak i x. To, co jsem napsal, je implementace vzorce, který jsi napsal. Problém asi je, že se snažíš implementovat něco, co nevíš, co je. Nastuduj si teorii - co je to regrese, její matematický popis, algoritmus a to implementuj.

Edit: Google toho najde dost.

hu

#16 vitamin
mam problem s tim, ze nevim jak se pocita

zkus dát do Google multiple regression. V češtině to nevypadalo na něco použitelného.

hu

#18 fnenks

S vypoctom ti nepomozem, ale toto vyzera nadejne:

http://stackoverflow.com/questions/2910361/efficient-multiple-linear-regression-in-c-sharp-net

Najtazsie bude asi spravyt inverznu maticu.

zkusím spíš napsat na matematické forum..

nasel jsem skvely kod ktery skoro funguje jak ma, je to psane v C# takze celkem pohodicka prepsat do C++

zdroj: http://www.codeproject.com/Articles/25335/An-Algorithm-for-Weighted-Linear-Regression

porovnal jsem s excelem vysledek a sedi presne, akurat jeden problem ze vzorec je Y = X1 * A1 + X2 * A2 + .. + Xn * An;
Ja potrebuju vzorec Y = X1 * A1 + X2 * A2 + .. + Xn * An + B;

Vzorec co zde píšeš je prostý součet součinů k sobě příslušnému parametru a proměnné, nazval bych to "parametrickou funkcí s n parametry (A1 .. An) a n proměnnými (X1 .. Xn)". Jak ten vzorec implementovat jsem tu naznačil už před 10 dny - vektor o n prvcích, jehož prvky jsou struktury, které mají dvě položky - k sobě příslušné A a X. 

Pokud to není prostý součet součinů jak vyplývá ze vzorce, který zde uvádíš, tak máš problém matematicky špatně popsaný a to vede k jeho chybné implementaci.

Regrese je pojem ze statistiky, zřejmě bude přesně matematicky popsaný.

hu

#23 hlucheucho
tohle nikdy nezvladnu, o linearni regresi nic nevim, nevim jak se pocita, tudiz mam problem s matematikou, ne s programovanim

Třída na mnohonásobnou regresi hotova. Kdyby někdo potřeboval tak je na adrese: https://github.com/Netix/Multiple-Linear-Regression