Login or Register to Ask a Question and Join Our Community

Sponsored Content
Full Discussion: Standard deviation in awk
Top Forums Shell Programming and Scripting Standard deviation in awk Post 302549171 by gd9629 on Monday 22nd of August 2011 11:19:10 AM
Old 08-22-2011
Hi Roboticus,

I've tried combining the codes you've given me into

Code:
function std_dev(data, count) {
    sum=0;
    for( x=1; x <= count; x++) {
        sum += data[x];
    }

    avg = sum/count;

    sumsq=0;
    for( x=1; x <= count; x++) {
        sumsq += (data[x] - avg)^2;
    }

    print sqrt(sumsq/count);
}


BEGIN {
    cnt = 0;
}


END {
    std_dev(array17, cnt); std_dev(array19, cnt); std_dev(array20, cnt);
}

{
   array17[cnt]=$17; array19[cnt]=$19; array20[cnt]=$20; cnt++;
   if (cnt==720) {
      print std_dev(array17, cnt) " " std_dev(array19, cnt)
       " " std_dev(array20, cnt);
   }
	
}

Unfortunately this just gives me a few numbers, the meaning of which I'm unsure lol

Code:
71.1077
0.178527
 
0.0675481
38.6887
0.825408
0.140382

How should I edit the script further to give the 3 separate columns of SD values?

Thanks!

---------- Post updated at 04:19 PM ---------- Previous update was at 09:58 AM ----------

Hmm I've tried this now

Code:
BEGIN {
    cnt = 0;
}


END {
    std_dev(array17, cnt); std_dev(array19, cnt); std_dev(array20, cnt);
}

{
   array17[cnt++] = $17; array19[cnt++] = $19; array20[cnt++] = $20;
   
   if (cnt==720) {
    print  std_dev(array17, cnt) " " std_dev(array19, cnt)
       " " std_dev(array20, cnt); 
	   cnt = 0; 
	  
   }
	
}

and get this as my results

Code:
272.963
1.14211
 
0.477465
251.055
1.05022
 
0.446297
232.296
1.00244
 
0.43713
222.966
0.970351
 
0.412753
219.27
0.957796
 
0.404489
216.511
0.948255

Any ideas?
Last edited by gd9629; 08-22-2011 at 11:02 AM..
 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Script for finding standard deviation

I have a CSV file that looks like 0,0,0,0,1,0,1,0,1,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0 10,11,7,0,4,12,2,3,7,0,11,3,12,4,0,5,5,4,5,0,8,6,12,0,9,3,3,0,2,7,8 19,11,7,0,4,14,16,10,8,2,13,7,15,6,0,76,6,4,10,0,18,10,17,1,11,3,3,0,9,9,8... (7 Replies)
Discussion started by: RJ17
7 Replies

2. Shell Programming and Scripting

Mean and Standard deviation

Hi all, I am new to shell scripting and wanna calculate the mean and standard deviation using shell programming. I have a file with letters that are repeating and their corresponding duration a 0.32 a 0.89 aa 0.34 aa 0.23 au 0.012 au 0.26... (4 Replies)
Discussion started by: lakshmikanth.pg
4 Replies

3. UNIX for Dummies Questions & Answers

Calculating the Standard Deviation for a column

Hi all, I want to calculate the standard deviation for a column (happens to be column 3). Does any know of simple awk script to do this? Thanks (1 Reply)
Discussion started by: kylle345
1 Replies

4. Shell Programming and Scripting

using awk to print average and standard deviation into a file

Hi I want to use awk to print avg and st deviation but it does not go into a file for column 1 only. I can do average and # of records but i cannot get st deviation. awk '{sum+=$1} END { print "Average = ",sum/NR}' thanks (1 Reply)
Discussion started by: phil_heath
1 Replies

5. Shell Programming and Scripting

AWK script for standard deviation / root mean square deviation

I have a file with say 50 columns, each containing a whole lot of data. Each column contains data from a separate simulation, but each simulation is related to the data in the last (REFERENCE) column $50 I need to calculate the RMS deviation for each data line, i.e. column 1 relative to... (12 Replies)
Discussion started by: chrisjorg
12 Replies

6. Shell Programming and Scripting

Finding standard deviation for all columns in a data file

