Computation of mean and standard deviation after supressing NaNs of an array
1 回表示 (過去 30 日間)
古いコメントを表示
nnmax = 4;
for nn=2:nnmax
D_filt(nn,:) = D_filt(nn-1,:);
[Dmax(nn,aa),Imax(nn,aa)] = max(D_filt(nn-1,:));
D_filt(nn,Imax(nn,aa)) = NaN ;
end
#Comments (loop indicates)
1. Scans from 0 to 360 degrees with 10 degree interval.
2. aa represents the same [0:10:360]
3. At each value of angle, the detector collects 10 data points
4. OBJECTIVE: To discard 4 data points out of the 10 and compute mean and SD
5. nn represnts the iteration number while supressing the bad data points
6. The output of the above loop displays
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 0
1 2 3 4 5 6 7 8 0 0
1 2 3 4 5 6 7 0 0 0
Instead of 4 rows, a single row matrix after supressing four bad points is desired, so that mean and SD can be computed with ease.
0 件のコメント
回答 (2 件)
Thomas Koelen
2015 年 4 月 8 日
編集済み: Thomas Koelen
2015 年 4 月 8 日
Hi Naga,
The term "average" usually encompasses several ways to measure what value best represents a sample. There are various measurements that are used and at times the term and measurement used depends on the situation.
A statistician or mathematician would use the terms mean and average to refer to the sum of all values divided by the total number of values, what you have called the average. This especially true if you have a list of numbers. In fact even in mathematics there are different "averages" or "means" and this one is more properly called the arithmetic mean.
I think the thing you are trying to accomplish is straight up using the function:
mean()
For the standard deviation use this function:
std()
Thomas Koelen
3 件のコメント
Thomas Koelen
2015 年 4 月 8 日
編集済み: Thomas Koelen
2015 年 4 月 8 日
a(isnan(a(:,1)),:)=[]
this should give you your vector without NaN, then you can just apply the mean() and the std() function.
have a look at this, a is your vector here.
Greg Dionne
2015 年 4 月 8 日
編集済み: Greg Dionne
2015 年 4 月 8 日
You can use mean(x,'omitnan') in R2015a. If you have the Statistics Toolbox, you can use nanmean.
Otherwise, try this:
n = sum(~isnan(x));
x(isnan(x)) = 0;
s = sum(x);
avg = s./n;
7 件のコメント
Greg Dionne
2015 年 4 月 9 日
In that case you could do:
sum(x(4,:),2) ./ sum(0~=x(4,:),2)
But it seems you are still confused as to why you are not getting NaN in your output like you are expecting. Perhaps if you posted your data (save it to a MAT file and post it to your question) and a quick blurb about what you want to see happen (i.e. what you expect the answer to be) we'll be able to point you in the right direction.
Good luck!
参考
カテゴリ
Help Center および File Exchange で Calendar についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!