how to access the individual elements of a matrix in python

There are powerful ways to do this all in one line, but by way of explanation, start by using split() to separate each line into individual data lists for further manipulation.

Each of which could be appended to an empty list, forming a multi-dimensional data set. Note that all the elements are strings and would have to be converted to numbers before the math.

my file is of the form

01 "\t" 10.19 "\t" 0.00 "\t" 10.65
02 "\t" 11.19 "\t" 10.12 "\t" 99.99


and i need to access the individual floating point numbers from it!
say for ex. the first no is 10.19.. i want to access this and add one to it.

Expand|Select|Wrap|Line Numbers

1. filename=open("half.transfac","r")
2. file_content=filename.readlines()
3. sam=""
4. for line in file_content:
5.     for char in line:
6.         if char=="\tchar\t\n":
7.             sam+=char
8.             print sam
9.  

for char accesss every digit and not the numbers{"10.19","0.00")etc.. how do i do this..help

my file is of the form

01 "\t" 10.19 "\t" 0.00 "\t" 10.65
02 "\t" 11.19 "\t" 10.12 "\t" 99.99


and i need to access the individual floating point numbers from it!
say for ex. the first no is 10.19.. i want to access this and add one to it.

Expand|Select|Wrap|Line Numbers

1. filename=open("half.transfac","r")
2. file_content=filename.readlines()
3. sam=""
4. for line in file_content:
5.     for char in line:
6.         if char=="\tchar\t\n":
7.             sam+=char
8.             print sam
9.  

for char accesss every digit and not the numbers{"10.19","0.00")etc.. how do i do this..help

Here is another way to access the numbers from a dictionary: You can perform mathematical operations on elements of the dictionary: To access individual characters: OR

THis is my file! AS u can see the first column gives the position,the second third ..etc are the part of the matrix.my problem here is that.I have to write a code to access the individual values.like if i say A[01] i should have the value 0.00 or say A[06] i should have a value 3.46. or C[04]=3.67 etc. and i have to add one to each of these elements.can you gimme the code. my code is not working.

NA bap

PO A C G T

01 0.00 3.67 0.00 0.00

02 0.00 0.00 3.67 0.00

03 0.00 0.00 0.00 3.67

04 0.00 3.67 0.00 0.00

05 3.67 0.00 0.00 0.00

06 3.46 0.00 0.22 0.00

07 0.00 0.00 3.67 0.00

08 0.00 0.00 0.00 3.67

09 0.00 0.00 0.00 3.67

10 0.00 3.67 0.00 0.00

11 3.67 0.00 0.00 0.00

12 3.67 0.00 0.00 0.00

13 0.00 0.00 3.67 0.00

14 0.00 0.00 0.00 3.67

15 0.00 0.00 3.67 0.00

16 0.00 3.67 0.00 0.00

//

//

NA bcd

Files operate much like lists, so I alway practice with a list, then go to a file:
Here's one way to look at such data in python:

<Moderator NOTE: Merged threads by OP on a single topic (after posting to the second)>

i have to consider the columns also:
like..this is my file
A T G C
01 1.00 2.00 3.00 4.00
02 2.00 3.00 4.00 5.00

now if i say A[01] i should have a value 1.00 or C[01]=4.00
those column names A,T,G,C. i am not able to format my query file properly hope it is understood

Files operate much like lists, so I alway practice with a list, then go to a file:
Here's one way to look at such data in python:

Expand|Select|Wrap|Line Numbers

1. rawdata = \
2. ['01 0.00 3.67 0.00 0.00',
3.  
4. '02 0.00 0.00 3.67 0.00',
5.  
6. '03 0.00 0.00 0.00 3.67',
7.  
8. '04 0.00 3.67 0.00 0.00',
9.  
10. '05 3.67 0.00 0.00 0.00',
11.  
12. '06 3.46 0.00 0.22 0.00',
13.  
14. '07 0.00 0.00 3.67 0.00',
15.  
16. '08 0.00 0.00 0.00 3.67',
17.  
18. '09 0.00 0.00 0.00 3.67',
19.  
20. '10 0.00 3.67 0.00 0.00',
21.  
22. '11 3.67 0.00 0.00 0.00',
23.  
24. '12 3.67 0.00 0.00 0.00',
25.  
26. '13 0.00 0.00 3.67 0.00',
27.  
28. '14 0.00 0.00 0.00 3.67',
29.  
30. '15 0.00 0.00 3.67 0.00']
31.  
32. datadictionary = {} # usually shorten the name to dd
33.  
34. for line in rawdata:
35.     items = line.split()
36.     key = items[0]
37.     datadictionary[key] = [float(item) for item in items[1:]]
38.  
39. print datadictionary['09'][3]

