0

マトリックスの列名がマトリックスの行名に含まれているかどうか、つまり、以下の例で colnames(abun) と rownames(x) に含まれているかどうかをテストしたいと思います

abun <- matrix(c(0.4,0,0.6,0.1,0.4,0.5), 
    nrow = 2, ncol = 3, byrow = TRUE, dimnames = list(c("x", "y"), 
    c("A","B","C")))

 abun
    A   B   C
x 0.4 0.0 0.6
y 0.1 0.4 0.5

x<-data.frame("Trait1" =c(1,1,0,1),
                    "Trait2"=c(1,1,1,1),
                    "Trait3" =c(1,1,0,1),
                    "Trait4" =c(1,0,1,1))
rownames(x)<-c("A","B","C","D") 

x
  Trait1 Trait2 Trait3 Trait4
A      1      1      1      1
B      1      1      1      0
C      0      1      0      1
D      1      1      1      1

更新: 私は関数を書いていますが、colnames(abun) が rownames(x) に含まれていない場合にエラー メッセージをスローしたいと考えています。私が試してみました：

if(colnames(abun) %in% rownames(x) = FALSE)
stop("species names in abun and x do not match")
4

2 に答える 2

1 
colnames(abun)[
colnames(abun) %in% rownames(x)
]

colnames(abun) %in% rownames(x)colnames(abun)は、 にどの要素が存在するかを示す true/false ベクトルを返しますrownames(x)。

于 2013-10-02T08:16:17.863 に答える
1

2 つのセットの交点を求めているだけですか?

intersect(c("A","B","C","D") ,
          c("A","B","C"))
[1] "A" "B" "C"

差を取得するには、次を使用しますsetdiff。

setdiff(c("A","B","C","D") ,
          c("A","B","C"))
[1] "D"
于 2013-10-02T08:26:31.173 に答える