[ncl-talk] plot lead/lag correlation with statistical significance

Dennis Shea shea at ucar.edu
Tue Nov 20 15:04:30 MST 2018


My opinion:

[1] You should remove the climatological annual cycle from each dataset.

*http://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml>
*http://www.ncl.ucar.edu/Document/Functions/Contributed/calcMonAnomTLL.shtml*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/calcMonAnomTLL.shtml>


[2] Typically, the 'swe' file contains a _FillValue (-999.0) when there is
no snow cover. I suggest setting this  to 0.0.

   sn= in1->swe(iStrt:iLast,{50:70},{20:140})  ;only specific region
averaged over total Eurasia
   sneu r=* wgt_areaave_Wrap*(sn,1.0,1.0,0);;;create area averaged snow
over Eurasia
   sneur = *where*(*ismissing*(sneur), 0.0, sneur)
   *printVarSummary*(sneur)
   p*rintMinMax(*sneur,0)
   print("---")
   ymd    := *cd_calendar*(sneur&time, -2)  ; yyyymmdd
   print(ymd+"  "+sneur)
   print("---")

[3] You *know* the climatology. You want to know the feedback of sst/sneur
anomalies on each other.
Use the anomalies in the correlation calculations.

[4] There are statistical issues. EG: the number of independent values to
be used in the testing.

Typically, successive monthly values of SST or SNE are not independent. If
(say) a January SNEUR is very large, then the February SNE will likely be
large also. So estimating the number of independent values is an issue.

 *
http://www.ncl.ucar.edu/Document/Functions/Built-in/equiv_sample_size.shtml
<http://www.ncl.ucar.edu/Document/Functions/Built-in/equiv_sample_size.shtml>*

Estimate for the SNEUR/SST anomalies separately. Use the smaller value.

[5] I am sure the IITM has staff members who know more about statistical
testing/inference than ncl-talk. I suggest talking with them.

Good luck

On Tue, Nov 20, 2018 at 6:13 AM Sujata Mandke <amin at tropmet.res.in> wrote:

> Dear NCL community,
> Greetings!
>  I have ploted lag/lead correlation (cc) between two area-averaged
> (Eurasian snow and East Pacific SST) monthly time series as xy plot
> (line).
> Further, I want to indicate which part of cc is statistically
> significant on this line using line markers.
>
> I had tested statistical significance of cc using “rtest”.
>  In the ncl script(attached),  I have plotted  line for
>  lead/lag correlation. How to mark part of this cc (line),
> which is statistically significant by using line markers.
>
> My questions are:
> (i) Is my statistical significance testing correct?
>
> (ii) how to indicate those points (by marker or any other way),
>  that are statistically significant on the   lead/lag cc line plot.
>
> I had extensively searched NCL-talk archives but
> did not find solution to my problem. My guess is that,
>  i had to use “where” function and then gsn_add_polymareker,
>  but do not know how to implement this in the ncl script.
> I am using NCL version 6.4.0 on linux machine.
>
> Any suggestion would be of great help.
> Many thanks in advance.
> With best regards
> Dr. Sujata Mandke
> Scientist, Indian Institute of Tropical Meteorology
> PUNE, INDIA_______________________________________________
> 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/20181120/f7dc0ed3/attachment.html>


More information about the ncl-talk mailing list