清空筛选条件
清空筛选条件

Calculate standard deviation for different groups

10 次查看（过去 30 天）
leonidas86
leonidas86 2018-7-23
评论： Guillaume 2018-7-23
I want to plot a Measurement with a bar figure. Therefore I want to create different groups and calculate the standard deviation for each group. I have a vector "Measurement" which is sorted according to increasing values. I get group IDs for the Measurement values by the use of the histc function.
[~,Measure_ID]=histc(Measurement_x,0:200:max);
In the next step I try to calculate the standard deviation with the accumarray function.
group_std = accumarray(Measure_ID,Measurement_y,[],@std);
The std command uses the follow mathematical formula:
The problem is that accumarray calculates the arithmetic average for each group separately. But I want to see the standard deviation regarding to the global arithmetic average. Is there a way to solve my issue?
Version: 2012b
  1 个评论
Jan
Jan 2018-7-23
编辑：Jan 2018-7-23
Do you mean:
std(Measurement_y)
? What exactly is "standard deviation regarding to the global arithmetic average"?

请先登录，再进行评论。

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

采纳的回答

Jan
Jan 2018-7-23
编辑：Jan 2018-7-23
If you want to set the x_bar values of the formula to mean(Measurement_y), why not implementing this as function?
function s = specialStd(x, m)
s = sqrt(sum((x - m) .^ 2) / (length(x) - 1));
end
and call it like:
m = mean(Measurement_x);
group_std = accumarray(Measure_ID, Measurement_y, [], @(x) specialStd(x, m));
  1 个评论
Guillaume
Guillaume 2018-7-23
specialStd, notexactlystd, I think we can agree that this thing shouldn't be called plain std!

请先登录，再进行评论。

更多回答（1 个）

Guillaume
Guillaume 2018-7-23
"The problem is that accumarray calculates the arithmetic average for each group separately".
It's not accumarray doing that, it's std, because that the definition of the standard eviation.
To do what you want, you'll have to implement your own standard deviation calculation:
whole_mean = mean(Mesurement_y);
group_notexactlystd = accumarray(Measure_ID, Measurement_y, [], @(v) sum((v - whole_mean).^2)/(numel(v)-1));

类别

在 Help Center 和 File Exchange 中查找有关 Surfaces, Volumes, and Polygons 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by