【发布时间】:2021-12-09 02:07:14
【问题描述】:
我正在尝试计算两个评估者和多个受试者之间的类内相关性。但是,我还希望能够确定评分者对每个主题的同意程度,以便我可以确定哪些领域的一致性高/一致性低,以及哪些方面可能会提高/降低分数。
我知道我无法计算 一个主题的两个评分者之间的 icc,但是找到每行(主题)上的值之间的百分比一致性将为我提供具体信息。
注意:这些值是连续整数,不是因数
我的数据目前看起来像:
subject rater1 rater2
1 val1 val2
2 val3 val4
3 val5 val6
4 val7 val8
我想改变另一列“协议”,以按行计算百分比协议。所以,它基本上看起来像这样:
subject rater1 rater2 agreement
1 val1 val2 agreement_val1 %
2 val3 val4 agreement_val2 %
3 val5 val6 agreement_val3 %
4 val7 val8 agreement_val4 %
我知道我会用它来变异:
df %>% mutate(agreement = ???)
我只是想弄清楚计算方法以及它如何应用于每行/观察/主题。谢谢大家
EDIT1:我基本上想在每行中使用类似的功能。
这不起作用,但具有这种性质:
df %>% mutate(agreement = agree(each row))
EDIT2:你们中的一些人希望我提供一个例子。
在此示例中,有两个评分者(rater1 和 rater2)。他们都在评估动物执行特定行为的时间(以秒为单位)。数据框如下所示:
behavior rater1 rater2
run 48 59
stand 23 91
sit 389 401
sleep 288 290
所以,“主题”是动物的行为。我可以运行 ICC() 或 icc() 测试,但这会给我一个基于两个评估者对所有行为的协议的值。但是,我希望能够看到两个评估者在每个人的行为中是如何同意/不同意的。希望这是有道理的。
【问题讨论】:
-
您研究过任何软件包吗?当我在谷歌搜索 R 包时,我倾向于查看 CRAN 或 Bioconductor 标签;为此,
"CRAN" "ICC"(引用必要)建议使用恰当命名的ICC包;虽然最近没有更新,但建议你看看。 -
您能否分享一个包含一些数字和预期结果的示例?我不确定您要执行哪种计算。
-
如果您提供带有数字和相应预期输出的实际数据,会更容易提供帮助。
-
@r2evans 感谢您提供信息!我现在正在查看 CRAN 文档。老实说,直到你提到它,我才发现 ICC 包,所以我很感激!
-
@denisafonin 刚刚提供了一个例子