regression coefficient

From: Xi Chang <xi.chang01_at_nyahnyahspammersnyahnyah>
Date: Mon Jan 07 2013 - 09:14:26 MST

I really need a help, I want to plot the regression coeff at 95% confidence
limit in 2D map.
---[snip]
   x = ispan(0,dimsizes(trmm&time)-1,1)*1.
   rc = regCoef(x,trmm)
   df = rc@nptxy-2 ; degrees of freedom
   tval = rc@tval ; t-statistic
   b = tval ; b must be same size as tval (and df)
   b = 0.5
   prob = (1 - betainc(df/(df+tval^2),df/2.0,b) )
   copy_VarCoords(rc,prob)
   rc95 = where(prob.lt.95,rc@_FillValue,rc)
--[snip]
 plot = gsn_csm_contour_map_ce(wks,rc,res)

but i got an error:
fatal:where: condition variable (parameter 1) dimension mismatch with
parameter 2

Anybody can give me some suggestions or maybe the codes are not correct (?).

Thanks,
Xi Chang.

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jan 7 09:14:37 2013

This archive was generated by hypermail 2.1.8 : Mon Jan 14 2013 - 14:31:38 MST