i have to consider the columns also:
like..this is my file
A T G C
01 1.00 2.00 3.00 4.00
02 2.00 3.00 4.00 5.00

now if i say A[01] i should have a value 1.00 or C[01]=4.00
those column names A,T,G,C. i am not able to format my query file properly hope it is understood

I hope you understand that you should be thinking "row 0n, column X, not the other way around. Rows enclose columns, so that is the first index that you deal with.

THis is my file! AS u can see the first column gives the position,the second third ..etc are the part of the matrix.my problem here is that.I have to write a code to access the individual values.like if i say A[01] i should have the value 0.00 or say A[06] i should have a value 3.46. or C[04]=3.67 etc. and i have to add one to each of these elements.can you gimme the code. my code is not working.

NA bap

PO A C G T

01 0.00 3.67 0.00 0.00

02 0.00 0.00 3.67 0.00

03 0.00 0.00 0.00 3.67

04 0.00 3.67 0.00 0.00

05 3.67 0.00 0.00 0.00

06 3.46 0.00 0.22 0.00

07 0.00 0.00 3.67 0.00

08 0.00 0.00 0.00 3.67

09 0.00 0.00 0.00 3.67

10 0.00 3.67 0.00 0.00

11 3.67 0.00 0.00 0.00

12 3.67 0.00 0.00 0.00

13 0.00 0.00 3.67 0.00

14 0.00 0.00 0.00 3.67

15 0.00 0.00 3.67 0.00

16 0.00 3.67 0.00 0.00

//

//

NA bcd

For simplicity, let us assume the data file looks like this:

PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00

Create a dictionary of dictionaries:

hey,
I don really understand what this line does.can u please explain?
values = [[float(s) for s in item] for item in [j[1:] for j in lineList]]
and

for i, item in enumerate(headerList):
dataDict[item] = {}
for key in lineDict:
dataDict[item][key] = lineDict[key]
I dont understand this either.
for this code,I am getting an error which says
value error:invalid literal for float():A

For simplicity, let us assume the data file looks like this:

PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00

Create a dictionary of dictionaries:

Expand|Select|Wrap|Line Numbers

1. fn = r'H:\TEMP\temsys\data7.txt'
2. lineList = [line.strip().split() for line in open(fn).readlines() if line != '\n']
3.  
4. headerList = lineList.pop(0)[1:]
5.  
6. # Key list
7. keys = [i[0] for i in lineList]
8. # Values list
9. values = [[float(s) for s in item] for item in [j[1:] for j in lineList]]
10.  
11. # Create a dictionary from keys and values
12. lineDict = dict(zip(keys, values))
13.  
14. dataDict = {}
15.  
16. for i, item in enumerate(headerList):
17.     dataDict[item] = {}
18.     for key in lineDict:
19.         dataDict[item][key] = lineDict[key][i]

Expand|Select|Wrap|Line Numbers

1. >>> dataDict['A']['05']
2. 3.6699999999999999
3. >>> globals().update(dataDict)
4. >>> A['05']
5. 3.6699999999999999
6. >>> 

hey,
I don really understand what this line does.can u please explain?
values = [[float(s) for s in item] for item in [j[1:] for j in lineList]]
and

for i, item in enumerate(headerList):
dataDict[item] = {}
for key in lineDict:
dataDict[item][key] = lineDict[key]
I dont understand this either.
for this code,I am getting an error which says
value error:invalid literal for float():A

The first question - The code is a list comprehension. Another way to write it would be: The result:

hey,
I don really understand what this line does.can u please explain?
values = [[float(s) for s in item] for item in [j[1:] for j in lineList]]
and

for i, item in enumerate(headerList):
dataDict[item] = {}
for key in lineDict:
dataDict[item][key] = lineDict[key]
I dont understand this either.
for this code,I am getting an error which says
value error:invalid literal for float():A

