¿Cómo comparar la repetibilidad (ICC) de los diferentes grupos?

Question

He calculado los valores de ICC para los dos grupos y quisiera comparar la CPI valores para determinar si los grupos difieren en su repetibilidad. En la literatura se han utilizado pruebas t para comparar la repetibilidad, pero no está claro para mí cómo hacer esto.

Por ejemplo, con los datos de prueba:

ID  gr  day behaviour
1   1   1   0.361
2   1   1   0.232
3   1   1   0.240
4   1   1   0.693
5   1   1   0.483
6   1   1   0.267
7   2   1   0.180
8   2   1   0.515
9   2   1   0.485
10  2   1   0.567
11  2   1   0.000
12  2   1   0.324
1   1   2   0.055
2   1   2   0.407
3   1   2   0.422
4   1   2   0.174
5   1   2   0.613
6   1   2   0.311
7   2   2   0.631
8   2   2   0.283
9   2   2   0.512
10  2   2   0.127
11  2   2   0.000
12  2   2   0.000

Puedo obtener la repetibilidad de las medidas para el grupo 1 y 2 de la siguiente manera:

library(ICC)
g1 <- ICCest(ID, behaviour, data=dummy[dummy$gr=="1",])
g2 <- ICCest(ID, behaviour, data=dummy[dummy$gr=="2",])

Pero, ¿cómo puedo determinar ahora si la repetibilidad del grupo 1 es diferente de la del grupo 2?

Answer 1

Dejando de lado cuestiones importantes acerca de la cuestión del estudio y los datos de demostración (y el tamaño de la muestra para obtener una estimación razonable de una CPI), la salida que está recibiendo de la ICCest función de los intervalos de confianza que se adjunta: como punto de partida para la comparación de grupos, usted podría considerar si existe una superposición entre cada intervalo de confianza y el otro grupo del cálculo del punto de la CPI.

En cualquier caso, la presentación de informes el cálculo del punto de la CPI y el intervalo de confianza para cada grupo va a ser más útil (y, por lo tanto, me gustaría recomendar la elaboración de informes de estos en cualquier instancia) que los informes sólo el punto de estimaciones y el resultado de algún tipo de prueba de hipótesis.

dummy <- structure(list(ID = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L,
                  11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L), 
           gr = c(1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 
                  1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L), 
           day = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 
                   2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), 
           behaviour = c(0.361, 0.232, 0.24, 0.693, 0.483, 0.267, 0.18, 0.515, 0.485,
                         0.567, 0, 0.324, 0.055, 0.407, 0.422, 0.174, 0.613, 0.311, 
                         0.631, 0.283, 0.512, 0.127, 0, 0)), 
           .Names = c("ID", "gr", "day", "behaviour"), 
          class = "data.frame", row.names = c(NA, -24L))

library(ICC)
ICCest(ID, behaviour, data=dummy[dummy$gr=="1",])
# First few lines of console output:
#$ICC
#[1] -0.1317788
#$LowerCI
#[1] -0.7728603
#$UpperCI
#[1] 0.6851783

ICCest(ID, behaviour, data=dummy[dummy$gr=="2",])
# First few lines of console output:
#$ICC
#[1] 0.1934523
#$LowerCI
#[1] -0.6036826
#$UpperCI
#[1] 0.8233986