來一份Python學習題

1、3*2*2 的輸出是多少？ (1分)8 % 4的輸出是多少？ (1分)32 + 32的輸出是什么? (1分)32 32的輸出是什么？(1分)Sheng Xin Bao Dian. find(x)和Sheng Xin BaoDian . find (X)的輸出分別是? (2分)一句話計算Sheng Xin Bao Dian字符串中n的數目？ (1分)寫出下面10段程序的輸出？(1分/段)aList = 1, 2, 3 bList = aList bList. append(4) aListaList = 1, 2, 3 cList = aList: cList. append(4) aLista

2、List = 1, 1, 2, 2, 3, 5, 4, 3 aSet = set (aList) aSet1, 2, 3 * 2i*2 for i in 1,2,3diet ( (i, i*2) for i in range (5)import re re. findall (，z I imageGPz/, ).joinvL Sheng , Xin , Bao , Dian)def sumNumber (a, b):return a + b sumNumber(2, 3)def sumNumber (a, b):return a + b print(sumNumber(2, 3)3.程序運行沒

3、有錯誤不代表你寫的程序完成了你的需 求，你要去查驗輸出結果是不是你想要的。2 5.關于程序調試0.在初寫程序時，可能會出現各種各樣的錯誤，常見的 有縮進不一致，變量名字拼寫錯誤，喪失冒號，文件 名未加引號等，這時要根據錯誤提示查看錯誤類型是 什么，出錯的是哪一行來定位錯誤。當然，有的時候 報錯的行自身不一定有錯，可能是其前面或后面的行 出現了錯誤。.用腦袋運行程序：當程序寫作完成后，自己嘗試對著 數據文件，一行一行的執(zhí)行程序，來看程序的運行是 否與自己想干的活一致，有沒有維漏。.當結果不符合預期時，要學會使用print來杳看每步 的操作是否正確，比方我讀入了字典，我就打印下字 典，看看讀入的是

4、不是我想要的，是否含有不該存在 的字符；或者在每個判斷句、函數調入的情況下打印 個字符，來跟蹤程序的運行軌跡.寫程序以下面列表中每個元素為key ,元素出現的次數為value , 構建一個字典，并遍歷字典按元素的ASCII碼順序輸出？ (5分) aList = a, b, c, a, d/e； A.對教案中腦筋急轉彎問題的解法進行優(yōu)化；問題是：現有100元 錢，需要買100個物品，其中鉛筆盒單價5元，筆單價3元，橡皮 單價0.5元，怎么組合可以把100元花完，同時三種物品的個數和 為100,請用編程解決。(3分).寫程序用高斯的計算方式計算1+2+3+.+100的加和。(3分).指出下面每個程

5、序運行時可能會出現的錯誤。(1分/段)aList = 1, 2,3 aDict = aDictaList = 1 b = aDictf a，if 1: print(z，Sheng xin bao dian great!”)32 + 32aList = 1, 2, 3 aList. add(4)aList = 1, 2, 3join(aList)int ( a)3/0for i in range(10) print (L)Sheng Xin * 3type 二 1. Python文件讀寫函數open的mode參數中r, w, a, t, b, x分別是什 么意思? （3分）. Python中如何

6、獲取當前所在的工作目錄？如何修改工作目錄？ （3 分）. Python中連接多個字符串的方法有哪些?優(yōu)缺點是什么? （3分）. print （-% （1/3）的輸出是什么? （2 分）.描述下語句import pandas as pd做了什么操作？ （2分）.教案中基因ENSEMBLE ID轉Gene Symbol程序用pandas實現 （GRCh38. idmap,ensm. id）o （5 分）. Jupyter中溉writefile, %run宏命令的用途是什么？ （2分）.找出 TP53 mRNA 序列中的 ORF （human_TP53_mRNA. fa）o （5 分）.列出大腸桿菌

7、基因組中限制性內切酶See!的切割位置（Ecoli. fa）0 （5 分）.計算data/testl. fa中每條序列的GC含量。（5分）.不使用pandas ,寫Python腳本處理Pandas教案中的TPM表 達矩陣的提取和合并？（ENCFF060LPA.tsv, ENCFF262OBL.tsv, ENCFF289HGQ.tsv, ENCFF673KYR.tsv） （8 分）.給定FASTA格式的文件（testl.fa和test2.fa）,寫一個程序cat. py讀入文件,并輸出到屏幕（2分）open（file）for. in loopprint()strip() function用到的知

8、識點13.給定FASTQ格式的文件(testl.fq),寫一個程序cat. py讀入文件，并輸出到屏幕(2分)0.同上.用到的知識點14.寫程序splitName.py,讀入test2.fa,并取原始序列名字第一個空 格前的名字為處理后的序列名字，輸出到屏幕(2分)0. split.字符串的索引.用到的知識點.輸出格式為：NM_001011874gcggcggcgggcgagcgggcgctggagtaggagctg.15.寫程序formatFasta. py,讀入test2.fa ,把每條FASTA序列連成 一行然后輸出(2分)0. join. strip.用到的知識點.輸出格式為：NM_0

9、01011874gcggcggcgggcTCCGCTGGCGTTCACC.CGGGGTCCGGAG16 .寫程序 formatFasta-2. py, 讀入test2.fa，把每條FASTA序列分 割成80個字母一行的序列（2分）0.字符串切片操作. range.用到的知識點.輸出格式為 NM_001011874 gcggcggcgc.（60 字母）.TCCGCTGACG #（每行80個字母） acgtgctacg.（60 個字母）.GCGTTCACCCACGTACGATG （最后一行可缺乏80個字母）17.寫程序sortFasta. py,讀入test2.fa,并取原始序列名字第一個空 格前

10、的名字為處理后的序列名字，排序后輸出（2分）0. sortdietaDictkey=aDictkey.append（value）用到的知識點.提取給定名字的序列（2分）0.用到的知識點. print fh, or fh.write（）.取模運算，4%2=0.寫程序grepFasta. py,提取中名字對應的test2.fa的序列，并輸出到屏幕。.寫程序grepFastq. pyz提取中名字對應 的testl.fq的序列，并輸出到文件。.寫程序 screenResult. py, 篩選 test.expr 中 foldChange 大于 2 的基因并且pa

11、dj小于0.05的基，可以輸出整行或只輸出基因名 字。（4分）0.邏輯與操作符and.文件中讀取的內容都為字符串，需要用int轉換為整 數，float轉換為浮點數.用到的知識點20.寫程序 transferMultipleColumToMatrix. py 將文件 （multipleColExpr.txt）中基因在多個組織中的表達數據轉換為矩陣 形式，并繪制熱圖。（6分）0. aDict key = . aDictkey key2 = value. if key not in aDict. aDict = ,ENSG00000000003, : ,A-431Z/:21.3, A-549” ,E

12、NSG00000000003M :,.用到的知識點.輸入格式（只需要前3列就可以）Gene Sample ValueUnit Abundance ENSG00000000003 A-431 21.3FPKM Medium ENSG00000000003 A-549 32.5FPKM Medium ENSG00000000003 AN3-CA 38.2FPKM Medium ENSG00000000003 BE WO 31.4FPKM Medium ENSGOOOOOOOOOO3 CACO-2 63.9FPKM High ENSGOOOOOOOOOO5 A-431 0.0 FPKMNot det

13、ected ENSG00000000005 A-549 0.0FPKM Not detected ENSG00000000005 AN3-CA0.0 FPKM Not detected ENSG00000000005BEWO 0.0 FPKM Not detectedENSG00000000005 CACO-2 0.0 FPKM Not detected.輸出格式 Name A-431 A-549 AN3-CA BEWOCACO-2 ENSG00000000460 25.2 14.2 10.6 24.414.2 ENSG00000000938 0.0 0.0 0.0 0.0 0.0ENSG00

14、000001084 19.1 155.1 24.4 12.6 23.5ENSG00000000457 2.8 3.4 3.8 5.8 2.921.寫程序 reverseComplementary. py 計算序列 ACGTACGTACGTCACGTCAGCTAGAC 的反向互補序列。（2 分）0. reverse. list（seq）.用到的知識點22.寫程序 collapsemiRNAreads. py 轉換 smRNA-Seq 的測序數據。（5分）0.輸入文件格式（mir.collapse, tab-分割的兩列文件，第 一列為序列，第二列為序列被測到的次數）ID_REFVALUE ACTG

15、CCCTAAGTGCTCCTTCTGGC 2 ATAAGGTGCATCTAGTGCAGATA 25 TGAGGTAGTAGTTTGTGCTGTTT 100 TCCTACGAGTTGCATGGATTC 41.輸出文件格式（mir.collapse.fa,名字的前3個字母為 樣品的特異標示，中間的數字表示第幾條序列，是序 列名字的唯一標示，第三局部是x加每個reads被測 到的次數。三局部用下劃線連起來作為fasta序列的 名字。）ESB_l_x2ACTGCCCTAAGTGCTCCTTCTGGC ESB_2_x25 ATAAGGTGCATCTAGTGCAGATA ESB_3_xl00 TGAGGT

16、AGTAGTTTGTGCTGTTT ESB_4_x4 TCCTACGAGTTGCATGGATTC23.簡化的短序列匹配程序（map.py）把short.fa中的序列比對到 ref.fa,輸出短序列匹配到ref.fa文件中哪些序列的哪些位置。（10 分）0. find.用到的知識點.輸出格式（輸出格式為bed格式,第一列為匹配到的染色體，第二列和第三列為匹配到染色體序列的起始終止位置（位置標記以0為起始，代表第一個位置；終止位置不包含在內，第一個例子中所示序列的位置 是（199,208（前閉后開，實際是chrl染色體第199- 206的序列，。起始）.第4列為短序列自身的序列.）。.附加要求：可以只匹配到給定的模板鏈，也可以考慮 匹配到模板鏈的互補鏈。這時第5列可以為短序列的 名字，第六列為鏈的信息，匹配到模板鏈為+, 匹配到互補鏈為-注意匹配到互補鏈時起始位置24.備注:也是從模板鏈的5端算起的。chrl 199 208 TGGCGTTCA chrl 207