我想要一些关于如何计算 R 中变量组合的 ICC 的帮助/示例示例。
例如，下面是示例数据集:- 
现在我想对所有 2 个基因组合(列是基因名称)和各自的 ICC 进行详尽的搜索，然后是 3 个基因组合，然后是 4 个基因组合。发布我想列出 ICC 值排序的组合
输出看起来像这样
我尝试使用下面的代码做同样的事情，但我正在努力完成这个详尽的搜索和组合事情
mat1 <- as.matrix(lmm_1[,1:10])
icc(mat1, model = c("t"),r0 = 0,conf.level = 0.95)
icc(lmer( YWHAZ~ ., data = mat1))
我的数据帧结构lmm_1如下，请帮忙
structure(list(`HKG1 (18S)` = c(140922.090379318, 114050.66128233,
231970.660810084), `HKG2 (GAPDH)` = c(5821646.94272285, 4054004.07527754,
8279934.32892469), `HKG3 (ACTB)` = c(513559.47728841, 682673.286801917,
890579.00053991), `HKG4 (HMBS)` = c(963904.825829829, 1298738.62059259,
2091567.68645239), `HKG5 (HPRT)` = c(5508533.08283376, 5516876.07150713,
10236502.8076606), `HKG6 (SDHA)` = c(115259.913117015, 127680.131645604,
260533.935040935), `HKG7 (PABPN1)` = c(1483557.685172, 1212117.22511327,
1823406.16115662), `HKG8 (TBP)` = c(1835256.87785655, 1633084.68373014,
2478000.07992693), `HKG9 (YWHAG` = c(5892446.87038693, 6039118.46057885,
11746762.2003624), `HKG10 (YWHAZ )` = c(250987.289730308, 250536.411475335,
814508.745534441), Group = structure(c(1L, 1L, 1L), .Label = c("1",
"2", "3", "4", "5", "6"), class = "factor"), ID = c(1, 2, 3)), row.names = c(NA,
3L), class = "data.frame")
最佳答案
这是一个可能有帮助的快速代码示例(第二个示例)。 https://www.kaggle.com/rtatman/intraclass-correlation-coefficients-example/
关于r - 计算R中所有列组合的类内系数，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/52134296/