How to print median values of matrix -awk? Post: 302999599

Sponsored Content

Top Forums Shell Programming and Scripting How to print median values of matrix -awk? Post 302999599 by RudiC on Friday 23rd of June 2017 04:54:16 AM

06-23-2017

Registered User

Try this for medians:

Code:

awk -F"\t" '
NR == 1
NR > 1          {for (i=2; i<=NF; i++) print $1, i, $i | "sort -k1,2 -k3bn > TMP"
                }

function PRMED()        {printf TFS "%s", MEDIAN
                         TFS = OFS
                         PRV2 = $2
                         CNT = 0
                        }

END             {while (1 == getline < "TMP")   {if ($1 != PRV1)        {PRMED()
                                                                         printf TRS "%s", $1
                                                                         TRS = ORS
                                                                         PRV1 = $1
                                                                        }
                                                 if ($2 != PRV2)        {PRMED()
                                                                        }
                                                 M[++CNT] = $3
                                                 CH       = int (CNT / 2)
                                                 MEDIAN   = CNT%2?M[CH+1]:(M[CH]+M[CH+1])/2
                                                }
                }
END             {PRMED()
                 printf ORS
                }
' OFS="\t" file
name 	s1	s2	s3	s4
g1	5	7	8	9
g2	8	8	8	8
g3	4	4	7.5	18

This User Gave Thanks to RudiC For This Post:

RudiC

View Public Profile for RudiC

Find all posts by RudiC

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Awk to print distinct col values

Hi Guys... I am newbie to awk and would like a solution to probably one of the simple practical questions. I have a test file that goes as: 1,2,3,4,5,6 7,2,3,8,7,6 9,3,5,6,7,3 8,3,1,1,1,1 4,4,2,2,2,2 I would like to know how AWK can get me the distinct values say for eg: on col2...

2. Shell Programming and Scripting

awk to median

hi! i have a file like the attachement. you can see on the last column, there is a marker from 1 to 64 for each time. I'd like to have the median for each marker: i want to get a median every 128 values the result is : for an hour and marker x, i have the median value thank you for...

3. Shell Programming and Scripting

Print a key with its all values using awk/others

input COL1 a1 b1 c1 d1 e1 f1 C1 10 10 10 100 100 1000 C2 20 20 200 200 200 2000 output C1 a1 10 1 C1 b1 10 1 C1 c1 10 1 C1 d1 100 2 C1 e1 100 2 C1 f1 1000 3 C2 ...

4. Shell Programming and Scripting

Help fixing awk code to print values from 2 files

Hi everyone, Please help on this: I have file1: <file title="Title 1 and 2"> <report> <title>Title 1</title> <number>No. 1234</number> <address>Address 1</address> <date>October 07, 2009</date> <description>Some text</description> </report> ...

5. UNIX for Advanced & Expert Users

Awk to print values of second file

Hello, I have a data file with 300,000 records in it, and another file which contains only the line numbers of roughly 13,000 records in the data file which have data integrity issues. I'm trying to find a way to print the original data by line number identified in the second file. How can I do...

6. Shell Programming and Scripting

Print minimum and maximum values using awk

Can I print the minimum and maximum values of values in first 4 columns ? input 3038669 3038743 3037800 3038400 m101c 3218627 3218709 3217600 3219800 m290 ............. output 3037800 3038743 m101c 3217600 3219800 m290

7. Shell Programming and Scripting

How to print in awk matching $1 values ,to $1,$4 example given.?

Hi Experts, I am trying to get the output from a matching pattern but unable to construct the awk command: file : aa bb cc 11 dd aa cc 33 cc 22 45 68 aa 33 44 44 dd aa cc 37 aa 33 44 67 I want the output to be : ( if $1 match to "aa" start of the line,then print $4 of that line, and...

8. Shell Programming and Scripting

awk print values between consecutive lines

I have a file in below format: file01.txt TERM TERM TERM ABC 12315 68.53 12042013 165144 ABC 12315 62.12 12042013 165145 ABC 12315 122.36 12052013 165146 ABC 12315 582.18 12052013 165147 ABC 12316 2.36 12052013 165141 ABC 12316 ...

9. Shell Programming and Scripting

awk print odd values

value=$(some command) for all in `echo $value` do awk checks each value (all) to see if it is a odd number. if so, prints the value done sounds easy enough but i've been unable to find anything on google.

10. Shell Programming and Scripting

Print values within groups of lines with awk

Hello to all, I'm trying to print the value corresponding to the words A, B, C, D, E. These words could appear sometimes and sometimes not inside each group of lines. Each group of lines begins with "ZYX". My issue with current code is that should print values for 3 groups and only is...

LEARN ABOUT DEBIAN

pgasetrealinitpercent

PGASetRealInitPercent(2)					      PGAPack						  PGASetRealInitPercent(2)

NAME

       PGASetRealInitPercent - sets the upper and lower bounds for randomly initializing real-valued genes.

DESCRIPTION

       For  each  gene	these bounds define an interval from which the initial allele value is selected uniformly randomly.  With this routine the
       user specifies a median value and a percent offset for each allele.

INPUT PARAMETERS

       ctx    - context variable
       median - an array containing the mean value of the interval
       percent
	      - an array containing the percent offset to add and subtract to the median to define the interval

OUTPUT PARAMETERS

       none

SYNOPSIS

       #include "pgapack.h"
       void  PGASetRealInitPercent(ctx, median, percent)
       PGAContext *ctx
       double *median
       double *percent

LOCATION

       real.c

EXAMPLE

       Set the initialization routines to select a value for each real-valued
       gene i uniformly randomly from the interval [i-v,i+v], where $v = i/2$.
       Assumes all strings are the same length.

       PGAContext *ctx;
       double *median, *percent;
       int i, stringlen;
       :
       stringlen = PGAGetStringLength(ctx);
       median  = (double *) malloc(stringlen*sizeof(double));
       percent = (double *) malloc(stringlen*sizeof(double));
       for(i=0;i<stringlen;i++) {
       median[i]  = (double) i;
       percent[i] = 0.5;
       }
       PGASetRealInitPercent(ctx, median, percent);

								     05/01/95						  PGASetRealInitPercent(2)