Lineaarne regressioon sõltuvalt diskreetsest faktorist
Sellelt jooniselt on näha, et seos pH ja läbipaistvuse vahel sõltub omakorda järvest. Otepää Pikajärve andmetel saame langustrendi - kui pH suureneb, siis läbipaistvus väheneb, Piigandi järve puhul sama väita oleks väga vale.
Statistilise järelduse tegemisel peab olema korrektne, üldistada ei saa seost kõikidele Eesti veekogudele, samuti tuleb olla kindel, et andmed on kogutud sesoonselt samal ajal, bioloogias on tavaliselt sesoonsus üks olulisi faktoreid, mis mõjutab paljude tunnuste väärtusi.
Kuidas faktorite mõju joonisel kajastada? Fikseerides teatavate faktorite nivood, saame illustreerida tunnustevahelisi seoseid (konkreetsemalt siis kahe erineva pideva arvtunnuse vahel)
Kõrvaloleva graafiku saaksime teha ka nii, et võtame mitte ainult 2 järve, vaid näiteks 2 erinevat tüüpi veekogud ja ühendame nende valimid, järve nime asemel kasutame järve tüüpi.
andmed=read.csv("2jarve.csv")
attach(andmed)
plot(PH, LP,col=NR,pch="*",main="Seos kahe järve andmetel",ylab="Läbipaistvus, m",xlab="pH")
1.järve joon
abline(lm(LP~PH,data=andmed[NR==1084,]),col="blue")
text(7,6,"Piigandi",col=4)
2.järve joon
abline(lm(LP~PH,data=andmed[NR==1078,]),col="red")
text(8.2,1.5,"Otepää Pikajärv",col="red")
Eelmine Järgmine
Statistilise järelduse tegemisel peab olema korrektne, üldistada ei saa seost kõikidele Eesti veekogudele, samuti tuleb olla kindel, et andmed on kogutud sesoonselt samal ajal, bioloogias on tavaliselt sesoonsus üks olulisi faktoreid, mis mõjutab paljude tunnuste väärtusi.
Kuidas faktorite mõju joonisel kajastada? Fikseerides teatavate faktorite nivood, saame illustreerida tunnustevahelisi seoseid (konkreetsemalt siis kahe erineva pideva arvtunnuse vahel)
Kõrvaloleva graafiku saaksime teha ka nii, et võtame mitte ainult 2 järve, vaid näiteks 2 erinevat tüüpi veekogud ja ühendame nende valimid, järve nime asemel kasutame järve tüüpi.
andmed=read.csv("2jarve.csv")
attach(andmed)
plot(PH, LP,col=NR,pch="*",main="Seos kahe järve andmetel",ylab="Läbipaistvus, m",xlab="pH")
1.järve joon
abline(lm(LP~PH,data=andmed[NR==1084,]),col="blue")
text(7,6,"Piigandi",col=4)
2.järve joon
abline(lm(LP~PH,data=andmed[NR==1078,]),col="red")
text(8.2,1.5,"Otepää Pikajärv",col="red")
Eelmine Järgmine