0

hi ,

my problem isthat i have two file s with this format

file 1 has two coloumns

protein id geneid

qqqq yyyy

tttt pppp

oooo llll


now i have one other file as cluster file as


cluster 1 : yyyy,pppp
cluster 2 : llll, yyyy,
.
.
.
.
cluster n : pppp,yyyy,llll

I want to map and find out cluster which belongs to qqqq instead of yyyy,tttt instaed of pppp...... and oooo instead of llll

3
Contributors
3
Replies
4
Views
7 Years
Discussion Span
Last Post by Beat_Slayer
0

And how are those clusters. I see that you want to search for the value on the fisrt column of the first file, I just don't know what would be a valid match on the second one.

A cluster that doesn't have some values, tat has some values uin some order....

What?

Edited by Beat_Slayer: n/a

0

Can you give exact math formula and/or example with few values and right result for that limited example.

0

This should do it. Yhere could be some tweaks depending on the data input.

proteins = open('protein.txt').readlines()
clusters = open('cluster.txt').readlines()

for protein in proteins:
    proteinid, geneid = protein.split(' ')
    geneid = geneid.rstrip('\n')
    for cluster in clusters:
        clusterid, genes = cluster.split(':')
        genes = genes.rstrip('\n')
        geneids = []
        geneids += genes.split(',')
        geneids = [x.lstrip() for x in geneids]
        if geneid in geneids:
            print protein.rstrip('\n')
            print cluster.rstrip('\n')
            print
This topic has been dead for over six months. Start a new discussion instead.
Have something to contribute to this discussion? Please be thoughtful, detailed and courteous, and be sure to adhere to our posting rules.