How to calculate the median for each day of the year in 91 years?

2 views (last 30 days)
Maria Mateus on 9 Oct 2013
Commented: Matt Kindig on 9 Oct 2013
I need to calculate the median value for each day of the year from a period of record that goes from Jan 1st 1915 to Dec 31 2006 (ie: median values for: 01-01-1915, 01-01-1016,...,01-01-2006).
The matrix is organized in the following order:
Year Month Day Streamflow_A1B Streamflow_B1 1915 01 01 1000 1010 1915 01 02 1100 1050 1915 01 03 2000 2100 ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... 2006 12 31 2000 2100
Any help with this will be very much appreciated. I am new at Matlab.
Thanks!

Matt Kindig on 9 Oct 2013
Edited: Matt Kindig on 9 Oct 2013
I think this is a good use of accumarray().
%some sample data, to illustrate approach
Data = [1915 1 1 1000 1010;
1915 1 2 1100 1050;
1915 1 3 2000 2100;
1918 1 1 500 800;
1918 1 2 1000 600;
2006 12 31 2000 2100;
2008 12 15 1800 1920];
%break into variables
Days = Data(:,2:3);
StreamA1B = Data(:,4);
StreamB1 = Data(:,5);
%get unique month/day combinations. Note, not all months will have the same
%number of days! Leap years can also complicate things a bit.
[MonthDays, ~, index] = unique( Days, 'rows');
%calculate median by day
StreamA1B_med = accumarray(index, StreamA1B, [], @median);
SreamB1_med = accumarray(index, StreamB1, [], @median);
Matt Kindig on 9 Oct 2013
No, this should give you all days of the month. If you look at the MonthDays variable, you'll see that all month/day combinations are covered.

Maria Mateus on 9 Oct 2013
Never mind. I got it. Thanks!

Categories

Find more on Resizing and Reshaping Matrices in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by