[ncl-talk] stat_dispersion error

[Dave Allured - NOAA Affiliate]

Please look at the function documentation carefully.  It says:

; STATX(22 to 27) require more than 1000 observations
; Otherwise, they are set to _FillValue

Output #19 says you have exactly 1000 values.  Therefore, _FillValue is the
expected result for #22 through #27.

--Dave


On Fri, Jul 15, 2016 at 4:05 PM, Lijun Diao <ljdiao at gmail.com> wrote:

> Dear NCL forum:
>
> I am trying to use *stat_dispersion* function to compute 99th percentile
> of data. At first, I tried the example 1 on the introduction page. However,
> as you can see in the following, the outputs 22 to 26 are missing values
> rather that the values shown on the web. My NCL is version 6.3.0. Is it due
> to my installation issue or anyone can compute the 99th percentile without
> problems?
>
> thanks for answering my question!
>
> ncl 4> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> ncl 5> statb = stat_dispersion(b, opt )
> (0)
> (0)        ===> Robust Dispersion Statistics <===
> (0)      [0]            Mean=1.96763
> (0)      [1]          StdDev=1.94921
> (0)      [2]             Min=0.00386701
> (0)      [3]          LowDec=0.217212
> (0)      [4]          LowOct=0.284452
> (0)      [5]          LowSex=0.383918
> (0)      [6]     LowQuartile=0.592704
> (0)      [7]          LowTri=0.783372
> (0)      [8]          Median=1.34768
> (0)      [9]         HighTri=2.21197
> (0)      [10]   HighQuartile=2.75091
> (0)      [11]        HighSex=3.55215
> (0)      [12]        HighOct=4.0649
> (0)      [13]        HighDec=4.39585
> (0)      [14]            Max=12.4941
> (0)      [15]          Range=12.4903
> (0)      [16]     Dispersion=6.40786
> (0)      [17]    RMS Anomaly=1.94824
> (0)      [18]      #   Total=1000
> (0)      [19]      #    Used=1000
> (0)      [20]      # Missing=0
> (0)      [21]      % Missing=0
> (0)      [22]     Lower 0.1%=9.96921e+36
> (0)      [23]     Lower 1.0%=9.96921e+36
> (0)      [24]     Lower 5.0%=9.96921e+36
> (0)      [25]     Upper 5.0%=9.96921e+36
> (0)      [26]     Upper 1.0%=9.96921e+36
> (0)      [27]     Upper 0.1%=9.96921e+36
> (0)      [28]       Skewness=1.85963
> (0)      [29]       Kurtosis=4.42162
> (0)
>
>
> -----------------------------------------------------------------------------
> Lijun Diao
> ldiao at uh.edu
> Research Assistant
> Department of Earth & Atmospheric Sciences
> University of Houston
> 4800 Calhoun RD., Bldg. S&R1-312
> Houston, TX 77204-5007
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160715/7164ef14/attachment.html