<div dir="ltr"><div><div><div><div><div><div><div><br>The question is difficult to answer.<br>You have not included information which ncl-talk constantly request of users.<br></div><div><br>[0] Always include the version of NCL you are using.<br></div><div><br></div><div>Always look at your data and the results of assorted calculations via printVarSummary and print<br></div><div><br></div>[1] printVarSummary(ts1) ; ?? (time) => [ntim] ?? or (time,lat,lon) ==> [ntim,nlat,mlon] ??<br></div><div> printVarSummary(ts2)<br><br></div><div>[2] ccr = escorc(ts2,ts1)<br></div><div> printVarSummary(ccr)<br></div><div><br>[3] Nx = num(.not.ismissing(ts1))<br> print(Nx) ; if ts1 is one-dimensional<br></div><div> printVarSummary(Nx) ; if ts1 & ts2 are three-dimensional<br></div><div><br></div> This give you the **total** number of non-missing points.<br></div> If the ts1 and ts2 arrays are three-dimensional, you would likely see a VERY LARGE number<br></div><div> Very large numbers result in everything being significant.<br></div><div><br></div>[4] If ts1 and ts2 are (time,lat,lon) and you have NCL versions 6.1.1 or<br> earlier, then you should reorder the arrays so time is the<br></div> rightmost dimensions . See escorc Example 5.<br><br></div> If you have NCL v6.2.1 and 3D arrays, you can use ccr = escorc_n(ts2,ts1,0)<br><br>====<br></div>Respond *only* to ncl-talk. Please, no personal salutation.<br><div><div><div><div><br></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 6, 2014 at 11:09 AM, Vanúcia Schumacher <span dir="ltr"><<a href="mailto:vanucia-schumacher@hotmail.com" target="_blank">vanucia-schumacher@hotmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr"><font color="#000000" face="Calibri,sans-serif"><span style="font-size:21px">I need to plot correlation and the significance level of 95%, which, I did r test and applied the example </span></font><h2 style="line-height:29.990400314331055px;font-size:21px;font-family:'Segoe UI Light','Segoe UI Web Light','Segoe UI Web Regular','Segoe UI','Segoe UI Symbol',HelveticaNeue-Light,'Helvetica Neue',Arial,sans-serif;font-weight:normal;margin:0px 0px 13px;padding-top:2px;color:rgb(68,68,68)">conOncon_4.ncl <span style="color:rgb(0,0,0);font-family:Calibri,sans-serif">of NCL, but I think there's something wrong with my script, because for all the data I use, plot area 100% the graph as significant, and I think there's something wrong, someone could help me on this?</span></h2><br><div style="line-height:21.299999237060547px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255)">....</div><div style="line-height:21.299999237060547px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255)"><br></div><div style="line-height:21.299999237060547px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255)"><div>; calculate cross correlations</div><div>;************************************************<span style="line-height:22.719999313354492px;font-size:12pt"> </span></div><div> ccr = escorc(ts2,ts1)</div><div> </div><div> copy_VarCoords(ts1,ccr)</div><div> </div><div>;*************************************************</div><div>; statistical significance of ccr</div><div>;*************************************************</div><div> </div><div> siglvl=0.05</div><div><br></div><div> Nx = num(.not.ismissing(ts1))</div><div><br></div><div> prob = rtest(ccr,Nx, 0) </div><div> </div><div> copy_VarCoords(ccr,prob)</div></div><div style="line-height:21.299999237060547px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255)"><br></div><div style="line-height:21.299999237060547px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255)">....</div><div style="line-height:21.299999237060547px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255)"><br></div><div style="line-height:21.299999237060547px;color:rgb(68,68,68);font-size:15px;background-color:rgb(255,255,255)"><div> res2@cnLevelSelectionMode = "ManualLevels" </div><div> res2@cnMinLevelValF = 0.00 ; set min contour level</div><div> res2@cnMaxLevelValF = 1.05 ; set max contour level</div><div> res2@cnLevelSpacingF = 0.05 ; set contour spacing</div><div><br></div><div>plot2 = gsn_csm_contour(wks,gsn_add_cyclic_point(prob(:,:)), res2) </div><div><br></div><div> plot2 = ShadeLtContour(plot2,<font style="line-height:normal" color="#ac193d"> 0.06, 17</font>) ; shade all areas less than the 0.05 contour level</div><div> overlay (plot, plot2)</div><div><br></div></div><div><font face="Comic Sans MS"><span style="line-height:17px">---</span></font></div><span style="line-height:17px;background-color:rgb(255,255,255)"><font style="font-size:12pt" face="Comic Sans MS" size="3"><font>Vanúcia Schumacher</font></font></span><div><span style="line-height:17px;background-color:rgb(255,255,255)"><font face="Comic Sans MS">Mestranda em Meteorologia - UFV<font style="font-size:12pt" size="3"><font><br style="line-height:17px">Meteorologista -UFPel</font><br></font></font></span></div><div><span style="line-height:17px;background-color:rgb(255,255,255)"><font face="Comic Sans MS"><font style="font-size:12pt" size="3"><font>Departamento de Meteorologia Agrícola - DEA</font></font></font></span></div><div><span style="line-height:17px;background-color:rgb(255,255,255)"><font style="font-size:12pt" face="Comic Sans MS" size="3"><font>Cel: (31) 9978 2522 </font></font></span></div><div><span style="line-height:17px;background-color:rgb(255,255,255)"><font style="font-size:12pt" face="Comic Sans MS" size="3"><font>DEA: (31) 3899 1890</font></font></span></div> </div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>