Sponsored Content
Top Forums Shell Programming and Scripting AWK script - extracting min and max values from selected lines Post 302595307 by grincz on Thursday 2nd of February 2012 04:57:56 PM
Old 02-02-2012
AWK script - extracting min and max values from selected lines

Hi guys!
I'm new to scripting and I need to write a script in awk.

Here is example of file on which I'm working

HTML Code:
 ATOM   4688  HG1 PRO A 322      18.080  59.680 137.020  1.00  0.00            
ATOM   4689  HG2 PRO A 322      18.850  61.220 137.010  1.00  0.00            
ATOM   4690  CD  PRO A 322      18.800  60.090 135.140  1.00  0.00            
ATOM   4691  HD1 PRO A 322      17.770  60.020 134.790  1.00  0.00            
ATOM   4692  HD2 PRO A 322      19.330  60.890 134.620  1.00  0.00            
ATOM   4693  C   PRO A 322      20.020  56.920 136.410  1.00  0.00            
ATOM   4694  O1  PRO A 322      18.780  56.610 136.300  1.00  0.00            
ATOM   4695  O2  PRO A 322      20.890  56.130 136.870  1.00  0.00            
ATOM   4696  C1   TB B 323      85.140  34.010  62.880  1.00  0.00            
ATOM   4697  C2   TB B 323      84.350  35.240  62.580  1.00  0.00            
ATOM   4698  C3   TB B 323      84.790  35.750  61.220  1.00  0.00            
ATOM   4699  C4   TB B 323      83.900  36.900  60.810  1.00  0.00            
ATOM   4700  O5   TB B 323      84.420  38.000  60.510  1.00  0.00            
ATOM   4701  O6   TB B 323      82.490  36.550  60.770  1.00  0.00            
ATOM   4702  C7   TB B 323      81.780  37.770  60.540  1.00  0.00            
ATOM   4703  C8   TB B 323      80.240  37.650  60.340  1.00  0.00
I would like to obtain min and max values of column 9 from lines having TB in 4th column. I was able to extract min and max values from whole file but couldn't restrict it to the lines with TB in 4th column. Any advice?

TY in advance!
 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

max values amd min values

Hello every one, I have following data ***CAMPAIGN 1998 CONTRIBUTIONS*** --------------------------------------------------------------------------- NAME PHONE Jan | Feb | Mar | Total Donated ... (12 Replies)
Discussion started by: devmiral
12 Replies

2. UNIX for Dummies Questions & Answers

extracting selected few lines through perl

How can I extract few lines(like 10 to 15, top 10 and last 10) from a file using perl. I do it with sed, head and tail in unix scripting. I am new to perl. Appreciate your help. (2 Replies)
Discussion started by: paruthiveeran
2 Replies

3. UNIX for Dummies Questions & Answers

Awk search for max and min field values

hi, i have an awk script and I managed to figure out how to search the max value but Im having difficulty in searching for the min field value. BEGIN {FS=","; max=0} NF == 7 {if (max < $6) max = $6;} END { print man, min} where $6 is the column of a field separated by a comma (3 Replies)
Discussion started by: Kirichiko
3 Replies

4. UNIX for Dummies Questions & Answers

Iterate a min/max awk script over time-series temperature data

I'm trying to iterate a UNIX awk script that returns min/max temperature data for each day from a monthly weather data file (01_weath.dat). The temperature data is held in $5. The temps are reported each minute so each day contains 1440 temperature enteries. The below code has gotten me as far as... (5 Replies)
Discussion started by: jgourley
5 Replies

5. Shell Programming and Scripting

trying to print selected fields of selected lines by AWK

I am trying to print 1st, 2nd, 13th and 14th fields of a file of line numbers from 29 to 10029. I dont know how to put this in one code. Currently I am removing the selected lines by awk 'NR==29,NR==10029' File1 > File2 and then doing awk '{print $1, $2, $13, $14}' File2 > File3 Can... (3 Replies)
Discussion started by: ananyob
3 Replies

6. UNIX for Dummies Questions & Answers

[Solved] Print a line using a max and a min values of different columns

Hi guys, I already search on the forum but i can't solve this on my own. I have a lot of files like this: And i need to print the line with the maximum value in last column but if the value is the same (2 in this exemple for the 3 last lines) i need get the line with the minimum value in... (4 Replies)
Discussion started by: MetaBolic0
4 Replies

7. Shell Programming and Scripting

Get the min avg and max with awk

aaa: 3 ms aaa: 2 ms aaa: 5 ms aaa: 10 ms .......... to get the 3 2 5 10 ...'s min avg and max something like min: 2 ms avg: 5 ms max: 10 ms (2 Replies)
Discussion started by: yanglei_fage
2 Replies

8. Shell Programming and Scripting

awk script to find min and max value

I need to find the max/min of columns 1 and 2 of a 2 column file what contains the special character ">". I know that this will find the max value of column 1. awk 'BEGIN {max = 0} {if ($1>max) max=$1} END {print max}' input.file But what if I needed to ignore special characters in the... (3 Replies)
Discussion started by: ncwxpanther
3 Replies

9. Shell Programming and Scripting

How to get min and max values using awk?

