"subr1" - 1 õppematerjal

Mitmene regressioonanalüüs ja mittelineaarne regressioonanalüüs

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])) text(20,4,valem1,pos=4) s1 = round(with(PD.KU, sqrt(sum((h-h.nsld)^2)/(length(h)-2))),2) subr1 = substitute(s[e]==s1,list(s1=s1)) text(20,2,subr1,pos=4) # mittelineaarne Mnsld.2<-nls(h~1.3+(d_k/(a+b*d_k))^3, PD.KU, start=list(a=1,b=0.4)) summary(Mnsld.2) curve(1.3+(x/(coefficients(summary(Mnsld.2))[1]+coefficients(summary(Mnsld.2))[2]*x))^3, from=0, to=35, col=3,add=T,lwd=2) rect(22,7,33,12,border=3,col=0) ac2 <- round(coefficients(summary(Mnsld.2))[1:2],3) valem2 <- substitute(h==1.3+(frac(d, b0+b1*d))^3, list(b0=ac2[1],b1=ac2[2])) text(23,10,valem2,pos=4)