Calculer le coefficient Intra de classe pour toutes les combinaisons de colonne dans R

voix
1

Je voudrais un peu d'aide / par exemple l'échantillon à la façon de calculer la CCI pour une combinaison de variables dans R.

Ex Voici ensemble de données échantillon: - ensemble

Maintenant, je veux faire une recherche exhaustive pour toutes les 2 combinaisons de gènes (colonnes sont des noms de gènes) et ICC respectifs, puis 3 combinaison de gènes, puis pour quatre. Post que je veux à la liste des valeurs ICC WRT triés combinaison

La sortie ressemblerait à quelque chose comme ça

Sortie

J'ai essayé de faire la même chose en utilisant le code ci-dessous, mais je me bats pour travailler cette chose de recherche et combinaison exhaustive

mat1 <- as.matrix(lmm_1[,1:10])
icc(mat1, model = c(t),r0 = 0,conf.level = 0.95)
icc(lmer( YWHAZ~ ., data = mat1))

Ma structure dataframe lmm_1 est comme ci-dessous, S'il vous plaît aider à ce

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)
Créé 02/09/2018 à 05:41
source utilisateur
Dans d'autres langues...                            


1 réponses

voix
1

Voici un exemple de code rapide qui pourrait aider (le deuxième exemple). https://www.kaggle.com/rtatman/intraclass-correlation-coefficients-example/

Créé 06/09/2018 à 00:17
source utilisateur

Cookies help us deliver our services. By using our services, you agree to our use of cookies. Learn more