6.2.1 Upper and Lower Smoothing

# Upper and Lower Smoothing streifen <- read.csv( " hamster.csv " ) x <- streifen[,1] y <- streifen[,2] a <- lowess(x,y,f=2/3) yhat <- a$y r <- y-yhat rplus <- lowess(x[sign(r)==1],r[sign(r)==1],f=2/3)$y rminus <- lowess(x[sign(r)==-1],r[sign(r)==-1],f=2/3)$y pdf("lowess5.pdf",width=6,height=6) par(mar=c(4,4,1,1)) plot(x,y,xlab="Anteil Winterschlaf [%]",ylab="Alter zum Todeszeitpunkt [Tage]", cex.lab=1.2) lines(a) lines(x[sign(r)==1],yhat[sign(r)==1]+rplus) lines(x[sign(r)==-1],yhat[sign(r)==-1]+rminus) dev.off()