Make multiple awk files into an executable Post: 302604972

Sponsored Content

Top Forums Shell Programming and Scripting Make multiple awk files into an executable Post 302604972 by jacobs.smith on Tuesday 6th of March 2012 11:37:34 AM

03-06-2012

Banned

Quote:

Originally Posted by birei

Give a try to next awk script:

Code:

$ cat 1.txt
chr1 14765298 14766727 def
chr1 16759093 16760238 def
chr1 16759236 16760238 def
chr1 20782516 20784428 him
chr1 20989962 20991078 her
chr2 31672150 31673532 abc
chr2 33157721 33158124 abc
chr3 34542283 34542962 abc
chr3 38248682 38251416 abc
chr4 58562053 58567653 abc
$ cat 2.txt
chr1 21438731 21439423 26.12
chr1 33939851 33940673 34.76
chr1 36779864 36780494 20.16
chr1 36817091 36817917 27.22
chr2 36977015 36977908 19.27
chr3 40475125 40475885 21.58
chr3 40483838 40484616 15.3
chr4 40502827 40503675 10.61
chr4 40532299 40533156 14.78
chr5 43593022 43594143 24.33
$ cat script.awk
BEGIN {
        if ( ARGC != 3 ) {
                print "Usage: awk -f script.awk <file1> <file2>"
                exit 0
        }
}

FNR == NR {
        f1_data[ FNR ] = $0
        next
}

FNR < NR {
        for ( i = 1; i <= length( f1_data ); i++ ) {
                split( f1_data[ i ], fields )
                if ( fields[1] != $1 ) {
                        next
                }
                substraction = fields[2] - $2
                if (substraction >= -1000000 && substraction <= 1000000 ) {
                        for ( j = 2; j <= length( fields ); j++ ) {
                                f1_line = (f1_line ? f1_line " " : "" ) fields[j]
                        }
                        printf "%s %s %d\n", $0, f1_line, substraction
                        f1_line = ""
                }
        }
}
$ awk -f script.awk 1.txt 2.txt 
chr1 21438731 21439423 26.12 20782516 20784428 him -656215
chr1 21438731 21439423 26.12 20989962 20991078 her -448769

Hi Birei,

Thanks for your time.

But, it is not producing any output. All I get is a blank output. I did exactly what you have written.

jacobs.smith

View Public Profile for jacobs.smith

Find all posts by jacobs.smith

10 More Discussions You Might Find Interesting

1. Solaris

How to make a script executable by all users?

I have a script in my home direcroty which upon execution gives the essential system information like memory,cpu etc and is currently owned by root:root. Now I want to see that every non root user will run this file and grab the reqired system info. I know this is some thing associated with chown...

2. Shell Programming and Scripting

Can I make "touch" create executable files by manipulating umask?

I'm getting to grips with this concept of the umask. What I thought was, setting umask uga+rwx would result in creating files with all permissions for everyone. Seems not to be the case though. Read and write bits get set, but not the execute bit. Is there some gap in my understanding, or is...

3. Shell Programming and Scripting

run and make an executable file

as i said before i'm a beginner in shell programming and i have two questions: how to run an executable file in shell scripts like for example let's say the file called "prog.exe", what's the shell command to run this file? also how can i make the shell file an executable file (if it is...

4. Shell Programming and Scripting

Unable to make script executable

Hello everybody, I'm unable to make my shell script an executable file. The details are as follows: PATH includes my $HOME/bin i.e. /rchome/rc1/bin HOME directory is /rchome/rc1 script name is prep_mig.sh permissions set are 755 It's executing if I give below command sh prep_mig.sh but...

5. Shell Programming and Scripting

How to make an editing script work for multiple files?

Hey everybody, I have a script for making a string substitution in a file. I am trying to modify it in order to make the same modifcation to multiples files. here is what I have so far. #!/bin/csh set p1="$1" shift set p2="$1" shift foreach x ($*) if ( { grep -w -c "$p1" $x } ) then mv...

6. UNIX for Dummies Questions & Answers

Using AWK: Extract data from multiple files and output to multiple new files

