myplot<-function(x,y,g,...)
#
#  Punktewolke von x und y, wobei die Zugehörigkeit zu den Gruppen von g mit verschiedenen
#    Farben gekennzeichnet wird. Dazu noch die Regressionsgerade
#
{
 if (class(g)=="numeric") g<-factor(g)    # Wandeln in factor, um die möglichen Werte von g zu erhalten
 plot(x,y,type="n",...)                   # leerer Rahmen
 c=0
 for (i in levels(g))                     # für jeden Wert von g           
 { c<-c+1 ; points(x[g==i],y[g==i],col=c) #   eine Punktewolke
   abline(lsfit(x[g==i],y[g==i]),col=c) } #   und Regressionsgeraden
                                          # abschließend eine Legende in die obere rechte Ecke
 legend("topright",levels(g),pch=rep(1,c),col=1:c)
}