在图中添加最小二乘线和LMS线

Question

我已经从软件包library（rounustbase）中加载了数据集andime2，并且我有兴趣与这些数据的对数合作。

library(robustbase)
x<-log(Animals2)
plot(x, main="Plot Animals2", col="darkgreen")

现在我需要添加最小二乘线，LMS线和线路 西格尔（Siegel），使用不同的颜色。最后，我还想展示两个估计的二维密度的图。第一个图应是密度的3D可视化（我试图使用命令persp），第二个图与命令image（应用于 密度估计）。

我有点困扰，也很感激。

Answer 1

假设mblm（）可以使用Siegel模型，则可以使用以下方式：

  library(robustbase)
  library(mblm)
  data(Animals2)
  x <- log(Animals2)
  plot(x, main="Plot Animals2", col="darkgreen")
  abline(lm(brain ~ body, data=x), col="red")
  abline(MASS::lqs(brain ~ body, data=x, method="lms"), col="blue")
  abline(mblm(brain ~ body, data=x, repeated=TRUE), col="black")
  legend("topleft", c("LM", "LMS", "Siegel"), 
         col = c("red", "blue", "black"), 
         lty = c(1,1,1), 
         inset=.01)

d2 <- MASS::kde2d(x$body, x$brain)
persp(d2)

image(d2, xlab="body", ylab="brain")

” -05-20 by reprex package （v2.0.1）