qrea.nas.orig <- apply(qrea,1,function(x)any(is.na(x)))   # die NA's merken, falls interpoliert wird (hier überflüssig)

qrea$tgl <- NULL
qrea$tgl2 <- NULL
qrea$tn <- NULL

#plot.ts(qrea,main=fn)
save(qrea,file=paste(fn,".dat",sep=""))

qrea.nas <- apply(qrea,1,function(x)any(is.na(x)))         # noch mal NA's merken

list(N=n,st=st,coef.hl=coef(hl.rlm),coef.hr=coef(hr.rlm),trig=trig,trigv=trigv,
na.hl=sum(is.na(qrea$hl)),
na.hr=sum(is.na(qrea$hr)),
na.vl=sum(is.na(qrea$vl)),
na.vr=sum(is.na(qrea$vr)),
nas.orig=cbind(beg=(1:length(qrea.nas.orig))[qrea.nas.orig & !c(F,qrea.nas.orig[-length(qrea.nas.orig)])],
  end=(1:length(qrea.nas.orig))[!c(qrea.nas.orig[-1],F) & qrea.nas.orig]),
nas=cbind(beg=(1:length(qrea.nas))[qrea.nas & !c(F,qrea.nas[-length(qrea.nas)])],
  end=(1:length(qrea.nas))[!c(qrea.nas[-1],F) & qrea.nas])
)