该博文更新已整理到新地址：

http://qinqianshan.com/r_code/r_draw/r-plotrix/R

当中的坐标中断一般都使用plotrix库中的axis.break(), gap.plot(), gap.barplot(), gap.boxplot()等几个函数来实现，例：

本博文已整理到新地址：http://qinqianshan.com/mothur-16s-double-end-sequencing/

老板不知道从哪弄了一堆肠道微生物的16s高通量数据丢给我，说了一通要求，然后让我分析。而我向来不善于拒绝，秉着“不会，硬着头皮也要弄出来”的原则，开始了回家前一段奇妙的mothur摸索之旅。

question:之前都是用qiime来处理数据，根据barcode来识别序列，给序列编号，同时去掉低质量的片段；而且不拼接。但这批来自苏州某测序公司的数据，既没有给我们barcode也没有给primer。所以有三个问题对我来说是新的：1，去低质量的片段；2双端的拼接；3给每条序列编号。

该博文的更新已整理到新地址：http://qinqianshan.com/biology/biodatabase/sliva/    

sliva rRNA数据库（http://www.arb-silva.de/）用来检查和比对RNA序列，既可以针对16S/18S,SSU，也可以针对23S/28S, LSU，包括了Bacteria, Archaea and Eukarya。同时也是ARB的官方指定数据库。

数据集

seqs      WS_bac            DSZ_bac           DSZ_arc            WS_arc

10         1            1            1            4.35505

20  &nb

1跨膜蛋白的定义

由双层脂类组成的生物膜结构是细胞的基本构成单元，然而，生物膜的大部分功能是由镶嵌在生物膜中的蛋白质来完成的，膜蛋白决定了生物膜的功能特性，因此，不同类型的生

物膜，其蛋白构成比例有很大的差别，以质量来算，从25%到75%

    在固氮菌的基因组中，固氮相关的nif( vnf/anf) 基因都是作为一个或几个基因簇存在的，称为固氮基因簇( nitrogen-fixationgene cluster) ， 是生物固氮的遗传基础，利用所有已知固氮菌nif 基因簇中的共有基因簇nifHDKENB发现新的可能的固氮菌。目前已知的固氮酶除了常见的由nifHDK

该博文已整理到新地址：

http://qinqianshan.com/python/py_module/matplotlib/

1，创建一个可以点的按钮

>>> def hello():

         print('hello there')

>>> from tkinter import *

>>> tk=Tk()

>>>

简单的输入：

在交互式环境的提示符>>>下，直接输入代码，按回车，就可以立刻得到代码执行结果。现在，试试输入100+200

该博客已整理到新地址：http://qinqianshan.com/ggplot2-legend图例的修改/

> head(msleep)

                        name      genus  vore        order conservation sleep_total sleep_rem sleep_cycle awake brainwt  bodywt

1                    Cheetah   Acinonyx carni    Carnivora           lc        12.1        NA          NA  11.9&nb

系统聚类法不仅需要度量个体与个体之间的距离，还要度量类与类之间的距离。类间距

离被度量出来之后，距离最小的两个小类将首先被合并成为一类。由类间距离定义的不同产

生了不同的系统聚类法。SPSS