Matematické Fórum

Nevíte-li si rady s jakýmkoliv matematickým problémem, toto místo je pro vás jako dělané.

Nástěnka
22. 8. 2021 (L) Přecházíme zpět na doménu forum.matweb.cz!
04.11.2016 (Jel.) Čtete, prosím, před vložení dotazu, děkuji!
23.10.2013 (Jel.) Zkuste před zadáním dotazu použít některý z online-nástrojů, konzultovat použití můžete v sekci CAS.

Nejste přihlášen(a). Přihlásit

#1 06. 06. 2015 16:35

lisakpodsity
Příspěvky: 150
Reputace:   
 

Statistika-kvadratická regrese

Dobrý den, toto bylo zadání ve zkoušce 4)    Pro data2.txt zjistěte a interpretujte nejdříve korelaci a následně zjistěte regresní vztah s využitím kvadratické regrese. 1. sloupec je závislá proměnná, druhý sloupec je nezávislou proměnnou. Výsledek regrese zapište ve formě regresní rovnice a interpretujte ho. Určete z něj kvalitu modelu.

Problém je v tom, že data nepocházejí z normálního rozdělení, ani u jedné proměnné. A jsou kladná i záporná, tudíž nejdou zlogaritmovat.

Nevíte co s tím ? Děkuji za nápady

Offline

 

#2 06. 06. 2015 16:39

Stýv
Vrchní cenzor
Příspěvky: 5710
Reputace:   215 
Web
 

Re: Statistika-kvadratická regrese

data nemají pocházet z normálního rozdělení. rezidua mají mít normíální rozdělení. a i když ho nemají, tak to není žádná velká tragédie

Offline

 

#3 06. 06. 2015 16:50 — Editoval lisakpodsity (06. 06. 2015 17:00)

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

Moc děkuji za bleskurychlou reakci.  Postupoval sem podle vzorce na kvadratickou regresi a vyšlo mi toto. 
Call:
lm(formula = (Var1) ~ Var2 + I(Var2^2), data = data2)

Coefficients:
(Intercept)         Var2    I(Var2^2) 
-2.177e+00   -5.845e-01   -2.616e-05 

Vubec nevím jestli to dává nějaký smysl ... A vubec nevím jak z toho mám dostat regresní rovnici

Offline

 

#4 06. 06. 2015 18:47

Eratosthenes
Příspěvky: 3111
Reputace:   140 
 

Re: Statistika-kvadratická regrese

ahoj ↑ lisakpodsity:,

pochybuju, že z toho někdo pozná, zda to dává smysl nebo ne, když nemá vstupní data. Já osobně vůbec nevím,  jaké máš "vzorce na kvadratickou regresi" a zápis

Call:
lm(formula = (Var1) ~ Var2 + I(Var2^2), data = data2)

je pro mě zajímavým čínským textem...


Budoucnost patří aluminiu.

Offline

 

#5 06. 06. 2015 19:24

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ Eratosthenes:
je to zápis z Erka.

Offline

 

#6 06. 06. 2015 19:24

Stýv
Vrchní cenzor
Příspěvky: 5710
Reputace:   215 
Web
 

Re: Statistika-kvadratická regrese

↑ Eratosthenes: to je výstup z Rka

↑ lisakpodsity: dává smysl ten výstup? ano. dává smysl něco psát do Rka, když netušíš, co vlastně děláš a proč, a nerozumíš pak ani výstupu? neřekl bych

Offline

 

#7 06. 06. 2015 22:09

Eratosthenes
Příspěvky: 3111
Reputace:   140 
 

Re: Statistika-kvadratická regrese

↑ Stýv:

No nevím sice ani to, co je Rko, ale pokládám za nesmírně pikantní ten detail, že zápisu

(Intercept)         Var2    I(Var2^2) 
-2.177e+00   -5.845e-01   -2.616e-05 

rozumím možná líp než ↑ lisakpodsity:, který by to Rko měl asi znát :-)


Budoucnost patří aluminiu.

Offline

 

#8 07. 06. 2015 07:50

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ Eratosthenes:
Ano,  tento fakt je pikantní :)  díky za tvoje cenné rady

Offline

 

#9 07. 06. 2015 12:21

