fill.md <-function(x,indices)
{                    # Auffüllen von x mit fehlenden Werten
                     # indices: Indices der fehlenden Positionen
                     # x      : Ergebnisvektor, bei dem fehlende Werte eliminiert sind
                     #
                     # Aufruf-Beispiel:
                     # indices <- which(is.na(mpg)|is.na(hp)|is.na(wt))
                     # res <- fillmd1(lm$residuals, indices])
                     #
   n<-length(indices)
   ii<-1             # nächste Position in indices
   
                     # Verschiebe x um eine Position nach hinten, wenn ein NA eingefügt werden muss
   while (ii<=n)
    {   lx<-length(x)
        if (indices[ii]<=lx) x[indices[ii]:lx]->x[(indices[ii]+1):(lx+1)]
        x[indices[ii]]<-NA
        ii<-ii+1}
   x 
}