안녕하세요
질문할게 자꾸 생깁니다. 초보자라 많은 선처 바랍니다.
R프로그램을 돌렸을때 funnel plot과 웹에서 하는 R 프로그램 funnel plot에 차이가 있습니다.
가로축 효과크기 눈금 수치가 서로 많이 달라 별도로 입력해야 값이 있는지 모르겠습니다.
예를 들면 r 값이 0. 5~0.9 으로 이루어져 있으면 0. 5와 1사이에 깔데기 모양이 있어야 하는데
여기서는 4-5사이에 있어요.
첨부물 올렸습니다.
안녕하세요
질문할게 자꾸 생깁니다. 초보자라 많은 선처 바랍니다.
R프로그램을 돌렸을때 funnel plot과 웹에서 하는 R 프로그램 funnel plot에 차이가 있습니다.
가로축 효과크기 눈금 수치가 서로 많이 달라 별도로 입력해야 값이 있는지 모르겠습니다.
예를 들면 r 값이 0. 5~0.9 으로 이루어져 있으면 0. 5와 1사이에 깔데기 모양이 있어야 하는데
여기서는 4-5사이에 있어요.
첨부물 올렸습니다.
안녕하세요 제 분야는 구조방정식쪽이라... meta분석에 대해 잘 모르지만... 학술적 자리가 아니기에 편하게 몇자 적어봅니다
일단 메타분석시 제일 먼저 확인해야 하는 것은 이질성 검증입니다
meta분석에서 이질성 검증을 위한 여러가지 통계적 방법이 제시되고 있지만, I^2 데이터 값이 98.6%입니다
25%미만이 이질성이 낮은 수준이고, 25~75%를 중간수준정도로 보고 있습니다
그리고 75% 이상을 높은 수준의 이질성을 보인다고 볼 수 있습니다
그런데 현재 98.6%이면 메타분석 자체가 불가하다고 생각해 볼 수 있습니다
물론 어느정도의 이질성은 변량효과모형으로 조정이 가능하지만..
육안으로 봐도 효과크기 및 특히 N의 차이가 심하게 보입니다
케이스중 4, 8, 13번이 유일하게 깔데 그래프 안으로 들어와있습니다 그리고 16번은 특이할 정도로 밖으로 벗어나 있습니다
16번 사례를 지우면 이전보다 조금 나은 그래프가 되는 걸 알 수 있습니다
그럼에도 90%가 넘는 이질성을 보이기에 사례들이 이질함을 갖는 원인을 찾아보시는게 맞습니다
그리고 원인이되는 그룹변수를 만들어 조절변수로 넣어주시면 될 것 같습니다
다만 조절변수로 metaANOVA를 확인할 경우 한 그룹당 최소 10개의 케이스가 있어야 합니다
추가적으로 첨부해 주신 그래프는 옵션이 좀 다르게 들어간 듯 합니다
웹메타에서 summary measure 를 COR가 아니라 ZCOR로 입력하시거나 4. Assign Variables에 other options에 backtransf=F를 입력하시면 Fisher's Z로 전환한 상관계수로 나오게 됩니다
그리고 6. Select Plots에서 show ALL study labels을 체크하시면 각 플랏에 케이스 번호가 나옵니다 필요하신 옵션으로 생각되어 말씀드립니다
서두에 말씀 드린것 처럼 저는 메타분석은 잘 모릅니다
저 또한 메타분석을 이제 막 공부하기 시작한 입장이라는 점을 말씀드립니다