how to put an if statement within for loop
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
Roisin Loughnane
el 4 de Oct. de 2017
Comentada: Jan
el 4 de Oct. de 2017
I compile climate data over many years by month within a for loop, and concatenate these to create a 3D array (of say January's over many years), and I want to take the element-wise mean to get the average over those years:
monthlydata = cell(1, numel(years)); %to hold the monthly data of each year. Preallocated.
years = 2003:2014;
for month = 1:12
for y = 1:numel(years)
monthlydata{y} = load(sprintf('sst_%d_%d.asc', years(y), month));
end
monthlydatac = cat(3, monthlydata{:}); %concatenate all cells along 3rd dimension
monthlymean = mean(monthlydatac, 3); %get mean along 3 dimensions
dlmwrite(sprintf('msst_%d.csv', month), monthlymean); %save
end
The problem is there is missing data for some months taking the form of numbers 250, 251, 252, 253, 254, 255. In some years observations occur in elements of the 3D array but in other years data is missing ( coded as 250:255) for that element. If this is the case I want to take the element-wise mean excluding these 'missing data' values. Would an 'if' statement work? Could anyone advise me on how to do this?
I am new to Matlab and programming and would be very greatly for your help.
4 comentarios
Jan
el 4 de Oct. de 2017
No, this does not help. What is "coded as one of those numbers 250:255"? How could we suggest an explicit solution based on the detail, that "some" elements of "some" years "may" have a value, that you "would like" to keep?
An exact, unequivocal, clear and unique definition is required, to implement this. We cannot do this for you. I assume the actual implementation will be very easy then.
Respuesta aceptada
Rik
el 4 de Oct. de 2017
Time for a pure guess: you want to ignore values in the 250-255 range and calculate the mean over the remaining data?
monthlydatac(monthlydatac>=250 & monthlydatac<=255)=NaN;
monthlymean = mean(monthlydatac, 3, 'omitnan');
2 comentarios
Más respuestas (0)
Ver también
Categorías
Más información sobre Multidimensional Arrays en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!