awk if then else statement


 
Thread Tools Search this Thread
Top Forums Shell Programming and Scripting awk if then else statement
# 1  
Old 04-02-2015
awk if then else statement

I am using awk as part of and if then else statement. I am trying to have the user enter a gene name and then a variant or variants and have a specific transcript assigned to the variants depending on the gene. Below is my code but the if then else statement is wrong. Basically, the gene name would be the condition that assigns a particular transcript to a variant.

The desired output would be, if the user enters GJB2 as the gene name and then c.283G>A for the variant there is a statement for that condition that saves a file with NM_004004.5:c.283G>A . However if MECP2 is entered as the gene name and c.48delC,c.255G>T are entered as the variants, the a file is saved:

Code:
 
NM_004992.3:c.48delC
NM_004992.3:c.255G>T

Thank you Smilie.

Code:
 syntax() {
    printf "\n\n"
    printf "Please enter the gene name  : "; read gene
    printf "Enter variant(s): "; IFS="," read -a variants
        
        [ -z "$gene" ] && printf "\n No ID supplied. Leaving match function." && sleep 2 && return
        [ "$gene" = "end" ] && printf "\n Leaving match function." && sleep 2 && return

		awk '{ 
		if(${gene} = "GJB2")
        for ((i=0; i<${#variants[@]}; i++))
              do printf "NM_004004.5:%s\n" "${variants}" >> c:/Users/cmccabe/Desktop/Python27/$gene.txt >> c:/Users/cmccabe/Desktop/Python27/out.txt;
	    done
		else if(${gene} = "MECP2")
        for ((i=0; i<${#variants[@]}; i++))
		      do printf "NM_004992.3:%s\n" "${variants}" >> c:/Users/cmccabe/Desktop/Python27/$gene.txt >> c:/Users/cmccabe/Desktop/Python27/out.txt;
        done
		else if(${gene} = "Phox2B")
        for ((i=0; i<${#variants[@]}; i++))
		      do printf "NM_003924.3:%s\n" "${variants}" >> c:/Users/cmccabe/Desktop/Python27/$gene.txt >> c:/Users/cmccabe/Desktop/Python27/out.txt 
	    }'
	name
}

# 2  
Old 04-02-2015
From the syntax of your read statement, I'm assuming you're using bash... But, in bash and in ksh, the way to reference then nth element of an array is ${array[$n]}; using ${array} just gives you the 0th element of the array.

Shell variables (whether they are scalars or arrays) can't be referenced by name in an awk script.

You can't write (overwrite or append) to two files using redirection in the shell or in awk. Either call printf twice (redirecting the output of each to one of your files), or use tee to make a copy in one file while directing the output to the other file.

Why would you invoke awk when there are no calculations being performed; why not just use shell if, for, and printf statements?

I'm surprised that you don' t want to print some kind of diagnostic if the entered gene name is not an empty string is not end, and is not one of your three expected gene names, but this seems to do something similar to what I think you were trying to do:
Code:
syntax() {
	printf "\n\n"
	printf "Please enter the gene name  : "; read gene
	[ -z "$gene" ] &&
		printf "\n No ID supplied. Leaving match function." &&
		sleep 2 &&
		return
	[ "$gene" = "end" ] &&
		printf "\n Leaving match function." &&
		sleep 2 &&
		return
	printf "Enter variant(s): "; IFS="," read -a variants

	code=""
	if [ "$gene" = "GJB2" ]
	then	code="NM_004004.5"
	elif [ "$gene" = "MECP2" ]
	then	code="NM_004992.3"
	elif [ "$gene" = "Phox2B" ]
	then	code="NM_003924.3"
	else	return
	fi
	for ((i = 0; i < ${#variants[@]}; i++))
	do	printf "%s:%s\n" "$code" "${variants[$i]}"
	done | tee -a c:/Users/cmccabe/Desktop/Python27/$gene.txt \
			>> c:/Users/cmccabe/Desktop/Python27/out.txt
}

These 2 Users Gave Thanks to Don Cragun For This Post:
# 3  
Old 04-03-2015
Thank you Smilie. I didn't even think about
Quote:
if the entered gene name is not an empty string is not end , and is not one of your three expected gene names
, but that is a good point that I will build in as well. Thanks.
Login or Register to Ask a Question

Previous Thread | Next Thread

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

awk statement help

There has to be a way to do this with awk or maybe I'm just focusing on the wrong tool and making this harder than it needs to be. I'm trying to do a file field lookup/join at a very large scale but the output changes has to change dramatically. I have an input file to do a field lookup from and... (3 Replies)
Discussion started by: brettcasper
3 Replies

2. Shell Programming and Scripting

Convert Update statement into Insert statement in UNIX using awk, sed....

Hi folks, I have a scenario to convert the update statements into insert statements using shell script (awk, sed...) or in database using regex. I have a bunch of update statements with all columns in a file which I need to convert into insert statements. UPDATE TABLE_A SET COL1=1 WHERE... (0 Replies)
Discussion started by: dev123
0 Replies

3. Shell Programming and Scripting

If statement in awk

I run my script "switch.sh" repeatedly (within 30 seconds). Each time script is triggered, script itself should kill all previous process. Here is my code: for pid in $(ps -fe | grep 'switch.sh' | grep -v grep | awk '{if ($2<$$) print $2}'); do sudo kill -9 $pid done sleep 30 ... (6 Replies)
Discussion started by: armatron
6 Replies

4. Shell Programming and Scripting

Help with awk statement

Hi I have a file with test test2 1000000657373 test1 test3 1000003849450 test2 test4 test5 100000837474 I cat the file and pipe it to an awk statement like so cat /tmp/file |awk '{if ($3 ~ "^*$" && $3 > 1024000000) print "/vol/"$1"/"$2;else if ($4 ~ "^*$" && $4 > 1024000000) print... (15 Replies)
Discussion started by: bombcan
15 Replies

5. Shell Programming and Scripting

awk if statement

hi, I have a l-column file of more than 10,000 lines with interspersed negative values. What I want to do is add a fixed number (360) everytime a negative value is encountered while leaving the positive ones as is. I need something that will read every line of the file and do the calculation... (3 Replies)
Discussion started by: ida1215
3 Replies

6. Shell Programming and Scripting

Help Regarding AWk and IF THEN ELSE Statement

Hi, I have a data file which contains record count. So doing wc -l rightfit_balancing_count.dat | awk '{print $1}'] gives me the record count stored in the file. Now, i want to send a mail from UNIX, if the record count is equal to 0,otherwise it should do nothing. Any help... (2 Replies)
Discussion started by: Shell_Learner
2 Replies

7. Shell Programming and Scripting

Awk 'if' statement help

Hello all, I'm very new to this (<5hrs!) please bear with me. My file looks like this 386259.448541 417069.155 154935.157 186.206 162 1 1 8 386259.448551 417068.53 154935.04 186.144 156 1 1 8 386259.448561 417067.911 154934.926 186.175 164 1 1 8 386259.450337 417086.643 154946.483 894.671... (4 Replies)
Discussion started by: rebeccab37
4 Replies

8. Shell Programming and Scripting

awk, if statement

Having a little trouble with awk and an if statement. I have a test setup which I am trying to only print the records which start with the month 03. Everything I tried, prints everything, even the 02 month 03/23/2010 12:47:51 ga2828 SUBMITTED FROM URL: test123.cgi show port count ... (2 Replies)
Discussion started by: numele
2 Replies

9. Shell Programming and Scripting

awk inside another awk statement

hi all, i have two files 1) a.txt one two three 2) abc "one" = 10 pqr "three" = 20 345 "two" = 0 this is what i want in third file (3 Replies)
Discussion started by: shishirkotkar
3 Replies

10. UNIX for Dummies Questions & Answers

if statement in awk

Hi Friends How do I do two things from one if statement inside awk? I want to run a script and create a new file from the same condition. awk '{ if ($2 == ""){print " "|"cd /local/test; ./script.ksh"}{cat > ran_true.txt}}' $IN_FILE Bolded are the two things I want to be done. Thanks (1 Reply)
Discussion started by: UNovIX
1 Replies
Login or Register to Ask a Question