r - R で凡例が正しくない理由

Question

Rのmsgpsパッケージを使用して、変更されたアイリスデータセットで回帰を試みています:

iris2 = iris
iris2$Species = as.numeric(iris2$Species)
str(iris2)
    'data.frame':   150 obs. of  5 variables:
     $ Sepal.Length: num  5.1 4.9 4.7 4.6 5 5.4 4.6 5 4.4 4.9 ...
     $ Sepal.Width : num  3.5 3 3.2 3.1 3.6 3.9 3.4 3.4 2.9 3.1 ...
     $ Petal.Length: num  1.4 1.4 1.3 1.5 1.4 1.7 1.4 1.5 1.4 1.5 ...
     $ Petal.Width : num  0.2 0.2 0.2 0.2 0.2 0.4 0.3 0.2 0.2 0.1 ...
     $ Species     : num  1 1 1 1 1 1 1 1 1 1 ...

library(msgps)
fit = msgps(as.matrix(iris2[1:4]), iris2$Species)
plot(fit, criteria='cp')
legend('topleft', names(iris2[1:4]), lty=1:4, col=1:4) 

伝説は明らかに間違っています。メインのグラフにはない赤い破線があります。メイン グラフには、凡例にはない水色の線があります。問題はどこだ？ご協力いただきありがとうございます。