Sponsored Content
Top Forums UNIX for Dummies Questions & Answers Interpolation if there is no exact match for value Post 302900893 by Akshay Hegde on Friday 9th of May 2014 10:43:22 AM
Old 05-09-2014
Code:
akshay@Aix:~/Desktop$ cat f1
BP         CM
752566 rs3094315
752721 rs3131972
753541 rs2073813
760300 rs11564776
12222222

Code:
akshay@Aix:~/Desktop$ cat f2
BP         CM
740857 1.984065
750235 2.009238
752566 2.012958
753269 2.013806
753541 2.014133
754745 2.014728
765948 2.019205
767038 2.019743
767070 2.019759
768448 2.020467
769551 2.020969
771521 2.021869

Code:
awk '

function find(Arr,key,arg,   low,upp,i){
			      low = upp = ""
			
			      for(i in Arr){
				            if(i<key)
						 {
						   low =  (low && low > i) ? low : i
						 }
					    if(i>key)
						 {   
					           upp =  (upp && upp < i) ? upp : i
						 }  	 
					     }
			    
			 low  =  (low != "") ? low : "-"
		         upp  =  (upp != "") ? upp : "-"
			 
			if(arg=="up"){return upp}
		   else if(arg=="dw"){return low}	 
}

function lin(x1,x3,y1,y3, x2){
				return (y1+(x2-x1)*(y3-y1)/(x3-x1))
			     }

 FNR==NR{ 
		A[$1] = $2
                next
        }
        {
		if($1 in A)
		{	
			# Exact match
			print $1,A[$1]
		} 
		else
		{
			# No match will interpolate
			x1=find(A,$1,"dw")
			x3=find(A,$1,"up")

		   # Whether we have upper and lower limit ?
		   if(x1 !="-" && x3 !="-")
			{
				# Yes interpolate
				y1=A[x1]
                                y3=A[x3] 
				print $1,lin(x1,x3,y1,y3, $1)
			}
			else
			{
				# No cannot interpolate so NaN
				print $1,"NaN"
			}
		}
	}
  ' f2 f1



Code:
akshay@Aix:~/Desktop$ ./interpolate 
BP CM
752566 2.012958
752721 2.01314
753541 2.014133
760300 2.01695
12222222 NaN


Last edited by Akshay Hegde; 05-09-2014 at 12:18 PM.. Reason: Bug fix
These 3 Users Gave Thanks to Akshay Hegde For This Post:
 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

perl exact match

How to emulate grep -o option in perl. I mean to print not all line, only the exact match. echo "2A2 BB" | perl -ne 'print if /2A2/' 2A2 BB I want to print only 2A2. (2 Replies)
Discussion started by: mirusnet
2 Replies

2. Shell Programming and Scripting

exact string match ; search and print match

I am trying to match a pattern exactly in a shell script. I have tried two methods awk '/\<mpath${CURR_MP}\>/{print $1 $2}' multipath perl -ne '/\bmpath${CURR_MP}\b/ and print' /var/tmp/multipath Both these methods require that I use the escape character. I am guessing that is why... (8 Replies)
Discussion started by: bash_in_my_head
8 Replies

3. Shell Programming and Scripting

Exact match question

Hi, I have a file like follows . . . White.Jack.is.going.home Black.Jack.is.going.home Red.Jack.is.going.home Jack.is.going.home . . . when I make: cat <file> | grep -w "Jack.is.going.home" it gives: White.Jack.is.going.home Black.Jack.is.going.home Red.Jack.is.going.home... (4 Replies)
Discussion started by: salih81
4 Replies

4. Solaris

grep exact match

Hi This time I'm trying to grep for an exact match e.g cat.dog.horse.cow.bird.pig horse.dog.pig pig.cat.horse.dog horse dog dog pig.dog pig.dog.bird how do I grep for dog only so that a wc -l would result 2 in above case. Thanks in advance ---------- Post updated at 06:33 AM... (4 Replies)
Discussion started by: rob171171
4 Replies

5. Shell Programming and Scripting

Exact match and #

Hi friends, i am using the following grep command for exact word match: >echo "sachin#tendulkar" | grep -iw "sachin" output: sachin#tendulkar as we can see in the above example that its throwinng the exact match(which is not the case as the keyword is sachin and string is... (6 Replies)
Discussion started by: neelmani
6 Replies

6. Shell Programming and Scripting

Match exact and append zero

file 11 2 12 6 13 7 114 6 011 7 if I'm searching for 11, output needed is output: 11 2 011 7 Code: awk '$1 ~ /^11$/' file I used the above to match exact, but it avoiding "011 7" line too, how to resolve this? (6 Replies)
Discussion started by: Roozo
6 Replies

7. Shell Programming and Scripting

Exact match using sed

I would like replace all the rows in a file if a row has an exact match to number say 21 in a tab delimited file. I want to delete the row only if it has 21 any of the rows but it should not delecte the row that has 542178 or 563421. I tried this sed '/\<21\>/d' ./inputfile > output.txt ... (7 Replies)
Discussion started by: Kanja
7 Replies

8. UNIX for Dummies Questions & Answers

Exact match question

Hi guys, I am using Centos 6.3. Actually I posted similar question but I still have some minor problem need be fixed. I have two files, file1:target: gi|57529786|ref|NM_001006513.1| mfe: -31.4 kcal/mol p-value: 0.006985 target: gi|403048743|ref|NM_001271159.1| mfe: -29.6 kcal/mol p-value:... (11 Replies)
Discussion started by: yuejian
11 Replies

9. Shell Programming and Scripting

Replacing exact match

Hi All, My Input file contains a 1000’s of lines in which I have to replace a a string to the other. Here the problem is, I have the lines in my Input as below. Cable Yes && !Pay TV && !ADS \noUE \Label="Cable Yes && !Pay TV && !ADS" I want to replace exactly the string Cable Yes &&... (37 Replies)
Discussion started by: am24
37 Replies

10. Shell Programming and Scripting

awk to update file based on partial match in field1 and exact match in field2

I am trying to create a cronjob that will run on startup that will look at a list.txt file to see if there is a later version of a database using database.txt as the source. The matching lines are written to output. $1 in database.txt will be in list.txt as a partial match. $2 of database.txt... (2 Replies)
Discussion started by: cmccabe
2 Replies
All times are GMT -4. The time now is 08:45 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy