接上回的故事,公司给了新一批的数据,但是批次效应比较重,需要去批次。去批次之后,fpkm出现了负值,但是counts在很多分析里不能直接用,所以需要counts转fpkm。目前网上关于人类基因组的已经比较多了,但是小鼠的我没咋找到,特别是biomart这个包我用的不熟练,绕了很多弯路,有点挠头。于是有了今天的文章,如何获取小鼠的基因长度。rm(list=ls())一键清空工作空间。listMarts()这个函数能够返回bioMRT可以连接的数据库列表,就是看看哪些数据库是能用的mart=useMart('ensembl')head(listDatasets(mart))dataset通过use