28087

data2<- read.tableCdrukarki.txt", header = FALSE) colnames(data2)<-c('czas¹,'liczbaDruk¹) czas<-data2$czas LDr<-data2$liczbaDruk

plot(czas,LDr)

hist(czas,col='red')

X<-czas Y<-LDr lm( Y - X)

Im.linear <- lm(Y - X)

Im.linear

summa ry (Im. linea r) residuals(Im.linear)

Im.linear.resids <- residuals(Im.linear) hist(lm.linear.resids, col='blue')

plot(X,Y)

lines(X,fitted(Im.linear)) plot(X,Y)

lines(X,fitted(Im.linear),col='blue') segments(X,fitted(lm.linear),X,Y,col='red')

qqno rm(resid(Im.linea r), col='red') pred.framę <- data.framę(X=l:160)

pp <- predict (Im .linear, int="p", newdata=p red. f ramę) pc <- predict (Im .linear, int=^,'c", newdata=pred .f ramę) plot(X,Y,ylim=range(Y, pp, na.rm=T)) pred.X <- pred.frame$X

matlines(pred.X, pc, lty=c(1,2,2), col="red") matlines(pred.X, pp, lty=c(1,3,3), col="blue")

p redict(Im.linea r,int="c") predict (Im.linear,int="p")

cor(X,Y, use="complete.obs") cor(data2,use="complete.obs‘') cor.test(X,Y)