v <- read.table("volt.dat", col.names = c("x", "y")) v0 <- read.table("volt0.dat", col.names = c("x", "y")) motif(); par(mfrow = c(2, 1)) plot(v0$x, v0$y, main = "volt0") plot(v$x, v$y, main = "volt (korrigiert)") plot(v0$x[-1], diff(v0$y), type = "b", main = "DIFF volt0") plot(v$x[-1], diff(v$y), type = "b", main = "DIFF volt (korrigiert)") attach("/u/maechler/s/goodies.Data") #-- for pl.ds -- s.v0 <- supsmu(v0$x, v0$y) s.v <- supsmu(v$x, v$y) pl.ds(s.v0$x, v0$y, s.v0$y, main = "Super smoother ( volt 0 )") l.v0 <- loess(y ~ x, data = v0) l.v <- loess(y ~ x, data = v) pl.ds(v$x, v$y, fitted(l.v), main = "loess [Default] ( volt 0 )")