极客教程R语言 如何计算类内相关系数
在这篇文章中，我们将讨论如何在R编程语言中计算类内相关系数。相关性是用来获得两个变量之间的关系。
- 如果数值为1，那么这种关系是正相关的
- 如果数值为-1，那么这种关系就是负相关关系
- 如果数值为0，那么这种关系就没有关联性。
类内相关系数用于了解不同的评分者是否能对元素进行可靠的评分。其范围从0开始，以1结束。
我们可以使用icc()方法，它在irr包中可用，代表评分者之间的可靠性。
语法 ：icc(data, model, type, unit)
其中 。
- data是输入的数据帧
- model是要使用的模型类型。有两种类型：oneway和twoway。
- type是要计算的评分者之间的关系。选项包括 “一致性 “或 “一致性”。
- 单位是用于分析的，是单个或平均数
返回 ：它将返回以下内容。
- model – 模型的类型 – 单向或双向的
- 类型–模型的类型，即一致性
- 对象–数据框架中每一列的数据总数
- 评分者–数据框架的列数
- ICC值–内部相关系数值
- Ftest – F值和P值
- CI – 相关值的置信区间。
例1：计算单向模型的ICC
在这个例子中，我们将创建一个有4列的数据框架，并计算一个单一单元的单向模型的ICC
# load the library
library(irr)
# create dataframe with 4 columns
data = data.frame(col1=c(1:10), col2=c(34:43),
col3=c(20:29), col4=c(56:65))
# calculate icc for oneway model
icc(data, model = "oneway",
type = "agreement",
unit = "single")
输出 。
在这个输出中，正如代码中提到的，模型是单向的，类型是协议。
由于数据框中有四列，所以Raters=4
在每一列中，行的数量是10，所以受试者=10
内部相关系数（ICC）为-0.304。
我们得到的F检验值是0,0692，P值是1。
最后，我们得到的95%的置信区间水平是在-0.332到-0.233之间。
例2：计算双向模型的ICC
在这个例子中，我们将创建一个有4列的数据框架，并计算一个单一单元的双向模型的icc。
# load the library
library(irr)
# create dataframe with 4 columns
data = data.frame(col1 = c(1:10),
col2 = c(34:43),
col3 = c(20:29),
col4 = c(56:65))
# calculate icc for twoway model
icc(data, model = "twoway",
type = "agreement",
unit = "single")
输出 。
在这个输出中，正如代码中提到的，模型是双向的，类型是协议。
由于数据框中有四列，所以Raters=4
在每一列中，行的数量是10，所以受试者=10
内部相关系数（ICC）为-0.0168。
我们得到的F检验值是-5.44指数，P值是1。
最后，我们得到的95%的置信区间水平是在0.001到0.00之间。
例3：计算一个平均单位的icc
在这个例子中，我们将创建一个有4列的数据框架，计算平均单位的双向模型的icc。
# load the library
library(irr)
# create dataframe with 4 columns
data = data.frame(col1 = c(1:10),
col2 = c(34:43),
col3 = c(20:29),
col4 = c(56:65))
# calculate icc for twoway model with average unit
icc(data, model = "twoway", type = "agreement",
unit = "average")
输出 。
在这个输出中，正如代码中提到的，模型是单向的，类型是协议。
由于数据框中有四列，所以Raters=4
在每一列中，行的数量是10，所以受试者=10
内部相关系数（ICC）为0.0639。
我们得到的F检验值是-5.44 exp，P值是1。
最后，我们得到的95%的置信区间水平是在0.005到0.257之间。
例4：计算单个单位的icc
在这个例子中，我们将创建一个有4列的数据框架，并计算单一单位的双向模型的icc。
# load the library
library(irr)
# create dataframe with 4 columns
data = data.frame(col1 = c(1:10),
col2 = c(34:43),
col3 = c(20:29),
col4 = c(56:65))
# calculate icc for twoway
# model with single unit
icc(data, model = "twoway", type = "agreement",
unit = "single")
输出 。
在这个输出中，正如代码中提到的，模型是双向的，类型是协议。
由于数据框中有四列，所以Raters=4。
在每一列中，行数为10，所以受试者=10
内部相关系数（ICC）为0.0168。
我们得到的F检验值是-5.44 exp，P值是1。
最后，我们得到的95%的置信区间水平是在0.001到0.00之间。