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

我能想到的解决方案是分别加1s和-1s（对于一致和不一致），但我不知道如何计算矩阵中的值。或者，有人有更好的方法来计算一致/不和谐，然后我很想知道。

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

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

1 个答案: