I don't know what's wrong with this code. I want to count from a blast report, how many hits correspond to a specific species (the keyword). So my idea is to loop through each line and when the keyword is found, add 1 to count and go to the next line, this because sometimes the name of the species is present more than once in the subject name.
#!/usr/bin/python
import sys
##usage: python filterbyword.py file keyword
file = open(sys.argv[1],'r')
keyword = sys.argv[2]
count = 0
for line in file:
while True:
if keyword not in line:
continue
else:
break
count = count + 1
print count
3 answers
The problem is, specifically, that your while loop is inside the for loop.
So if the keyword is not in the line, then your script keeps iterating through the while loop, and it appears to "hang".
All you have to do is take out the while loop and test directly:
for line in file:
if keyword in line:
count = count + 1
The test should only be done once per line, but you can verify this by trying it out with test input that contains multiple instances of a keyword on a line, and seeing if the final count matches what you expect.
•
written
by
Alex Reynolds
What is happening with this script ? If it not looping again for next search, try to useseek()method. Something like file.seek(0)
•
written
by
GouthamAtla
Why would you seek to the start of a file on each iteration?
You are incrementing count regardless of whether keyword is in line or not. Maybe you want:
#!/usr/bin/python
import sys
##usage: python filterbyword.py file keyword
file = open(sys.argv[1],'r')
keyword = sys.argv[2]
count = 0
for line in file:
while True:
if keyword not in line:
continue
else:
count += 1
break
print count
Or even simpler:
for line in file:
if keyword in line:
count+=1
print count
The simpler one is what I was using but then I got confused about whether there would be an addition for each time the word appears in the line. I don't want that, I want to have one count for each line with the word, not a count for each time the word is in the file.
Log in to answer this question.
-
Parallelize microRNA targets
written by aranyak111 •
-
Import files in newick format with python script
written by ArusjakGevorgyan •
-
Removing gap using array
written by skjobs1234 •
-
Arranging header in nr database via Python (Unable to input multiple files)
written by AsoInfo
-
Extracting fasta sequences with identical header
written by seta
-
Python IndexError: list index out of range
written by larafrommer •
-
Extracting Accession Numbers With Ncbixml.Parse
written by jed.tressler •
-
How To Extract Desire Genes Blast Xml Result From A Big Blast Xml File
written by KJ Lim •
-
Converting Python Code To R
written by ruchiksy •
-
Find Sequences In Fasta File From List Of Id'S
written by arronslacey
what's wrong with `grep` ?
I was gonna ask that (along with why not awk for more advanced grepping), but maybe OP wishes to add this functionality as a module to existing code?
OMG you are right! How didn't I think about grep. Shame on me!!