[ncl-talk] correlation and the significance

[Vanúcia Schumacher]

I need to plot correlation and the significance level of 95%, which, I did r test and applied the example  conOncon_4.ncl‏ 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?
....
; calculate cross correlations;************************************************   ccr    =   escorc(ts2,ts1)                    copy_VarCoords(ts1,ccr) ;*************************************************;  statistical significance of ccr;*************************************************   siglvl=0.05
  Nx    = num(.not.ismissing(ts1))
  prob  = rtest(ccr,Nx, 0)                      copy_VarCoords(ccr,prob)
....
 res2 at cnLevelSelectionMode = "ManualLevels"  res2 at cnMinLevelValF      = 0.00        ; set min contour level res2 at cnMaxLevelValF      = 1.05        ; set max contour level  res2 at cnLevelSpacingF     = 0.05        ; set contour spacing
plot2   = gsn_csm_contour(wks,gsn_add_cyclic_point(prob(:,:)), res2) 
  plot2   = ShadeLtContour(plot2, 0.06, 17)  ; shade all areas less than the  0.05 contour level  overlay (plot, plot2)
---Vanúcia SchumacherMestranda em Meteorologia - UFVMeteorologista -UFPel
Departamento de Meteorologia Agrícola - DEACel: (31) 9978 2522 DEA: (31) 3899 1890 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141006/f42a4c22/attachment.html