Bối cảnh:
Tôi muốn tiến hành hồi quy meta bằng các nghiên cứu có (1) một số kết quả / cấu trúc (= đa biến) và (2) nhiều kích cỡ hiệu ứng cho mỗi kết quả này vì các biện pháp khác nhau. Đây là một kế hoạch hy vọng giải thích nó tốt nhất:
- Nghiên cứu 1, Kết quả A, Hiệu ứng kích thước 1
- Nghiên cứu 1, Kết quả A, Hiệu ứng kích thước 2
- Nghiên cứu 1, Kết quả B, Kích thước hiệu ứng 3
- Nghiên cứu 2, Kết quả A, Hiệu ứng kích thước 4
- Nghiên cứu 2, Kết quả C, Hiệu ứng kích thước 5
- Nghiên cứu 2, Kết quả C, Hiệu ứng kích thước 6
- ...
Các nghiên cứu so sánh phương tiện của hai nhóm về kết quả và kích cỡ hiệu ứng khác nhau là Hedge g.
Một ví dụ thực tế sẽ là "Bộ nhớ làm việc" có thể được chia thành các kết quả khác nhau (Baddeley, 1974), ví dụ "Vòng lặp âm vị học", "Bảng vẽ phác thảo trực quan" hoặc "Điều hành trung tâm".
Ví dụ: Nghiên cứu 1 đánh giá "Vòng lặp âm vị" (Kết quả A) với hai biện pháp khác nhau (= Kích thước hiệu ứng 1 và 2) và "Điều hành trung tâm" (Kết quả B) với một biện pháp (= Kích thước hiệu ứng 3).
Vấn đề:
Một cách tiếp cận đa biến thích hợp đòi hỏi phải biết mọi mối tương quan giữa kích thước hiệu ứng và kết quả để ước tính hiệp phương sai. Tuy nhiên tôi không biết (1) mối tương quan giữa các kích cỡ hiệu ứng khác nhau trong cùng một nghiên cứu và (2) mối tương quan giữa các kết quả của các nghiên cứu khác nhau. Tôi có thể ước tính chúng hoặc cố gắng tìm ít nhất một vài mối tương quan để làm việc, nhưng điều đó có nghĩa là rất nhiều tìm kiếm tài liệu bổ sung mà tôi muốn tránh.
Giải pháp (cho đến nay):
Tôi đã bắt gặp một vài phương pháp xử lý các vấn đề tương tự.
Ước tính phương sai mạnh mẽ (Hedges, 2010) là một cách tiếp cận hay để đối phó với nhiều kích cỡ hiệu ứng. Tuy nhiên tôi vẫn phải đoán một mối tương quan và thực hiện phân tích độ nhạy và dường như không thể so sánh một số kết quả (nghĩa là chỉ có hồi quy meta đơn biến).
Cách tiếp cận đa cấp của Van den Noorgate (2014) đầy hứa hẹn vì không cần ước tính bất kỳ mối tương quan nào bằng cách cho phép thay đổi giữa kích thước hiệu ứng và giữa kích thước hiệu ứng trong các nghiên cứu. Một phân tích meta đa biến đa cấp (= kết quả khác nhau và nhiều kích cỡ hiệu ứng như trong sơ đồ trên) và hồi quy meta đa biến đơn biến (= nhiều kích thước hiệu ứng nhưng không có sự khác biệt giữa các kết quả) được mô tả.
Sử dụng gói metafor trong R Tôi tự hỏi liệu tôi có thể kết hợp cả hai cách tiếp cận đa cấp và thực hiện hồi quy meta đa biến đa cấp hay không. Ví dụ về phân tích meta đa cấp và hồi quy meta đa biến bằng metafor được đưa ra ở đây http://www.metafor-project.org/doku.php/analyses:konstantopoulos2011 (đa cấp) và tại đây http: //www.metafor- project.org/doku.php/analyses:berkey1998 (đa biến). (Xin lưu ý rằng ví dụ đa cấp được liên kết ở trên thực sự mô tả một cách tiếp cận để đối phó với các phụ thuộc phân cấp (ví dụ: các nghiên cứu được thực hiện bởi cùng một phòng thí nghiệm nghiên cứu). Thay vào đó tôi sử dụng phương pháp đa cấp được mô tả bởi Van den Noorgate.)
Các biến:
ES: Kích thước hiệu ứng (Hedge g)
VI: Phương sai kích thước hiệu ứng
Pub_Year: Năm xuất bản với tư cách là công cụ dự đoán trong hồi quy meta
ES_ID: Mỗi kích thước hiệu ứng đều có một ID duy nhất, bất kể chúng thuộc về nghiên cứu hay kết quả nào.
Outcome_ID: Các kết quả tương tự có cùng một ID (ví dụ: "Vòng lặp âm vị" = 1, "Điều hành trung tâm" = 2), bất kể chúng thuộc về nghiên cứu nào.
Research_ID: Kích thước hiệu ứng của cùng một nghiên cứu có cùng ID (ví dụ: Kích thước hiệu ứng của nghiên cứu 1 = 1, Kích thước hiệu ứng của nghiên cứu 2 = 2), bất kể kết quả của chúng thuộc về kết quả nào.
Mã R trong metafor cho phân tích meta đa biến đa cấp:
rma.mv (ES, VI, mods = ~ Outcome_ID -1, Random = list (~ 1 | Study_ID, ~ 1 | ES_ID), data = data.set)
- mods = ~ Outcome_ID -1 gọi cho cách tiếp cận đa biến và liệt kê kích thước hiệu ứng trung bình cho mọi kết quả.
- Random = list (~ 1 | Study_ID, ~ 1 | ES_ID) là cách tiếp cận đa cấp được mô tả bởi Van den Noorgate. Nó cho phép thay đổi ngẫu nhiên giữa các kích thước hiệu ứng trong các nghiên cứu (~ 1 | Study_ID) và giữa các kích thước hiệu ứng (~ 1 | ES_ID). Bạn cũng có thể tiến hành phân tích này bằng gói metaSEM. Kết quả giống hệt nhau.
Mã R trong metafor cho hồi quy meta đa biến đa cấp:
rma.mv (ES, VI, mods = ~ Outcome_ID + Kết quả: Tôi (Pub_Year-mean (Pub_Year)) -1, Random = list (~ 1 | Study_ID, ~ 1 | ES_ID), dữ liệu = data.set)
- mods = ~ Outcome_ID + Kết quả: Tôi (Pub_Year-mean (Pub_Year)) -1 bây giờ yêu cầu hồi quy meta đa biến với Năm xuất bản tập trung quanh ý nghĩa như một công cụ dự đoán.
Sử dụng tùy chọn profile () trong metafor, Sơ đồ khả năng hồ sơ trông ổn. Tuy nhiên, tôi vẫn tự hỏi nếu tôi không định lượng quá mức mô hình hoặc nếu có gì đó sai khi kết hợp các đối số mod và ngẫu nhiên theo cách này.
Mong ý kiến của bạn, đề xuất, ý tưởng, cách tiếp cận khác, mọi thứ ;-) Cảm ơn!
Cập nhật, trả lời câu trả lời của Wolfgang:
Trước hết: Cảm ơn rất nhiều về câu trả lời chi tiết của bạn và các liên kết bổ sung mà bạn đã cung cấp. Tôi không biết về danh sách gửi thư mô hình hỗn hợp R-sig. Vì vậy, cảm ơn! Tôi đánh giá cao điều đó rất nhiều.
Hãy để tôi thử tổng hợp mọi thứ và điều chỉnh nó phù hợp với tình huống của tôi để xem tôi có hiểu mọi thứ ngay tại đây không. Tôi có thể làm những điều sau đây:
Có được mối tương quan: Thật không may là mối tương quan không được báo cáo. Ban đầu, phân tích tổng hợp bao gồm hơn 50 nghiên cứu. Gần một nửa các nghiên cứu đã mất hoặc không báo cáo dữ liệu. Mọi tác giả của những nghiên cứu này đều được liên lạc và tôi nhận được 4 phản hồi trong số 26 yêu cầu (sau 2 tháng chờ đợi). Nhưng đó là một vấn đề báo cáo chung không được thảo luận ở đây.
Nếu tôi dự đoán sơ bộ về tất cả các mối tương quan tôi có thể:
Tiến hành phân tích tổng hợp đa biến và hồi quy meta như trong Berkey et al. (1998) ví dụ và làm một phân tích độ nhạy.
Sử dụng mô hình phân tích tổng hợp đa biến được trang bị này và làm việc với hàm mạnh (). Tuy nhiên, không có hồi quy meta dựa trên hàm mạnh () dường như có thể có trong metafor. Và hàm mạnh mẽ () được mô tả trong blog của James Pustejovsky chỉ hoạt động với các hồi quy meta đơn biến. Vì vậy, nếu tôi hiểu đúng, các ước tính của hàm mạnh () ít nhiều sẽ xác nhận các ước tính của mô hình đã được trang bị của tôi (?).
Trực tiếp đi cho các phương pháp mạnh mẽ và sử dụng gói robumeta. Tuy nhiên không có phân tích tổng hợp đa biến nào là có thể. Tôi tìm thấy một mã SASđể xử lý vấn đề này. Nhưng mã được phát triển 3 năm trước và dường như nó chưa bao giờ thực sự được thảo luận. Cuối cùng, khi sử dụng robumeta, tôi phải tóm tắt rất nhiều kết quả khác nhau thành một phân tích tổng hợp lớn hoặc tôi phải tiến hành một số phân tích tổng hợp đơn biến cho mỗi kết quả mà tôi muốn tránh.Nếu tôi không muốn đoán bất kỳ mối tương quan nào, tôi có thể đi theo phương pháp đa cấp như được mô tả bởi Van den Noorgate bằng metafor, metaSEM hoặc SAS. Tuy nhiên, có một số hạn chế khi sử dụng phương pháp này so với cách tiếp cận đa biến dựa trên các mối tương quan. Ngoài ra, tôi không chắc liệu có thể hồi quy meta đa biến đa cấp hay không. Gói metaSEM chỉ mô tả phân tích meta đa biến đa cấp hoặc hồi quy meta đa biến đơn biến.
Thật không may, tôi không quen với việc sử dụng các phương pháp lấy mẫu lại trong phân tích tổng hợp. Tôi đã nghiên cứu các ví dụ của bạn nhưng tôi không chắc nó có thể giúp tôi giải quyết vấn đề tương quan / đa biến thế như thế nào. Ý bạn là tôi nên cố gắng ước tính các mối tương quan bằng cách sử dụng bootstrapping? Và nếu vậy, tôi không chắc giá trị nào sẽ tương quan vì số lượng phương tiện hoặc kích cỡ hiệu ứng trong và giữa các nghiên cứu khác nhau.
Việc đơn giản hóa mô hình được mô tả bởi Riley và các đồng nghiệp nghe có vẻ thú vị. Tôi ghi nhớ nó, mặc dù tôi muốn làm việc với một trong những phương pháp được mô tả ở trên.
robust()
trên các rma.mv
mô hình mà bạn đã viết trong câu hỏi của bạn (với biến Study_ID
là phân cụm). Vì các mô hình này giả định các lỗi lấy mẫu độc lập, chúng bị sai. Cách tiếp cận mạnh mẽ sẽ cho bạn kết quả phù hợp ngay cả khi phần ma trận phương sai hiệp phương sai của mô hình bị sai.