NaN Suite

Descriptive Statistics for N-D matrices ignoring NaNs.

Description

The collection of functions in the NaN Suite enhances the functionality of the original 'nan'-functions in the Statistics Toolbox for N-D matrices. In addition, the new function nansem.m calculates the standard error of the mean ignoring NaNs.

The NaN Suite functions are used similar to the original functions (e.g. mean, std, var), that is they take a matrix as the first input, and the dimension along which the operation should be computed as the second input. For example, nanmean(A,4) calculates the mean of A along the 4th dimension of A ignoring NaNs.

• MATLAB release: MATLAB 6.5 (R13)
• Other Files: nansuite/nanmax.m, ; nansuite/nanmean.m, ; nansuite/nanmedian.m, ; nansuite/nanmin.m, ; nansuite/nansem.m, ; nansuite/nanstd.m, ; nansuite/nansum.m, ; nansuite/nanvar.m

[Julie Golomb]

so useful -- thank you!

[Can Ozan Tan]

pretty handy

[Daphne Weihs]

Thanks!
This happened to be just what I needed.

[Rudolf von Steiger]

 

[dgfgf ddf]

thanks a lot

[Uri Merhav]

Great minds think alike, I see. I was just about to write out this code myself when I realized someone out there must also be using nan as a non-measurement, and want it ignored while computing average and standard deviation.

Good job.

[Santosh Verma]

nice

[anony mouse]

genius!

[Pieter Kuppens]

Isn't life as simple as this:
nanmean(a) = mean(a(find(~isnan(a))))

[David Redish]

Thank you for fixing this. We have been bitten regularly by Mathworks' strange decision to define nanxyz(all-nans) as 0, which propogates nans incorrectly. This suite has helped us tremendously.

[C D]

This is a great suite of functions. I did notice the max/min toolbox functions in M7 filtered NaN's from the calculations. However, the other functions did not. So, that's where the NaN suite came in handy.

[miguel Ampudia]

is there any function to do quintiles ignoring NaN? Thanks!

[Andrei Terebilo]

Little bug here (otherwise great, thanks):

>> nanmedian([0 0 0 0 0 1])
ans =
   NaN

Just modify line 71:
i = find(y==0);
to
i = find(s==0);