Làm thế nào để so sánh độ lặp lại (ICC) của các nhóm khác nhau?

Tôi đã tính các giá trị ICC cho hai nhóm và bây giờ muốn so sánh các giá trị ICC để xác định xem các nhóm có khác nhau về độ lặp lại hay không. Trong tài liệu người ta chỉ đơn giản sử dụng các bài kiểm tra t để so sánh độ lặp lại nhưng tôi không rõ làm thế nào để làm điều này.

Ví dụ: với dữ liệu giả:

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

Tôi có thể nhận được các biện pháp lặp lại cho nhóm 1 và 2 như sau:

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

Nhưng làm thế nào bây giờ tôi có thể xác định xem độ lặp lại của nhóm1 có khác với nhóm2 không?

Bỏ qua các vấn đề quan trọng về câu hỏi nghiên cứu và dữ liệu demo (và cỡ mẫu để có ước tính hợp lý về ICC), đầu ra bạn nhận được từ ICCesthàm có các khoảng tin cậy được đính kèm: như một điểm khởi đầu để so sánh các nhóm, bạn có thể xem xét liệu có sự chồng chéo giữa mỗi khoảng tin cậy và ước tính điểm của nhóm khác về ICC.

Ở mọi mức độ, báo cáo ước tính điểm của ICC và khoảng tin cậy cho mỗi nhóm sẽ hữu ích hơn (và do đó tôi khuyên bạn nên báo cáo những điều này trong mọi trường hợp) hơn là chỉ báo cáo ước tính điểm và kết quả của một số loại kiểm định giả thuyết.

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