Lineaarne regressioonijoon koos usalduspiiridega
Antud joonisel on lisatud regressioonsirgele ka läbipaistvuse keskväärtuse usalduspiirid ja prognoosipiirid sõltuvalt pH tasemest. Prognoosipiirid näitavad missugused läbipaistvuse väärtused saame 95 protsendil juhtudest, kui me teeme uue katse.
Näiteks jooniselt loeme välja, et pH=8.8 korral peaks Otepää Pikajärves meie andmetele tuginedes olema läbipaistvus vahemikus 0.3-0.9. Keskmine läbipaistvus sama pH juures aga 0.5 - 0.78
Regressioonijoone headust näitab see, et enamus vaatlusi jääb usalduspiiride vahele ja selle lähedusse.
Andmed failis 2jarve.csv
Eraldada tuleb Otepää Pikajärve andmed.Otepaa=subset(andmed,NR==1078)
attach(Otepaa)
Kõigepealt teeme mudeli
mudel1=lm(LP~PH)
X=c(8,9.3) # prognooside arvutamiseks anname ette miinimum ja maksimumväärtuse pH jaoks
F=predict(mudel1,data.frame(PH=X)) #prognoosid
CI=predict(mudel1,data.frame(PH=X),interval="conf") #usalduspiirid
PI=predict(mudel1,data.frame(PH=X),interval="pred") #prognoosipiirid
plot(PH,LP,main="Otepää Pikajärv",ylab="Labipaistvus, m",xlab="pH",ylim=c(0,1.2),xlim=c(8,9.3),pch="*")
lines(X,F,lty=1,col="black") #joonte lisamine, enne peab olema hajuvusdiagramm tellitud
lines(X,CI[,2],lty=2,col="blue");lines(X,CI[,3],lty=2,col="blue") #usalduspiir sinine
lines(X,PI[,2],lty=2,col="red");lines(X,PI[,3],lty=2,col="red") #prognoosipiir punane
Eelmine Järgmine
Näiteks jooniselt loeme välja, et pH=8.8 korral peaks Otepää Pikajärves meie andmetele tuginedes olema läbipaistvus vahemikus 0.3-0.9. Keskmine läbipaistvus sama pH juures aga 0.5 - 0.78
Regressioonijoone headust näitab see, et enamus vaatlusi jääb usalduspiiride vahele ja selle lähedusse.
Andmed failis 2jarve.csv
Eraldada tuleb Otepää Pikajärve andmed.Otepaa=subset(andmed,NR==1078)
attach(Otepaa)
Kõigepealt teeme mudeli
mudel1=lm(LP~PH)
X=c(8,9.3) # prognooside arvutamiseks anname ette miinimum ja maksimumväärtuse pH jaoks
F=predict(mudel1,data.frame(PH=X)) #prognoosid
CI=predict(mudel1,data.frame(PH=X),interval="conf") #usalduspiirid
PI=predict(mudel1,data.frame(PH=X),interval="pred") #prognoosipiirid
plot(PH,LP,main="Otepää Pikajärv",ylab="Labipaistvus, m",xlab="pH",ylim=c(0,1.2),xlim=c(8,9.3),pch="*")
lines(X,F,lty=1,col="black") #joonte lisamine, enne peab olema hajuvusdiagramm tellitud
lines(X,CI[,2],lty=2,col="blue");lines(X,CI[,3],lty=2,col="blue") #usalduspiir sinine
lines(X,PI[,2],lty=2,col="red");lines(X,PI[,3],lty=2,col="red") #prognoosipiir punane
Eelmine Järgmine