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

【问题标题】：How to find the number of discordant and concordant pairs in R?如何在R中找到不一致和一致对的数量？
【发布时间】：2015-07-16 06:36:34
【问题描述】：

我正在尝试在临床试验中查找不一致和一致对的数量，并且遇到了提供函数 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（用于一致和不一致）相加，但我不知道如何计算矩阵中的值。或者，是否有人有更好的方法来计算一致/不一致，然后我很想知道。

【问题讨论】：

你想要的输出是什么样的？
我想计算一致和不一致对。总体目的是进行功效计算以获得样本量。

标签： r matrix statistics

【解决方案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

【讨论】：