我计划运行一些线性混合效果模型(对我来说是一种新方法)。我读到一个应该报告ICC(类间相关系数)。
我下载了几个包,但是计算失败。
ICC(DF, missing = T)
data.frame(x.s, subs = rep(paste("S", 1:n)错误。Obs, sep = "), nj)):参数表示不同的行数:898、2245警告信息:在stack.data.frame(x)中:非向量列将被忽略
这是我的数据:
DF <- structure(list(ID = c("SR6", "YLG19", "YLG19", "SR5", "SR2",
"TG5", "FB7", "SR9", "KBU15", "FB5"), sub_group = structure(c(2L,
2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L), .Label = c("European Bullhead",
"Salmonids"), class = "factor"), taxa = c("salmo.trutta", "oncorhynchus.mykiss",
"oncorhynchus.mykiss", "salmo.trutta", "salmo.trutta", "salmo.trutta",
"cottus.gobio", "cottus.gobio", "cottus.gobio", "cottus.gobio"
), sampling.site = c("oberer.seebach.ritrodat", "ybbs.lunz.grossau",
"ybbs.lunz.grossau", "oberer.seebach.ritrodat", "oberer.seebach.ritrodat",
"tagles.unten", "faltlbach", "oberer.seebach.ritrodat", "kothbergbach.unten",
"faltlbach"), body_weight_g = c(4L, 8L, 8L, 20L, 26L, 42L, 6L,
10L, 4L, 6L), PUFA = structure(c(3L, 4L, 2L, 3L, 2L, 1L, 3L,
1L, 4L, 3L), .Label = c("SDA", "EPA", "ARA", "DHA"), class = "factor"),
organ = structure(c(2L, 3L, 3L, 3L, 4L, 4L, 4L, 3L, 1L, 3L
), .Label = c("Brain", "Eyes", "Liver", "Muscles"), class = "factor"),
isotopic_value = c(-36.7301983, -39.5973755, -40.549113,
-35.6261828, -36.4038883, -46.085506, -39.0796303, NA, -41.6335499,
-41.484535)), row.names = c(289L, 488L, 487L, 280L, 242L,
367L, 52L, 308L, 189L, 19L), class = "data.frame")
这是我的LMM:
isotopic_value ~ organ + body_weight_g + (1 | ID)
我做错了什么?欢呼,Nadine
cols_ID <- sapply(1:ncol(DF),function(x) is.numeric(DF[,x]))
ICC(DF[,which(cols_ID)])