R - 拟合平滑曲线,PDF显示

时间:2014-11-12 22:14:44

标签: r curve-fitting

我在一张PDF中有几个散点图,每个图显示Gage A(感兴趣的量具)和参考量具之间的流量值。目的是找到与Gage A具有最佳相关性的参考量具。相关性可以是直线或曲线。

RefGage.list =与Gage A比较的参考量表列表

Tables.list =包含Gage A流量值的表名称列表和参考量表的并发流量值。例如,如果有10个参考量具,Tables.list将有10个文件名,每个参考量具有一个表/文件。

以下是我创建PDF的代码:

pdf(file=Plot.file, width=8.5, height=11)
par(mfrow=c(4,4), oma=c(1,1,6,1))
Test <- lapply(1:length(RefGage.list), function(z) {
   a <- as.data.table(read.csv(Tables.list[z], header = TRUE))
   plot(GageAFlow ~ RefGageFlow, data = a, sub = paste("R2 =", summary(lm(GageAFlow ~ RefGageFlow, a))$adj.r.squared))
   abline(lm(GageAFlow ~ RefGageFlow, a), col = "red")
  })
  title("Streamflow", line=3, cex.main=1.1, outer=TRUE)
  dev.off()

我想做的是创建像Excel这样的东西,我可以显示拟合的直线和拟合的平滑曲线(多项式),然后显示两条线的R2。上面的代码仅用于显示直线(abline),并且我在副标题中显示了R2。不是最好的地方,但我不知道其他任何方式。 (对不起,还在学习R)。我尝试了黄土,但曲线不是那么简单和平滑。我使用了&#34; lm($ PRflow~poly($ Indexflow,2))&#34;但是生成的线条呈锯齿状,看起来像一条移动平均线。请帮帮忙?

0 个答案:

没有答案