Jozef3
Příspěvky: 276
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ lisakpodsity:
Smysl tento výstup samozřejmě dává. Znamená to, že regresní funkce má předpis $y = -2,177 - 0,5845 \cdot x - 0,00002616  \cdot x^{2}$.
Vidíte ješte někde problém?

Offline

 

#10 07. 06. 2015 14:20

Eratosthenes
Příspěvky: 3111
Reputace:   140 
 

Re: Statistika-kvadratická regrese

ahoj ↑ Jozef3:,

no, trochu problém vidím. Podle toho, jak dopadla parabola, mně trochu uniká smysl kvadratické regrese a ke svým "cenným radám" pro ↑ lisakpodsity: bych přidal ještě jednu: kvadratický člen směle škrtni a na kvadratickou regresi zapoměň - lineární docela stačí. Její koeficienty zaokrouhli na tolik platných míst, kolik mají vstupní data (pochybuji, že jich bylo šest...).


Budoucnost patří aluminiu.

Offline

 

#11 07. 06. 2015 14:39 — Editoval lisakpodsity (07. 06. 2015 14:41)

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ Eratosthenes:
Děkuji za Vaše rychlé odpovědi. Rád bych provedl lineární regresi, ale vím že vstupní předpoklad je normalita dat, kterou tu nesplňuji, tudíž ta kvadratická. Ještě by mě zajímalo, jak z něj určím kvalitu modelu...

Offline

 

#12 07. 06. 2015 15:06

Jozef3
Příspěvky: 276
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ lisakpodsity:
Normalita dat není žádným vstupním předpokladem pro lineární regresi. Normalita je v regresní analýze dobrá akorát k tomu, že pak Vám vyjdou přesně p-hodnoty a intervaly spolehlivosti pro regresní parametry. Bez splnění předpokladu normality jsou tyto p-hodnoty a intervaly spolehlivosti pouze asymptoticky správné.

Kvalitu modelu můžete posuzovat podle několika kritérií. Může Vás třeba zajímat, jestli jsou splněny předpoklady homoskedasticity nebo nekorelovanosti reziduí (což už vskutku jsou předpoklady lineární regrese), jestli nejsou v datech odlehlá pozorování atd. Toto posoudíte nejlépe pomocí obrázků (což takhle použít v Rku funkci plotLM z balíčku mffSM?). Nebo Vás může zajímat předpovědní kvalita modelu, kterou posoudíte např. podle koeficientu determinace (to R-squared ve výpisu funkce lm).

P.S: To co zde nazýváte "kvadratickou regresí" je klasická lineární regrese. Jde totiž o to, že regresní funkce je lineární v regresních parametrech, a to, že do modelu dáte regresor v druhé mocnině na tom nic nemění.

Offline

 

#13 07. 06. 2015 15:15 — Editoval lisakpodsity (07. 06. 2015 15:31)

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ Jozef3:
Mnohokrát děkuji za Vaši obsáhlou a rychlou odpověď. Přecházím z programu Statistica, takže ještě si uplně s Erkem netykám. Mohl byste mi poradit,  jak docílím tohoto  koeficientu determinace (to R-squared ve výpisu funkce lm). Nikde jsem to bohužel nenašel....
Multiple R-squared:  0.5132,    Adjusted R-squared:  0.5125   nalezeno. Koeficient není příliš vysoký, zhruba polovina rozptylu je vysvětlena regresí...

Offline

 

#14 07. 06. 2015 15:33

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ lisakpodsity:
//forum.matweb.cz/upload3/img/2015-06/84003_Rplot.jpeg


Vyšel mi z toho tento graf :-)

Offline

 

#15 07. 06. 2015 15:50

Eratosthenes
Příspěvky: 3111
Reputace:   140 
 

Re: Statistika-kvadratická regrese

↑ lisakpodsity:

no, já tedy nevím, statistik nejsem, ale tady bych asi regresi úplně zahodil. Připadá mi to jako naprostý chaos...


Budoucnost patří aluminiu.

Offline

 

#16 07. 06. 2015 16:01

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ Eratosthenes: nevím co je to za graf,  ale qq plot s regresní přímkou vypadá takhle,  to už celkem smysl dává :)


//forum.matweb.cz/upload3/img/2015-06/85707_Rplot01.jpeg

