This page was created by the IDL library routine mk_html_help2.

Last modified: Mon Jun 10 10:37:11 2013.


Directory Listing of Routines


Routine Descriptions

DIFFERENCE_TEST

[Next Routine] [List of Routines]
Name:
  difference_test
Purpose:
  平均値検定を行うプログラム。
  χ二乗検定によって正規分布との適合を検定する。
  どちらも正規分布に従う場合はWelch検定
  どちらかが正規分布に従わない場合はマンホイットニー検定のみを用いる。
Syntax:
  difference_test,x,y,result,sl=**
Keywords:
  result:検定結果を'0':ウェルチ検定使用-判定は異,'1':ウェルチ検定使用-判定は同,'2':マンホイットニー-異,'3':マンホイットニー-同、で返す
  sl:有意水準。指定しない場合はsl=0.05で検定。
  test_sel:行う検定を指定。'2'はマンホイットニー検定、'1'はウェルチ検定、'0'はχ二乗検定によって正規分布との適合を検定
  2012/12/13 改変 by 濵口

(See iugonet/tools/statistical_package/difference_test.pro)


MANN_WHITNEY_TEST

[Previous Routine] [Next Routine] [List of Routines]
NAME:
  Mann_Whitney_test
  
PURPOSE:
   2データ間の分布の差を検定する。 母集団の分布は問わない 
   
KEYWORDS:
  sl:significance level 
     If you don't set sl, sl = 0.05 
  mv:missing value
     If you don't set mv, mv = only NaN 

結果は分布が同じと判定されれば「3」、異なっていれば「2」を返す。
z:検定統計量
Z0:有意水準における標準正規分布の臨界値
c:コメント

CODE:
R. Hamaguchi, 13/02/2012.

MODIFICATIONS:
A. Shinbori, 01/05/2013.

ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/mann_whitney_test.pro)


NORMALITY_TEST

[Previous Routine] [Next Routine] [List of Routines]
name
 normality_test
purpose
 カイ二乗適合度検定を用いて、正規分布に従っているか検定する
syntax
 result=normality_test(x,sl=sl,mv=mv)
Keywords:
  x:検定を行う観測データ
  sl:有意水準
     入力しなければ5%で検定
  mv:欠損値
     入力しなければNaNのみを欠損値として処理;
結果は、正規分布に従っている場合は0、従っていない場合は1を返す

CODE:
R. Hamaguchi, 13/02/2012.

MODIFICATIONS:
A. Shinbori, 01/05/2013.

ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/normality_test.pro)


S_TRANS

[Previous Routine] [Next Routine] [List of Routines]

NAME:
s_trans

PURPOSE:
 Calculate the local power spectrum for given time-series data and return the structure
 which contains {st:, ph: freq:, time:} where st is the S Transform.

SYNTAX:
 results = s_trans(time-series)

KEYWOARDS:
 \help               explains all the keywords and parameters
 \verbose            flags errors and size
 \samplingrate       if set returns array of frequency
 \maxfreq
 \minfreq
 \freqsmplingrate
 \piecewisenumber    divides the time series, and passes back array
 \power              returns the power spectrum
 \abs                returns the absolute value spectrum
 \Rremoveedge        removes the edge with a 5% taper, and takes out least-sqares fit parabola

 Added a "mask out edges" keyword, which STs a line (ie st of edges) and thresholds the returned st matrix.
 The value of masked edges is the percent at which to  make the threshold default = 5 %.

 Added an EXAMPLE keyword, will display a time series and the amplitude of the ST to the current graphics device
 WARNING, will call !P.multi=0

===================== Modified hilbert transform ================================
;;;;;;;;;;;;;;;;;;;;
+
 NAME:
  hilbert

 PURPOSE:
  Return a series that has all periodic terms shifted by 90 degrees.

CODE:
 A. Shinbori, 04/06/2012.

MODIFICATIONS:
 
 
ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/s_trans.pro)


TREND_TEST

[Previous Routine] [Next Routine] [List of Routines]

