Perform Time Series Analysis And Forecasting Using R Programming Language

This article illustrates how to perform time-series analysis and forecasting using the R programming language.
time series analysis and forecasting with R

This article illustrates how to perform time-series analysis and forecasting using the R programming language. Time series analysis refers to an important statistical technique for studying the trends and characteristics of collecting data points indexed in chronological order. On the other hand, time series forecasting involves the task of getting insights from recorded time series data and making future predictions based on them. 

If you are unfamiliar with R and its basic concepts, check out this article before proceeding. We have already covered several articles on time series analysis and time series forecasting but with Pythonic code.

Practical implementation using R

The code here has been implemented using RStudio IDE (version 1.2.1335). You can download RStudio from here. Step-wise explanation of the code is as follows:

Time-Series Analysis

  1. Install rmeta R package.


Check if the package has been installed by displaying the whole list of packages   

installed in the library.



 package ‘rmeta’ successfully unpacked and MD5 sums checked
 The downloaded binary packages are in
  1. Read the data containing age of various kings when they died. The dataset is available here.
 kings_data <- scan("",skip=3)


R data1
  1. Store the kings’ data into a time series object for performing time series analysis.
 kings_ts_data <- ts(kings_data) 


R time series1
  1. Similarly, read the New York City (NYC)’s data containing a monthly record of a number of births in the city. (Raw data is available here).
 NYC_data <- scan("")


R data2
  1. Store NYC’s data into an R time series object. ‘frequency’ parameter of the ts() function should be set to 12 for month-wise data. Also, we specify start=c(1950,1)) as a parameter to ts() indicating that the data should start from the year 1950 and there should 1 sample per year for each month.
 NYC_ts_data <- ts(NYC_data, frequency=12, start=c(1950,1))

Sample condensed output:

R time series2
  1. Read the data containing monthly records of a souvenir shop in Australia. Find the data here
 shop_data <- scan("")


R data3
  1. Convert the souvenir shop’s data into a time series object;
 shop_ts_data <- ts(shop_data, frequency=12, start=c(1980,1))


R time series3
  1. Plot the time-series version of all the three datasets (kings’ data, NYC’s data and souvenir shop’s data).



time series plot 1



time series plot 2



time series plot 3
  1. Time series decomposition is a process of decomposing the time series data into components viz. A trend component and an irregular component. A seasonal data additionally has a seasonal component. 

We first smoothen the kings’ time series data using SMA() function of the TTR package (where SMA stands for ‘simple moving average’) for getting the trend component.

 kings_ts_SMA2 <- SMA(kings_ts_data,n=2) 

(‘n’ parameter here specifies the order of SMA)



Similarly, we can change the order to say 8 and observe the change in the trend.

 kings_ts_SMA8 <- SMA(kings_ts_data,n=8)


It can be observed that incrementing the order of SMA smoothens the plot more i.e. reduces fluctuations in the trend.

  1. For seasonal data such as that of births in NYC, the decomposition can be carried out using the decompose() function since it also has a seasonal component apart from the trend and irregular ones.

NYC_ts_comp <- decompose(NYC_ts_data)

Now, we can separate obtain each of the three components of the decomposed NYC data.

Seasonal component:


Sample condensed output:

R time series decomposition 1

Trend component:


Sample condensed output:

R time series decomposition 2

Irregular component:


Sample condensed output:

R time series decomposition 3

Plot all the components of the NYC time series in a single plot.


Sample condensed output:

R time deries decomposition plot
  1. The seasonal time series data of NYC can be seasonally adjusted by subtracting its seasonal component from the original time series.
 NYC_ts_adjusted <- NYC_ts_data - NYC_ts_comp$seasonal


R time series adjusted

Time-Series Forecasting

  1. Read the data containing records of annual precepitation (in inches) in London (data is available here). 
 rain_data <- scan("",skip=1)


time series forecasting data
  1. Create a time series object of the precipitation data for the year starting from 1815 and plot that time-series.
 rain_ts_data <- ts(rain_data, start=c(1815))


time series forecasting plot
  1. We can forecast the rainfall in future using simple exponential smoothing technique since the rain data has no seasonality and is an additive time-series.

 The HolWinters() functions with ‘beta’ and ‘gamma’ parameters set to False can be used for simple exponential smoothing of the data. It estimates the value of smoothing factor ‘alpha’ whose value lies between 0 and 1. Larger the value of alpha, more is the importance/weightage given to the historic data for making predictions.

 rain_fc <- HoltWinters(rain_ts_data, beta=FALSE, gamma=FALSE)


forecast result

The forecasts made by the HoltWinters() function can be obtained from the named variable called ‘filleted’ of the rain_fc variable. It is in this variable that the HoltWinters() function stores its output.



The complete output will show that HolWinters() makes predictions for that period only which is covered by our original time-series.

  1.  Plot the forecasted as well as original time series in a single plot.



The ‘red’ line in the above plot shows the forecasted time series which is, as can be 

seen, has much less fluctuations than the original time series (shown in black).

Download our Mobile App

Nikita Shiledarbaxi
A zealous learner aspiring to advance in the domain of AI/ML. Eager to grasp emerging techniques to get insights from data and hence explore realistic Data Science applications as well.

Subscribe to our newsletter

Join our editors every weekday evening as they steer you through the most significant news of the day.
Your newsletter subscriptions are subject to AIM Privacy Policy and Terms and Conditions.

Our Upcoming Events

15th June | Bangalore

Future Ready | Lead the AI Era Summit

15th June | Online

Building LLM powered applications using LangChain

17th June | Online

Mastering LangChain: A Hands-on Workshop for Building Generative AI Applications

20th June | Bangalore

Women in Data Science (WiDS) by Intuit India

Jun 23, 2023 | Bangalore

MachineCon 2023 India

26th June | Online

Accelerating inference for every workload with TensorRT

MachineCon 2023 USA

Jul 21, 2023 | New York

Cypher 2023

Oct 11-13, 2023 | Bangalore

3 Ways to Join our Community

Telegram group

Discover special offers, top stories, upcoming events, and more.

Discord Server

Stay Connected with a larger ecosystem of data science and ML Professionals

Subscribe to our Daily newsletter

Get our daily awesome stories & videos in your inbox

Can Apple Save Meta?

The iPhone kicked off the smartphone revolution and saved countless companies. Could the Pro Reality headset do the same for Meta?