[ncl-talk] Problem in plotting t-test figure

[Mary Haley]

Did you look at our examples page?

http://www.ncl.ucar.edu/Applications/

There are hundreds of examples showing you how to create different types of
graphics.

In particular, see example climo_1.ncl at:

http://www.ncl.ucar.edu/Applications/climo.shtml

--Mary


On Tue, Oct 20, 2015 at 8:52 PM, Steven Kong <kongksk at gmail.com> wrote:

> hi all,
>
> Current i would like to a plot a t-test output, a significant change of
> daily mean surface temperature between 01 Janaury and 02 January. Could any
> one help me?
>
> ;;;;;;;;;;;;;;;;;;;;;daily mean temperature at 01 January;;;;;;;;;;;;;;
>
> temp1=a->T2(0,:,:)
> temp1=0
>   do i = 0,23
>     temp1 = temp1 + a->T2(i,:,:)
>   end do
>
>   mean1 = (temp1)/24
>
> ;;;;;;;;;;;;;;;;;;;;;daily mean temperature at 02 January;;;;;;;;;;;;;;
>
> temp2=b->T2(0,:,:)
> temp2=0
>   do i = 0,23
>     temp2 = temp2 + b->T2(i,:,:)
>   end do
>
>   mean2 = (temp2)/24
>
> ;-------------T-Test: Area with 95% significant
> difference----------------------------------------------
>
> r_ave1 = avg(mean1)
> r_var1 = variance(mean1)
>
> r_ave2 = avg(mean2)
> r_var2 = variance(mean1)
>
> probt = 100.*(1. -ttest(r_ave1,r_var1,24,r_ave2,r_var2,24, True, False)) ;
> all areas whose value > 95 are significant
> aveX = where(probt.lt.95.,r_ave1 at _FillValue, r_ave1)
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151026/db7a1319/attachment.html