Offline

 

#17 07. 06. 2015 16:31

Jozef3
Příspěvky: 276
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ lisakpodsity:
No vidíte. Koeficient determinace je to Multiple R-squared a vzhledem k vysokému počtu pozorování bych vůbec neřekl, že hodnota 0,5132 je nízká - to je celkem slušné. Z Vámi uváděného prvního grafu se zdá, že předpoklad homoskedasticity (konstantního rozptylu reziduí) je celkem uspokojivě splněn. Naopak z QQ-grafu je vidět, že s předpokladem normality reziduí nebude všechno v pořádku.
Ještě byste měl pomocí funkce plotLM získat jeden graf. Je tomu tak? Pokud ano, ukažte nám ho.

Offline

 

#18 07. 06. 2015 16:46

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ Jozef3:
//forum.matweb.cz/upload3/img/2015-06/88371_Rplot.jpeg
ano ještě tento graf

Offline

 

#19 07. 06. 2015 16:54

Jozef3
Příspěvky: 276
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ lisakpodsity:
Tento jsem zrovna nemyslel, ale budiž. Tady vidíte, že Váš model má akorát jedno vlivné pozorování, a to je pozorování s pořadovým číslem 1 224.
Myslel jsem graf Scale-Location, kde na x-ové ose je Fitted values a na y-ové ose odmocnina z abs. hodnoty Standardized residuals.

Offline

 

#20 07. 06. 2015 17:02

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ Jozef3:
//forum.matweb.cz/upload3/img/2015-06/89185_Rplot01.jpeg
děkuji.  takže jestli to tedy chápu správně, zadání bylo   
4)    Pro data2.txt zjistěte a interpretujte nejdříve korelaci a následně zjistěte regresní vztah s využitím kvadratické regrese. 1. sloupec je závislá proměnná, druhý sloupec je nezávislou proměnnou. Výsledek regrese zapište ve formě regresní rovnice a interpretujte ho. Určete z něj kvalitu modelu

toto je regresní rovnice: y=-2,177 - 0,5845x - 0,00002615x^2
ale nev9m jak bzch to tedz intepretoval..

Offline

 

#21 07. 06. 2015 17:20

Jozef3
Příspěvky: 276
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ lisakpodsity:
Jak byste to interpretoval?
Třeba takto: 1) Zřejmě jste zjistil, že kvadratický člen v regresním modelu je nevýznamný (což si ostatně formálně ověříte, když se podíváte na p-hodnotu u členu I(Var2^2) ve výstupu funkce lm).
2) Zjistil jste, že pokud se x zvětší o jednotku, tak očekávaná hodnota y se zmenší o přibližně 0,5845.
3) Pokud je x nulové, tak očekávaná hodnota y je -2,7177.

To jsou asi tři nejzajímavější výstupy z Vaší regresní analýzy.

Offline

 

#22 07. 06. 2015 18:24

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ Jozef3:
Díky a kvalita modelu je dle R-squared poměrně dobrá :-)
Ještě bych potřeboval poradit ve vedlejším vláknu...

Offline

 

#23 07. 06. 2015 19:37

Jozef3
Příspěvky: 276
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ lisakpodsity:
Myslíte s tou analýzou rozptylu jednoduchého třídění? To je snadné, tam jste jednoznačně zamítl hypotézu o rovnosti středních hodnot.

Offline

 

#24 07. 06. 2015 19:39

lisakpodsity
Příspěvky: 150
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ Jozef3: Děkuji, ano podle nízkého p value.
//forum.matweb.cz/upload3/img/2015-06/98707_gr.png
tady je k tomu i ten graf. Jenže jak odpovědět na tuto otázku ? Zjistěte, zda existuje rozdíl mezi hodnotami sloupce Z v závislosti na vytvořených skupinách.

Neexistuje tedy rozdíl ?

Offline

 

#25 07. 06. 2015 20:52

Jozef3
Příspěvky: 276
Reputace:   
 

Re: Statistika-kvadratická regrese

↑ lisakpodsity:
Existuje, hypotézu o rovnosti středních hodnot jednotlivých skupinek jste zamítl.

Offline

 

Zápatí

Powered by PunBB
© Copyright 2002–2005 Rickard Andersson