不同种属之间 同源基因的批量转换(R语言)

具体代码如下

#安装homologene这个R包
install.packages('homologene')
#加载homologene这个R包
library(homologene)
#我的列表是人类的基因列表
options(stringsAsFactors = F)

#读取代转换基因列表(如下图

genelist1<-read.table("human.txt",header = T,sep = "\t")

 #将上述读入的数据框格式转换为字符串的列表格

genelist2 <- genelist1$symbol


#使用homologene函进行转换
#@genelist2是要转换的基因
#@inTax是输入的基因列表所属的物种号,10090是小鼠
#@outTax是要转换成的物种号,9606是人
convt <- homologene(genelist2, inTax = 9606, outTax = 10090)

 #输出
write.table(convt,file = "convt.csv",sep = ",",row.names = F)

#输出homologene支持的物种号
homologene::taxData
#小鼠转换为human
mouse2human(genelist2)
#人源转换为小鼠
human2mouse(genelist2)

参考:如何快速查找物种间对应的同源基因 - 知乎

  • 4
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值