r - 将x = y行添加到散点图

我正在使用汽车包装中的散点图函数来生成散点图。
我希望能够在绘图中生成应该为x = y的引用线。
我尝试使用abline，它确实添加了一行，但它不是x = y行。有人可以帮忙吗？

我的代码如下:

scatterplot(phenos$P1~pheno$P0, data=pheno,spread=FALSE,ylab="6 month timepoint", xlab="Baseline Timepoint", jitter=list(x=1, y=1))
abline(0,1)

谢谢。

最佳答案

实际上，这是相当困难的/骇人听闻的，因为scatterplot()在内部使用layout，这使得很难控制图形驱动程序当前正在使用的子图。 (更新:这比我想象的要难-设置par("mfg")一定是偶然地起作用了。)

组成数据(更新:使用均值x和y不等于零且彼此不相等的数据，因为它更清楚地说明了单纯使用abline()的困难)

set.seed(1)
d <- data.frame(x=rnorm(10,mean=10),y=rnorm(10,mean=12))
library(car)

尝试我的旧策略(该策略实际上不起作用，或者只能以不可预测的方式起作用):

scatterplot(y~x,data=d,reset.par=FALSE)
k <- 1
for (i in 1:2) {
   for (j in 1:2) {
      par(mfg=c(i,j,2,2))
        abline(0,1,lwd=3,col=k)
        k <- k+1
  }

}

根据我的操作方式，我会收到警告和错误，也可能会得到虚假的答案。我是否在函数内执行scatterplot()似乎很重要... ??

第二次尝试，更加保守:从头开始重新构建布局。

 scatterplot(y~x,data=d)
 uu <- par("usr")
 ## mimic layout frolm car:::scatterplot.default.  Would be different if we were drawing only one
 ## of x-boxes or y-boxes
 layout(matrix(c(1, 0, 3, 2), 2, 2), widths = c(5, 95),
        heights = c(95, 5))
 oldmar <- par(mar=rep(0,4))  ## zero out margins so we can plot in sub-boxes without errors
 ## now skip through the first two sub-plots
 par(new=TRUE); plot.new(); par(new=TRUE); plot.new()
 par(oldmar)  ## reset margins
 ## blank plot with user limits set and 'interior' axis calculations
 plot(0:1,0:1,xlab="",ylab="",xlim=uu[1:2],ylim=uu[3:4],xaxs="i",yaxs="i")
 ## add annotation
 abline(a=0,b=1,col=4,lwd=3)

考虑到此解决方案的工作量和脆弱性，实际上最好是破解scatterplot以有选择地允许另外指定abline()或向维护者索要该功能……