5 Robuste Schätzung linearer Trends

# robuste Regressionsgerade

x <- c(0.5,0.9,1,1.2,2.3,2.6,2.8,3.3,3.9,4.2,4.4,6.5)
y <- c(0.5,0.9,0.7,1.3,1.5,2.0,2.5,2.1,2.6,3.4,3.2,1.1)

pdf("gerade1.pdf",width=5,height=5)
par(mar=c(4,4,1,1))
plot(x,y,cex.lab=1.3)
abline(lsfit(x,y)$coef)
dev.off()

# robuste Regressionsgerade

x <- c(0.5,0.9,1,1.2,2.3,2.6,2.8,3.3,3.9,4.2,4.4,6.5)
y <- c(0.5,0.9,0.7,1.3,1.5,2.0,2.5,2.1,2.6,3.4,3.2,1.1)

pdf("gerade2.pdf",width=5,height=5)
par(mar=c(4,4,1,1))
plot(x,y,cex.lab=1.3)
abline(line(x,y)$coef)
dev.off()