MATLAB Answers

How to convert a time serie with daily data to a time serie with weekly data by taking the average of every 7 values of the daily vector?

2 views (last 30 days)
I have a vector with daily data. I need to average every 7 values of it to obtain a new vector with weekly data.

  2 Comments

Fernanda Suarez Jaimes
Fernanda Suarez Jaimes on 10 Mar 2020
Sure!
%Time Vector
xtime1= [1:1:6000]'; %daily
xtime2= [1:1:42000]'; %weekly
%Time Series 1- daily entries
%comand randn generats 42000 random numbers that are normaly distributed.
rng('default')
time_series1=randn(42000,1);
%generating time series with 6000 entries log-normal distributed
%the code for this distribution is based on the code available at
%mathworks.com
rng('default'); % So that numbers can be repeated
time_series2 = lognrnd(0,0.25, [6000 1]); %generating time series with mu set to zero and sigma 0.25

Sign in to comment.

Accepted Answer

Ameer Hamza
Ameer Hamza on 11 Mar 2020
As you mentioned in your comment, you have a vector of type double. In that case, you can average 7 elements as follow:
x = rand(1000,1); % random data
x = padarray(x, [ceil(numel(x)/7)*7 - numel(x) 0], 0, 'post'); % make sure that number of elements are multiple of 7
mean_x = mean(reshape(x, 7, []), 1)';

  3 Comments

Fernanda Suarez Jaimes
Fernanda Suarez Jaimes on 12 Mar 2020
Is the second x a new vector? or am I re arranging the first x vector?
Is this code correct?
rng('default')
time_series1 = normrnd(0,1,42000,1); %daily time series multiple of 7, 42000
time_series3 = padarray (time_series1, [ceil (numel (time_series1) / 7) * 7 - numel (time_series1) 0], 0, 'post' ); %
mean_ts3 = mean (reshape (time_series3, 7, []), 1);
Ameer Hamza
Ameer Hamza on 12 Mar 2020
Fernanda, the second x does not need to be the same, your code is correct. But you will get an error because you are giving extra spaces between function name and the input arguments. This is generally not a problem, but it will cause an error when done inside [ ] brackets. So change the third line
time_series3 = padarray (time_series1, [ceil(numel (time_series1) / 7) * 7 - numel(time_series1) 0], 0, 'post' );

Sign in to comment.

More Answers (1)

Stephan
Stephan on 10 Mar 2020
Edited: Stephan on 10 Mar 2020
The usual way to do this is by using the retime function.