在线图中显示黄土平滑趋势线，＃ggplot2。

Question

如何在情节中显示黄土平滑的趋势线？请帮助处理警告消息：“删除了包含非有限值的19行（stat_smooth）”。

我的数据：

yrcnt<-read.table(header = TRUE, text = "year outcome pop rate pred.SC
1  1995    2306  87592001 2.632660 0.9626214
2  1996    2221  87628543 2.534562 0.9599941
3  1997    2202  81872629 2.689544 0.9573667
4  1998    2316  88200076 2.625848 0.9547394
5  1999    2456  96200312 2.553006 0.9521121
6  2000    2526  99565063 2.537035 0.9494848
7  2001    2511  95951330 2.616952 0.9468575
8  2002    2537  96976191 2.616106 0.9442302
9  2003    2618 101673130 2.574918 0.9416028
10 2004    2644 104554479 2.528825 0.9389755
11 2005    2594 100522055 2.580528 0.9363482
12 2006    2620 105787278 2.476668 0.9337209
13 2007    2722 108946407 2.498476 0.9310936
14 2008    2788 112200567 2.484836 0.9284663
15 2009    2706 104491560 2.589683 0.9258389
16 2010    2773 108651896 2.552187 0.9232116
17 2011    2764 109632577 2.521148 0.9205843
18 2012    2694 107594922 2.503836 0.9179570
19 2013    2673 107553219 2.485281 0.9153297")

http://tutorials.iq.harvard.edu/R/Rgraphics/Rgraphics.html 

我的代码：

    p1 <- ggplot(yrcnt, aes(y = log(rate), x = year))
    yrcnt$pred.SC <- predict(lm(year ~ log(rate), data = yrcnt))
    p1 + geom_line(aes(color = rate)) +geom_line(aes(y = pred.SC))
    p1 + geom_line(aes(color = rate)) + geom_smooth()
    p4 <- p1 + geom_line(aes(color = rate)) + geom_smooth(color="red")
    p4 + scale_x_continuous(name = "Years",limits = c(1995, 2013),breaks = 1994:2014) +
      scale_y_continuous(name = "Pancreatic Cancer Hospitalization Rate, 1995-2013",limits = c(2.4, 2.7),breaks = seq(2.4, 2.7, by = 0.1)) +
      ggtitle("Long Term Trend in Pancreatic Cancer Hospitalizations")
>`geom_smooth()` using method = 'loess'
p4=Base plot with trendline

无法合并到基础图= p4

的比例图

Answer 1

在对scale_y_continuous()的函数调用中，删除参数

                       limits = c(2.4, 2.7), 
                       breaks = seq(2.4, 2.7, 
                                    by = 0.1)) 

因为y轴的真实极限在0.9和1之间，并且您将它们设置在2.4和2.7之间的范围内。我不知道你是否需要费率或日志（费率）。

另一种选择是

library('ggplot2')
p1 <- ggplot(yrcnt, aes(y = rate, x = year))

######### lm() args flipped, then 
######### wrapped in exp() function.
yrcnt$pred.SC <- exp(predict(lm( log(rate) ~ year, data = yrcnt))) 

p1 + geom_line(aes(color = rate)) +geom_line(aes(y = pred.SC))

p4 <- p1 + geom_line(aes(color = rate)) + geom_smooth(color="red", method="loess")
p4 + scale_x_continuous(name = "Years",limits = c(1995, 2013),breaks = 1994:2014) +
        scale_y_continuous(name = "Pancreatic Cancer Hospitalization Rate, 1995-2013",limits = c(2.4, 2.7),breaks = seq(2.4, 2.7, by = 0.1)) +
        theme(legend.position ="none") +
        ggtitle("Long Term Trend in Pancreatic Cancer Hospitalizations")