我在一张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;但是生成的线条呈锯齿状,看起来像一条移动平均线。请帮帮忙?