comp.lang.idl-pvwave archive: archive » Re: Incorrect behavior of /NAN

Home » Public Forums » archive » Re: Incorrect behavior of /NAN

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: Incorrect behavior of /NAN [message #30025 is a reply to message #30024]

Thu, 04 April 2002 10:19

Paul van Delst
Messages: 364
Registered: March 1997

Senior Member

Kenneth Bowman wrote:
>
> I tried to incite some interest in this last week, but had no takers,
> so I'll try again. ;-)
>
> I believe that the behavior of the TOTAL function is incorrect in
> the case where all of the data are NaN and the /NAN keyword is set.
>
> Here are several possible cases:
>
> IDL> a = !values.f_nan
> IDL> print, total([a, 1.0])
> NaN
> IDL> print, total([a, 1.0], /nan)
> 1.00000
> IDL> print, total([a, a])
> NaN
> IDL> print, total([a, a], /nan)
> 0.00000
>
> I believe that the last case is incorrect.
>
> The documentation for /NAN says "Elements with the value NaN are treated
> as missing data." In the last case there are no valid data, so how can
> their sum be zero?
>
> I think this is probably a simple error in the TOTAL algorithm.
> Rather than computing the "sum of all the non-NaN values", it is
> probably doing "set sum to zero and then add all non-Nan values".

Seems reasonable to me. It only sums the data not flagged as missing. If all the data is
missing, the sum of nothing is, well, nothing.

> For comparison, the MEAN function does not behave this way.
>
> IDL> print, mean([a, 1.0])
> NaN
> IDL> print, mean([a, 1.0], /nan)
> 1.00000
> IDL> print, mean([a, a])
> NaN
> IDL> print, mean([a, a], /nan)
> NaN

I think the last result *here* is wrong :o) Shouldn't it be 1.0? What is MEAN in this case
but:

total([a,a],/NAN) / n_elements(non-nan values only)

Both the numerator and denominator are 0.

> Before I submit something (at least a question) to RSI, does anyone
> have comments?

Yes. Check your data *before* you TOTAL or MEAN them. :o)

paulv

--
Paul van Delst Religious and cultural
CIMSS @ NOAA/NCEP purity is a fundamentalist
Ph: (301)763-8000 x7274 fantasy
Fax:(301)763-8545 V.S.Naipaul

Report message to a moderator

[Message index]

		Re: Incorrect behavior of /NAN By: Paul van Delst on Thu, 04 April 2002 11:38
		Re: Incorrect behavior of /NAN By: Jaco van Gorkom on Thu, 04 April 2002 10:38
		Re: Incorrect behavior of /NAN By: Paul van Delst on Thu, 04 April 2002 10:19
		Re: Incorrect behavior of /NAN By: Craig Markwardt on Thu, 04 April 2002 09:32
		Re: Incorrect behavior of /NAN By: K. Bowman on Fri, 05 April 2002 10:20
		Re: Incorrect behavior of /NAN By: Mark Hadfield on Thu, 04 April 2002 13:06
		Re: Incorrect behavior of /NAN By: Craig Markwardt on Thu, 04 April 2002 18:26

Previous Topic:	Need help identifying left most zero bit
Next Topic:	Interfacing the functions

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sun Nov 30 07:42:00 PST 2025

Total time taken to generate the page: 1.11936 seconds