找到约 15 条结果

文章 什么是GSEA?

...影响。 这与之前讲述的GO富集分析不同。GO富集分析是先筛选差异基因,再判断差异基因在哪些注释的通路存在富集;这涉及到阈值的设定,存在一定主观性并且只能用于表达变化较大的基因,即我们定义的显著差异基因。而GSE...

文章 IF=4.3|棉花GT64基因家族的泛基因组分析及对黄萎病的表达响应

...并从Pfam数据库下载GT64的隐马尔科夫模型PF09258。通过HMMER搜索及Pfam、SMART、CDD数据库二次确认,最终共鉴定出39个GT64基因,分布在8个棉花基因组中。 对39个GT64基因家族成员进行氨基酸长度、蛋白质分子量、等电点、染色体位置...

文章 提取miRNA前后500bp序列

...9703021        29703094        1e-14   83.8 然后对比对结果筛选,尽量选择完全匹配,并且匹配长度最长的。筛选后的blast.out结果如下: bta-miR-34c     Chr20   100.00  113     0       0       1       113     39385825        39385713   ...

文章 Plant Communications | 浙大团队解析比克氏棉全基因组进化规律及其色素腺体和棉酚形成机制

...操》、《群体研究合辑》,课程中方法均提炼自高分期刊发表权威文献,提供全部分析代码及示例数据。习得之后即可对自己的群体测序数据进行分析,分析所得图片可直接用于发表。该系列课程在网易云课堂、腾讯课堂累计...

文章 菠萝WRKY基因家族分析

...进行了详细回答,感兴趣的可以到 www.omicsclass.com 网站上搜索关键词“tandem”就可以得到详细答案,我这里就不详细介绍了。今天想说的是一篇刚刚发表的关于基因家族分析的文章,其正好利用了MCScanX分析了菠萝中的WRKY基因家...

文章 OPLS-DA分析

...与Y相关和不相关的两类信息,通过去除不相关的差异来筛选差异变量。 通过OPLS-DA分析,每个代谢物可以得出一个VIP值,即变量重要性投影(VariableImportance inProjection,VIP),VIP值越大,代表该物质对于区分两组所具有的贡献越...

文章 层次聚类方法详解

...有的 cluster 之间的距离来计算合并后的新的 cluster 和剩余 cluster 之间的距离,这种计算方法将由下一部分中的 Lance-Williams 方法给出)。 Centroid/UPGMC 方法给每一个 cluster 计算一个质心,两个 cluster 之间的距离即为对应的...

文章 Pacbio三代全长转录组(Iso-Seq ) 采用CD-HIT合并转录本

Pacbio三代全长转录组采用isoseq3 进行了转录本的筛选,聚类和校正之后得到高质量的转录本序列,可以对序列再进行一定的合并操作。 采用cDNA_Cupcaake 文档中提供的cd-hit方案,运行命令如下: cd-hit-est -i <input> -o <output> ...

问题 基因家族解析报错

... : Starting MSA/Trees Species tree: Using 969 orthogroups with minimum of 81.2% of species having single-copy genes in any orthogroup Inferring multiple sequence alignments for species tree ------------------------------------------------------- 2025-06-07 14:57:14 : Done 0 of 969 2025-06-07 ...

文章 组学大讲堂云服务器conda环境配置:

组学大讲堂云服务器conda 经安装, 大家如果需要使用 conda 安装服务器上没有的软件可以运行以下命令配置一下: 1.添加conda命令到环境变量PATH中: echo 'export PATH=$PATH:/share/biosoft/conda/miniconda3/condabin/' >>~/.bashrc source ~/.b...

文章 perl或者python代码:批量修改进化树中基因ID或者物种ID名字 nwk文件

...修改,这里可以利用perl或者python 再结合正则表达式批量搜索替换,给出示例代码: python实现替换: import re,os tree="(chicken,((mouse,rat),(chimp,human)));" names=re.findall(r"[a-zA-Z0-9_]+",tree) f=open("newtree.nwk","w") for i in names: if not re.match(r...

文章 在家学习基因家族分析技能!-在线直播,2020年11月12日开班!

...10个分析技能,您也能写一篇基因家族的文章啦! 1、搜索基因家族成员蛋白保守结构域,鉴定基因家族成员; 2、预测基因家族成员的顺势作用元件; 3、构建系统进化发育树; 4、基因家族成员MOTIF分析及绘图; 5、...

文章 实验室linux生物信息服务器搭建CentOS7

...录,通常您自己创建的文件,都保存在这里,您最好给它分配一个分区 /usr----应用程序目录。大部分的软件都安装在这里,如果您计划安装许多软件,建议也给它分配一个分区 /var----如果您要作一些服务器方面的应用,可以考...

文章 GenomicsDBImport 的多样本变异检测

...一次建库用的参数:# --genomicsdb-workspace-path $my_database \#向存在的库添加新数据的参数:#--genomicsdb-update-workspace-path $my_database \$gatk --java-options "-Xmx$Xmx -Xms$Xms" \GenomicsDBImport \--genomicsdb-workspace-path $my_database \--intervals $interval_list \--b...

文章 R及R包的更新、移除、卸载

...load the namespace when a package is being detached. 如果要移除所有加载的包,可参考https://stackoverflow.com/questions/7505547/detach-all-packages-while-working-in-r 中的回答: lapply(names(sessionInfo()$loadedOnly), require, character.only = TRUE)invisible(lapply(paste0('p...