NAME:
  trend_test
PURPOSE:
  データの傾向性を検定。
SYNTAX:
  trend_test,y
KEYWORDS:
    y:データ
    Z:統計検定量
    sl:有意水準指定用 (Default is 0.05)

CODE:
 R. Hamaguchi, 17/01/2013.

MODIFICATIONS:
 A. Shinbori, 01/05/2013.
 
ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/trend_test.pro)


UCHANGE_POINT_CHECKER

[Previous Routine] [Next Routine] [List of Routines]

NAME:
uchange_point_checker
  
PURPOSE:
Identify the change point of time-serise data.
  
SYNTAX:
uchange_point_checker,vname1
  

CODE:
R. Hamaguchi, 17/01/2013.

MODIFICATIONS:
A. Shinbori, 08/05/2013.
 
ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/uchange_point_checker.pro)


UCROSS_COR

[Previous Routine] [Next Routine] [List of Routines]
NAME:
ucross_cor

PURPOSE:
  Calculate the cross correlation between two time-serise data 
  and perform a test of correlation. 

CALLING SEQUENCE:
 ucross_cor, vname1,vname2
 
INPUT:
 vname1 = first tplot variable name
 vname2 = second tplot variable name

OUTPUT:
 cross_cor = cross correlation between two time-serise data
 x_cor = auto correlation of the first tplot variable
 y_cor = auto correlation of the second tplot variable
 
KEYWORDS:
  sl  = Significant level of correlation test.
  low_lag = lag time of left direction
  high_lag = lag time of right direction
  
EXAMPLE:
  ucross_cor, tplot1, tplot2

CODE:
R. Hamaguchi, 15/01/2013.

MODIFICATIONS:
A. Shinbori, 31/01/2013.

ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/ucross_cor.pro)


UDATA_INTERPOLATION

[Previous Routine] [Next Routine] [List of Routines]
NAME:
udata_interpolation

PURPOSE:
 Perform the data interpolation for two time-serise data. 

CALLING SEQUENCE:
 udata_interpolation, vname1,vname2
 
INPUT:
 vname1 = first tplot variable name
 vname2 = second tplot variable name
 
KEYWORDS:
  st_time0  = start time of data interpolation
  ed_time0 = end time of data interpolation
  reverse = exchange of two data sets
    If set, the data interplation is performed on the basis of time range of data2.
  set_interval = set up time interval to interplate the two data sets.
  If not set, the time interval is determined on the basis of the minimum interval
  in the two data sets. 
  
EXAMPLE:
   udata_interpolation, tplot1, tplot2

CODE:
R. Hamaguchi, 15/01/2013.

MODIFICATIONS:
A. Shinbori, 20/02/2013.

ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/udata_interpolation.pro)


USPEC_COH

[Previous Routine] [Next Routine] [List of Routines]
NAME:
uspec_coh

PURPOSE:
  Calculate the coherence and phase difference between two time-serise data. 

CALLING SEQUENCE:
 uspec_coh, vname1,vname2
 
INPUT:
 vname1 = first tplot variable name
 vname2 = second tplot variable name

OUTPUT:
 main_period = the period given the maximum coherence.
 
KEYWORDS:
  deltat  = time interval of observation data
  width = the width of window function
    Default is 10.
  window = type of window function (hanning, boxcar, gaussian, triangle)
    Default is 'hanning' window.
  xsize = plot size of x direction. Default is 1400. 
  ysize = plot size of y direction. Default is 700.
  
EXAMPLE:
   uspec_coh, tplot1, tplot2

CODE:
R. Hamaguchi, 15/01/2013.

MODIFICATIONS:
A. Shinbori, 31/01/2013.

ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/uspec_coh.pro)


USTRANS_PWRSPC

[Previous Routine] [Next Routine] [List of Routines]
PROCEDURE: 
 ustrans_pwrspc, varname, newname = newname
 