The second question - This is 'headerList': The values in 'headerList' will be 'keys' in 'dataDict'. 'dataDict' will be the main dictionary, and the values will be subdictionaries. A dictionary key is associated with a value - in this case the value will be another dictionary. Variable 'keys' contain the subdictionary keys: 'lineDict' is a temporary dictionary created to make it easier to compile the data in the form you wanted: Using 'enumerate' on 'headerList', Python gives me these values: Here's an interactive example showing what is happening inside the nested for loop: I hope this helps you understand what is happening.

for this code,I am getting an error which says
value error:invalid literal for float():A

Check variable 'lineList'. It should look like this:

The header list is supposed to take the first four letters right? like A,T,G,C.it is not happening
and my linelist has the entire file(I mean with the A,T,G,C)

Check variable 'lineList'. It should look like this:

Expand|Select|Wrap|Line Numbers

1. >>> for line in lineList:
2. ...     print line
3. ...     
4. ['01', '0.00', '3.67', '0.00', '0.00']
5. ['02', '0.00', '0.00', '3.67', '0.00']
6. ['03', '0.00', '0.00', '0.00', '3.67']
7. ['04', '0.00', '3.67', '0.00', '0.00']
8. ['05', '3.67', '0.00', '0.00', '0.00']
9. ['06', '3.46', '0.00', '0.22', '0.00']
10. ['07', '0.00', '0.00', '3.67', '0.00']
11. ['08', '0.00', '0.00', '0.00', '3.67']
12. ['09', '0.00', '0.00', '0.00', '3.67']
13. ['10', '0.00', '3.67', '0.00', '0.00']
14. ['11', '3.67', '0.00', '0.00', '0.00']
15. ['12', '3.67', '0.00', '0.00', '0.00']
16. ['13', '0.00', '0.00', '3.67', '0.00']
17. ['14', '0.00', '0.00', '0.00', '3.67']
18. ['15', '0.00', '0.00', '3.67', '0.00']
19. ['16', '0.00', '3.67', '0.00', '0.00']
20. >>> 

fn = 'half.txt'
fn_=open("half.txt","r")
file_content=fn_.readlines()

for line in file_content:
linelist=line.strip().split()
print linelist

headerList = linelist.pop(0)[1:]
print headerList

when i do this my code has a headerlist which has
0
1
2
3
4
5
6
7
8
9
0
1
2
3
4
5
6
WHY IS THIS HAppening?
and what does the strip( ) do?

headerList = linelist.pop(0)[1:]
WHAT DOES THIS LINE DO?

Check variable 'lineList'. It should look like this:

Expand|Select|Wrap|Line Numbers

1. >>> for line in lineList:
2. ...     print line
3. ...     
4. ['01', '0.00', '3.67', '0.00', '0.00']
5. ['02', '0.00', '0.00', '3.67', '0.00']
6. ['03', '0.00', '0.00', '0.00', '3.67']
7. ['04', '0.00', '3.67', '0.00', '0.00']
8. ['05', '3.67', '0.00', '0.00', '0.00']
9. ['06', '3.46', '0.00', '0.22', '0.00']
10. ['07', '0.00', '0.00', '3.67', '0.00']
11. ['08', '0.00', '0.00', '0.00', '3.67']
12. ['09', '0.00', '0.00', '0.00', '3.67']
13. ['10', '0.00', '3.67', '0.00', '0.00']
14. ['11', '3.67', '0.00', '0.00', '0.00']
15. ['12', '3.67', '0.00', '0.00', '0.00']
16. ['13', '0.00', '0.00', '3.67', '0.00']
17. ['14', '0.00', '0.00', '0.00', '3.67']
18. ['15', '0.00', '0.00', '3.67', '0.00']
19. ['16', '0.00', '3.67', '0.00', '0.00']
20. >>> 

<snip>what does the strip( ) do?

strip() removes whitespace (or any characters you tell it to) from BOTH ends of a string. Like this:

headerList = linelist.pop(0)[1:]
WHAT DOES THIS LINE DO?

strip() removes whitespace (or any characters you tell it to) from BOTH ends of a string. Like this:

Expand|Select|Wrap|Line Numbers

1. >>> s = "   had_whitespace   "
2. >>> t = s.strip()
3. >>> print repr(t)  # repr() prints the string representation of a varialbe 
4. 'had_whitespace'
5. >>> 

headerList = linelist.pop(0)[1:]
WHAT DOES THIS LINE DO?

Since I didn't write that, I'll just have to do my best to describe it:
Pop(0) takes the zeroth element from linelist (which looks like another list). The result is stored in an unseen, temporary variable. [1:] takes all but the zeroth element from that temporary variable and that (shortened list) is stored in headerList. You can check it out by doing something like
Hope that helps.

thanks for that . but could you give ur mail id..because I am not able to paste my exact file here and am stuck with this problem.plz help!

Since I didn't write that, I'll just have to do my best to describe it:
Pop(0) takes the zeroth element from linelist (which looks like another list). The result is stored in an unseen, temporary variable. [1:] takes all but the zeroth element from that temporary variable and that (shortened list) is stored in headerList. You can check it out by doing something like

Expand|Select|Wrap|Line Numbers

1. temp = linelist.pop(0)
2. print temp
3. headerList = temp[1:]   # this is called a "slice" from the temp list.
4. print headerList

Hope that helps.

thanks for that . but could you give ur mail id..because I am not able to paste my exact file here and am stuck with this problem.plz help!

Paste what you can here and I'll format it for you.

Paste what you can here and I'll format it for you.

NA bap
PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00
//
//
NA bcd
PO A C G T
01 42.55 8.75 145.86 8.14
02 0.14 0.53 204.64 0.00
03 126.83 78.02 0.11 0.34
04 0.21 0.17 0.00 204.92
05 0.00 12.38 0.43 192.50
06 174.48 0.95 1.32 28.56
07 79.53 4.70 100.44 20.64
//
//
NA bin
PO A C G T
01 0.45 8.27 0.00 11.39
02 0.00 0.00 10.02 10.09
03 5.80 1.39 0.00 12.93
04 12.33 5.18 2.60 0.00
05 12.43 0.00 0.00 7.68
06 18.55 0.00 1.57 0.00
07 0.05 0.58 0.00 19.48
08 20.11 0.00 0.00 0.00
09 20.06 0.05 0.00 0.00
10 20.11 0.00 0.00 0.00
11 0.00 15.33 0.00 4.78
12 20.06 0.05 0.00 0.00
13 14.99 0.35 4.78 0.00
14 13.64 2.42 3.37 0.68
15 5.03 0.00 15.08 0.00
16 7.23 0.45 10.94 1.49
//
//
I want to add one to each and every element of the file(of course except the first column giving the position)
and want to access the elements like A[01]=,c[02]=,etc..

NA bap
PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00
//
//
NA bcd
PO A C G T
01 42.55 8.75 145.86 8.14
02 0.14 0.53 204.64 0.00
03 126.83 78.02 0.11 0.34
04 0.21 0.17 0.00 204.92
05 0.00 12.38 0.43 192.50
06 174.48 0.95 1.32 28.56
07 79.53 4.70 100.44 20.64
//
//
NA bin
PO A C G T
01 0.45 8.27 0.00 11.39
02 0.00 0.00 10.02 10.09
03 5.80 1.39 0.00 12.93
04 12.33 5.18 2.60 0.00
05 12.43 0.00 0.00 7.68
06 18.55 0.00 1.57 0.00
07 0.05 0.58 0.00 19.48
08 20.11 0.00 0.00 0.00
09 20.06 0.05 0.00 0.00
10 20.11 0.00 0.00 0.00
11 0.00 15.33 0.00 4.78
12 20.06 0.05 0.00 0.00
13 14.99 0.35 4.78 0.00
14 13.64 2.42 3.37 0.68
15 5.03 0.00 15.08 0.00
16 7.23 0.45 10.94 1.49
//
//
I want to add one to each and every element of the file(of course except the first column giving the position)
and want to access the elements like A[01]=,c[02]=,etc..

Lacking your working program pasted here, I'm afraid that the best I can tell you is it should be as simple as: or even in a loop through your data structure.

when i execute i get an error that says
keyerror: '01'

<removed quote=bvdet to make this post appear.>

'headerList': The values in 'headerList' will be 'keys' in 'dataDict'. 'dataDict' will be the main dictionary, and the values will be subdictionaries. A dictionary key is associated with a value - in this case the value will be another dictionary. Variable 'keys' contain the subdictionary keys: 'lineDict' is a temporary dictionary created to make it easier to compile the data in the form you wanted: Using 'enumerate' on 'headerList', Python gives me these values: Here's an interactive example showing what is happening inside the nested for loop: I hope this helps you understand what is happening.

Expand|Select|Wrap|Line Numbers

