Skip to main content

Column subtractor - Python

I am still struggling with the after effects of Linux reinstallation on my system.
Two major problems:
1) Firefox and openoffice is still not in English. Its giving me hard time.
2) The functions in openoffice  version of excel is not working.

I am yet to sort that out. In the mean time, I desperately needed it to do simple task. subtracting values of one column from one in another column and write the difference in next column.

So, I had input in this form:

------------------


0    2.179139    0.951816
1    2.314531    1.006480
2    2.191692    0.951938
3    2.079704    0.952747
4    2.092761    0.939545
5    2.362056    0.984791
6    2.153595    0.957934
7    2.213173    0.951768
8    2.175687    0.953091

--------------------

I needed to subtract  column 3 data from column 2 data. So, I tried this:

---------------------------
#!/usr/bin/env python

fo=open('name.txt' +'_output','w')
with open('name.txt') as fd:
 for line in fd:
   columns=line.split()
   columns=map(float,columns)
   fo.write ("%s \t %s \n" %(line.strip(), columns[1] - columns[2]))

-----------------------------

It creates an output with difference written in next column :


---------------------------------
0    2.179139    0.951816      1.227323
1    2.314531    1.006480      1.308051
2    2.191692    0.951938      1.239754
3    2.079704    0.952747      1.126957
4    2.092761    0.939545      1.153216
5    2.362056    0.984791      1.377265
6    2.153595    0.957934      1.195661
7    2.213173    0.951768      1.261405
8    2.175687    0.953091      1.222596

---------------------------------

It works fine, only problem is it ends with an error:

IndexError: list index out of range

Although it doesn't effect the result.
So, for the moment I got the result I wanted but I will come back and take care of this small issue.








Comments

Popular posts from this blog

Material Studio - selective molecule movement

So, working with Material Studio is awesome. You get to handle bigger, complicated and even unimaginable system. But I am facing problems in very minute functioning. I recently faced two problems, but managed to resolve on my own. Somehow, tutorials doesn't help in resolving them and you end up doing some exercise on your own to get the solution. Among these, the problem of bringing two molecular structures in one window and then adjusting one molecule by moving it around, while other structure remains untouched, as it is. Actually, this should not be a problem to do. There is a simple way of doing this. Select the molecule to be moved and then there is a symbol (button) to select, which allows you to selectively move that molecule only. This is the same way it can be done in Material studio too.  But among all the dozens of symbols, the new person handling MS, kind of gets lost. So, what he does is check the tutorials, for any help. But the tutorials don't talk about suc...

How would COVID bahave in August - few projections

Entering August with COVID by side   Entering in the month of August, COVID is still here by our side. Situation seems grim, worse than earlier, as number of active cases as well as the count of deceased patients are increasing day by day. Yesterday we crossed the mark of 50,000 cases in a day. Closely competing with Brazil with US counts also in range now. On the other side, few other numbers like reduced mortality rate to 2.18%, recovery rate of 64.6% and slightly increased doubling rate of 21 days can reduce the anxiety. My projections of 16.75 lakhs for India might be crossed by coupe of thousands of more cases after today’s count is added. This number would reach to a huge count of 28.5 lakhs by August 15. Maharashtra just doesn’t like competition. Just like what Delhi tried few days back, Andhra Pradesh also managed to cross Maharashtra’s per day count for2 consecutive days. Seems like not a good idea. Maharashtra immediately went to the next zone. The state which w...

Rigorous method to create lib file

Sometimes you totally fail to go ahead with the prepin and frcmod files you have generated from the newly generated PDB file even after you follow the procedure given in previous two posts. And you keep wondering? In such a case you need a more detailed procedure. First take a look at the pdb file. For example: Here is one such pdb file generated from the gabedit that troubled me for some time. --------------- HEADER    PROTEIN COMPND    UNNAMED AUTHOR    GENERATED BY GABEDIT 2.3.8 at Tue Jul 26 16:53:48 2011 ATOM  1      C   C     1       -2.907  -0.997  -0.385  1.00  300.00          C 0.0000   ATOM  2      N   N     1       -2.060  -0.981  -1.189  1.00  300.00 ...