Sponsored Content
Top Forums Shell Programming and Scripting AWK script for standard deviation / root mean square deviation Post 302589723 by chrisjorg on Thursday 12th of January 2012 12:03:39 PM
Old 01-12-2012
Code:
2.91187 4.25656 7.3225   ..... until column 50 
3.4187   2.67656 6.3225
3.54117 6.27656 4.3225
5.61187 6.27656 2.3225
....          ...           ....

The output should just be a column of numbers where each entry represents the calculated root-mean-square deviation (RMSD) of each entry of the 50 columns, i.e. (here I am just giving you random numbers)
Code:
RMSD
4.31185 
3.4185   
2.64115
4.71183 
....

The error on running the script as it is is
Code:
awk -f rmsd merge.pmf 
awk: syntax error at source line 6 source file rmsd
 context is
          END  >>> 
 <<< 
awk: bailing out at source line 7

Moderator's Comments:
Mod Comment Use code tags please, check your PMs.

Last edited by zaxxon; 01-12-2012 at 01:59 PM.. Reason: code tags
 

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

Standard deviation in awk

Hi all, I need to find the standard deviation of each column of a dataset below for each hour. The data is given in 5 second intervals as shown below DATE TIME FRAC_DAYS_SINCE_JAN1 FRAC_HRS_SINCE_JAN1 EPOCH_TIME ... (11 Replies)
Discussion started by: gd9629
11 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
g_rms(1)					 GROMACS suite, VERSION 4.5.4-dev-20110404-bc5695c					  g_rms(1)

NAME
g_rms - calculates rmsd's with a reference structure and rmsd matrices VERSION 4.5.4-dev-20110404-bc5695c SYNOPSIS
g_rms -s topol.tpr -f traj.xtc -f2 traj.xtc -n index.ndx -o rmsd.xvg -mir rmsdmir.xvg -a avgrp.xvg -dist rmsd-dist.xvg -m rmsd.xpm -bin rmsd.dat -bm bond.xpm -[no]h -[no]version -nice int -b time -e time -dt time -tu enum -[no]w -xvg enum -what enum -[no]pbc -fit enum -prev int -[no]split -skip int -skip2 int -max real -min real -bmax real -bmin real -[no]mw -nlevels int -ng int DESCRIPTION
g_rms compares two structures by computing the root mean square deviation (RMSD), the size-independent rho similarity parameter ( rho) or the scaled rho ( rhosc), see Maiorov & Crippen, Proteins 22, 273 (1995). This is selected by -what. Each structure from a trajectory ( -f) is compared to a reference structure. The reference structure is taken from the structure file ( -s). With option -mir also a comparison with the mirror image of the reference structure is calculated. This is useful as a reference for 'significant' values, see Maiorov & Crippen, Proteins 22, 273 (1995). Option -prev produces the comparison with a previous frame the specified number of frames ago. Option -m produces a matrix in .xpm format of comparison values of each structure in the trajectory with respect to each other structure. This file can be visualized with for instance xv and can be converted to postscript with xpm2ps. Option -fit controls the least-squares fitting of the structures on top of each other: complete fit (rotation and translation), transla- tion only, or no fitting at all. Option -mw controls whether mass weighting is done or not. If you select the option (default) and supply a valid .tpr file masses will be taken from there, otherwise the masses will be deduced from the atommass.dat file in GMXLIB. This is fine for proteins, but not neces- sarily for other molecules. A default mass of 12.011 amu (carbon) is assigned to unknown atoms. You can check whether this happend by turn- ing on the -debug flag and inspecting the log file. With -f2, the 'other structures' are taken from a second trajectory, this generates a comparison matrix of one trajectory versus the other. Option -bin does a binary dump of the comparison matrix. Option -bm produces a matrix of average bond angle deviations analogously to the -m option. Only bonds between atoms in the comparison group are considered. FILES
-s topol.tpr Input Structure+mass(db): tpr tpb tpa gro g96 pdb -f traj.xtc Input Trajectory: xtc trr trj gro g96 pdb cpt -f2 traj.xtc Input, Opt. Trajectory: xtc trr trj gro g96 pdb cpt -n index.ndx Input, Opt. Index file -o rmsd.xvg Output xvgr/xmgr file -mir rmsdmir.xvg Output, Opt. xvgr/xmgr file -a avgrp.xvg Output, Opt. xvgr/xmgr file -dist rmsd-dist.xvg Output, Opt. xvgr/xmgr file -m rmsd.xpm Output, Opt. X PixMap compatible matrix file -bin rmsd.dat Output, Opt. Generic data file -bm bond.xpm Output, Opt. X PixMap compatible matrix file OTHER OPTIONS
-[no]hno Print help info and quit -[no]versionno Print version info and quit -nice int 19 Set the nicelevel -b time 0 First frame (ps) to read from trajectory -e time 0 Last frame (ps) to read from trajectory -dt time 0 Only use frame when t MOD dt = first time (ps) -tu enum ps Time unit: fs, ps, ns, us, ms or s -[no]wno View output .xvg, .xpm, .eps and .pdb files -xvg enum xmgrace xvg plot formatting: xmgrace, xmgr or none -what enum rmsd Structural difference measure: rmsd, rho or rhosc -[no]pbcyes PBC check -fit enum rot+trans Fit to reference structure: rot+trans, translation or none -prev int 0 Compare with previous frame -[no]splitno Split graph where time is zero -skip int 1 Only write every nr-th frame to matrix -skip2 int 1 Only write every nr-th frame to matrix -max real -1 Maximum level in comparison matrix -min real -1 Minimum level in comparison matrix -bmax real -1 Maximum level in bond angle matrix -bmin real -1 Minimum level in bond angle matrix -[no]mwyes Use mass weighting for superposition -nlevels int 80 Number of levels in the matrices -ng int 1 Number of groups to compute RMS between SEE ALSO
gromacs(7) More information about GROMACS is available at <http://www.gromacs.org/>. Mon 4 Apr 2011 g_rms(1)
All times are GMT -4. The time now is 01:37 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy