r - 如何找到R中不一致和一致对的数量？

Question

我正在尝试查找临床试验中不一致和一致对的数量，并且遇到了提供 ConDis.matrix 函数的“asbio”库。 (http://artax.karlin.mff.cuni.cz/r-help/library/asbio/html/ConDis.matrix.html)

他们给出的示例数据集是:

crab<-data.frame(gill.wt=c(159,179,100,45,384,230,100,320,80,220,320,210),
body.wt=c(14.4,15.2,11.3,2.5,22.7,14.9,1.41,15.81,4.19,15.39,17.25,9.52))
attach(crab)
crabm<-ConDis.matrix(gill.wt,body.wt)
crabm

结果如下:

   1  2  3  4  5  6  7  8  9 10 11 12
1  NA NA NA NA NA NA NA NA NA NA NA NA
2   1 NA NA NA NA NA NA NA NA NA NA NA
3   1  1 NA NA NA NA NA NA NA NA NA NA
4   1  1  1 NA NA NA NA NA NA NA NA NA
5   1  1  1  1 NA NA NA NA NA NA NA NA
6   1 -1  1  1  1 NA NA NA NA NA NA NA
7   1  1  0 -1  1  1 NA NA NA NA NA NA
8   1  1  1  1  1  1  1 NA NA NA NA NA
9   1  1  1  1  1  1 -1  1 NA NA NA NA
10  1  1  1  1  1 -1  1  1  1 NA NA NA
11  1  1  1  1  1  1  1  0  1  1 NA NA
12 -1 -1 -1  1  1  1  1  1  1  1  1 NA

我能想到的解决方案是将 1 和 -1 分别相加(表示一致和不一致)，但我不知道如何计算矩阵中的值。或者有人有更好的方法来计算一致/不一致，那么我很想知道。

Answer 1

您找到的解决方案是

sum(crabm == 1, na.rm = TRUE)
[1] 57
sum(crabm == -1, na.rm = TRUE)
[1] 7

你可以尝试(C...一致，D...不一致对):

library(DescTools)
tab <- table(crab$gill.wt, crab$body.wt)
ConDisPairs(tab)[c("C","D")]

$C
[1] 57

$D
[1] 7