Map multiple annotations using pandas 


A simple pandas solution to map multiple annotations for a protein. 


A protein or gene file will have annotations curated by different methods. Most frequently, biologists will encounter more than one annotation for a single protein. It is a task in itself to pick the right annotation. 


One of the simple ways is to consolidate them and pick the right ones after enough evidence is known. 


Coding in matlab or other languages may require more number of lines to achieve the same output. Here, a simple 'groupby' of pandas can do produce the same outputin seconds !


# Map multiple annotation for protein ID and join with a delimiter
'''
 Sample input
 A_xx  Annotation1
 A_xx  Annotation2
 B_xx  Annotation1 

 Sample output
 A_xx Annotation1, Annotation2
 B_xx Annotation1

'''
import pandas as pd
data = pd.read_csv('input_file.txt', delimiter='\t')
dfc = data.groupby(['PROTID'])['ANNOTATION'].apply(", ".join)
print dfc
dfc.to_csv('Output.txt', sep='\t', encoding='utf-8')
Location: Szeged, Hungary

Comments

Post a Comment

Popular posts from this blog

Pick Matching lines with list of keywords

#Simple code to find the occurrence of list of search terms in single line in a huge file. Search_terms = [ 'A' , 'B' , 'C' , 'D' ] with open ( 'BigFile.txt' , 'r' ) as infile : entries = infile . read () each_line = entries . splitlines () new_list = [] for ix , row in enumerate ( each_line ): element = row . split ( "\t" ) if set ( Search_terms ) == set ( element ): #Use <= if you want no strict option new_list . append ( element ) print ix + 1 #List the matching line number ! thefile = open ( 'Output.txt' , 'w' ) for item in new_list : print >> thefile , item
Read more

Presence-absence Matrix to Fasta format

Convert the binary matrix to fasta format with this simple code in Python !. Recommended for larger file sizes ! Sample Input: Sp1 1 1 1 1 Sp2 0 0 0 0 Sp3 1 1 0 0 Sample Output: >Sp1 ['1', '1', '1', '1'] >Sp2 ['0', '0', '0', '0'] >Sp3 ['1', '1', '0', '0'] Workaround the output file in any text file to remove [,'  ] to generate final output as: >Sp1 1111 >Sp2 0000 >Sp3 1100 __author__ = 'Arun Prasanna' ''' This is a simple python code to convert binary matrix into fasta format. There  are many public softwares available. But each one has size limits (<=2MB).  This code processes 61 x 51000 character matrix in less than 10 seconds in  python 2.7! ''' with open ( 'infile_matrix-cp' , 'r' ) as infile : entries = infile . read () . strip () each_line = entries . ...
Read more