Re: ttest......

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Wed, 02 Sep 2009 11:16:23 -0600

It looks to me like your data are (lat,lon,time)

where 'time' is size 4.

As Adam said, s1,s2 are scalars containing the number
of statistically independent observations.
Hence, s1=s2=4

Adam Phillips wrote:
> Hello,
> Please take a look at the ttest documentation on the ncl website:
> http://www.ncl.ucar.edu/Document/Functions/Built-in/ttest.shtml
> s1,s2 are scalars containing the number of statistically independent
> observations.
>
> There are numerous examples that show how to use the function located at
> the bottom of the page.
> Best regards,
> Adam
>
> Sabeerali(sebi) wrote:
>> Hello support,
>>
>> I want to do the ttest for composite..when I searched in the NCL
>> tutorial I found the argument needed for ttest is ave1,ave2,
>> var1,var2,s1,s2, iflag,tval_opt
>> How I can find the arguments s1 and s2. What is mean by this?
>> The dimension of the ave1,ave2,var1,var2 are nlatXnlon..
>>
>> Any help would be appreciable
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>> begin
>> ;**************************************************************
>> in2=addfile("mn_olr_1983-05.nc <http://mn_olr_1983-05.nc>","r")
>> olr=in2->olr
>> olrflt=short2flt(olr)
>> mn_clim_olr=clmMonTLL(olrflt)
>> ;printVarSummary(mn_clim_olr)
>> Mon_olr_anom=calcMonAnomTLL(olrflt,mn_clim_olr)
>> Mon_olr_anom@_FillValue=32767
>> Mon_olr_anom_at_missing_value=32767
>>
>> ;**********************************************************************
>> ; significant test
>> ;*********************************************************************
>> four_month_runavg=runave_Wrap(Mon_olr_anom(lat|:,lon|:,time|:),4,0)
>> JJAS_olr_anom_1983_05=four_month_runavg(:,:,6::12)
>>
>> JJAS_olr_anom_4yrs =
>> (/JJAS_olr_anom_1983_05(:,:,0),JJAS_olr_anom_1983_05(:,:,7),JJAS_olr_anom_1983_05(:,:,16),JJAS_olr_anom_1983_05(:,:,21)/)
>>
>> ave1=dim_avg_n(JJAS_olr_anom_4yrs,0)
>> var1=dim_variance_n(JJAS_olr_anom_4yrs,0)
>> ave2=dim_avg_Wrap(JJAS_olr_anom_1983_05)
>> var2=dim_variance_Wrap(JJAS_olr_anom_1983_05)
>>
>> end
>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> 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 02 2009 - 11:16:23 MDT

This archive was generated by hypermail 2.2.0 : Thu Sep 03 2009 - 07:55:08 MDT