1. file_=open("half1.txt","r")
2. file_content=file_.readlines()
3. linelist=[line.strip().split() for line in file_content if line != '\n']
4.  
5. headerlist=linelist.pop(0)[1:]
6.  
7. keys=[i[0] for i in linelist]
8.  
9.  
10. values=[[float(s) for s in item] for item in [j[1:] for j in linelist]]
11.  
12. linedict =dict(zip(keys,values))
13.  
14. datadict={}
15. for i,item in enumerate(headerlist):
16.     datadict[item]={}
17.     for key in linedict:
18.         print key
19.         datadict[item][key]=linedict[key][i]
20.         print datadict['A']['01']
21.  

when i execute i get an error that says
keyerror: '01'

this is my program. i am getting an error that says
key error:'01' and i do want to add one to all the elements
(i mean the float values)
can you see the error?

Expand|Select|Wrap|Line Numbers

1. file_=open("half1.txt","r")
2. file_content=file_.readlines()
3. linelist=[line.strip().split() for line in file_content if line != '\n']
4.  
5. headerlist=linelist.pop(0)[1:]
6. print headerlist
7.  
8. keys=[i[0] for i in linelist]
9.  
10.  
11. values=[[float(s) for s in item] for item in [j[1:] for j in linelist]]
12.  
13. linedict =dict(zip(keys,values))
14.  
15. datadict={}
16. for i,item in enumerate(headerlist):
17.     datadict[item]={}
18.     for key in linedict:
19.         print key
20.         datadict[item][key]=linedict[key][i]
21.         print datadict['A']['01']

this is my program. i am getting an error that says
key error:'01' and i do want to add one to all the elements
(i mean the float values)
can you see the error?

I've answered above.

You must learn to use code tags. Instructions are on the right hand side of the page when you are making your post (or reply). There is also very much helpful information in the How to ask a question section of our Posting Guidelines. Thanks

Lacking your working program pasted here, I'm afraid that the best I can tell you is it should be as simple as:

Expand|Select|Wrap|Line Numbers

1. temp = A[01]
2. print temp
3. A[01] = temp + 1
4. print A[01]

or even

Expand|Select|Wrap|Line Numbers

1. print A[01]
2. A[01] += 1
3. print A[01]

in a loop through your data structure.

It looks like you have three data sets in the file. Which one do you want to work on?

Notice that I simplified my example code by eliminating some of the lines in the data file. It appears that you have not adjusted the code to account for that.

I have modified the first part of my example code to read the first data set in the OP file: This will add a set amount to every element in the data:

so here is the code that you suggested for creating dictionaries for a file(matrix)
now i have gotta do something like..
supposing the sequence entered is something like "ATATTA".. so A is in the first position so A[01]*T[02]*A[03]*T[04]*T[05]*A[06]=??
the sequence is going to be entered by the user everytime(so it will keep changing)
how do i do this?? what changes should i do??hope I am clear!!
THIS is the file containing the matrix
NA bap
PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00
//
//
NA bcd
PO A C G T
01 42.55 8.75 145.86 8.14
02 0.14 0.53 204.64 0.00
03 126.83 78.02 0.11 0.34
04 0.21 0.17 0.00 204.92
05 0.00 12.38 0.43 192.50
06 174.48 0.95 1.32 28.56
07 79.53 4.70 100.44 20.64
//
//
NA bin
PO A C G T
01 0.45 8.27 0.00 11.39
02 0.00 0.00 10.02 10.09
03 5.80 1.39 0.00 12.93
04 12.33 5.18 2.60 0.00
05 12.43 0.00 0.00 7.68
06 18.55 0.00 1.57 0.00
07 0.05 0.58 0.00 19.48
08 20.11 0.00 0.00 0.00
09 20.06 0.05 0.00 0.00
10 20.11 0.00 0.00 0.00
11 0.00 15.33 0.00 4.78
12 20.06 0.05 0.00 0.00
13 14.99 0.35 4.78 0.00
14 13.64 2.42 3.37 0.68
15 5.03 0.00 15.08 0.00
16 7.23 0.45 10.94 1.49
//
//