Hi, I need your kind help to get min and max values from file based on value in $5 . File1 SP12.3 stc 2240806 2240808 + ID1_N003 ID2_N003T0 SP12.3 sto 2241682 2241684 + ID1_N003 ID2_N003T0 SP12.3 XE 2239943 2240011 + ID1_N003 ID2_N003T0 SP12.3 XE 2240077 2241254 + ID1_N003 ... (12 Replies)
Discussion started by: redse171
12 Replies

10. Shell Programming and Scripting

awk Sort 2d histogram output from min(X,Y) to max(X,Y)

I've got Gnuplot-format 2D histogram data output which looks as follows. 6.5 -1.25 10.2804 6.5404 -1.25 10.4907 6.58081 -1.25 10.8087 6.62121 -1.25 10.4686 6.66162 -1.25 10.506 6.70202 -1.25 10.3084 6.74242 -1.25 9.68256 6.78283 -1.25 9.41229 6.82323 -1.25 9.43078 6.86364 -1.25 9.62408... (1 Reply)
Discussion started by: chrisjorg
1 Replies
RODS(1) 						      General Commands Manual							   RODS(1)

NAME
rods - Raster3D preprocessor for ball-and-stick models SYNOPSIS
rods [-h] [-b] [-radius R] [-Bcolor Bmin Bmax] Rods reads a file describing atom colours and/or a PDB coordinate file and produces a file containing Raster3D descriptor records. The file produced by rods may be fed directly to render or it may be combined with descriptor files produced by other Raster3D utilities. EXAMPLES
To describe a simple bonds-only model coloured by residue type: cat mycolours.pdb protein.pdb | rods | render > mypicture.png To render the same molecule as ball-and-stick: cat mycolours.pdb protein.pdb | rods -b | render > mypicture.png OPTIONS
-h Suppress header records in output. By default rods will produce an output file which starts with header records containing a default set of scaling and processing options. The -h flag will suppress these header records so that the output file contains only object descrip- tors. This option is useful for producing files which describe only part of a scene, and which are to be later combined with descriptor files produced by other programs. -b By default rods will describe bonds only; the -b flag will cause it to include spheres at the atom positions also, yielding a ball-and- stick representation. -radius R By default rods will draw bonds as cylinders with a 0.2A radius. The radius option allows you to change this cylindrical radius. -Bcolor Bmin Bmax Assign colors based on B values rather than from atom or residue types. Atoms with B <= Bmin will be colored dark blue; atoms with B >= Bmax will be colored light red; atoms with Bmin < B < Bmax will be assigned colors shading smoothly through the spectrum from blue to red. DESCRIPTION
The input to rods consists of a single text file containing colour information and atomic coordinates in PDB data bank format. Coordinates are output as Raster3D descriptor records with colours and sphere radii assigned according to the COLO records described below. Ball-and- stick figures have atoms drawn at 0.2 * VanderWaals radius, connected by rods with a default 0.2A cylindrical radius. Bonds are drawn for atoms which lie closer to each other than 0.6 * (sum of VanderWaals radii). By default the output file contains a set of header records as required by the render program. The header is constructed to include a TMAT matrix corresponding to the transformation matrix contained in file setup.matrix (if it exists), or to the Eulerian angles contained in file setup.angles (if it exists). Colours are assigned to atoms using a matching process, using COLOUR records prepended to the input PDB file. If no COLOUR records are present in the input file, atoms will receive default CPK colors (C=grey, O=red, N=blue, S=yellow, P=green, other=magenta). Raster3D uses a pseudo-PDB record type with the same basic layout as the above but with COLO in the first 4 columns: Columns 1 - 4 COLO 7 - 30 Mask (described below) 31 - 38 Red component 39 - 46 Green component 47 - 54 Blue component 55 - 60 van der Waals radius in Angstroms 61 - 80 Comments Note that the Red, Green, and Blue components are in the same positions as the X, Y, and Z components of an ATOM or HETA record, and the van der Waals radius goes in place of the Occupancy. The Red, Green, and Blue components must all be in the range 0 to 1. The Mask field is used in the matching process as follows. First the program reads in and stores all the ATOM, HETA, and COLO records in input order. Then it goes through each stored ATOM/HETA record in turn, and searches for a COLO record that matches the ATOM/HETA record in all of columns 7 through 30. The first such COLO record to be found determines the colour and radius of the atom. In order that one COLO record can provide colour and radius specifications for more than one atom (e.g., based on residue or atom type, or any other criterion for which labels can be given somewhere in columns 7 through 30), the "#" symbol is used as a wildcard. I.e. a # in a COLO record matches any character in the corresponding column in an ATOM or HETA record. All other characters must match literally to count as a match. Note that the very last COLO record in the input should have # symbols in all of columns 7 through 30 in order to pro- vide a colour for any atom whose ATOM/HETA record fails to match any previous COLO record. This idea of matching masks for colour specifi- cations is due to Colin Broughton. ENVIRONMENT
The files setup.matrix and setup.angles, if they exist, affect the header records produced by rods. SOURCE
web URL: http://www.bmsc.washington.edu/raster3d/raster3d.html contact: Ethan A Merritt University of Washington, Seattle WA 98195 merritt@u.washington.edu SEE ALSO
render(l), ribbon(l), balls(l) AUTHORS
Ethan A Merritt Raster3D 8 May 1999 RODS(1)
All times are GMT -4. The time now is 04:04 PM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy