I’ve got a NetAtmo weather station. One can download the measurements from its web interface as a CSV file. I wanted to give time series analysis with the extraction of seasonal components (‘decomposition’) a try, so I thought it would be a good opportunity to use the temperature measurements of my weather station. The data is available.

To make things visually a little easier, I only tried this with 14 days of temperature measurements, including all measurements from November 1st, 2015 till November 14th, 2015.

The raw data looks like this (on the x-axis, there is a running number of measurements):

data <- read.table("temps.csv", header = T, sep = "t")
plot(data$temp, type = “l”, bty = “n”, ylab = “Temperature in °C”)

Where does the 285 come from? The NetAtmo station takes one measurement every 5 minutes. So, in one hour, there should be 60 / 5 = 12 measurements and in one day there should be 60 / 5 * 24 = 288 measurements. However, if we check table(data$date), we see that most of the days, the station made 3 measurements less per day, so I use this as my frequency.

In the last step, we decompose and plot the time series.

plot(stl(data$temp, s.window = “periodic”))

The plot is divided into the raw data (again), the seasonal component as extracted by the function stl, the overall trend of the time series and the differences between trend and data. Nice.