Hi All, I want someone to modify the below script from this forum so that it can be used for all columns in the file( instead of only printing column 3 mean and standard deviation values). I don't know how to loop around all the columns. ... (3 Replies)
Discussion started by: ks_reddy
3 Replies

7. Shell Programming and Scripting

calculating row-wise standard deviation using awk

Hi, I have a file containing 100,000 rows-by-120 columns and I need to compute for the standard deviation for each row. Any idea on how to calculate row-wise standard deviation using awk? My sample data looks like this: input data: 23 35 12 25 16 17 18 19 29 12 12 26 15 14 15 23 12 12... (2 Replies)
Discussion started by: ida1215
2 Replies

8. Shell Programming and Scripting

Computing average and standard deviation from multiple text files

Hello there, I found an elegant solution to computing average values from multiple text files awk '{for (i=1;i<=NF;i++){if ($i!~"n/a"){a+=$i}else{b++}}}END{for (i=1;i<=FNR;i++){for (j=1;j<=NF;j++){printf (a/(3-b))((b>0)?"~"b" ":" ")};printf "\n"}}' file1 file2 file3 I tried to modify... (2 Replies)
Discussion started by: charmmilein
2 Replies

9. Shell Programming and Scripting

Output mean and standard deviation of a row

I have a file that looks that this: 820 890 530 1650 1600 1800 1850 1900 2270 1640 2300 1670 2080 2200 2350 1150 1630 2210 I would like to output the mean and standard deviation of each row so that my final output would look like this 820 890 530 746.667 155.849 1650 1600 1800... (5 Replies)
Discussion started by: kayak
5 Replies

10. Shell Programming and Scripting

SMA (Single Moving Average) and Standard Deviation

Hello Team, I am using the following awk script to calculate the SMA (Single Moving Average) for an specific period but now I would like to include the standard deviation output. Could you please help me to modify this awk shell script awk -F, -v points=5 ' { a = $2; ... (4 Replies)
Discussion started by: csierra
4 Replies

LEARN ABOUT DEBIAN

net::gpsd3::return::gst

Net::GPSD3::Return::GST(3pm)				User Contributed Perl Documentation			      Net::GPSD3::Return::GST(3pm)

NAME

       Net::GPSD3::Return::GST - Net::GPSD3 Return GST Object

SYNOPSIS

DESCRIPTION

       Provides a Perl object interface to the GST (GPS Pseudorange Noise Statistics) object returned by the GPSD daemon.

	 {
	   'class' => 'GST',
	   'device' => '/dev/cuaU0',
	   'tag' => '0x0130',
	   'time' => '1970-01-01T00:00:00.00Z',
	   'lat' => '0',
	   'lon' => '0',
	   'alt' => '0.002',
	   'rms' => '0',
	   'orient' => '0',
	   'major' => '0'
	   'minor' => '-0',
	 }

METHODS PROPERTIES

   class
       Returns the object class

   string
       Returns the JSON string

   parent
       Returns the parent Net::GPSD3 object

   device
   tag
   time
       Returns a unix epoch time

   timestamp
       Returns a W3C formated date

   datetime
       Returns a DateTime object

   rms
       Total RMS standard deviation of ranges inputs to the navigation solution

   major
       Standard deviation (meters) of semi-major axis of error ellipse

   minor
       Standard deviation (meters) of semi-minor axis of error ellipse

   orient
       Orientation of semi-major axis of error ellipse (true north degrees)

   lat
       Standard deviation (meters) of latitude error

   lon
       Standard deviation (meters) of longitude error

   alt
       Standard deviation (meters) of altitude error

BUGS

       Log on RT and Send to gpsd-dev email list

SUPPORT

       DavisNetworks.com supports all Perl applications including this package.

       Try gpsd-dev email list

AUTHOR

	 Michael R. Davis
	 CPAN ID: MRDVT
	 STOP, LLC
	 domain=>michaelrdavis,tld=>com,account=>perl
	 http://www.stopllc.com/

COPYRIGHT

       This program is free software licensed under the...

	 The BSD License

       The full text of the license can be found in the LICENSE file included with this module.

SEE ALSO

       Net::GPSD3, Net::GPSD3::Return::Unknown::Timestamp

perl v5.14.2							    2011-11-09					      Net::GPSD3::Return::GST(3pm)
All times are GMT -4. The time now is 11:30 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy