Question

我试图在一个图中绘制两个locfit模型但是我无法得到第二个locfit模型来绘制置信区间。我创建了两个locfit模型：

1_fit = locfit(Y~Time,data=data_1)
2_fit = locfit(Y~Time,data=data_2)

每个都可以使用以下内容自行绘制，并使用以下内容进行绘制：

plot(1_fit,band="local",type = "l", xlab = "Time", ylab = "Y-Axis",ylim=c(0,22), 
  col = "red",lwd = 5,font=3,main="Local Poly Fit 1",cex.lab=1.5, cex.axis=1.5, 
  cex.main=1.5, cex.sub=1.5)

然而，当我尝试使用以下方法将另一个locfit模型绘制到图中时：

lines(2_fit,col="blue")

我只能添加locfit线而不是置信区间。我试图这样做：

lines(2_fit,band="local",col="blue")

但是我收到了这条消息而没有置信区间：

Warning message: In plot.xy(xy.coords(x, y), type = type, ...) : "band" is not a graphical parameter

我也考虑过使用lines.locfit，但没有运气，因为R只是说它无法找到功能线.locfit。

我有一个方法是使用以下方法将两个图放在同一个窗口中：

par(mfrow=c(2,1))

但是我想避免这种情况，因为如果它们在同一个情节中，它会使这些情节更具可比性。

Answer 1

答案由Richard Telford在评论中提供。以下代码能够实现我所需要的：

    plot(1_fit,band="local",type = "l", xlab = "Time", ylab = "Y-Axis",ylim=c(0,22), xlim=c(0,12), col = "red",lwd = 5,font=3,main="Local Poly Fit 1",cex.lab=1.5, cex.axis=1.5, cex.main=1.5, cex.sub=1.5)`

    par(new = TRUE)

    plot(2_fit,band="local",type = "l", xlab = "Time", ylab = "Y-Axis",ylim=c(0,22), xlim=c(0,12),col = "blue",lwd = 5,font=3,main="Local Poly Fit 1",cex.lab=1.5, cex.axis=1.5, cex.main=1.5, cex.sub=1.5)`

我需要确保ylim和xlim以及main，ylab和xlab相同。同样来自理查德的旁注，1_fit不是一个合法的名称，我在这里使用它只是作为占位符名称，但似乎很好的传递知识。

R绘制具有置信区间的多个locfit模型

1 个答案: