Re: confidence interval/level

Hi Dennis,
   Thanks for your help. I followed your instruction, but I have got the error as follw,

fatal:wherefunc: dimension sizes of parameter 0 and parameter 1 do not match
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 59 in file spectrum_ana.ncl

I think that we should use

 phase=where(spec@coher_probability(1).le.spec@coher,phase@_FillValue,phase)

not what you said. But I am still not sure we should use spec@coher or sqrt(spec@coher)? Because the spec@coher is the squared of coherence.

   Thanks.

2013-09-18



朱学明



发件人：Dennis Shea <shea@ucar.edu>
发送时间：2013-09-18 07:29
主题：Re: [ncl-talk] confidence interval/level
收件人："朱学明"<557xiaozhu@163.com>
抄送："ncl-talk@ucar.edu USERS"<ncl-talk@ucar.edu>

    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 Wed Sep 18 07:46:35 2013