so here is the code that you suggested for creating dictionaries for a file(matrix)
now i have gotta do something like..
supposing the sequence entered is something like "ATATTA".. so A is in the first position so A[01]*T[02]*A[03]*T[04]*T[05]*A[06]=??
how do i do this?? what changes should i do??
THIS is the file containing the matrix
NA bap
PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00
//
//
NA bcd
PO A C G T
01 42.55 8.75 145.86 8.14
02 0.14 0.53 204.64 0.00
03 126.83 78.02 0.11 0.34
04 0.21 0.17 0.00 204.92
05 0.00 12.38 0.43 192.50
06 174.48 0.95 1.32 28.56
07 79.53 4.70 100.44 20.64
//
//
NA bin
PO A C G T
01 0.45 8.27 0.00 11.39
02 0.00 0.00 10.02 10.09
03 5.80 1.39 0.00 12.93
04 12.33 5.18 2.60 0.00
05 12.43 0.00 0.00 7.68
06 18.55 0.00 1.57 0.00
07 0.05 0.58 0.00 19.48
08 20.11 0.00 0.00 0.00
09 20.06 0.05 0.00 0.00
10 20.11 0.00 0.00 0.00
11 0.00 15.33 0.00 4.78
12 20.06 0.05 0.00 0.00
13 14.99 0.35 4.78 0.00
14 13.64 2.42 3.37 0.68
15 5.03 0.00 15.08 0.00
16 7.23 0.45 10.94 1.49
//
//

<Mod NOTE: merged duplicate threads - again>

HERE IS MY code!!plz help how to proceed??

Thanks for using code tags.
First thing to do is try it. It looks like it should work.

Right after that, READ the rest of the POSTING GUIDELINES. You have committed several infractions and must consider yourself warned.

I am sorry about that. I will follow the rules hence forth.but i dont know how to proceed with my problem further.

I am sorry about that. I will follow the rules hence forth.but i dont know how to proceed with my problem further.

Ok. All is forgiven. Here's the thing. I can't give you working code.
I can give you examples of what I see going on here so that you can try things on your own. It looks to me like you could multiply the specified elements (after your matrix has been created) like this:

Thanks a lot for all the help! I still have a small doubt.what does this line do?
I mean this is supposed to refer to the line which doesnt start with "PO' right? when i say print line. It is printing
PO A C G T
but this is not what it means right?

Thanks a lot for all the help! I still have a small doubt.what does this line do?

Expand|Select|Wrap|Line Numbers

1.  
2. while not  line.startswith('PO'):
3.  
4.  
5.     line=f.next()
6.     print line
7.  

I mean this is supposed to refer to the line which doesnt start with "PO' right? when i say print line. It is printing
PO A C G T
but this is not what it means right?

f.next() merely advances one line in the file. The object is to advance to the first line that starts with "PO".

I asked you a question in an earlier post about your data. It looks like there are three data sets in your data file. The example code I posted only parses the first data set.

okay,then why do we say,
while not?? it should be while..right??
and yes there are actually a lot of data sets in my file around 15 to twenty first am trying to make it work for one.then i will do it for the rest(sorry for not answering I am really going crazy with the program)i owe u an apology
And i have one more doubt:
[code]
datadict={}
for i,item in enumerate(headerlist):
datadict[item]={}
for key in linedict:
print key
# when i say a print key here, it is printing the keys twice and in an unorderly #fashion.why is this happening?
datadict[item][key]=linedict[key][i]
for keymain in datadict:
for keysub in datadict[keymain]:
datadict[keymain][keysub]+=1.0
print datadict
looking fwd for ur reply!
cheers!

okay,then why do we say,
while not?? it should be while..right??
and yes there are actually a lot of data sets in my file around 15 to twenty first am trying to make it work for one.then i will do it for the rest(sorry for not answering I am really going crazy with the program)i owe u an apology
And i have one more doubt:
[code]
datadict={}
for i,item in enumerate(headerlist):
datadict[item]={}
for key in linedict:
print key
# when i say a print key here, it is printing the keys twice and in an unorderly #fashion.why is this happening?
datadict[item][key]=linedict[key][i]
for keymain in datadict:
for keysub in datadict[keymain]:
datadict[keymain][keysub]+=1.0
print datadict
looking fwd for ur reply!
cheers!

I should be 'while not'. We want to skip the lines until the line starts with 'PO'. You can use this same method to advance into later data sets.

No problem about not answering.

Dictionaries are unordered collections of data. You can print in an orderly fashion like this:

my file:
my file:
NA bap
PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00
//
//



This is my code:
I have written the problems in comment form in the code!

can you execute my code and see why this is happening? I have been stuck with this for a day now.
waiting for your reply
cheers!

my file:
my file:
NA bap
PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00
//
//



