首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往
您找到你想要的搜索结果了吗?
是的
没有找到

用Matlab筛选mirbase,一种基于miRBase数据库的无参的miRNA数据分析方法与流程

本发明涉及转录组测序领域,具体涉及一种在miRBase数据库中无本物种参考miRNA数据的miRNA测序的数据分析方法。...技术实现要素: 为了克服现有技术所存在的上述缺陷,本发明的目的在于提供一种基于miRBase数据库的无参的miRNA数据分析方法。...为了实现本发明的目的之一,所采用的技术方案是:一种基于miRBase数据库的无参的miRNA数据分析方法,包括如下步骤: 步骤一,文件准备步骤: 准备并读取config文件,读取后生成相应的shell脚本...过滤序列用于后续分析; 步骤三,sRNA分类注释步骤: 将去重后的序列与Rfam数据库进行blast比对,筛选出碱基错配数小于2的结果,注释出其中的非编码RNA序列, 将其余的小RNA序列与miRBase...作为示例的,将其余的小RNA序列与miRBase数据库中该物种的miRNA成熟体序列进行Blast比对,筛选出碱基错配数小于2的结果,注释为已知的miRNA序列,同时计算测到的miRNA表达量,进行表达模式分析

56530

【R语言】获取最新的人的所有miRNA的ID号

#加载mirbase.rds文件,里面保存了人的所有miRNA的成熟体ID和miRNA名字 load("mirbase.rds") 其实,前面小编就用视频给大家介绍过,如何使用Excel来提取人的所有的...接下来小编就给大家讲讲如何使用R来从miRBase数据库中下载人的最新的miRNA注释信息,然后使用R来出来提取所有的miRNA的ID号。对miRBase这个数据库还不了解的小伙伴,请猛戳下面链接。...☞miRBase数据库介绍及miRNA数据下载 #miRBase数据库中人的miRNA的注释文件 link="https://www.mirbase.org/ftp/CURRENT/genomes/hsa.gff3...=",fixed=F),"[",2:3)))) #将miRNA的ID号和名字保存到mirbase.rds中 saveRDS(human_mirs,file="mirbase.rds") #读取mirbase.rds...中的内容,可以赋给任意变量名 mirbase=readRDS("mirbase.rds") #查看前几行 head(mirbase) 这段代码中用了saveRDS和readRDS这样一对函数来保存和读取数据

1.1K40

不同版本miRNA名字相互转换

miRBase数据库介绍及miRNA数据下载 如果我们要使用这套数据,做差异表达分析,然后再去预测miRNA的靶基因的时候,问题就来了。...现在绝大多数的miRNA靶基因预测软件使用的都是最新版本的miRbase数据库中miRNA的名字,绝大多数的miRNA名字都会带-3p或者-5p。...我们需要使用的是miRbase converter这个工具。 把你需要转换的miRNA的名字贴到下面的对话框中,选择相应的miRbase版本号,点击submit,就大功告成了。...参考资料: ☞GEO数据库数据检索方法(一) ☞GEO数据库数据检索方法(二) ☞GEO数据库数据页面介绍(三) ☞零代码差异表达分析工具:GEO2R ☞GEO芯片数据差异表达分析 ☞miRBase

40930

miRNA注释包:miRBaseVersions.db

这个包是基于miRbase数据库的。 由于数据库不断的增长和变化,miRNA的名称可能在不同的版本中有不同的名称,甚至不再被列为有效的miRNA。...键类型“MIMAT”是主表,包含所有支持的miRBase版本的所有记录。以前缀“VW-MIMAT”开头的键类型被称为SQL views。...你可以理解为不同的版本,例如,键类型“VW-MIMAT-22.0”是来自“MIMAT”表的SQL views,它只保存来自miRBase版本22.0的记录。...我们可以看到,结果返回了不同miRBase版本中加入的所有miRNA名称。参数columns = "*"表示返回所有列。...当然,除了用注释包以外,我们可从数据库找那个下载所有的miRNA信息文件,下面是地址: ftp://mirbase.org/pub/mirbase/CURRENT/ ? 下载后打开文件是这样的。

1.9K10
领券