MATLAB: How to calculate yearly average value for large data set with uneven time spacing

data analysis

I want to calculate yearly average value of 'd18O PDB' from the following data set (1984-1707):
<ftp://ftp.ncdc.noaa.gov/pub/data/paleo/coral/east_pacific/secas-10yr-iso.txt>
I have downloaded the text file into 2723 times 2 array and changed ('floored') the year to a whole number but I don't know how to go from there. Since the time spacing is uneven, it is hard to figure out how to solve the average value efficiently.

Best Answer

I couldn’t get the ftp call to work (I would not object to your sharing that since you obviously did), so I just copied the text file.
fidi = fopen('secas-10yr-iso.txt', 'rt');
Data = textscan(fidi, '%f%f%f%f', 'HeaderLines',4, 'CollectOutput',1);
fclose(fidi);
data = cell2mat(Data);
d18O_PDB = data(:,4);
Date_Data = [data(:,2) data(:,4)];
Y1 = Date_Data(end,1);                                          % Start Year
years = floor(Date_Data(:,1)-Date_Data(end,1))+1;               % Create Indices By Year
Yu = unique(years);                                             % Unique Year Indices
mean_d18O_PDB = accumarray(years, Date_Data(:,2), [], @mean);   % Calculate Annual Means
Output = flipud([Yu+floor(Y1)  mean_d18O_PDB]);                 % Matrix Of Means By Year
The flipud call is necessary because the first year is the first index, and flipping the matrix puts the years in the order they are in the original file.
EDIT: Added the ‘Date_Data’ (2723 x 2) array to correspond with the array in the Question.
Related Question