Plotting abline on time series with quarterly entries

[dmd560]

I plotted GDP data that is recorded quarterly by turning it into a time series object.

library(ggplot2)
library(ggfortify)

Y <- ts(GDP$GDP, start(1947,1), frequency = 4)

Y <- ts(GDP$GDP, start(1947,1), frequency = 4)
autoplot(GDP$GDP) +
ggtitle("Time Plot: Quarterly GDP") +
ylab("GDP")

I can't seem to figure out how to plot an abline onto this graph.

[FJCC]

Can you please include code to produce a small subset of the GDP data frame?

[dmd560]

Is this what you mean?

head(GDP)

A tibble: 6 x 2

DATE GDP
chr dbl
1 1/1/1947 2033.
2 4/1/1947 2028.
3 7/1/1947 2023.
4 10/1/1947 2055.
5 1/1/1948 2086.
6 4/1/1948 2120.

[nirgrahamuk]

try

dput(head(GDP))

[dmd560]

This is what it returns

dput(head(GDP))
structure(list(DATE = c("1/1/1947", "4/1/1947", "7/1/1947", "10/1/1947",
"1/1/1948", "4/1/1948"), GDP = c(2033.061, 2027.639, 2023.452,
2055.103, 2086.017, 2120.45)), class = c("tbl_df", "tbl", "data.frame"
), row.names = c(NA, -6L))

[FJCC]

I could not get autoplot to work with the GDP data frame. I have never used autoplot but it seems you must define a method for that kind of object. I did plot the time series Y and found that the x values are large negative numbers because the data are plotted as being in the year 1.

GDP <- structure(list(DATE = c("1/1/1947", "4/1/1947", "7/1/1947", "10/1/1947",
                        "1/1/1948", "4/1/1948"), 
               GDP = c(2033.061, 2027.639, 2023.452, 2055.103, 2086.017, 2120.45)), 
          class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, -6L))

library(ggplot2)
library(ggfortify)

Y <- ts(GDP$GDP, start(1947,1), frequency = 4)

PLT <- autoplot(Y) +
         ggtitle("Time Plot: Quarterly GDP") +
         ylab("GDP") 

#The dates are in year 1
PLT$data$Index
#> [1] "0001-01-01" "0001-04-01" "0001-07-01" "0001-10-01" "0002-01-01"
#> [6] "0002-04-01"

#Their numeric value is ~ -719162, the number of days between Jan 1, 1970 and Jan 1, 1
as.numeric(PLT$data$Index)
#> [1] -719162 -719072 -718981 -718889 -718797 -718707

autoplot(Y) +
  ggtitle("Time Plot: Quarterly GDP") +
  ylab("GDP") + geom_abline(slope = 1/-719162, intercept = 2035)

^{Created on 2020-02-21 by the reprex package (v0.3.0)}

That may not be very helpful because you seem to be able to plot the GDP data frame. Try storing the result of autoplot in a variable PLT as I did above and look at the value of

as.numeric(PLT$data$Index)

[dmd560]

Yeah, I was getting the same results. The only way I came up to figure it out was to not use the data as a time series object and replace the date with an order like this:

GDP["order"] <- seq(1,290)
# run regression
lr <- lm(log(GDP)~order, data = GDP)
# scatter plot
p <- ggplot(GDP, aes(x=order, y=log(GDP)))+geom_point()
# fitted line plot
p + geom_abline(intercept = lr$coefficients[1], slope = lr$coefficients[2])

Thank you for your help

[dromano]

You don't have to replace dates by indices if you need them:

library(tidyverse)

GDP <- 
  structure(
    list(
      DATE = c("1/1/1947", "4/1/1947", "7/1/1947", "10/1/1947","1/1/1948", "4/1/1948"),
      GDP = c(2033.061, 2027.639, 2023.452, 2055.103, 2086.017, 2120.45)
      ), 
    class = c("tbl_df", "tbl", "data.frame"), 
    row.names = c(NA, -6L)
    )

GDP %>% 
  ggplot(aes(DATE, GDP)) +
  geom_point() +
  geom_abline(slope = 20, intercept = 2025)

^{Created on 2020-02-22 by the reprex package (v0.3.0)}