1

我正在编写一个相当长的函数，它要求所有 colnames(abun) 都存在于 rownames(x) 中，反之则不然。如果不满足要求，我已经设计了它，以便函数抛出错误消息。除了错误消息，我还想告诉用户哪些 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","E")))

abun
    A   B   E
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               



if(any(colnames(abun) %in% rownames(x) != TRUE))
stop("The following species names in abun are missing trait information")
4

2 回答 2

2

回到你之前的问题

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

这应该返回您需要的值。

于 2013-10-02T08:56:08.157 回答
1

像这样的东西？

if(any(colnames(abun) %in% rownames(x) != TRUE))
stop("The following species names in abun are missing trait information:",
     paste(setdiff(colnames(abun), rownames(x)), collapse=" "))

感谢@Hadley 建议 setdiff！

于 2013-10-02T08:56:40.093 回答