par(mar=c(4.5,4.5,1,1)) plot(NULL,xlim=c(0,40),ylim=c(0,25),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 points(h~d_k,data=subset(PD.KU),lwd=1) with(subset(PD., pl=="KU"),rug(d_k)) 1. Sirge h=a+b*d M1 <- lm(h~d_k, data=PD.KU) summary(M1) D<-0:40 M1.pred <- predict(M1,newdata=data.frame(d_k=D)) lines(D,M1.pred, col="red") coefficients(M1)[1] coefficients(M1)[2] # dobavit' p-value v tablicu v vide * summary(M1)$adj.r.squared summary(M1)$sigma # sqrt(sum(M1$residuals^2)/(length(M1$residuals)-2)) AIC(M1) > coefficients(M1)[1] (Intercept) 7.758512 > coefficients(M1)[2] # dobavit' p-value v tablicu v vide * d_k 0.5027412 > summary(M1)$adj.r.squared [1] 0.7619479 > summary(M1)$sigma # sqrt(sum(M1$residuals^2)/(length(M1$residuals)-2)) [1] 1.823029 > AIC(M1) [1] 177.6237 2. Ruutparabool h=a+b*d+c*d^2 M2 <- lm(h~d_k + I(d_k^2), data=PD.KU)
, . 1841 . (kak i vse ostalnoe etot avtomobil towe ne idealen i u nego est mnowestvo nedostatkov naprimer on o4en dolgo zarewaetsa i na nem mowna o4en malo proehat) · . · , « » . · . · . 2000- · . · . . Akamulatorq derwat malo energii toest nado 4awe zarewat mawinu , (, ) . . (No i takwe u etoi mowini mnogo preimuwestv .ewe mowna dobavit odnu vew 4to eta mawina dorogaja , , ne smotra na to 4to mnogie piwut ili govorat 4to oni dewevle 4em obq4naja toplevnaja mawina ) · . · . · . · , . · . . (V Estonii Tartu mq mowem uvidet Taksi gde ezdijut tolko na elektromobile) -- . , , . : Toyota RAV4 EV, ZENN, ZAP Xebra, General Motors EV1, Chevrolet Volt, Volvo C30