setwd("D:/生信代码复现/1.新TCGA/1.新TCGA")#💚💚💚💚💚💚💚设置工作目录,改成自己的工作目录
#读取hsa.gff3的内容，跳过#开始的行
mir=read.table("hsa.gff3",comment.char = "#",sep="\t",stringsAsFactors = F)
#第三列为miRNA的行包含成熟体信息，具体在第九列
mature=mir[mir$V3=="miRNA",9]
#根据;Alias=,;Name=,;Derives_from=来拆分第九列的内容
#提取拆分开的向量中的第二和三个元素，MIMAT0027618  hsa-miR-6859-5p
#转置之后，强制转换成数据框，去除重复
human_mirs=data.frame(unique(t(sapply(strsplit(mature,";.*?=",fixed=F),"[",2:3))))
#将miRNA的ID号和名字保存到mirbase.rds中
saveRDS(human_mirs,file="mirbase.rds")

#读取mirbase.rds中的内容，可以赋给任意变量名
mirbase=readRDS("mirbase.rds")

#查看前几行
head(mirbase)