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

He calculado los valores de ICC para dos grupos y ahora me gustaría comparar los valores de ICC para determinar si los grupos difieren en su repetibilidad. En la literatura, la gente simplemente ha utilizado las pruebas t para comparar la repetibilidad, pero no me queda claro cómo hacerlo.

Por ejemplo, con los datos ficticios:

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 las medidas de repetibilidad 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 grupo1 es diferente del grupo2?

Dejando de lado los problemas materiales sobre la pregunta de estudio y los datos de demostración (y el tamaño de la muestra para obtener una estimación razonable de un ICC), el resultado que obtiene de la ICCestfunción tiene intervalos de confianza adjuntos: como punto de partida para comparar grupos, puede considerar si existe una superposición entre cada intervalo de confianza y la estimación puntual del otro grupo de la CPI.

En cualquier caso, informar la estimación puntual de la CPI y el intervalo de confianza para cada grupo será más útil (y, por lo tanto, recomendaría informar esto en cualquier caso) que informar solo las estimaciones puntuales y el resultado de algún tipo de prueba de hipotesis.

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