t <- seq(1:36)
t
y <- c(1.2, 0.7, 2.0, 2.2, 3.2, 2.7, 3.8, 4.3, 3.6, 2.9, 2.9, 7.7, 
       1.4, 0.7, 2.2, 2.3, 3.5, 2.4, 4.0, 4.1, 3.3, 2.7, 3.1, 6.9, 
       1.3, 0.8, 1.9, 2.1, 2.9, 2.5, 4.1, 4.0, 3.2, 2.6, 3.0, 7.6)
y
summary(y)
plot(NULL, xlim=c(-1, 37), ylim=c(-3, 10), xlab="Zeit", ylab="")
points(t,y)

lines(t,y)

Trend <- lm(y~t)
G <- fitted.values(Trend)
G
abline(Trend, col=134)

Trendbereinigt <- c(y - G)
Trendbereinigt

points(t, Trendbereinigt, col=50)
lines(t, Trendbereinigt, col=50)

Ss <- 1:12
for (i in 1:12)
{Ss[i] <- 1/3*(   Trendbereinigt[i + 12*0] 
                + Trendbereinigt[i + 12*1] 
                + Trendbereinigt[i + 12*2]
                )
}
Ss

Null <- rep(0, 36)
Null

S <- Ss + Null 
S

points(t, S, col=25)
lines(t, S, col=25)

R <- Trendbereinigt - S
R

windows()
plot(t, R)