Moving mean/std based on number of NaN

조회 수: 1 (최근 30일)
Dave
Dave 2021년 8월 28일
답변: Walter Roberson 2021년 8월 28일
Hello,
Is there a way to compute movstd but replace with NaN if the window includes 2 or more NaN?
A=[4, 8,6,3,NaN,NaN,-1,9,4,5]'
This window is 4
M=movstd(A,4,'omitnan','Endpoints','discard')
M is
2.217355783
2.516611478
2.121320344
2.828427125
7.071067812
5
4.11298756
With the condition to replace with NaN if the window includes 2 or more NaN
M should be
2.217355783
2.516611478
NaN
NaN
NaN
5
4.11298756

채택된 답변

Matt J
Matt J 2021년 8월 28일
편집: Matt J 2021년 8월 28일
A=[4, 8,6,3,NaN,NaN,-1,9,4,5]';
M=movstd(A,4,'omitnan','Endpoints','discard');
idx=movsum(isnan(A),4,'omitnan','Endpoints','discard')>=2;
M(idx)=nan
M = 7×1
2.2174 2.5166 NaN NaN NaN 5.0000 4.1130

추가 답변 (1개)

Walter Roberson
Walter Roberson 2021년 8월 28일
A=[4, 8,6,3,NaN,NaN,-1,9,4,5]'
A = 10×1
4 8 6 3 NaN NaN -1 9 4 5
window = 4;
M=movstd(A,window,'omitnan','Endpoints','discard')
M = 7×1
2.2174 2.5166 2.1213 2.8284 7.0711 5.0000 4.1130
An = isnan(A.');
starts = strfind([false, An], [false,true,true])
starts = 5
stops = strfind([An,false], [true,true,false])+1
stops = 6
M(unique(cell2mat(arrayfun(@(start,stop) start-window+2:stop-1, starts, stops, 'uniform', 0)))) = nan
M = 7×1
2.2174 2.5166 NaN NaN NaN 5.0000 4.1130

카테고리

Help CenterFile Exchange에서 Descriptive Statistics에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by