PURPOSE:
 Calculate the local power spectrum for given time-series data and return the structure
 which contains {st:, ph: freq:, time:} where st is the S Transform and
   store the tplot variable.
                   
INPUT:
 varname = variable passed to get_data, example - uts_mag_ccnv

KEYWORDS:
 newname: set output variable name
 \help               explains all the keywords and parameters
 \verbose            flags errors and size
 \samplingrate       if set returns array of frequency
 \maxfreq            maximum frequency performing the S transform
 \minfreq            minimum frequency performing the S transform
 \freqsmplingrate    frequency interval 
 \power              returns the power spectrum
 \abs                returns the absolute value spectrum
 \Rremoveedge        removes the edge with a 5% taper, and takes out least-sqares fit parabola
 
 
HISTORY:
 14-aug-2012, Atsuki Shinbori
$LastChangedBy$
$LastChangedDate$
$LastChangedRevision$
$URL$

(See iugonet/tools/statistical_package/ustrans_pwrspc.pro)


UTREND_TEST

[Previous Routine] [Next Routine] [List of Routines]

NAME:
utrend_test
  
PURPOSE:
Test whether there is a significant trend of time-serise data or not.
  
SYNTAX:
utrend_test,vname1,sl=sl
  
KEYWORDS:
sl:significant level (Default is 0.05)

CODE:
R. Hamaguchi, 17/01/2013.

MODIFICATIONS:
A. Shinbori, 01/05/2013.
 
ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/utrend_test.pro)


WELCH_TEST

[Previous Routine] [Next Routine] [List of Routines]
NAME:
  welch_test
PURPOSE:
  2つのデータ間の平均値が同一かを検定する。仮定は両分布とも正規分布かつ異分散。
SYNTAX:
  result=Welch_test(x,y,sl=**,mv=**)
KEYWORDS:
  sl:有意水準。設定しなければ5%で検定。
  mv:欠損値。設定しなければNaNのみを処理する。
結果は同一ならば「1」を、異なれば「0」を返す。
z:検定統計量
t_sl:有意水準におけるt分布の臨界値

CODE:
 R. Hamaguchi, 13/12/2012.

MODIFICATIONS:
 A. Shinbori, 01/05/2013.
 
ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/welch_test.pro)


  UDIFFERENCE_TEST

[Previous Routine] [List of Routines]
NAME:
  udifference_test

PURPOSE:
  平均値検定を行うプログラム。
  χ二乗検定によって正規分布との適合を検定する。
  どちらも正規分布に従う場合はWelch検定
  どちらかが正規分布に従わない場合はマンホイットニー検定のみを用いる。
  
SYNTAX:
  difference_test,vname1,vname2,sl=**
  
KEYWORDS:
  result:検定結果を'0':ウェルチ検定使用-判定は異,'1':ウェルチ検定使用-判定は同,'2':マンホイットニー-異,'3':マンホイットニー-同、で返す
  sl:有意水準。指定しない場合はsl=0.05で検定。
  test_sel:行う検定を指定。'2'はマンホイットニー検定、'1'はウェルチ検定、'0'はχ二乗検定によって正規分布との適合を検定
  
CODE:
R. Hamaguchi, 13/02/2012.

MODIFICATIONS:
A. Shinbori, 01/05/2013.

ACKNOWLEDGEMENT:
 $LastChangedBy:  $
 $LastChangedDate:  $
 $LastChangedRevision:  $
 $URL $

(See iugonet/tools/statistical_package/udifference_test.pro)


href="mddb/overlay_map_obs.pro">iugonet/tools/mddb/overlay_map_obs.pro)


-serise data. SYNTAX: uchange_point_checker,vname1 CODE: R. Hamaguchi, 17/01/2013. MODIFICATIONS: A. Shinbori, 08/05/2013. ACKNOWLEDGEMENT: $LastChangedBy: $ $LastChangedDate: $ $LastChangedRevision: $ $URL $

dBy: $ $LastChangedDate: S LxHdRevision: $ $URL $

(See iugonet/tools/statistical_package/udifference_test.pro)