我正在尝试使用 R 中的 msgps 包对修改后的 iris 数据集进行回归:
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)
传说显然是错误的。它有一条红色的虚线,主图中没有。主图有一条浅蓝色线,图例中没有。问题出在哪里?谢谢你的帮助。