Hi, I'd like to process multiple files. For example: file1.txt file2.txt file3.txt Each file contains several lines of data. I want to extract a piece of data and output it to a new file. file1.txt ----> newfile1.txt file2.txt ----> newfile2.txt file3.txt ----> newfile3.txt Here is...

7. Shell Programming and Scripting

Make multiple files of equal length

I have 150 files with 4 columns each but variable row lengths that I need to combine by column. I do not have any common column. I want to use "paste " command in unix to do it but before that I have to get all my files to be of equal length. Is there a way using awk or sed to fill up n no. of...

8. SuSE

How To make bootable USB with multiple ISO Files?

Hi All, I would need your assistance to make a bootable USB with SUSE LINUX Enterprise Server I have already downloaded relevant OS (Trail Version) packages @ 1) SLES-11-SP4-DVD-i586-GM-DVD1 2) SLES-11-SP4-DVD-i586-GM-DVD2 when I tried to open these packages with PowerISO one of the...

9. Shell Programming and Scripting

Bash script make itself executable

Is there a way to make this make itself executable? Thanks. :-) cat > somescript.sh << \EOF #!/bin/bash block_count=$(sudo tune2fs -l /dev/sda1 | awk '/^Block count:/ {print $NF}') reserved_block_count=$(sudo tune2fs -l /dev/sda1 | awk '/^Reserved block count:/ {print $NF}') perl -e...

10. Shell Programming and Scripting

Xargs to call python executable to process multiple bam files

I am running the below loop that to process the 3 bam files (which isn't always the case). A .py executable is then called using | xargs sh to further process. If I just run it with echo the output is fine and expected, however when | xargs sh is added I get the error. I tried adding | xargs...

LEARN ABOUT DEBIAN

numm.one-bit-instrument

one-bit-instrument(7)						  Numm Tutorials					     one-bit-instrument(7)

NAME

       one bit instrument - how to make a gnarly synth with numm-run

SYNOPSIS

       numm-run FILE

DESCRIPTION

       In  this  tutorial we will learn how to live-code a one-bit synthesizer controlled by mouse position.  It is intended as a gentle introduc-
       tion to development with numm-run.

       To get started, create a text file with the following method stubs:

	      def audio_out(a):
		  pass
	      def video_out(a):
		  pass

       Save the file as onebit.py, and then launch it with numm-run:

	      numm-run onebit.py

       You should see a blank window appear.  We will now make some sound and light by changing the value of a	in  the  audio_out  and  video_out
       functions:

	      def audio_out(a):
		  a[::100] = 2**15
	      def video_out(a):
		  a.flat[::100] = 255

       Save  the  file,  and you should see and hear the sketch update.  This is using numpy's array-indexing to turn every hundredth audio sample
       and pixel-color on.  The audio sample rate is by default 44100Hz, so it produces a series of clicks that will be perceived as a 441Hz tone.
       We can turn this into an instrument by connecting mouse motion to frequency:

	      period = 100
	      def audio_out(a):
		  a[::period] = 2**15
	      def video_out(a):
		  a.flat[::period] = 255
	      def mouse_in(type,px,py,button):
		  global period
		  period = px*1000

       Finally, let's use the keyboard to record and jump to notes.  The first time you press a key, it records the period, and subsequant depres-
       sions play the saved period:

	      period = 100
	      record = {}
	      def audio_out(a):
		  a[::period] = 2**15
	      def video_out(a):
		  a.flat[::period] = 255
	      def mouse_in(type,px,py,button):
		  global period
		  period = px*1000
	      def keyboard_in(type,key):
		  global period
		  if record.has_key(key):
		      period = record[key]
		  elif record.has_key(key):
		      record[] = period

SEE ALSO

       numm-run(1), numm.getting-started(7), numm.spectral-analysis(7)

numm								   February 2012					     one-bit-instrument(7)

10 More Discussions You Might Find Interesting

1. Solaris

How to make a script executable by all users?

Discussion started by: chrs0302

2. Shell Programming and Scripting

Can I make "touch" create executable files by manipulating umask?

Discussion started by: tphyahoo

3. Shell Programming and Scripting

run and make an executable file

Discussion started by: _-_shadow_-_

4. Shell Programming and Scripting

Unable to make script executable

Discussion started by: jitu.keshwani