Tính toán độ tin cậy giữa các nhà nghiên cứu trong R với số lượng xếp hạng thay đổi?


9

Wikipedia gợi ý rằng một cách để xem xét độ tin cậy giữa các nhà nghiên cứu là sử dụng mô hình hiệu ứng ngẫu nhiên để tính toán tương quan nội hàm . Ví dụ về tương quan nội hàm nói về việc nhìn vào

σα2σα2+σϵ2

từ một mô hình

Yij=μ+αi+ϵij

"Trong đó Y ij là quan sát thứ j trong nhóm thứ i , là một ý nghĩa tổng thể không quan sát được, α i là một hiệu ứng ngẫu nhiên không quan sát được chia sẻ bởi tất cả các giá trị trong nhóm i và ε ij là một thuật ngữ tiếng ồn không quan sát được."

Đây là một mô hình hấp dẫn đặc biệt bởi vì trong dữ liệu của tôi, không có người đánh giá nào đánh giá tất cả mọi thứ (mặc dù hầu hết đã xếp hạng 20+) và mọi thứ được đánh giá một số lần thay đổi (thường là 3-4).

Câu hỏi số 0: "Nhóm i" trong ví dụ đó ("nhóm i") có phải là một nhóm các thứ được xếp hạng không?

Câu hỏi số 1: Nếu tôi đang tìm kiếm độ tin cậy giữa các nhà nghiên cứu, tôi không cần một mô hình hiệu ứng ngẫu nhiên với hai thuật ngữ, một cho người đánh giá và một cho điều được đánh giá? Rốt cuộc, cả hai đều có thể thay đổi.

Câu hỏi số 2: Làm thế nào tôi thể hiện tốt nhất mô hình này trong R?

Có vẻ như câu hỏi này có một đề xuất đẹp mắt:

lmer(measurement ~ 1 + (1 | subject) + (1 | site), mydata)

Tôi đã xem xét một vài câu hỏi và cú pháp của tham số "ngẫu nhiên" cho lme là không rõ ràng đối với tôi. Tôi đọc trang trợ giúp cho lme , nhưng mô tả cho "ngẫu nhiên" là không thể hiểu được đối với tôi mà không có ví dụ.

Câu hỏi này có phần tương tự như một dài danh sách các câu hỏi , với này gần nhất. Tuy nhiên, hầu hết không đề cập đến R một cách chi tiết.


Mô hình hiệu ứng hỗn hợp và hiệu ứng ngẫu nhiên được mã hóa theo cùng một cách trong R. Xem ncbi.nlm.nih.gov/pmc/articles/PMC3402032 để biết thêm thông tin về tuto!
noé

Câu trả lời:


6

Mô hình bạn tham chiếu trong câu hỏi của bạn được gọi là "mô hình một chiều". Nó giả định rằng các hiệu ứng hàng ngẫu nhiên là nguồn phương sai duy nhất có hệ thống. Trong trường hợp độ tin cậy giữa các nhà nghiên cứu, các hàng tương ứng với các đối tượng đo lường (ví dụ: các đối tượng).

Mô hình một chiều : trong đó là giá trị trung bình của tất cả các đối tượng, là hiệu ứng hàng và là hiệu ứng dư.

xij=μ+ri+wij
μriwij

Tuy nhiên, cũng có "mô hình hai chiều". Chúng giả định rằng có phương sai liên quan đến hiệu ứng hàng ngẫu nhiên cũng như hiệu ứng cột ngẫu nhiên hoặc cố định. Trong trường hợp độ tin cậy giữa các nhà nghiên cứu, các cột tương ứng với các nguồn đo lường (ví dụ: các bộ đo).

Mô hình hai chiều : trong đó là trung bình cho tất cả các đối tượng, là hiệu ứng hàng, là hiệu ứng cột, là hiệu ứng tương tác và là hiệu ứng dư. Sự khác biệt giữa hai mô hình này là bao gồm hoặc loại trừ hiệu ứng tương tác.

xij=μ+ri+cj+rcij+eij
xij=μ+ri+cj+eij
μricjrcijeij

