qrea$tgl <- c(rep(0,100),qrea$tg[1:(length(qrea$tg)-100)]) # Verschiebung um 100 Abtastungen
qrea$tgl2 <- c(rep(0,200),qrea$tg[1:(length(qrea$tg)-200)]) # Verschiebung um 200 Abtastungen
qrea$tn <- 1:n

hl.rlm <- rlm( hl ~ I(hpos[tgl])+factor(tn>trig[trigv==31][1]) ,
 subset(qrea,tgl==tgl2&tgl>0&tgl<8),method="MM",init=c(16650,338,0))
hr.rlm <- rlm( hr ~ I(hpos[tgl])+factor(tn>trig[trigv==31][1]) ,
 subset(qrea,tgl==tgl2&tgl>7&tgl<15),method="MM",init=c(14650,338,0))