This is my code:

Expand|Select|Wrap|Line Numbers

1. f=open("deeps1.txt","r")
2. line=f.next()
3. while not line.startswith('PO'):
4.     line=f.next()
5.  
6. headerlist=line.strip().split()[1:]
7. linelist=[]
8.  
9.  
10. line=f.next().strip()
11. while not line.startswith('/'):
12.     if line != '':
13.         linelist.append(line.strip().split())
14.     line=f.next().strip()
15.  
16. keys=[i[0] for i in linelist]
17. values=[[float(s) for s in item] for item in [j[1:] for j in linelist]]
18. array={}
19. linedict=dict(zip(keys,values))
20. keys = linedict.keys()
21. keys.sort()
22. for key in keys:
23.     array=[key,linedict[key]]
24.  
25. datadict={}
26. datadict1={}
27. for i,item in enumerate(headerlist):
28.     datadict[item]={}
29.     print datadict[item]#this one returns empty dictionary
30.     for key_ in linedict:
31.         print item# all items are getting printed
32.         datadict[item][key_]=linedict[key_][i]
33.         # but for the print statement below its saying key error:'C'
34.        print datadict['C']['01']
35.  

I have written the problems in comment form in the code!

can you execute my code and see why this is happening? I have been stuck with this for a day now.
waiting for your reply
cheers!

For that to work, you'll have to print every new addition to the datadict.

my file:
my file:
NA bap
PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00
//
//



This is my code:

Expand|Select|Wrap|Line Numbers

1. f=open("deeps1.txt","r")
2. line=f.next()
3. while not line.startswith('PO'):
4.     line=f.next()
5.  
6. headerlist=line.strip().split()[1:]
7. linelist=[]
8.  
9.  
10. line=f.next().strip()
11. while not line.startswith('/'):
12.     if line != '':
13.         linelist.append(line.strip().split())
14.     line=f.next().strip()
15.  
16. keys=[i[0] for i in linelist]
17. values=[[float(s) for s in item] for item in [j[1:] for j in linelist]]
18. array={}
19. linedict=dict(zip(keys,values))
20. keys = linedict.keys()
21. keys.sort()
22. for key in keys:
23.     array=[key,linedict[key]]
24.  
25. datadict={}
26. datadict1={}
27. for i,item in enumerate(headerlist):
28.     datadict[item]={}
29.     print datadict[item]#this one returns empty dictionary
30.     for key_ in linedict:
31.         print item# all items are getting printed
32.         datadict[item][key_]=linedict[key_][i]
33.         # but for the print statement below its saying key error:'C'
34.        print datadict['C']['01']
35.  

I have written the problems in comment form in the code!

can you execute my code and see why this is happening? I have been stuck with this for a day now.
waiting for your reply
cheers!

I made a few minor changes to your code:

