r - 无法在 R 中的条形图上添加线性趋势线

Question

我创建了一个条形图barplot，然后我想显示线性趋势。我使用abline但图中没有显示线性趋势线。我想知道问题是什么。谢谢。

set.seed(100)
Mydata=rnorm(65)
Year=1950:2014
barplot(Mydata)
fit=lm(Mydata~Year)
abline(fit)

正如@G5W 指出的那样，fit=lm(Mydata~I(Year-1950)). 但新的问题是趋势线太“长”。如第二张图所示，趋势线超出了条形图。有什么建议吗？

Answer 1

如果你可以使用ggplot：

library(ggplot2)    

df <- data.frame(Mydata, Year)

ggplot(df, aes(x = Year, y = Mydata)) + 
  geom_bar(stat = "identity") + 
  geom_smooth(method = "lm")

Answer 2

要扩展 @bouncyball 的评论，如果需要，可以使用更高的线宽 ( lwd) 值来模仿。barplot

plot(Year, Mydata, type = 'h',lwd=5,col = "grey")
abline(fit, lty =2)

编辑

先复制这个函数

barplot2 <- function(x, y, lty = 1, lwd = 1, col = "grey", border = "black"){
w = ((max(x) - min(x))/length(x)) * 0.75
plot(x, y, type = 'p', pch = NA, yaxt = "n", xaxt = "n", xlab = "", ylab = "")
for (i in 1:length(x)){
x1  = x[i] - w/2
x2 = x[i] + w/2
y1 = 0
y2 = y[i]
polygon(x = c(x1,x2,x2,x1), y = c(y1,y1,y2,y2), lty = lty, lwd = lwd, col = col, border = border)
}
}

然后制作条形图

barplot2(Year,Mydata)

然后添加ablineclip从plotrix库

ablineclip(fit, x1 = min(Year), x2 = max(Year), y1 = min(Mydata), y2 = max(Mydata))