std for a growing data range

Question

mick strife 2014-4-25

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127206-std-for-a-growing-data-range

回答： mick strife 2014-4-26

采纳的回答： Roger Stafford

在 MATLAB Online 中打开

Hello :),

i am looking for an effective way to calculate following problem:

for i=1:1000
    tmp=data(1:i,:);
    rsStd(i)=std(tmp(:));
end

"data" is a matrix with k rows and j columns. What i dont want is to calculate the std for every column. What i would need is the std for a growing data range. By each loop step the data range for the calculation of the std grows as you can see in the code snippet.

Does anyone know probably a faster solution and could give me a hint? I am thinking something with vectorization but currently i have no idea. Many thanks!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Roger Stafford 2014-4-25

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127206-std-for-a-growing-data-range#answer_134621

编辑：Roger Stafford 2014-4-25

在 MATLAB Online 中打开

 n = (1:size(data,1))'*size(data,2);
 sq = cumsum(sum(data.^2,2));
 sm = cumsum(sum(data,2));
 rsStd = sqrt((sq-sm.^2./n)./(n-1)); % <-- Corrected

(I'm assuming you want the unbiased version of standard deviation.)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Jos (10584) 2014-4-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/127206-std-for-a-growing-data-range#answer_134622

编辑：Jos (10584) 2014-4-25

在 MATLAB Online 中打开

Here is the completely vectored code for the running sample standard deviation across the rows of a matrix X:

X = ceil(10*rand(13,4)) % example data
N = cumsum(ones(size(X)),1)
S1 = cumsum(X,1) 
S2 = cumsum(X.^2,1) 
R = sqrt((N.*S2 - S1.^2) ./ (N.*(N-1)))
% this can be a one-liner if you insist ..
% check
k = 4 ; [R(k,:) ; std(X(1:k,:))] % the same