Boxplot and mean for selected range of temporal datasets

Question

Adi Purwandana 2024-11-5，11:11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2164085-boxplot-and-mean-for-selected-range-of-temporal-datasets

评论： Voss 2024-11-7，1:25

data_w.mat

Hello there,

I have a datasets containing 3 parameters, let's say z_i as depth, HOUR as time (in numeric format), and EP as the parameter I want to check its variability over depth and time. My intention is to get the EP median, min-max, and mean for selected range of z_i for each time, and then plot them all in a box plot and mean (in one plot), with x-axis as time and y-plot is the EP (better plotted in log-10 scale as the difference is very small).

See my file attached and the illustration of the flow process as I described above below:

Thank you!

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Mathieu NOE 2024-11-5，12:37

hello

have you started something ? what issue are you facing ?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Shishir Reddy 2024-11-6，9:13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2164085-boxplot-and-mean-for-selected-range-of-temporal-datasets#answer_1541680

在 MATLAB Online 中打开

Hi Adi

As per my understanding, you would like to filter a dataset based on a specified depth range (‘z_i’) and calculate statistical measures (median, min, max, mean) of the parameter ‘EP’ for each hour (‘HOUR’).

From the ‘data_w.mat’ file, I see that ‘EP’ and ‘HOUR’ are matrices with 6 rows (possibly representing different observations or repetitions) and 9991 columns , while ‘z_i’ is a vector with 9991 elements, representing depth.

You can process this data in MATLAB as follows –

1. Load the data

load('data_w.mat'); % This loads EP, HOUR, and z_i

2. Filter the data

filtered_indices = (z_i >= z_min) & (z_i <= z_max); % z_min, z_max, depends on the selected range
filtered_EP = EP(:, filtered_indices);
filtered_HOUR = HOUR(:, filtered_indices);
filtered_EP = filtered_EP(:);% Flattening the matrices
filtered_HOUR = filtered_HOUR(:);

3. Compute the statistics

% Assumption – ‘unique_hours’, ‘medians’, ‘mins’, ‘maxs’, ‘means’ are declared.
for i = 1:length(unique_hours)
    hour_data = filtered_EP(filtered_HOUR == unique_hours(i));    
    medians(i) = median(hour_data);    
    mins(i) = min(hour_data);    
    maxs(i) = max(hour_data);    
    means(i) = mean(hour_data);    
end

4. Plotting

figure;
hold on;
boxplot(filtered_EP, filtered_HOUR, 'Colors', [0.7 0.7 0.7], 'Symbol', '');
plot(unique_hours + 1, means, 'ro-', 'LineWidth', 1.5, 'DisplayName', 'Mean');
set(gca, 'YScale', 'log');

For more information regarding the ‘boxplot’ function, kindly refer to the following documentation - https://www.mathworks.com/help/stats/boxplot.html

I hope this helps.

8 个评论
显示 6更早的评论隐藏 6更早的评论

Voss 2024-11-6，16:06

编辑：Voss 2024-11-6，22:52

在 MATLAB Online 中打开

data_w.mat

load('data_w.mat'); % This loads EP, HOUR, and z_i

z_min = min(z_i);

z_max = max(z_i);

filtered_indices = (z_i >= z_min) & (z_i <= z_max); % z_min, z_max, depends on the selected range

filtered_EP = EP(:, filtered_indices);

filtered_HOUR = HOUR(:, filtered_indices);

filtered_EP = filtered_EP(:);% Flattening the matrices

filtered_HOUR = filtered_HOUR(:);

unique_hours = unique(filtered_HOUR);

unique_hours(isnan(unique_hours)) = [];

NH = numel(unique_hours);

medians = zeros(1,NH);

mins = zeros(1,NH);

maxs = zeros(1,NH);

means = zeros(1,NH);

exp_log_means = zeros(1,NH);

for i = 1:NH

hour_data = filtered_EP(filtered_HOUR == unique_hours(i));

medians(i) = median(hour_data);

mins(i) = min(hour_data);

maxs(i) = max(hour_data);

means(i) = mean(hour_data,'omitnan');

exp_log_means(i) = exp(mean(log(hour_data),'omitnan'));

end

figure;

hold on;

boxplot(filtered_EP, filtered_HOUR);

plot(1:NH, means, 'go-', 'LineWidth', 1.5);

plot(1:NH, exp_log_means, 'ko-', 'LineWidth', 1.5);

set(gca, 'YScale', 'log');

[hh,mm,ss] = hms(hours(unique_hours));

xticklabels(string(datetime(0,0,0,hh,mm,round(ss),'Format','HH:mm:ss')))

Adi Purwandana 2024-11-7，0:59

在 MATLAB Online 中打开

Solved by adding this lines:

datetimeDate = datetime(unique_hours, 'ConvertFrom', 'datenum','Format','dd/MM/yy HH:mm');
xticklabels(string(datetimeDate))

Voss 2024-11-7，1:25

Ah, so HOUR is actually days.

请先登录，再进行评论。

Boxplot and mean for selected range of temporal datasets

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

回答（1 个）

8 个评论
显示 6更早的评论隐藏 6更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Boxplot and mean for selected range of temporal datasets

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

回答（1 个）

8 个评论 显示 6更早的评论隐藏 6更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

8 个评论
显示 6更早的评论隐藏 6更早的评论