如何添加一个控件线条-不同的地块会有所不同？

Question

我有一个庞大的数据集，给出了许多治疗的多个参数的数据，其中一个是对照。

我已经写了一个脚本，允许我生成治疗$参数的多个散点图，现在我希望添加一个空白线(h=“控件”)，但我一直收到错误。

“对照”是许多处理方法之一，在数据集中位于“处理方法”列下。我是不是把它合并错了？

我试过but (h=treatment==“控制”)，但没有成功。另外，任何关于如何思考这些问题的建议-我不想把它们称为翻译问题-将非常感谢，并将有助于我的自我故障排除之旅。

编辑：

例如:对于数据(“USArrests”)，我的脚本已经为所有列(谋杀、袭击等)生成了散点图；对于我的数据框架，有几个数据点/参数)。如果我试图用USArrests做类似的事情，我会想要一个加州$参数的基准测试行。

提前感谢！

Answer 1

如果使用abline(h="control")或abline(h=treatment=="control")，它将不起作用，因为第一个是尝试绘制字符，第二个是布尔值。(h=..)需要数字输入。

由于您知道来自哪一行，因此可以执行以下操作：

library(datasets)
wh = which(rownames(USArrests)=="California")

如果使用配对，则可以像这样绘制：

plotfun = function(x,y){
points(x,y,pch=18);abline(h=y[wh],lty=8,col="blue")
}
pairs(USArrests,upper.panel=NULL,panel=plotfun)

​

​

如果你有一个for循环：

PAIRS = combn(colnames(USArrests),2)
par(mfrow=c(2,3))
for(i in 1:ncol(PAIRS)){
plot(USArrests[,PAIRS[,i]])
INTERCEPT = USArrests[wh,PAIRS[2,i]]
abline(h=INTERCEPT,lty=8,col="blue")
}

​

​