Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
An Introduction to Python – Part III Dr. Nancy Warter-Perez Overview 2-D Lists List comprehensions Zip File I/O Split Functions Programming Workshop #3 Introduction to Python – Part III 2 Python List Comprehensions Precise way to create a list Consists of an expression followed by a for clause, then zero or more for or if clauses Ex: >>> [str(round(355/113.0, i)) for i in range(1,6)] ['3.1', '3.14', '3.142', '3.1416', '3.14159'] Ex: replace all occurrences of G or C in a string of amino acids with a 1 and A and T with a 0 >>> x = "acactgacct" >>> y = [int(i=='c' or i=='g') for i in x] >>> y [0, 1, 0, 1, 0, 1, 0, 1, 1, 0] Introduction to Python – Part III 3 Creating 2-D Lists To create a 2-D list L, with C columns and R rows initialized to 0: L = [[]] #empty 2-Dlist L = [[0 for col in range(C)] for row in range(R)] To assign the value 5 to the element at the 2nd row and 3rd column of L L[2][3] = 5 Introduction to Python – Part III 4 Zip – for parallel traversals Visit multiple sequences in parallel Ex: >>> L1 = [1,2,3] >>> L2 = [5,6,7] >>> zip(L1, L2) [(1,5), (2,6), (3,7)] Ex: >>> for(x,y) in zip(L1, L2): … print x, y, '--', x+y 1 5 -- 6 2 6 -- 8 3 7 -- 10 Introduction to Python – Part III 5 More on Zip Zip more than two arguments and any type of sequence Ex: >>> T1, T2, T3 = (1,2,3),(4,5,6),(7,8) >>> T3 (7,8) >>> zip(T1, T2, T3) [(1,4,7),(2,5,8)] -- truncates to shortest sequence Introduction to Python – Part III 6 Dictionary Construction with zip Ex: >>> keys = ['a', 'b', 'd'] >>> vals = [1.8, 2.5, -3.5] >>> hydro = dict(zip(keys,vals)) >>> hydro {'a': 1.8, 'b': 2.5, 'd': -3.5} Introduction to Python – Part III 7 File I/O To open a file myfile = open('pathname', <mode>) modes: 'r' = read 'w' = write Ex: infile = open("D:\\Docs\\test.txt", 'r') Ex: outfile = open("out.txt", 'w') – in same directory Introduction to Python – Part III 8 Common input file operations Operation input = open ('file', 'r') S = input.read() S = input.read(N) S = input.readline() L = input.readlines() Interpretation open input file read entire file into string S Read N bytes (N>= 1) Read next line Read entire file into list of line strings Introduction to Python – Part III 9 Common output file operations Operation Interpretation output = open('file', 'w') create output file output.write(S) Write string S into file output.writelines(L) Write all line strings in list L into file output.close() Manual close (good habit) Introduction to Python – Part III 10 Extracting data from string – split String.split([sep, [maxsplit]]) - Return a list of the words of the string s. If the optional argument sep is absent or None, the words are separated by arbitrary strings of whitespace characters (space, tab, newline, return, formfeed). If the argument sep is present and not None, it specifies a string to be used as the word separator. The optional argument maxsplit defaults to 0. If it is nonzero, at most maxsplit number of splits occur, and the remainder of the string is returned as the final element of the list (thus, the list will have at most maxsplit+1 elements). Introduction to Python – Part III 11 Split Ex: >>> x = "a,b,c,d" >>> x.split(',')['a', 'b', 'c', 'd'] >>> x.split(',',2)['a', 'b', 'c,d'] Ex: >>> y = "5 33 a 4" >>> y.split()['5', '33', 'a', '4'] Introduction to Python – Part III 12 Functions Function definition def adder(a, b, c): return a+b+c Function calls adder(1, 2, 3) -> 6 Introduction to Python – Part III 13 Functions – Polymorphism >>>def fn2(c): … a=c*3 … return a >>> print fn2(5) 15 >>> print fn2(1.5) 4.5 >>> print fn2([1,2,3]) [1,2,3,1,2,3,1,2,3] >>> print fn2("Hi") HiHiHi Introduction to Python – Part III 14 Functions - Recursion def fn_Rec(x): if x == []: return fn_Rec(x[1:]) print x[0], y = [1,2,3,4] fn_Rec(y) >>> 4 3 2 1 Introduction to Python – Part III 15 Programming Workshop #3 Create a text file called "test1.txt" with the following data: # Sample data 12345 # more data 6 7 8 9 10 Create another text filed called "test2.txt" with the following data: # More test data # With more header info ABCDEFG Write a script to do the following 1. Prompt the user for a filename 2. Open the file 3. Read the file into a list of strings. 4. If the line does not begin with a '#' print the line to the screen. Test your script on test1.txt and test2.txt. Introduction to Python – Part III 16 Programming Homework #2P Write a program to prompt the user for a scoring matrix file name and read the data into a dictionary Download a representative set of PAM and Blossum Scoring Matrix Files Scoring matrices should be downloaded from ftp://ftp.ncbi.nih.gov/blast/matrices/ Due Date: Thursday, November 8th (must submit to get an extension) One Week Automatic Extension: Thursday, November 15th Introduction to Python – Part III 17 Example Scoring Matrix File Introduction to Python – Part III 18 Algorithm for Homework 2P Step 1: Create an empty list (of dictionaries) Step 2: Prompt the user for the scoring matrix file name Step 3: Open the file and read the contents as a list of strings. Ignore the comment lines Step 4: When you reach a line that doesn’t start with '#' read in the amino acid symbols and split them into your keys for your dictionary Step 5: Read in the rest of the lines one at a time. For each line: Step 5a. Slice off the first character (amino acid). Step 5b. For the rest of the string split into individual numbers and convert to a list of integers (use a list comprehension). This is your data for your dictionary. Step 5c. Zip the keys and data together and convert into a dictionary. Step 5d. Add the dictionary to the list of dictionaries Step 6: After you’ve read all lines, create the dictionary of dictionaries by zipping the keys and the list of dictionaries and convert into a dictionary. Introduction to Python – Part III 19