我有这个代码用于绘制x,y并显示回归线和拟合统计数据:
for (i in unique(Data$SeasAlltxt)){
print (i)
subdata1 <- subset(Data, SeasAlltxt == i)
for (j in unique(Data$ALSCIDtxtall)){
subdata2 <- subset(subdata1, ALSCIDtxtall == j)
plot(subdata2$PeakToGone, subdata2$NO3_AVG, xlim = c(0, 70))
abline(fit <- lm(NO3_AVG~PeakToGone, data = subdata2))
rmse <- round(sqrt(mean(resid(fit)^2)), 2)
r2 <- round(summary(fit)$r.squared, 3)
ycoord <- max(subdata2$NO3_AVG)
eqn <- bquote(~~ r^2 == .(r2) * "," ~~ RMSE == .(rmse))
text(40, ycoord, eqn, pos = 4)
}
}
这适用于大多数情节,并在情节上显示文字。但是,有些情节不显示文字。这些图的轴包含文本的原点坐标(40,ycoord)的值,因此我对文本没有出现的原因感到困惑。我尝试将这些坐标更改为绘图的中心,并且文本仍未显示在相同的图上。
有什么想法?