首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

R」从gtf文件抽取基因idname

参考文章http://www.bioinfo-scrounger.com/archives/342计算FPKM值,发现计算完每个基因下所有外显子总长度后,记录都是ENSEMBL gene id,而我需要是...奇怪是GenomicFeatures既然把GTF文件读取进去了还抽取基因id了,但它就是不提供抽gene symbol功能。...最后还是自己动手,完成了6万个gene feature转换。 整个提取操作包装为函数了,输入可以是文件名或已经导入gtf文件数据框(最好还是文件吧)。由data.table包支持,速度杠杠!...message("Treat input as file") input = data.table::fread(input, header = FALSE) } else { data.table....*" gene_id = sub(pattern_id, "\\1", input[[9]]) gene_name = sub(pattern_name, "\\1", input[[9

4.6K50
您找到你想要的搜索结果了吗?
是的
没有找到

SQL 日期时间类型

date:日历日期,包括年(四位),月日。 time: 一天时间,包括小时,分秒。可以用变量time(p)来表示秒小数点后数字位数(默认是0)。 ...如果指定with timezone,则时区信息也会被存储 日期时间类型值可按如下方式说明: date:‘2018-01-17’ time:‘10:14:00’ timestamp:‘2018-01-...17 10:14:00.45’ 日期类型必须按照如上年月日格式顺序指定。...timetimestamp秒部分可能会有小数部分。...我们可以利用cast e as t形式表达式来讲一个字符串(或字符串表达式)e转换成类型t,其中t是date,time,timestamp一种。字符串必须符合正确格式,像本段开头说那样。

3.1K60

Java时间日期处理

一、日期处理 旧版本 Date 在Java 1.0,对日期时间支持只能依赖java.util.Date类。它在易用性上许多问题,下面就谈谈这个类缺点。 缺点一:易用性较差。...Java 8 日期时间:LocalDateLocalTime LocalDate类实例是一个不 可变对象,它只提供了简单日期,并不含当天时间信息。另外,它也不附带任何与时区相关信息。...日期时间种类都不包含时区信息。...时区处理是新版日期时间API新增 加重要功能,使用新版日期时间API时区处理被极大地简化了。跟其他日期时间类一 样,ZoneId类也是无法修改。...每个特定 ZoneId对象都由一个地区ID标识,比如: ZoneId romeZone = ZoneId.of("Europe/Rome"); 地区ID都为“{区域}/{城市}”格式,这些地区集合设定都由英特网编号分配机构

2.7K40

R RStudio 安装及 R Profile 配置 & 爬取 CRAN 上所有 R名称、发布日期标题

,事实上只要用 3.0.0 之后版本都差距不是很大。我电脑上 R 版本是 R version 3.6.2 (2019-12-12)。...安装常用一些 R 包 可以运行下面的几句命令快速安装一些 R 包: # 安装 devtools install.packages("devtools") # tidyverse 系列 R 包 install.packages...# 这是设定 R 基础绘图系统绘图主题 basetheme::basetheme("brutal") Mac 系统参考: # 把日期环境设置成英文 suppressMessages(suppressWarnings...一个小练习:爬取 CRAN 上所有 R名称、发布日期标题 可以爬清华镜像源列表:https://mirrors.tuna.tsinghua.edu.cn/CRAN/web/packages.../available_packages_by_date.html 显然,这是个表格,很容易爬取: library(rvest) library(tidyverse) # lubridate 是处理日期一个

3.5K40

了解 HTML ID 类之间区别。

每当我们决定学习新事物时,我们都会面临各种各样困难。理解我们想要学习概念是很重要。今天,我们将学习两个在成为程序员或开发人员时每天都会遇到常用概念。那就是 ID CLASS 概念。...对它们理解获取护照一样重要,如果您希望出国旅行的话。首先,什么是 ID?用简单英语来说,ID 代表身份证件。每个人都有某种身份,可以被识别。...在上面解释身份证明文件类比,当两个或更多人拥有完全相同身份名称、文件号、出生日期等时,这意味着某些地方出现了问题,需要重新检查更正。编程世界也是同样概念。...对于 Class 来说,不同的人可以拥有完全相同特征,比如姓名、号码、出生日期等,一切仍然都会很好。用我自己的话来描述,我会将类描述为松散灵活。它们不在乎具体性。...看一下当您编写代码时,类 ID 是如何在 HTML 写入示例。

11010

R语言数据分析利器data.table包 —— 数据框结构处理精讲

将一个R对象转化为data.tableR可以时矢量,列表,data.frame等,keep.rownames决定是否保留行名或者列表名,默认FALSE,如果TRUE,将行名存在"rn"行,keep.rownames...="id",行名保存在"id"行。...比如此例取出DT X 列为"a"行,"a"进行merge。on参数第一列必须是DT第一列 DT[....n列,.N(总列数,直接在j输入.N取最后一列),:=(直接在data.table上添加列,没有copy过程,所以快,有需要的话注意备份),.SD输出子集,.SD[n]输出子集第n列,DT[,....,采取data.table链接符合表达式 DT[v>1, sum(y), by=v] #对v列进行分组后,取各组v>1行出来,各组分别对定义y求和 DT[, .N, by=x] #用

5.6K20

r」dplyr 里 join 与 base 里 merge 存在差异

今天在使用连接操作时发现:虽然都是合并操作函数,dplyr 包里 *_join() 基础包里面的 merge() 存在差异,不同数据结构,结果也会存在偏差。...r4 r1 r3 r2 #> 1 S1 S2 S2 S1 S1 #> 2 S2 S1 S1 S2 S2 看起来似乎有点不可理喻,但实际上上面我构造数据集是有点特别的:前 2 个子集第 3 个子集是没有可以连接...本质上是 data.table 体格泛型函数不支持类似基础包操作。 如何编写代码支持对上述数据集连接操作?...: 构造两个集合 to_join be_join,to_join 初始化为数据集第一个子集,而 be_join 为其他子集。...如果 be_join 不为空,进行如下循环: 如果存在,则将这个子集 to_join 按共同列合并 如果不存在,使用循环位移一位,将当前 be_join 第 2 个子集移动为 第 1 个。

1.5K30

R语言ggtree:将进化树序列id改成物种名称

通常我们会使用比对好fasta文件构建进化树,fasta文件中大于号后内容就是最终进化树上文字标签。如果拿到进化树文件后你想替换掉其中一些内容,那该怎么办呢?...本篇推文介绍一下使用R语言ggtree包实现这个目的 这个问题是来源于公众号一位读者提问 ?...大家可以关注我公众号 小明数据分析笔记本 留言相关问题,如果我恰巧会的话,我会抽出时间介绍对应解决办法 首先你已经有了构建好进化树文件 (Synergus:0.1976902387,(((((Periclistus...image.png 第一列x就是进化树中原本序列名称 第二列y是想要替换成id名称 读入进化树文件 library(treeio) tree<-read.newick("ggtree_practice_aligned.fasta.treefile...image.png 把这个新进化树写出到文件里 write.tree(tree1@phylo,file = "pra.nwk") 这样就达成目的了 这里导出进化树文件没有了最初支持率信息,我们再通过一行代码给他加上就好了

2.5K10

R︱高效数据操作——data.table包(实战心得、dplyr对比、key灵活用法、数据合并)

(参考来源:R高效数据处理包dplyrdata.table,你选哪个?) ?...data.table包提供了一个非常简洁通用格式:DT[i,j,by],可以理解为:对于数据集DT,选取子集行i,通过by分组计算j。...(ID)] 三种数据筛选方式,dplyr包、base基础包、data.table包。其中,dplyr是select语句,data.table要注意.()表达方式。...data.table,还有一个比较特立独行函数: 使用:=引用来添加或更新一列(参考:R语言data.table速查手册) DT[, c("V1","V2") := list(round(exp(V1...参考文献: 些许案例,代码参考自以下博客,感谢你们辛勤: 1、R语言data.table简介 2、超高性能数据处理包data.table 3、R语言data.table速查手册 4、R高效数据处理包

7.7K43

Rstackunstack函数

我们用R做数据处理时候,经常要对数据格式进行变换。例如将数据框(dataframe)转换成列表(list),或者反过来将列表转换成数据框。...那么今天小编就给大家介绍一对R函数来实现这样功能。 这一对函数就叫做stackunstack。从字面意思上来看就是堆叠去堆叠,就像下面这张图展示这样。...那么R里面这两个函数具体可以实现什么样功能呢?下面这张图可以帮助大家来理解。unstack就是根据数据框第二列分组信息,将第一列数据划分到各个组,是一个去堆叠过程。...一、unstack 下面我们来看几个具体例子 例如现在我们手上有一个数据框,里面的数据来自PlantGrowth 我们可以先看看PlantGrowth 内容,第一列是重量,第二列是不同处理方式...df = PlantGrowth unstacked_df = unstack(df) unstacked_df 结果如下,因为这里ctrl,trt1trt2样本刚好都是10个,所以这里结果看上去还像是一个数据框

5.1K30

Rgrepgrepl函数

在日常数据分析过程,我们经常需要在一个字符串或者字符串向量查找是否包含我们要找东西,或者向量那几个元素包含我们要查找内容。...这个时候我们会用到R中最常用两个函数,grepgrepl。...其实grep这个函数也并非是R所特有的,在linux模式匹配也用grep这个函数,前面我就给大家简单介绍过☞Linux xargs grep zgrep命令。...我们先来看看grepgrepl这两个函数用法。 这两个函数最大区别在于grep返回找到位置,grepl返回是否包含要查找内容。接下来我们结合具体例子来讲解。...☞讨论学习Rgrepl函数 参考资料: ☞Linux xargs grep zgrep命令 ☞讨论学习Rgrepl函数

2.3K10

NetT-sql日期函数操作

net日期函数代码:  代码 数据库日期函数: 函数 参数/功能 GetDate( ) 返回系统目前日期与时间 DateDiff (interval,date1,date2) 以interval...指定方式,返回date2 与date1两个日期之间差值 date2-date1 DateAdd (interval,number,date) 以interval指定方式,加上number之后日期...DatePart (interval,date) 返回日期date,interval指定部分所对应整数值 DateName (interval,date) 返回日期date,interval指定部分所对应字符串名称...参数 interval设定值如下: 值 缩 写(Sql Server) (Access ASP) 说明 Year Yy yyyy 年 1753 ~ 9999 Quarter Qq q 季...1 ~ 4 Month Mm m 月1 ~ 12 Day of year Dy y 一年日数,一年第几日 1-366 Day Dd d 日,1-31 Weekday Dw w 一周日数,一周第几日

1.3K60

一行代码对日期插值

在分析时,我们为了获得完整时间序列就需要“插入”那些丢失日期。 举一个例子: ? 这个数据集中有5行观测,2组分类(id等于12)。...我们看到每个id对应date都是有缺失,例如从2001-01-09直接跳到了2001-01-12,当中少了10号11号。 如何只用一行代码就高效优美地把这些缺失日期补上呢?...在merge过程,我们指定iddate变量必须匹配,也即on = .(id, date)语句作用: # 把CJ函数merge回原始数据集 dt[CJ, on = ....例如,在我们样例数据集sampleid=1观测对应日期最小值为01-08,最大值为01-14,而我们希望填充这两个日期“之间”所有值。...思路情况1类似,我们先构造CJ数据集,只不过在这里我们seq函数起讫点不再是固定值,而是每个id对应日期最大值与最小值: # 建立完整日期序列 # 注意minmax函数作用 CJ <- dt

1.4K30

PHPSession ID实现原理分析实例解析

session 工作机制: 为每个访问者创建一个唯一 id (UID),并基于这个 UID 来存储变量。UID 存储在 cookie ,亦或通过 URL 进行传导。...+ 当前时间(微妙)+ PHP自带随机数生产器) 从以上hash_func(*)数据采样值内容分析,多个用户在同一台服务器时所生产PHPSESSIONID重复概率极低。...生成session文件名规则即为sess_PHPSESSID,session文件存在session.save_path。)...rand(100,999); var_dump($_SESSION); Session实例问题 现有系统A,B; 假设A系统是可以独立运行web系统,即可以浏览器直接处理session, B系统是基于...cookie,首先session是一个只要活动就不会过期东西,只要开启cookie,每一次会话,session_id都不会改变,我们可以根据session_id来判断用户是否是正常登陆,防止用户伪造

4.4K10
领券