sorry,
that doesnt seem to do anything.
:( or should i add it somewhere else?
waiting for ur reply
cheers!

CODE TAGS! This code: give me this output:

thank you sir!that works for me too..but what exactly i want when i type
datadict['C']["01'] is 3.67(according to my file that is printed below).But it is printing the entire thing
this is what my file says:
NA bap
PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00
//
//
so how do I change the code accordingly?
waiting,
cheers!

thank you sir!that works for me too..but what exactly i want when i type
datadict['C']["01'] is 3.67(according to my file that is printed below).But it is printing the entire thing
this is what my file says:
NA bap
PO A C G T
01 0.00 3.67 0.00 0.00
02 0.00 0.00 3.67 0.00
03 0.00 0.00 0.00 3.67
04 0.00 3.67 0.00 0.00
05 3.67 0.00 0.00 0.00
06 3.46 0.00 0.22 0.00
07 0.00 0.00 3.67 0.00
08 0.00 0.00 0.00 3.67
09 0.00 0.00 0.00 3.67
10 0.00 3.67 0.00 0.00
11 3.67 0.00 0.00 0.00
12 3.67 0.00 0.00 0.00
13 0.00 0.00 3.67 0.00
14 0.00 0.00 0.00 3.67
15 0.00 0.00 3.67 0.00
16 0.00 3.67 0.00 0.00
//
//
so how do I change the code accordingly?

Expand|Select|Wrap|Line Numbers

1. f=open("deeps1.txt","r")
2. line=f.next()
3. while not line.startswith('PO'):
4.     line=f.next()
5.  
6. headerlist=line.strip().split()[1:]
7. linelist=[]
8.  
9.  
10. line=f.next().strip()
11. while not line.startswith('/'):
12.     if line != '':
13.         linelist.append(line.strip().split())
14.     line=f.next().strip()
15.  
16. keys=[i[0] for i in linelist]
17. values=[[float(s) for s in item] for item in [j[1:] for j in linelist]]
18. array={}
19. linedict=dict(zip(keys,values))
20. keys = linedict.keys()
21. keys.sort()
22. for key in keys:
23.     array=[key,linedict[key]]
24.  
25. datadict={}
26. datadict1={}
27. for i,item in enumerate(headerlist):
28.     datadict[item]={}
29.     print item
30.     for key_ in linedict:
31.         datadict[item][key_]=linedict[key_][i]
32.         print linedict[key_][i] 
33. print datadict['C']['01']
34.  

waiting,
cheers!

I tested your code exactly as you posted it: There are other print statements in the code.

Thanks for all the help,I wouldnt have come this far without the help of all u ppl.
there is a new problem
this is the code.as I already posted i want to find something like
A[01]*T[02]*A[03]
But the problem I am facing is that the datadict has keys like "01","02" but the in the loop of seq i have 1,2,3,4. and i cant start from zero whatsoever. how can i make the looping of my seq to 01,02 etc..if i say
for i in range('01',len(seq)):
its taking it as a string!
waiting for ur reply,
cheers!!

oops am sorry!
The one you suggested is working:
this is the code:
This seems to do for the first letter only "A"
waiting for your reply,
cheers!

sorry!!
its working:)
cheers!!

I am doing a simple task here,and I am getting error.cant understand why that is happening!
I am trying to find out a score for a sequence after creating all the dictionaries:
if the seq="ACGT"
value of A and T is 0.3
VALUE of C and G is 0.2
score=val(A)*val(c)*val(G)*val(T)
so it should be score=0.3*0.2*0.2*0.3=3.6
My error is mentioned in comment form in the last lines of the code
this is my code:
waiting 4 ur reply
cheers!

I am doing a simple task here,and I am getting error.cant understand why that is happening!
I am trying to find out a score for a sequence after creating all the dictionaries:
if the seq="ACGT"
value of A and T is 0.3
VALUE of C and G is 0.2
score=val(A)*val(c)*val(G)*val(T)
so it should be score=0.3*0.2*0.2*0.3=3.6
My error is mentioned in comment form in the last lines of the code
this is my code:

Expand|Select|Wrap|Line Numbers

1. f=open("deeps1.txt","r")
2. line=f.next()
3. while not line.startswith('PO'):
4.     line=f.next()
5.  
6. headerlist=line.strip().split()[1:]
7. linelist=[]
8.  
9.  
10. line=f.next().strip()
11. while not line.startswith('/'):
12.     if line != '':
13.         linelist.append(line.strip().split())
14.     line=f.next().strip()
15.  
16. keys=[i[0] for i in linelist]
17. values=[[float(s) for s in item] for item in [j[1:] for j in linelist]]
18. array={}
19. linedict=dict(zip(keys,values))
20. keys = linedict.keys()
21. keys.sort()
22. for key in keys:
23.     array=[key,linedict[key]]
24.  
25. datadict={}
26. datadict1={}
27. for i,item in enumerate(headerlist):
28.     datadict[item]={}
29.     for key_ in linedict:
30.         datadict[item][key_]=linedict[key_][i]
31.  
32. for keymain in datadict:
33.     for keysub in datadict[keymain]:
34.         datadict[keymain][keysub]+=1.0
35. #print datadict['T']['16']
36. seq="CGTCAG"
37.  
38. res=1
39. for i in range(0,len(seq)):
40.     key=seq[i]
41.     res*=datadict[key]["%02d"%(i+1)]
42.     #print res
43.     score=1
44.     value={"A":"0.3","T":"0.3","C":"0.2","G":"0.2"}
45.     for it in value:
46.         for item in seq:
47.             if it==key:
48.                 score=score*value[it]
49.                 print score# I get an error that says TypeError: can't multiply sequence by non-int of type 'str'
50.  
51.  
52. print res
53.  

waiting 4 ur reply
cheers!

You are receiving the error because the values in dictionary 'value' are strings. Either define them as numbers (e.g. "A":0.3,"T":0.3,....) or convert to float in the calculation: