Python学习之路和隐藏特征

针对生信领域的零基础爱好者及生信分析中遇到的种种问题,生信领域知名公众号“生信宝典”团队组织了中科院系统项目经验丰富的一线科研人员开展系列培训活动。

小学生都学Python了,你还不知道怎么开始文中介绍了Python的应用广泛,功能强大,提供了Python的在线学习视频和资料等。

学习程序语言不是一件难事,也不是一件简单事。为什么编程这么难中翻译了一篇编程学习的心路历程。

(图例“编程信心与能力”:纵轴为信心值,横轴为能力水平,虚线从左至右依次分割出手牵手蜜月期、混沌悬崖、绝望沙漠、令人兴奋的上升期四个阶段,第5条虚线标志着工作准备就绪)

在我看来,学习编程是学以致用,学习方式是硬着头皮去写。

读文档是蜜月期,读读就过去,谁都会。

写程序是混沌悬崖,知道是这么回事,就是写不出来;

调程序是绝望沙漠,怎么看自己写的都对,就是编译器不开眼;

程序正确了就是兴奋期,万里长征迈出又一步。

这些过程没有捷径,但可以加速,就是我们的零基础Python编程班,应用Python处理生物信息数据和作图

广告做完了,来点Python的特性。

Python编程有其特有的方式 phyonic,下面是python编程一些有意思的用法和技巧,总结出来,以飨大家。

1.原位替换

a = 5
b = 6
c = 7
a, b = b,a
print("a is",a)
print("b is",b)
a is 6
b is 5
a, (b,c) = c, (a,b)
print("a is",a)
print("b is",b)
print("c is",c)
a is 7
b is 6
c is 5
first, *middle_all, last = (1,2,3,4,5,6)
middle_all
[2, 3, 4, 5]

2.链似比较

x = 5
print("1 < x < 10 is", 1 < x < 10)
print("10 > x <= 9", 10 > x <= 9)
1 < x < 10 is True
10 > x <= 9 True

3.列表索引,反序

a = [1,2,3,4,5]
a[::2]
[1, 3, 5]
a[::-1]
[5, 4, 3, 2, 1]

4.列表解析、字典解析、元组解析

# 获得系列坐标点
a = ((i,j) for i in range(3) for j in range(2))
a
<generator object <genexpr> at 0x7fa8d04d7fc0>
for i in a:
    print(i)
(0, 0)
(0, 1)
(1, 0)
(1, 1)
(2, 0)
(2, 1)
str1 = "I love sheng xin bao dian"
print([i for i in str1.split() if i.endswith('n')])
['xin', 'dian']
a = {i:i*2 for i in range(5)}
a
{0: 0, 1: 2, 2: 4, 3: 6, 4: 8}
[(x, y) for x in range(4) if x % 2 == 1 for y in range(4)]
[(1, 0), (1, 1), (1, 2), (1, 3), (3, 0), (3, 1), (3, 2), (3, 3)]

5.集合操作

a = set(["sheng", "xin", "bao","dian","best","tutotials"])
b = set(["hong", "ji", "yin","zu","best","tutotials"])
a | b  # union
{'bao', 'best', 'dian', 'hong', 'ji', 'sheng', 'tutotials', 'xin', 'yin', 'zu'}
a & b # intersection
{'best', 'tutotials'}
a ^ b # Symmetric Difference
{'bao', 'dian', 'hong', 'ji', 'sheng', 'xin', 'yin', 'zu'}

6.Negative round

print("round整数:",str(round(1234.5678, -2)))

print("round小数:",str(round(1234.5678, 2)))
round整数: 1200.0
round小数: 1234.57

7.多行字符串的优雅嵌套

# \可以,但是第二行需要起头
system_command = "s-plot pheatmap -f matrix \
-t heatmap -a TRUE"
print(system_command)
s-plot pheatmap -f matrix -t heatmap -a TRUE
# 字符串中包含换行符
# 切第二行要起头,不然会有较多空格
system_command = """s-plot pheatmap -f matrix
-t heatmap -a TRUE""" 
print(system_command)
print(system_command.replace('\n', ' '))
s-plot pheatmap -f matrix 
-t heatmap -a TRUE
s-plot pheatmap -f matrix  -t heatmap -a TRUE
# 类元组的写法,既可以跨行,又可以自由格式
# 需要注意2点
   # 类元组,无逗号
   # 字符串连接时不会自动加空格,空格需要保存在字符串里面 
system_command = ("s-plot pheatmap -f matrix "
                  "-t heatmap -a TRUE")
print(system_command)
s-plot pheatmap -f matrix -t heatmap -a TRUE
# 多一步join;
system_command = ["s-plot pheatmap -f matrix",
                  "-t heatmap -a TRUE"]
print(' '.join(system_command))
s-plot pheatmap -f matrix -t heatmap -a TRUE

8.zip转换两个列表为字典

keyL = [1,2,3]
valueL = ['a','b','v']
for i ,j in zip(keyL, valueL):
    print(i,j)
1 a
2 b
3 v
import pprint
pprint.pprint(dict(zip(keyL, valueL)))
{1: 'a', 2: 'b', 3: 'v'}
dict([(i,j) for i ,j in zip(keyL, valueL)])
{1: 'a', 2: 'b', 3: 'v'}

9.enumerate索引列表 (不再使用len)

a = ['s','x','b','d']

# Preferred way
for index, item in enumerate(a):
    print(index,item)

print("\n")

# Old way
for i in range(len(a)):
    print(i,a[i])
    
0 s
1 x
2 b
3 d


0 s
1 x
2 b
3 d

10.矩阵转置

a = [(1,2), (3,4), (5,6)]
b = zip(*a)
for i in b:
    print(i)
(1, 3, 5)
(2, 4, 6)

11.sum的另一用法,二维数组秒变1维

aList = [[1, 2, 3], [4, 5], [6], [7, 8, 9]]
sum(aList, [])
[1, 2, 3, 4, 5, 6, 7, 8, 9]

12.解释正则表达式

import re
re.compile("^[a-z]*$", re.DEBUG)
AT AT_BEGINNING
MAX_REPEAT 0 MAXREPEAT
  IN
    RANGE (97, 122)
AT AT_END
re.compile("^[a-z][0-9]+$", re.DEBUG)
AT AT_BEGINNING
IN
  RANGE (97, 122)
MAX_REPEAT 1 MAXREPEAT
  IN
    RANGE (48, 57)
AT AT_END
re.compile("^[a-z]([0-9]+)$", re.DEBUG)
AT AT_BEGINNING
IN
  RANGE (97, 122)
SUBPATTERN 1 0 0
  MAX_REPEAT 1 MAXREPEAT
    IN
      RANGE (48, 57)
AT AT_END

13.for..else;若for循环中未执行break,则else会被执行

found = False
for i in range(2,5):
    if i == 1:
        found = True
        break
if not found: 
    print("i was never 1")
i was never 1
for i in range(2,5):
    if i == 1:
        break
else:
    print("i was never 1")
i was never 1

14.启动网络服务器,用于文件预览或传输

# run in commang line
# python -m http.server 8000
Serving HTTP on 0.0.0.0 port 8000 (http://0.0.0.0:8000/) ...

15.python打印九九乘法表 (列表解析)

print('\n'.join([' '.join(['%s*%s=%-2s' % (y,x,x*y) for y in range(1,x+1)]) for x in range(1,10)]))
1*1=1 
1*2=2  2*2=4 
1*3=3  2*3=6  3*3=9 
1*4=4  2*4=8  3*4=12 4*4=16
1*5=5  2*5=10 3*5=15 4*5=20 5*5=25
1*6=6  2*6=12 3*6=18 4*6=24 5*6=30 6*6=36
1*7=7  2*7=14 3*7=21 4*7=28 5*7=35 6*7=42 7*7=49
1*8=8  2*8=16 3*8=24 4*8=32 5*8=40 6*8=48 7*8=56 8*8=64
1*9=9  2*9=18 3*9=27 4*9=36 5*9=45 6*9=54 7*9=63 8*9=72 9*9=81

16.map, filter

# 过滤大于4的元素
a = [3,4,5]
[i for i in a if i<4]
[3]
list(filter(lambda x: x<4, a))
[5]
# 每个元素加2
[i+2 for i in a]
[5, 6, 7]
map(lambda x: x+2, a)
<map at 0x7fa8d0502358>
list(map(lambda x: x+2, a))
[5, 6, 7]

算2的1000次方的各位数之和

sum(map(int, str(2**1000)))
1366

17.python打印心型

print('\n'.join([''.join([('Love'[(x-y) % len('Love')] if ((x*0.05)**2+(y*0.1)**2-1)**3-(x*0.05)**2*(y*0.1)**3 <= 0 else ' ') for x in range(-30, 30)]) for y in range(30, -30, -1)]))
                veLoveLov           veLoveLov               
            eLoveLoveLoveLove   eLoveLoveLoveLove           
          veLoveLoveLoveLoveLoveLoveLoveLoveLoveLov         
         veLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveL        
        veLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLov       
        eLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLove       
        LoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveL       
        oveLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLo       
        veLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLov       
        eLoveLoveLoveLoveLoveLoveLoveLoveLoveLoveLove       
         oveLoveLoveLoveLoveLoveLoveLoveLoveLoveLove        
          eLoveLoveLoveLoveLoveLoveLoveLoveLoveLove         
          LoveLoveLoveLoveLoveLoveLoveLoveLoveLoveL         
            eLoveLoveLoveLoveLoveLoveLoveLoveLove           
             oveLoveLoveLoveLoveLoveLoveLoveLove            
              eLoveLoveLoveLoveLoveLoveLoveLove             
                veLoveLoveLoveLoveLoveLoveLov               
                  oveLoveLoveLoveLoveLoveLo                 
                    LoveLoveLoveLoveLoveL                   
                       LoveLoveLoveLov                      
                          LoveLoveL                         
                             Lov                            
                              v                             
  1. python打印曼德勃罗集合
print('\n'.join([''.join(['*'if abs((lambda a:lambda z,c,n:a(a,z,c,n))(lambda s,z,c,n:z if n==0else s(s,z*z+c,c,n-1))(0,0.02*x+0.05j*y,40))<2 else' 'for x in range(-80,20)])for y in range(-20,20)]))
                                                                                *                   
                                                                                                    
                                                                                                    
                                                                          **                        
                                                                     ***********                    
                                                                    ************                    
                                                                      *********                     
                                                            *    * ************  * *                
                                                    ****** * *************************** *          
                                                    *************************************** *****   
                                                     *******************************************    
                                               *** ******************************************** *   
                                                **************************************************  
                                             *******************************************************
                      *        *              ***************************************************** 
                     **** ******* *          *******************************************************
                      *****************     ******************************************************* 
                   *********************** *********************************************************
                   *********************** ******************************************************** 
             **** ********************************************************************************  
*********************************************************************************************       
             **** ********************************************************************************  
                   *********************** ******************************************************** 
                   *********************** *********************************************************
                      *****************     ******************************************************* 
                     **** ******* *          *******************************************************
                      *        *              ***************************************************** 
                                             *******************************************************
                                                **************************************************  
                                               *** ******************************************** *   
                                                     *******************************************    
                                                    *************************************** *****   
                                                    ****** * *************************** *          
                                                            *    * ************  * *                
                                                                      *********                     
                                                                    ************                    
                                                                     ***********                    
                                                                          **                        

19.打开一个关于Python漫画的网站

import antigravity
CHENTONG
版权声明:本文为博主原创文章,转载请注明出处。
alipay.png WeChatPay.png

CHENTONG

CHENTONG
积微,月不胜日,时不胜月,岁不胜时。凡人好敖慢小事,大事至,然后兴之务之。如是,则常不胜夫敦比于小事者矣!何也?小事之至也数,其悬日也博,其为积也大。大事之至也希,其悬日也浅,其为积也小。故善日者王,善时者霸,补漏者危,大荒者亡!故,王者敬日,霸者敬时,仅存之国危而后戚之。亡国至亡而后知亡,至死而后知死,亡国之祸败,不可胜悔也。霸者之善著也,可以时托也。王者之功名,不可胜日志也。财物货宝以大为重,政教功名者反是,能积微者速成。诗曰:德如毛,民鲜能克举之。此之谓也。

生信宝典文章集锦

生信的作用越来越大,想学的人越来越多,不管是为了以后发展,还是为了解决眼下的问题。但生信学习不是一朝一夕就可以完成的事情,也许你可以很短时间学会一个交互式软件的操作,却不能看完程序教学视频后就直接写程序。也许你可以跟着一个测序分析流程完成操作,但不懂得背后的原理,不知道什么...… Continue reading

转录组专题分析开课啦

Published on October 06, 2018