不积小流 , 无以成江海 !分享 http://blog.sciencenet.cn/u/xiongchaoliang

博文

按标题搜索
用R做heatmap图(附代码)
2014-12-3 22:35
众所周知,R是统计分析大师,拥有大量开源的包,同时R也是一个绘图大师,只不过它的各种参数让小白们头疼。今天就给大家介绍一下如何用R做芯片表达数据的heatmap图(PS:当然可以用于其他数据)。 代码如下图:我主要用了heatmap.2函数,来自gplots包,该包不需要特异加载。(下面的代码中还贴了用limma包筛选差异表达 ...
个人分类: 【技术-可视化】|10118 次阅读|没有评论
fasta与fastq格式文件解读
2014-12-3 22:32
1 、 FASTA 文件的格式 在生物信息学中,FASTA格式(又称为Pearson格式)是一种基于文本的、用于表示核苷酸序列或氨基酸序列的格式。在这种格式中碱基对或氨基酸用单个字母来表示,且允许在序列前添加序列名及注释。 FASTA 文件以序列表示和序列作为一个基本单元,各行记录信息如下: 第一行是由大于号 ...
个人分类: 【测序概述】|15826 次阅读|没有评论
cytocape软件的使用
2014-12-3 22:28
在此为大家介绍一个开放源码的生物资讯软件 –Cytoscape ,它可以建构可视化的分子交互作用网络,并可将已有的基因表达信息 (gene expression profiles) 整合进此网络中,轻易观察分子间 ( 蛋白质 — 蛋白质 或 蛋白质 —DNA…) 的关联性。 Cytoscape 是 Institute ...
个人分类: 【技术-软件】|6968 次阅读|没有评论
高通量测序数据的质控工具---fastx_toolkit软件使用说明
2014-12-3 22:22
高通量测序数据下机后的原始 fastq 文件,包含 4 行,其中一行为质量值,另外一行则为对应序列,我们都了解高通量的数据处理首先要进行质量控制,这些过程包括去接头、过滤低质量 reads 、去除低质量的 3 ’ 和 5 ’ 端,去除 N 较多的 reads 等,而针对高通量测序数据的质控软件也有很多,在这里给 ...
个人分类: 【技术-软件】|15392 次阅读|没有评论
利用Tophat和Cufflinks进行RNA-Seq的数据分析流程(入门简介)
2014-12-3 22:15
其 基本流程图如下图 : 1、 主要软件及其版本 T ophat v2.0.8b 、bowtie2、cufflinks v2.1.1 、samtools 0.1.7 、python 2.6.5 、R 3.0 2、分析过程文件 说明 1 )tophat分析: 将每组样本的reads与参考基因组进行比对 分析结果文件:accepted_hits.bam、deletions.bed ...
个人分类: 【转录组-mRNA分析】|10747 次阅读|没有评论
使用UCSC和IGV查看reads在基因组上分布情况
2014-12-2 20:39
如果想要查看reads在基因组上的情况,可选用以下两种方法(当然方法很多,在此仅列举以下两种): 1.UCSC 2.IGV(Integrative Genomics Viewer) 以上两种方法均可以以多种格式的文件作为inFile,下面以BAM格式的inFile为例: (一)UCSC(适用于查看文件不太大的inFile中reads在基因组上的分布情况) 首先从R ...
个人分类: 【技术-可视化】|45764 次阅读|没有评论
查看linux下已设置好的环境变量所在路径的方法
2014-12-2 20:00
要想查看linux下已设置好的环境变量所在路径可用以下方法: $ which + 环境变量的名字
个人分类: 【技术-脚本】|2752 次阅读|没有评论
从哪里可以获取chromsize文件?
2014-12-2 17:01
通过UCSC可以获取chromsize文件,但是要注意版本(下面以hg19为例): UCSC---Downloads---Human---hg19---Annotation database---chromInfo.txt.gz
个人分类: 【基因组分析】|3988 次阅读|没有评论
bam 转 bedgraph
2014-11-30 00:19
有两种方法可以将bam转成bedgraph 方法1: samtools view my.bam my.sam wiggles my.sam my.wiggles.bedgraph 方法2: samtools sort my.bam my.sorted.bam genomeCoverageBed -bga -ibam my.sorted.bam -g chrom.sizes my.bedtools.bedgraph
个人分类: 【技术-脚本】|10722 次阅读|没有评论
bed 转 bam/sam ?
2014-11-30 00:17
bed如何转 bam/sam??? 可用bedtools的bedToBam: $ bedToBam Program: bedToBam (v2.12.0) Author: Aaron Quinlan (aaronquinlan@gmail.com) Summary: Converts feature records to BAM format. Usage: bedToBam -i bed/gff/vcf -g genome Options: ...
个人分类: 【技术-脚本】|6959 次阅读|没有评论

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-9-20 20:50

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部