Mitmene regressioonanalüüs ja mittelineaarne regressioonanalüüs
library(car)
vif(Mdhhv)
vif(Mdhv)
vif(Mdh)
anova(Md,Mdh)
anova(Md,Mdhv)
anova(Md,Mdhhv)
anova(Mdhv,Mdhhv)
anova(Mdh,Mdhhv)
par(mfcol=c(1,1), las=1)
plot(Md)
# 2: Näslundi funkts.
PD.KU<-subset(PD.1,pl=="KU")
par(mar=c(4.7,4.7,1,1))
plot(h~d_k, PD.KU, xlim=c(0,35),ylim=c(0,20),xlab="Diameeter, cm", ylab="Kõrgus, m")
abline(v=seq(0,40,10),lty=3,col="grey75")
abline(h=seq(0,25,5),lty=3,col="grey75") # abijooned
#lineariseerides
Mnsld.1 <- lm(I(d_k/(h-1.3)^(1/3))~d_k, PD.KU)
summary(Mnsld.1)
PD.KU$h.nsld<- 1.3+(PD.KU$d_k/
(Mnsld.1$coefficients[1]+Mnsld.1$coefficients[2]*PD.KU$d_k))^3
with(PD.KU, sqrt(sum((h-h.nsld)^2)/(length(h)-2)))
curve(1.3+(x/(Mnsld.1$coefficients[1]+Mnsld.1$coefficients[2]*x))^3, from=0, to=35,
col="red",add=T,lwd=2,lty=2)
rect(19,1,30,6,border=2,lty=2,col=0)
ac1 <- round(Mnsld.1$coefficients,3)
valem1 <- substitute(h==1.3+(frac(d, b0+b1*d))^3, list(b0=ac1[1],b1=ac1[2]))
annaabi.ee 
