-1

我正在尝试使用 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) 

在此处输入图像描述

传说显然是错误的。它有一条红色的虚线,主图中没有。主图有一条浅蓝色线,图例中没有。问题出在哪里?谢谢你的帮助。

4

1 回答 1

-1

我不确定,但我的猜测是您的代码legend(... lty=1:4, col=1:4)只会更改图例中的颜色和线型,而不是绘图本身。我也会尝试将该代码放入plot()函数中。这是为了确保 R 知道collty值在绘图函数和图例函数中同步。

于 2015-06-22T17:02:15.263 回答