找到约 15 条结果

文章 怎样画一张曼哈顿图CMplot

...图我们也是可以画出来的。 在 如何绘制SNP密度图 文章中，我们利用R中的 CMplot 包绘制了SNP密度图，其实 CMplot 包还可以绘制曼哈顿图哦！这里就来说说具体怎么操作吧。 CMplot 包链接：https://github.com/YinLiLin/R-CMplot 安装并...

文章 8月相约南京--群体遗传进化GWAS分析、T2T基因组、泛基因家族、宏基因组等主题培训！

...效的生信学习方式，是实现弯道超车的最佳选择： ① 集中时间培训，有效治疗“拖延症”，激发体内“自律因子”，学习更加高效； ② 讲师现场手把手指导，及时解决学习中遇到的难题； ③ 和同领域同行进行交流，掌握科...

文章 GEO数据挖掘生物信息文章解读（直肠癌）

...据库全称GENE EXPRESSION OMNIBUS，是由美国国立生物技术信息中心NCBI创建并维护的基因表达数据库。它创建于2000年，收录了世界各国研究机构提交的基因表达数据（主要包括芯片表达数据，也包含一些高通量测序表达数据）。这里...

文章 在Windows系统中安装并使用seqkit工具

...OWS\system32.将下载的seqkit.exe软件复制到本地C:\Windows\System32中 运行软件： 1. 打开windows系统的Windows power shell页面。输入命令：seqkit version。确认是否下载安装成功： 2. 之后就可以使用seqkit工具处理生信文件了。 seqkit工具的命...

文章 获取关系矩阵（WGCNA)

在进行WGCNA分析的过程（学习WGCNA）中需要基于表达矩阵转换关系矩阵，结合power值构建邻接矩阵，并由此构建TOM矩阵最终构建网络。在代码实现的过程中往往： 1、计算power值 2、基于power直接利用adjacency()由表达矩阵--邻接矩...

文章 iTRAQ/TMT蛋白质组学分析基本知识点

...不同来源的样品的蛋白鉴定。 3、蛋白质组学相关的文章中，选择的差异倍数经常有1.2倍、1.3倍、1.5倍、2倍等待的，这有没有一个统一的标准？ 目前蛋白组学分析大多数是按照差异倍数2倍或1.5倍进行分析的，但是我们要根据...

文章 blast XML结果解析，获得更多blast比对信息 perl

...，包括比对的序列信息，长度信息等等；接下来我们可以用bioperl包来解析，获取自己想要的信息； #!/usr/bin/perluse strict;use Bio::SearchIO;use Data::Dumper;my $fi     = $ARGV[0];my $out = $ARGV[1];my $format;$format||="blastxml";my $searchio     = new B...

文章 快速生成iTOL配置文件美化你的进化树

...带来了不便，这里介绍一个工具帮助大家可以从表格文件中快速生成itol需要的配置文件美化进化树。下面给出例子： . 示例数据 1. 进化树数据 abunt-tree.nwk： ((OTU_47:0.101892516,OTU_46:0.205918075)0.999:0.12637554249999997,((OTU_7:0.079452354,(OT...

文章 SCI文章投稿各阶段的电子邮件模板合集

...得尊重，因为第一次的投稿不知道具体负责的编辑，只能用通用的Editors）: On behalf of my co-authors, we thank you very much for giving us an opportunity to revise our manuscript, we appreciate editor and reviewers very much for their positive and constructive comments and ...

文章 hmmsearch搜索蛋白保守结构域结果说明

如果是用hmmsearch命令查询我们序列中的保守结构域： hmmsearch --domtblout hmm_WRKY.txt --cut_tc WRKY.hmm Arabidopsis_thaliana.TAIR10.pep.all.fa target length是我们输入序列总长度；query length 输入的查询序列也就是蛋白结构域的总长度；筛选的时...

文章 GEO数据库之芯片数据下载

...的数据信息，它只能是单一平台的数据，往往，GSE 和GDS中会包含多个GSM数据； 一个GPL Accession，则对应一个platform信息。 R包安装与加载 GEOquery ## try http:// if https:// URLs are not supportedsource("https://bioconductor.org/biocLite.R")biocLite("GEOq...

文章 农艺性状表型数据与环境关联互作分析R代码

...tGUI")library(GGEBiplotGUI)library("GGEBiplots")library("reshape2")library(vegan)library(corrplot)library("psych")library(Cairo)setwd("D:/potato/")data1<-read.table("张北data.txt",header = F,row.names = 1,comment.char = "",sep = "\t",encoding = "UTF-8")data.t<-as.data.frame(t(data1))data.l&l...

问题 转录组KEGG富集分析出现问题。

代码本中参考代码： cd $workdir/ mkdir 6.enrich cd 6.enrich #GO,KEGG 富集分析以下脚本只支持模式物种见表格：https://www.omicsclass.com/article/1244 #非模式物种参考：http://guangchuangyu.github.io/cn/2017/07/clusterprofiler-maize/ Rscript $scriptdir/enrichKEG...

文章 你了解SAM和BAM文件吗

...是列之间由Tab隔开，每一字段具体含义参考下图： 其中： 1. QNAME 表示reads名称； 2. FLAG：表示比对的结果，由数字表示，不同的数值含义不同，其列表如下： 中文解释： 1 ： 代表这个序列采用的是PE双端测序2： 代表...

问题 16S 差异比较使用compare_errorbar.r的时候报错

我用这个代码是他的前一步不报错，但是这一步报错，代表我的分组信息这些都是准确的 这个代码不报错 for i in {2..6};do     Rscript $scriptdir/meta_stat.r -i ../05.export_data/taxa_summary/absolute_abundance/feature_table_tax_L$i.txt \ -m $metadata -...