Невозможно добавить линейную линию тренда на баррель в 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) 

EDIT

Первая копия этой функции

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) 
} 
} 

Затем сделать barplot

barplot2(Year,Mydata) 

Затем добавьте ablineclip из plotrix библиотеки

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