[ncl-talk] Regarding trend_manken function

Ipsita Putatunda ipsita.putatunda at gmail.com
Sun Jan 3 10:47:42 MST 2021


Thanks a lot Dennis for clarifying my doubts.

But what is the difference between pc(0) and PC? And why we are doing
(1-pc(0)).
As per the function tren_manken documentation pc(0) is giving the
probability to reject the null Hypothesis. Then what is PC actually?

Thanks in advance,
Ipsita

On Sun, Jan 3, 2021, 10:17 PM Dennis Shea <shea at ucar.edu> wrote:

> From Example 1 which compares NCL's *trend_manken*
> <http://www.ncl.ucar.edu/Document/Functions/Built-in/trend_manken.shtml>
> with results an *R* function
>
>      c  = (/-1,2,-3,4,-5,6,-7,8,-9,10,-11,12/)
>      pc = *trend_manken*(c, False, 0)     ; pc(0)=*0.2683*,   pc(1)=1.0
>
> For comparison, on the last example R (after loading the 'Kendall'
> library) returns:
>
>     < *library(Kendall)*
>     > q  <- c(-1,2,-3,4,-5,6,-7,8,-9,10,-11,12)
>     > x  <- c(1,2,3,4,5,6,7,8,9,10,11,12)
>     > qx <-Kendall(x,y)
>     > qx
> R returns
>       tau = 0.0909, *2-sided *pvalue =*0.7317*    (Note: 1-0.7317= *0.2683* which matches NCL)
> ====
>
> So two-sided
> ====
>
> Run the above example in Python.  That should give you a baseline comparison.
> *******************************************************************************
> Another R example:   *https://www.statology.org/mann-kendall-trend-test-r/* <https://www.statology.org/mann-kendall-trend-test-r/>
> *******************************************************************************
>
> prc = (/31.69,29.77,31.70,33.06,31.31,32.72,31.18,29.90,29.17,31.48,28.11,32.61 \
>        ,31.31,30.96,28.40,30.68,33.67,28.65,30.62,30.21,28.79,30.92,30.92,28.13 \
>        ,30.51,27.63,34.80,32.10,33.86,32.33,25.69,30.60,32.85,30.31,27.71,30.34 \
>        ,29.14,33.41,33.51,29.90,32.69,32.34,35.01,33.05,31.15,36.36,29.83,33.70 \
>        ,29.81,32.41,35.90,37.45,30.39,31.15,35.75,31.14,30.06,32.40,28.44,36.38 \
>        ,31.73,31.27,28.51,26.01,31.27,35.57,30.85,33.35,35.82,31.78,34.25,31.43 \
>        ,35.97,33.87,28.94,34.62,31.06,38.84,32.25,35.86,32.93,32.69,34.39,33.97 \
>        ,32.15,40.16,36.32/)
>
>  pc = *trend_manken(*prc, False, 0)
>  print(pc)
>  print("============")
>  *PC = 1-pc(0)*
>  print("PC=1-pc(0)="+PC)
>
> Variable: pc
> Type: float
> Total Size: 8 bytes
>             2 values
> Number of Dimensions: 1
> Dimensions and sizes:	[2]
> Coordinates:
> (0)	0.9997079
> (1)	0.03999999
> (0)	============
> (0)	*PC=1-pc(0)=0.000292063   <=== matches R*
>
>
> On Sat, Jan 2, 2021 at 11:45 PM Ipsita Putatunda via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Dear ncl users,
>> I have gone through the function trend_manken, but I have few doubts as
>> follows-
>> The probability we are getting from this function is two-tailed or
>> one-tailed?
>> Using this function I am getting one of my trend probability 0.97 whereas
>> using python I am getting the value 0.022 which is giving a
>> significantly strong trend.
>> Any help in this regard would be highly appreciated.
>>
>> Thanking you,
>> Ipsita
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at mailman.ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210103/95cd6bc6/attachment.html>


More information about the ncl-talk mailing list