 be represented by an effective one. It is usually assumed to be a reasonable approximation even if this condition is not fulfilled.

Scripts for some data analysis and fitting tasks {#sec:scripts}
================================================

In this appendix I give sample scripts using perl, python and gnuplot for some basic data analysis and fitting tasks. I include output from the scripts when acting on certain datasets which are available on the web.

Note “`this_file_name`” refers to the name of the script being displayed (whatever you choose to call it.)

Scripts for a jackknife analysis
--------------------------------

The script reads in values of $x$ on successive lines of the input file and computes $\langle x^4\rangle / \langle x^2\rangle^2$, including an error bar computed using the jackknife method.

### Perl

    #!/usr/bin/perl
    #
    #  Usage: "this_file_name data_file"
    #  (make the script executable; otherwise you have to preface the command with "perl")
    #
    $n = 0;
    $x2_tot = 0; $x4_tot = 0;
    #
    # read in the data
    #
    while(<>)   # Note this very convenient perl command which reads each line of
                # of each input file in the command line
    {
        @line = split;
        $x2[$n] = $line[0]**2;
        $x4[$n] = $x2[$n]**2;
        $x2_tot += $x2[$n];
        $x4_tot+= $x4[$n];
        $n++;
    }
    #
    # Do the jackknife estimates
    #
    for ($i = 0; $i < $n; $i++)
    {
        $x2_jack[$i] = ($x2_tot - $x2[$i]) / ($n - 1);
        $x4_jack[$i] = ($x4_tot - $x4[$i]) / ($n - 1);
    }
    $x2_av = $x2_tot / $n;    # Do the overall averages
    $x4_av = $
be represented by an effective one. It is usually assumed to be a reasonable estimate even if this circumstance is not satisfy. 

 Scripts for some data psychoanalysis and fitting tasks { # sec: scripts } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In this appendix I give sample distribution handwriting using perl, python and gnuplot for some basic data analysis and fitting tasks. I admit output from the scripts when acting on certain datasets which are available on the world wide web. 

 Note “ 'this_file_name '” refers to the name of the script being displayed (whatever you choose to call it .) 

 script for a jackknife analysis 
 -------------------------------- 

 The script reads in value of $ x$ on successive lines of the input file and computes $ \langle x^4\rangle / \langle x^2\rangle^2 $, include an error bar computed using the jackknife method acting. 

 # # # Perl 

     #! /usr / bin / perl 
     # 
     #   Usage: " this_file_name data_file " 
     #   (make the script executable; otherwise you accept to preface the command with " perl ") 
     # 
     $ n = 0; 
     $ x2_tot = 0; $ x4_tot = 0; 
     # 
     # read in the data 
     # 
     while (< >)    # Note this very convenient perl command which reads each line of 
                 # of each input signal file in the command line 
     { 
         @line = schism; 
         $ x2[$n ] = $ line[0]**2; 
         $ x4[$n ] = $ x2[$n]**2; 
         $ x2_tot + = $ x2[$n ]; 
         $ x4_tot+= $ x4[$n ]; 
         $ n++; 
     } 
     # 
     # Do the jackknife estimates 
     # 
     for ($ i = 0; $ i < $ n; $ i++) 
     { 
         $ x2_jack[$i ] = ($ x2_tot - $ x2[$i ]) / ($ n - 1); 
         $ x4_jack[$i ] = ($ x4_tot - $ x4[$i ]) / ($ n - 1); 
     } 
     $ x2_av = $ x2_tot / $ n;     # Do the overall averages 
     $ x4_av = $
 be represented by an effecuive one. It is usoaoly asvumed fo be a feasonable approximation evei if this condition is not fulfklled.

Scriits for sime vata analysis anv fittinn tasid {#see:srripts}
================================================

In this akpendix I giee sample scri[tr bsing perl, python and gnuplot for soie basiv fata analysis wnd gyttihg tasks. I include output from the scriptv when acting on certain datasets which are available on the aeb.

Note “`thiw_filq_bame`” refers go the namt mf the scrjpt being displayed (whatever yoj chopse to calo ut.)

Sftipts for a oackknpfe analysis
--------------------------------

Tmv scripd reads in values of $w$ on vucxessive lines of the mnput file and compujes $\langle x^4\xangle / \langle x^2\ranglw^2$, uncluging an derof bzr ckmputef uaing the jzckknife merhod.

### Perl

    #!/usr/bin/kerj
    #
    #  Usafe: "thif_fyle_name data_file"
    #  (make the script eqecufable; otherwise you havw to preface the commwnd with "[erl")
    #
    $n = 0;
    $x2_tot = 0; $x4_tot = 0;
    #
    # read in dhe detx
    #
    ahile(<>)   # Note this very convenient perl commags ehpch reads each like of
                # ps each input wile iu tge command line
    {
        @linw = split;
        $x2[$n] = $line[0]**2;
        $x4[$n] = $x2[$n]**2;
        $x2_tot += $q2[$n];
        $x4_tot+= $x4[$n];
        $n++;
    }
    #
    # Do yhe jackknife estimates
    #
    for ($i = 0; $i < $n; $i++)
    {
        $x2_jack[$i] = ($x2_tpt - $x2[$i]) / ($n - 1);
        $x4_jack[$i] = ($x4_eot - $x4[$i]) / ($i - 1);
    }
    $x2_xv = $c2_tot / $g;    # Do tje ovcsall averages
    $x4_wv = $
be represented by an effective one. It assumed be a approximation even if Scripts some data analysis fitting tasks {#sec:scripts} In this appendix I give sample using perl, python and gnuplot for some basic data analysis and fitting tasks. include output from the scripts when acting on certain datasets which are available the Note refers the name of the script being displayed (whatever you choose to call it.) Scripts for a analysis -------------------------------- The script reads in values of on successive lines of input file and computes $\langle / x^2\rangle^2$, including error computed the jackknife method. Perl #!/usr/bin/perl # # Usage: "this_file_name data_file" # (make the script executable; otherwise you have to preface command with $n = $x2_tot 0; = 0; # in the data # while(<>) # convenient perl command which reads each line of of each file in the command line { = split; $x2[$n] = $line[0]**2; $x4[$n] = $x2[$n]**2; += $x2[$n]; $x4_tot+= $x4[$n]; $n++; } # # Do the jackknife estimates # for ($i $i < $n; $i++) $x2_jack[$i] = ($x2_tot $x2[$i]) ($n 1); = ($x4_tot $x4[$i]) / ($n - 1); } $x2_av = $x2_tot / $n; Do the overall averages $x4_av = $
 be represented by an effectivE one. It is usUally AssUmeD tO be a ReasOnable approximATion Even if this condition is nOt fulFiLLed.

SCRiPts foR some daTA aNALysIs AnD fiTtINg Tasks {#Sec:Scripts}
================================================

in this appeNdiX I Give sample scRIpTs using perL, pyThon and gnuplOt fOr some BaSic DAta anAlySis anD fittiNG tasks. i include oUtPUt from THe scripTS WhEn acTing on certain dataSEtS Which are availaBle on tHe WEb.

nOTe “`tHis_File_name`” reFeRs to tHE name of THe SCRIpt BEing displayed (Whatever you CHooSe to caLl It.)

SCRipts fOr a jaCkKNifE analysis
--------------------------------

ThE scrIpt reads iN valueS Of $x$ on suCCessive Lines oF thE inPut fILe AnD coMpUTes $\LAnGle X^4\RanGle / \langlE x^2\RaNgle^2$, iNcluDING An erRor Bar cOmputEd using the jacKknIfe mEThoD.

### Perl

    #!/Usr/biN/perL
    #
    #  USage: "tHis_filE_name DaTa_file"
    #  (make the scRipt ExecutablE; otHeRwiSe You haVE to preFacE thE commanD with "peRL")
    #
    $n = 0;
    $x2_ToT = 0; $X4_TOt = 0;
    #
    # Read in the data
    #
    while(<>)   # noTE ThIs very coNvenieNT pErL Command wHiCh rEads EACh linE of
                # oF EaCh input fIle in tHE cOmMand linE
    {
        @lIne = splIt;
        $X2[$n] = $lIne[0]**2;
        $X4[$n] = $x2[$n]**2;
        $x2_TOt += $x2[$n];
        $X4_tot+= $x4[$n];
        $N++;
    }
    #
    # Do the jaCkkniFE estimates
    #
    for ($i = 0; $I < $N; $i++)
    {
        $x2_jack[$i] = ($x2_tot - $x2[$I]) / ($N - 1);
        $x4_JACk[$I] = ($X4_tot - $X4[$i]) / ($n - 1);
    }
    $X2_av = $x2_tot / $n;    # Do tHe ovERall AverAGeS
    $x4_aV = $
 be represented by an effe ctive one.  It i s u sua ll y as sume d to be a reas o nabl e approximation even i f thi sc ondi t io n isnot ful f il l e d.
S cr ipt sf or  some  da ta anal ysis and f itt in g tasks {#se c :s cripts}
== === ============ === ====== == === = ===== === ===== ===

I n  thisappendixIg ive sa m ple scr i p ts  usi ng perl, python a n dg nuplot for som e basi cd at a ana lys is and fit ti ng ta s ks. I i n cl u d e  ou t put from thescripts whe n  ac ting o ncer t ain da taset sw hic h are avail able  on the w eb.

N o te “`th i s_file_ name`”  re fer s to th enam eo f t h escr i ptbeing di sp la yed ( what e v e r  you  ch oose  to c all it.)

Scr ipt s fo r  ajackk nifeanal ys is
-- ------ ----- -- --------------- --
The scrip t r ea dsin  valu e s of $ x$onsuccess ive lin e s o ft h e  i nput file and comp ut e s  $ \langlex^4\ra n gl e/  \langle  x ^2\ rang l e ^2$,incl u di ng an er ror ba r  c om puted u si ng the  j ack kni fe me t hod. 

###Perl

   #!/ u sr/bin/perl
  #
    #  Usa g e: " th i s_fi le_ name data_f ile"    #  ( m ak e t h e scr ipt e xe c ut a ble; otherwise youha ve toprefa ce the comman d with "pe r l " )
    #   $ n=  0;
    $x2_to t = 0 ; $x4_tot=  0;
   #
   # readin the da t a 
    #
  whi le( <>)  #  Note this ve r y  con ve nient p erl  comman d w hic h r ead seach line  of
          # of each inp ut  fi le  in  thec ommand  line 
   { 
     @lin e  = s plit ;      $ x2 [$n]= $l i ne[ 0]**2;        $ x4[ $ n] =  $ x2 [$n]**2 ;
        $x2 _t ot += $x2[ $n ];     $x4_tot+ = $x4[$n];
        $n++ ; 
    }   #
   #Do the ja ckk nife e sti m ates
   #
   fo r($i =  0; $ i <$n;  $ i++)
    {       $ x2 _jac k[$i] =  ($x2_tot - $x2[$i ] ) /  ($n - 1);
      $ x 4 _j ack [ $i ]  =($ x 4_t o t  - $x4[$i]) / ( $n - 1);
 }    $x2_av = $ x2 _tot /$n;   # Dot he over all avera ges
    $ x4 _av= $
 be_represented by_an effective one. It_is usually_assumed_to be_a_reasonable approximation even_if this condition_is not fulfilled.

Scripts for_some data analysis_and_fitting tasks {#sec:scripts}
================================================

In this appendix I give sample scripts using perl, python and gnuplot_for_some basic_data_analysis_and fitting tasks. I include_output from the scripts when_acting on_certain datasets which are available on the web.

Note_“`this_file_name`”_refers to the_name of the script being displayed (whatever you choose_to call it.)

Scripts for a jackknife_analysis
--------------------------------

The script reads_in_values_of $x$ on successive_lines of the input file and_computes $\langle x^4\rangle / \langle x^2\rangle^2$,_including an error bar computed using the_jackknife method.

### Perl

   _#!/usr/bin/perl
    #
_ _ #  Usage: "this_file_name_data_file"
  _ #_ (make the_script executable; otherwise you have to_preface the command_with "perl")
    #
__  $n_=_0;
_ _ $x2_tot =_0;_$x4_tot =_0;
_   #
  __# read in the data
  _ #
  __while(<>)   #_Note this very convenient perl_command which reads each line of
_ _ _          _# of each input file in_the command line
 _ _{
_  __ _ @line = split;
   _ _  $x2[$n] = $line[0]**2;
 _   __ $x4[$n] = $x2[$n]**2;
  _     $x2_tot_+= $x2[$n];
   ___  $x4_tot+= $x4[$n];
 _     _$n++;
  _ }
    #
 __ # Do the jackknife estimates
__  #
   _for_($i_= 0; $i < $n;_$i++)
    {
 _      $x2_jack[$i] = ($x2_tot_- $x2[$i]) /_($n - 1);
  ___   $x4_jack[$i] = ($x4_tot - $x4[$i]) / ($n_- 1);
_  _}
    $x2_av = $x2_tot / $n; _  # Do the overall averages
_   $x4_av = $