Regrese ve VB - 2. díl
 x   TIP: Přetáhni ikonu na hlavní panel pro připnutí webu
Reklama

Regrese ve VB - 2. dílRegrese ve VB - 2. díl

 

Regrese ve VB - 2. díl

Google       Google       5. 3. 2007       12 713×

V minulém díle jsme skončili u regrese lineární, ale lineární průběh není jediným průběhem, který potřebujeme zkoumat. Mezi základní průběhy patří ještě průběh logaritmický a exponenciální. První z nich bude předmětem tohoto dílu.

Reklama
Reklama

Logaritmická regrese slouží k výpočtu přibližných výsledků z několika zadaných hodnot, ty musí být, stejně jako při regresi lineární, dvě. Tedy dvě dvojice hodnot. Logaritmický průběh je doufám všem známý, jinak byste to asi ani nečetli. Takže se pustím rovnou do zpracování kódu. K tomu ale potřebujeme nějaké vzorečky.

Při psaní programu se ale setkáváme s problémem, tím je přirozený logaritmus, pro něj totiž VB nemá funkci, proto využijeme následujícího vzorce, kde je eulerovo číslo.

Nyní tedy již ke kódu. Prvně si opět ukážeme, jak řešit příklad z přímo zadaných hodnot a potom si ukážeme, jak vytvořit dynamičtější program.

Private Sub Command1_Click()
    e = 2.718281828
    
    xa = 1 'vzorky x
    xb = 2
    
    ya = 3 'vzorky y
    yb = 4
    
    n = 2 'počet vzorků
    
    'výpočet regresního koeficientu B
    B = ((n * (((Log(xa) / Log(e)) * ya) + ((Log(xb) / Log(e)) * yb))) - (((Log(xa) / Log(e)) + (Log(xb) / Log(e))) * (ya + yb))) / ((n * (((Log(xa) / Log(e)) ^ 2) + ((Log(xb) / Log(e)) ^ 2))) - (((Log(xa) / Log(e)) + (Log(xb) / Log(e))) ^ 2))
    Print B

    'výpočet regresního koeficientu A
    A = ((ya + yb) - (B * ((Log(xa) / Log(e)) + (Log(xb) / Log(e))))) / (n)
    Print A
    
    x = 3 'hodnota kterou chceme dopočítat
    
    y = A + (B * (Log(x) / Log(e))) 'závěrečný vzorec pro výpočet y
    Print y
End Sub

Vzoreček je dosti dlouhý, proto je možná vhodnější si ho rozepsat na části. Viz další kód.

 Private Sub Command1_Click()
    e = 2.718281828
    
    xa = 1 'vzorky x
    xb = 2
    
    ya = 3 'vzorky y
    yb = 4
    
    n = 2 'počet vzorků
    
    part1 = n * (((Log(xa) / Log(e)) * ya) + ((Log(xb) / Log(e)) * yb))
    part2 = ((Log(xa) / Log(e)) + (Log(xb) / Log(e))) * (ya + yb)
    part3 = n * (((Log(xa) / Log(e)) ^ 2) + ((Log(xb) / Log(e)) ^ 2))
    part4 = ((Log(xa) / Log(e)) + (Log(xb) / Log(e))) ^ 2
    
    'výpočet regresního koeficientu B
    B = (part1 - part2) / (part3 - part4)
    Print B

    
    'výpočet regresního koeficientu A
    A = ((ya + yb) - (B * ((Log(xa) / Log(e)) + (Log(xb) / Log(e))))) / (n)
    Print A
    
    x = 3 'hodnota kterou chceme dopočítat
    
    y = A + (B * (Log(x) / Log(e))) 'závěrečný vzorec pro výpočet y
    Print y
End Sub

Nyní se vrhneme na program dynamický, postup bude stejný jako pro regresi lineární, pouze vzorečky se budou lišit.

Private Sub Command1_Click()
    Const e = 2.718281828
    Dim data(256, 1) As Single
    Dim sumlnxy As Single
    Dim sumlnx As Single
    Dim sumy As Single
    Dim sumlnx2 As Single
    Dim n As Byte
    Dim x As Single
    Dim y As Single
    Dim A As Single
    Dim B As Single
    
    n = InputBox("Zadej počet vzorů", "Požadavek na data", "2")
    temp = n - 1
    
    For i = 0 To temp
    x = InputBox("Zadej hodnotu x dat pro n = " & i + 1, "Požadavek na data", "0")
    y = InputBox("Zadej hodnotu y dat pro n = " & i + 1, "Požadavek na data", "0")
    data(i, 0) = x
    data(i, 1) = y
    Next i
    
    sumxy = 0
    sumx = 0
    sumy = 0
    sumx2 = 0
    
    For i = 0 To temp
    sumlnxy = sumlnxy + (Log(data(i, 0)) / Log(e)) * data(i, 1)
    sumlnx = sumlnx + (Log(data(i, 0)) / Log(e))
    sumy = sumy + data(i, 1)
    sumlnx2 = sumlnx2 + ((Log(data(i, 0)) / Log(e)) * (Log(data(i, 0)) / Log(e)))
    Next i
    
    B = ((n * sumlnxy) - (sumlnx * sumy)) / ((n * sumlnx2) - (sumlnx * sumlnx))
    Print B
    A = (sumy - (B * sumlnx)) / n
    Print A
    
    Min = InputBox("Zadej interval grafu hodnotu min: ", "Požadavek na data", "0")
    Max = InputBox("Zadej interval grafu hodnotu max: ", "Požadavek na data", "0")
    
    For x = Min To Max
        y = A + (B * (Log(x) / (Log(e))))
        Print "x = " & x & ", y = " & y
    Next x
