zoukankan      html  css  js  c++  java
  • R %in% ? merge() match()

    in

    > rc_24[1:5,]
               X24h.B3.1 X24h.B3.2 X24h.D.1 X24h.D.2
    B3-1000001      1298      1176     1347     2165
    B3-1000002       602       593      560     1030
    B3-1000003       272       286      296      432
    B3-1000004      4700      3895     4233     6764
    B3-1000005       102       112      124      232
    

    in错误

    t1 <- rc_24[1:5,]
    t2 <- rc_24[3:6,]
    a = t1[,t1[,3]%in%t2[,3]]
    a
    

    (subscript) logical subscript too long

    交换t1t2位置

    t1 <- rc_24[1:5,]
    t2 <- rc_24[3:6,]
    a = t1[,t2[,3]%in%t1[,3]]
    a
    
    > a
               X24h.B3.1 X24h.B3.2 X24h.D.1
    B3-1000001      1298      1176     1347
    B3-1000002       602       593      560
    B3-1000003       272       286      296
    B3-1000004      4700      3895     4233
    B3-1000005       102       112      124
    

    match

    c=match(t1[,3],t2[,3])
    c
    t1[c]
    
    > c
    [1] NA NA  1  2  3
    

    match数据框

    c=match(t1[,3],t2)
    c
    

    c
    [1] NA NA 9 10 11

    merge把基因名加到差异分析后

    anno <- read.table("annotation.table.xls",sep = ' ',header = F,fill = T)
    an = anno[,c(1,13)]
    an(gene <- an)V13
    an(pp <- an)V1
    an <-an[,-c(1,2)]
    head(an)

              pp        gene
    1 B3-1000001        dnaA
    2 B3-1000002 DPO3B, dnaN
    3 B3-1000003        recF
    4 B3-1000004        gyrB
    5 B3-1000005           _
    6 B3-1000006      K07024
    
    rc3max=as.data.frame(rc3max)
    rc3max$pp <- rownames(rc3max)
    rc39_a = merge(rc3max,an,by='pp')
    head(rc39_a)
    

    head(rc39_a)
    pp X3h.B3.1 X3h.B3.2 X24h.B3.1 X24h.B3.2 gene
    1 B3-1000001 3940 2795 1298 1176 dnaA
    2 B3-1000002 2968 1730 602 593 DPO3B, dnaN
    3 B3-1000003 828 525 272 286 recF
    4 B3-1000004 8008 7042 4700 3895 gyrB
    5 B3-1000005 55 68 102 112 _
    6 B3-1000006 401 424 230 219 K07024

    match

  • 相关阅读:
    Windows安装node.js
    02ARM体系结构
    01ARM介绍
    01软件过程模型
    创建一个JavaWeb工程
    python之面向对象01
    python学生管理系统
    python文件、文件夹的相关操作
    python文件的定位读写
    python文件的写入与备份
  • 原文地址:https://www.cnblogs.com/impw/p/13888455.html
Copyright © 2011-2022 走看看