Đưa ra mô hình hai chiều, bạn có thể tính toán một trong bốn hệ số ICC: ICC thống nhất điểm số duy nhất (C, 1), ICC thống nhất điểm trung bình (C, k), thỏa thuận điểm số ICC (A, 1) hoặc thỏa thuận điểm trung bình ICC (A, k). Các ICC điểm đơn áp dụng cho các phép đo đơn (ví dụ: các chỉ số riêng lẻ), trong khi các ICC điểm trung bình áp dụng cho các phép đo trung bình (ví dụ: giá trị trung bình của tất cả các phép đo ). Các ICC nhất quán loại trừ phương sai cột khỏi phương sai mẫu số (ví dụ: cho phép các bộ đo thay đổi theo phương tiện riêng của chúng), trong khi các ICC thỏa thuận bao gồm phương sai cột trong phương sai mẫu số (ví dụ: yêu cầu các bộ đo thay đổi theo cùng một giá trị trung bình).xijx¯i

Dưới đây là các định nghĩa nếu bạn giả sử một hiệu ứng cột ngẫu nhiên:

Định nghĩa ICC hiệu ứng ngẫu nhiên hai chiều (có hoặc không có hiệu ứng tương tác) :

ICC(C,1)=σr2σr2+(σrc2+σe2) or σr2σr2+σe2
ICC(C,k)=σr2σr2+(σrc2+σe2)/k or σr2σr2+σe2/k
ICC(A,1)=σr2σr2+(σc2+σrc2+σe2) or σr2σr2+(σc2+σe2)
ICC(A,k)=σr2σr2+(σc2+σrc2+σe2)/k or σr2σr2+(σc2+σe2)/k

Bạn cũng có thể ước tính các giá trị này bằng cách sử dụng bình phương trung bình từ ANOVA:

Ước tính ICC hai chiều :

ICC(C,1)=MSRMSEMSR+(k1)MSE
ICC(C,k)=MSRMSEMSR
ICC(A,k)=MSR-MSE
ICC(A,1)=MSRMSEMSR+(k1)MSE+k/n(MSCMSE)
ICC(A,k)=MSRMSEMSR+(MSCMSE)/n

Bạn có thể tính các hệ số này trong R bằng cách sử dụng gói ir :

icc(ratings, model = c("oneway", "twoway"),
type = c("consistency", "agreement"),
unit = c("single", "average"), r0 = 0, conf.level = 0.95)

Người giới thiệu

McGraw, KO, & Wong, SP (1996). Hình thành các suy luận về một số hệ số tương quan nội hàm. Phương pháp tâm lý, 1 (1), 30 bóng46.

Cây bụi, PE, & Fleiss, JL (1979). Tương quan nội bộ: Sử dụng trong việc đánh giá độ tin cậy của người đánh giá. Bản tin tâm lý, 86 (2), 420


Cảm ơn câu trả lời tuyệt vời! Trong mô hình hai chiều trong icc trong R, làm thế nào để chúng ta biểu diễn lựa chọn ngẫu nhiên các bộ đo trên mỗi hàng? Ý tôi là, hãy tưởng tượng chúng ta có một nhóm gồm 100 người và mỗi môn được khoảng 5-10 người trong số họ. Một kịch bản như vậy có thể được xử lý bởi gói icc?
michal

Mỗi rater nên có cột riêng trong ma trận bạn cung cấp cho hàm icc. Mặt khác, tính toán là giống nhau cho các mô hình hiệu ứng ngẫu nhiên và hỗn hợp - sự khác biệt chính là trong diễn giải (mức độ tổng quát của các kết quả có thể được xem xét).
Jeffrey Girard

Cảm ơn câu trả lời! Tôi đang cố gắng làm điều đó, có hầu hết NA trong các ô (và chỉ có một vài giá trị với số thực trên mỗi cột, trong đó một người đánh giá cụ thể đã xếp hạng một chủ đề tương ứng với một hàng). Tuy nhiên, trong đầu ra tôi nhận được một văn bản nói rằng không có đối tượng nào được ghi lại (ví dụ: Chủ đề = 0 Raters = 9). Có lẽ nó có nghĩa là bất cứ nơi nào ít nhất một NA được tìm thấy, toàn bộ hàng được lọc ra? Nhưng sau đó, làm thế nào tôi có thể biểu thị xếp hạng bị thiếu từ một người đánh giá?
michal

Hmm có thể là một hạn chế của chức năng icc cụ thể này. Tôi có một kịch bản MATLAB có thể xử lý tình huống này. Bạn có tình cờ truy cập MATLAB không?
Jeffrey Girard

1
Vâng, hãy xem trang web của tôi: mrelabilities.jmgirard.com
Jeffrey Girard
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.