End Sub

A ještě bych upozornil na to, že nemůžete hledat logaritmy čísla nula a čísel menších a taky na to že, příště toto naše povídání zakončíme exponenciální regresí a pak se přesuneme ke tvorbě grafů.

×Odeslání článku na tvůj Kindle

Zadej svůj Kindle e-mail a my ti pošleme článek na tvůj Kindle.
Musíš mít povolený příjem obsahu do svého Kindle z naší e-mailové adresy kindle@programujte.com.

E-mailová adresa (např. novak@kindle.com):

TIP: Pokud chceš dostávat naše články každé ráno do svého Kindle, koukni do sekce Články do Kindle.

2 názory  —  2 nové  
Hlasování bylo ukončeno    
0 hlasů
Google
(fotka) Jiří ChytilAutor programuje ve VB, zajímá se o elektrotechniku, studuje na SOŠ Elektrotechnické - obor číslicová technika.
Web    

Nové články

Obrázek ke článku Malware KONNI se úspěšně skrýval 3 roky. Odhalil ho bezpečnostní tým Cisco Talos

Malware KONNI se úspěšně skrýval 3 roky. Odhalil ho bezpečnostní tým Cisco Talos

Bezpečnostní tým Cisco Talos odhalil celkem 4 kampaně dosud neobjeveného malwaru, který dostal jméno KONNI. Ten se dokázal úspěšně maskovat od roku 2014. Zpočátku se malware zaměřoval pouze na krádeže citlivých dat. Za 3 roky se ale několikrát vyvinul, přičemž jeho současná verze umožňuje útočníkovi z infikovaného počítače nejenom krást data, ale i mapovat stisky na klávesnici, pořizovat screenshoty obrazovky či v zařízení spustit libovolný kód. Pro odvedení pozornosti oběti zasílali útočníci v příloze také obrázek, zprávu a výhružkách severokorejského režimu či kontakty na členy mezinárodních organizací.

Reklama
Reklama
Obrázek ke článku Pouze jedna z deseti lokálních firem ví o pokutách plynoucích z GDPR

Pouze jedna z deseti lokálních firem ví o pokutách plynoucích z GDPR

Trend Micro, celosvětový lídr v oblasti bezpečnostních řešení a VMware, přední světový dodavatel cloudové infrastruktury a řešení pro podnikovou mobilitu, oznámily výsledky výzkumu mezi českými a slovenskými manažery zodpovědnými za ochranu osobních údajů, který zjišťoval, jak jsou připraveni na nové nařízení o ochraně osobních údajů (GDPR). Většina firem v České republice a na Slovensku nad 100 zaměstnanců je již s novým nařízením GDPR obeznámena. Výzkum provedený ve spolupráci s agenturou Ipsos ukázal, že téměř 8 firem z 10 o nařízení ví, přičemž jeho znalost je o něco vyšší na Slovensku (89 %) než v České republice (69 %).

Obrázek ke článku Vyděračský software Locky se vrací, tváří se jako potvrzení platby, odhalil tým Cisco Talos

Vyděračský software Locky se vrací, tváří se jako potvrzení platby, odhalil tým Cisco Talos

Jeden z nejznámějších ransomwarů, Locky, se vrací. Po většinu roku 2016 patřil mezi nejrozšířenější vyděračské softwary. Ke svému šíření využíval emailové kampaně s infikovanými přílohami. Ransomware Locky byl rozesílán prostřednictvím botnetu (internetový robot zasílající spamy) Necurs. Jeho aktivita na konci roku 2016 téměř upadla a spolu s ní i šíření ransomwaru Locky. Před několika týdny se Necurs opět probudil a začal posílat spamy nabízející výhodný nákup akcií. Dne 21. dubna zaznamenal bezpečnostní tým Cisco Talos první velkou kampaň ransomwaru Locky prostřednictvím botnetu Necurs za posledních několik měsíců.

Obrázek ke článku Dovozci baterií mění logistiku, letadlo nahrazuje námořní doprava

Dovozci baterií mění logistiku, letadlo nahrazuje námořní doprava

Dovozci baterií do mobilů či notebooků upouštějí od letecké přepravy zboží. V letošním roce plánují dovézt až 80 % produktů lodí. Přitom před 5 lety byla většina baterií do mobilních přístrojů dovezených do České republiky přepravována letadlem. Za proměnou způsobu transportu akumulátorů stojí zpřísnění pravidel pro leteckou přepravu, která přinášejí vyšší náklady i náročnou agendu.

loadingtransparent (function() { var po = document.createElement('script'); po.type = 'text/javascript'; po.async = true; po.src = 'https://apis.google.com/js/plusone.js'; var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(po, s); })();
Hostujeme u Českého hostingu       ISSN 1801-1586       ⇡ Nahoru Webtea.cz logo © 20032017 Programujte.com
Zasadilo a pěstuje Webtea.cz, šéfredaktor Lukáš Churý