Re: confidence interval/level

    spec = specxy_anal(x,y,d,sm,pct)
    PHASE = spec@phase
    PHASE@_FillValue = 1e20
    PHASE = where(spec@coher_probability(1).le.0.95,PHASE@_FillValue,PHASE)

=========
See: http://www.ncl.ucar.edu/Applications/xy.shtml#ex25

The right (R) plot should

  resR = True ; plot mods desired
  resR@xyMarkLineModes = "Markers" ; choose which have markers
  resR@xyMarkers = 16 ; choose type of marker
  resR@xyMarkerColor = "red" ; Marker color
  resR@xyMarkerSizeF = 0.005 ; Marker size (default 0.01)

   plot = gsn_csm_xy2(wks,spec@frq,spec@coher,PHASE,resL,resR)

On 9/17/13 4:04 PM, 朱学明 wrote:
> Hi all,
> I have tried to use the function specxy_anal get the coherance and phase of two timeseries data based on the example script spec_2.ncl. But I just want to get the coherence and phase satisfy the 0.95 confidence level, the picture likes that in the attachment. How can I do that?
> It looks like that the output of function specx_ci is not what I want.
>
> Any help would be appreciated. Thanks.
>
> 2013-09-18
>
>
>
> 朱学明
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Sep 17 17:29:52 2013