{"id":543,"date":"2018-07-14T09:56:07","date_gmt":"2018-07-14T00:56:07","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/how-to-draw-funnel-plot-for-checking-publication-bias\/"},"modified":"2024-10-14T08:40:50","modified_gmt":"2024-10-13T23:40:50","slug":"how-to-draw-funnel-plot-for-checking-publication-bias","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/how-to-draw-funnel-plot-for-checking-publication-bias\/","title":{"rendered":"R \u3067\u30e1\u30bf\u30a2\u30ca\u30ea\u30b7\u30b9\u306e\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u3092\u66f8\u304f\u65b9\u6cd5"},"content":{"rendered":"\n

\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u3068\u306f\u4f55\u304b\uff1f<\/p>\n\n\n\n

\u51fa\u7248\u30d0\u30a4\u30a2\u30b9\u3068\u306f\u4f55\u304b\uff1f<\/p>\n\n\n\n

\u30d5\u30a9\u30ec\u30b9\u30c8\u30d7\u30ed\u30c3\u30c8\u3068\u306e\u9055\u3044\u306f\u4f55\u304b\uff1f<\/p>\n\n\n\n\n\n\n\n

\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u306e\u524d\u306b\u51fa\u7248\u30d0\u30a4\u30a2\u30b9\u3068\u306f\uff1f<\/h2>\n\n\n\n

\u51fa\u7248\u30d0\u30a4\u30a2\u30b9\u3068\u306f\u3001\u516c\u8868\u3055\u308c\u3066\u3044\u308b\u8ad6\u6587\u306b\u504f\u308a\u304c\u3042\u308b\u3053\u3068\u3060\u3002<\/p>\n\n\n\n

\u7d71\u8a08\u5b66\u7684\u6709\u610f\u306b\u306a\u3063\u305f\u7d50\u679c\u304c\u516c\u8868\u3055\u308c\u3084\u3059\u3044\u3002<\/p>\n\n\n\n

\u307e\u305f\u3001\u901a\u8aac\u3092\u8986\u3059\u30bb\u30f3\u30bb\u30fc\u30b7\u30e7\u30ca\u30eb\u306a\u5185\u5bb9\u304c\u516c\u8868\u3055\u308c\u3084\u3059\u3044\u3002<\/p>\n\n\n\n

\u7d71\u8a08\u5b66\u7684\u6709\u610f\u3067\u306a\u304b\u3063\u305f\u7d50\u679c\u306f\u516c\u8868\u3055\u308c\u306b\u304f\u3044\u3002<\/p>\n\n\n\n

\u901a\u8aac\u901a\u308a\u3067\u8a71\u984c\u6027\u306b\u4e4f\u3057\u3044\u7d50\u679c\u306f\u516c\u8868\u3055\u308c\u306b\u304f\u3044\u3002<\/p>\n\n\n\n

\u3053\u3046\u3044\u3046\u504f\u308a\u304c\u8ad6\u6587\u516c\u8868\u306e\u4e16\u754c\u306b\u306f\u5b58\u5728\u3059\u308b\u3002<\/p>\n\n\n\n

\u30e1\u30bf\u30a2\u30ca\u30ea\u30b7\u30b9\u306b\u306f\u3001\u7279\u306b\u30b5\u30f3\u30d7\u30eb\u30b5\u30a4\u30ba\u304c\u5c0f\u3055\u304b\u3063\u305f\u305f\u3081\u306b\u7d71\u8a08\u5b66\u7684\u6709\u610f\u306b\u306a\u3089\u306a\u304b\u3063\u305f\u8a66\u9a13\u7d50\u679c\u304c\u91cd\u8981\u3067\u3042\u308b\u3002<\/p>\n\n\n\n

\u30b5\u30f3\u30d7\u30eb\u30b5\u30a4\u30ba\u304c\u5c0f\u3055\u304b\u3063\u305f\u305f\u3081\u306b\u7d50\u8ad6\u304c\u5f97\u3089\u308c\u306a\u304b\u3063\u305f\u547d\u984c\u306b\u7b54\u3048\u308b\u306e\u304c\u30e1\u30bf\u30a2\u30ca\u30ea\u30b7\u30b9\u306e\u4f7f\u547d\u306e\u4e00\u3064\u3067\u3042\u308b\u3002<\/p>\n\n\n\n

\u3053\u3046\u3044\u3046\u7406\u7531\u3067\u3001\u30e1\u30bf\u30a2\u30ca\u30ea\u30b7\u30b9\u3092\u5b9f\u65bd\u3059\u308b\u969b\u306b\u306f\u3001\u51fa\u7248\u30d0\u30a4\u30a2\u30b9\u306e\u7a0b\u5ea6\u3092\u78ba\u8a8d\u3057\u3066\u304a\u304f\u5fc5\u8981\u304c\u3042\u308b\u3002<\/p>\n\n\n\n

\u305d\u306e\u305f\u3081\u306e\u65b9\u6cd5\u306e\u4e00\u3064\u304c\u3001\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u3060\u3002<\/p>\n\n\n\n

\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u3068\u306f\uff1f\u3000\u51fa\u7248\u30d0\u30a4\u30a2\u30b9\u3092\u30c1\u30a7\u30c3\u30af\u3059\u308b<\/h2>\n\n\n\n

\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u306e\u524d\u306b\u30d5\u30a1\u30f3\u30cd\u30eb\u3068\u306f\u4f55\u304b\uff1f<\/p>\n\n\n\n

\u30d5\u30a1\u30f3\u30cd\u30eb\u3068\u306f\u6f0f\u6597\uff08\u308d\u3046\u3068\u3001\u3058\u3087\u3046\u3054\uff09\u306e\u3053\u3068\u3060\u3002<\/p>\n\n\n\n

\u6db2\u4f53\u3092\u53e3\u306e\u72ed\u3044\u74f6\u306a\u3069\u306b\u5165\u308c\u308b\u3068\u304d\u306b\u4f7f\u3046\u9006\u4e09\u89d2\u5f62\u306e\u5f62\u306e\u3082\u306e\u3060\u3002<\/p>\n\n\n\n

\u500b\u3005\u306e\u7814\u7a76\u7d50\u679c\u306e\u6563\u3089\u3070\u308a\u65b9\u304c\u3001\u6f0f\u6597\u3092\u88cf\u8fd4\u3057\u305f\u3088\u3046\u306a\u5f62\u306b\u306a\u308c\u3070\u7406\u60f3\u3067\u3001\u51fa\u7248\u30d0\u30a4\u30a2\u30b9\u306f\u5c0f\u3055\u3044\u3068\u5224\u65ad\u3067\u304d\u308b\u3002<\/p>\n\n\n\n

\uff1e\uff1e\u3082\u3046\u7d71\u8a08\u3067\u60a9\u3080\u306e\u306f\u7d42\u308f\u308a\u306b\u3057\u307e\u305b\u3093\u304b\uff1f\u00a0<\/a><\/strong><\/span><\/p>\r\n\"\"<\/a>\r\n

\u21911\u4e07\u4eba\u4ee5\u4e0a\u306e\u533b\u7642\u5f93\u4e8b\u8005\u304c\u8cfc\u8aad\u4e2d<\/span><\/strong><\/span><\/p><\/div>

R\u3067\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u3092\u66f8\u3044\u3066\u307f\u308b<\/h2>\n\n\n\n

R\u3067\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u3092\u66f8\u304f\u65b9\u6cd5\u3092\u89e3\u8aac\u3059\u308b\u3002<\/p>\n\n\n\n

metafor\u30d1\u30c3\u30b1\u30fc\u30b8\u306efunnel()\u3092\u4f7f\u3046\u3002<\/p>\n\n\n\n

\u307e\u305a\u6700\u521d\u306e\u4e00\u56de\u3060\u3051metafor\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3059\u308b\u3002<\/p>\n\n\n\n

install.packages<\/span>(<\/span>\"metafor\"<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
\u30c7\u30fc\u30bf\u306f\u4ee5\u4e0b\u306e\u901a\u308a\u3002<\/p>\n\n\n\n
a <-<\/span> c<\/span>(<\/span>3<\/span>,<\/span>7<\/span>,<\/span>5<\/span>,<\/span>102<\/span>,<\/span>28<\/span>,<\/span>4<\/span>,<\/span>98<\/span>,<\/span>60<\/span>,<\/span>25<\/span>,<\/span>138<\/span>,<\/span>64<\/span>,<\/span>45<\/span>,<\/span>9<\/span>,<\/span>57<\/span>,<\/span>25<\/span>,<\/span>65<\/span>,<\/span>17<\/span>)<\/span>\nn1 <-<\/span> c<\/span>(<\/span>38<\/span>,<\/span>114<\/span>,<\/span>69<\/span>,<\/span>1533<\/span>,<\/span>355<\/span>,<\/span>59<\/span>,<\/span>945<\/span>,<\/span>632<\/span>,<\/span>278<\/span>,<\/span>1916<\/span>,<\/span>873<\/span>,<\/span>263<\/span>,<\/span>291<\/span>,<\/span>858<\/span>,<\/span>154<\/span>,<\/span>1195<\/span>,<\/span>298<\/span>)<\/span>\nc <-<\/span> c<\/span>(<\/span>3<\/span>,<\/span>14<\/span>,<\/span>11<\/span>,<\/span>127<\/span>,<\/span>27<\/span>,<\/span>6<\/span>,<\/span>152<\/span>,<\/span>48<\/span>,<\/span>37<\/span>,<\/span>188<\/span>,<\/span>52<\/span>,<\/span>47<\/span>,<\/span>16<\/span>,<\/span>45<\/span>,<\/span>31<\/span>,<\/span>62<\/span>,<\/span>34<\/span>)<\/span>\nn0 <-<\/span> c<\/span>(<\/span>39<\/span>,<\/span>116<\/span>,<\/span>93<\/span>,<\/span>1520<\/span>,<\/span>365<\/span>,<\/span>52<\/span>,<\/span>939<\/span>,<\/span>471<\/span>,<\/span>282<\/span>,<\/span>1921<\/span>,<\/span>583<\/span>,<\/span>266<\/span>,<\/span>293<\/span>,<\/span>883<\/span>,<\/span>147<\/span>,<\/span>1200<\/span>,<\/span>309<\/span>)<\/span>\ndat <-<\/span> data.frame<\/span>(<\/span>a,<\/span>n1,<\/span>c,<\/span>n0)<\/span>\n<\/code><\/pre>\n\n\n\n
metafor\u30d1\u30c3\u30b1\u30fc\u30b8\u306eescalc()\u3092\u4f7f\u3046\u3002<\/p>\n\n\n\n
measure=\u3067\u6307\u6a19\u3092\u6307\u5b9a\u3002\u3053\u306e\u30c7\u30fc\u30bf\u306f\u30aa\u30c3\u30ba\u6bd4\u3002<\/p>\n\n\n\n
\u63a8\u5b9a\u5024yi \u3068\u5206\u6563vi\u304c\u8a08\u7b97\u3055\u308c\u308b\u3002<\/p>\n\n\n\n
rma.uni()\u3067\u3001method=”REML”\u3092\u6307\u5b9a\u3057\u3001\u5236\u9650\u4ed8\u304d\u6700\u5c24\u63a8\u5b9a\u91cf\uff08REML\uff09\u3067\u7d71\u5408\u30aa\u30c3\u30ba\u6bd4\u306895%\u4fe1\u983c\u533a\u9593\u3092\u6c42\u3081\u308b\u3002<\/p>\n\n\n\n
library<\/span>(<\/span>metafor)<\/span>\ndat.escalc <-<\/span> escalc<\/span>(<\/span>measure=<\/span>\"OR\"<\/span>,<\/span> ai=<\/span>a,<\/span> n1i=<\/span>n1,<\/span> ci=<\/span>c,<\/span> n2i=<\/span>n0,<\/span> data=<\/span>dat)<\/span>\nres.reml <-<\/span> rma.uni<\/span>(<\/span>yi,<\/span> vi,<\/span> method=<\/span>\"REML\"<\/span>,<\/span> data=<\/span>dat.escalc)<\/span>\n<\/code><\/pre>\n\n\n\n
\u305d\u3057\u3066\u3001\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u3092\u66f8\u304f\u3002<\/p>\n\n\n\n
funnel<\/span>(<\/span>res.reml)<\/span>\n<\/code><\/pre>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
X\u8ef8\u306b\u63a8\u5b9a\u5024\u3001Y\u8ef8\u306b\u6a19\u6e96\u8aa4\u5dee\u3092\u53d6\u3063\u3066\u3001\u30d7\u30ed\u30c3\u30c8\u3055\u308c\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
\u4e2d\u592e\u306e\u7dda\u306f\u7d71\u5408\u30aa\u30c3\u30ba\u6bd4\u306e0.79\u3002<\/p>\n\n\n\n
\u56f3\u5185\u306e\u53f3\u4e0b\u3001\u63a8\u5b9a\u5024\u3082\u6a19\u6e96\u8aa4\u5dee\u3082\u5927\u304d\u3081\u306e\u7814\u7a76\u304c\u306a\u3044\u3053\u3068\u304c\u308f\u304b\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u51fa\u7248\u30d0\u30a4\u30a2\u30b9\u306e\u5b58\u5728\u3092\u7591\u308f\u305b\u308b\u3002<\/p>\n\n\n\n
\u305d\u306e\u9650\u754c\u3092\u77e5\u3063\u305f\u3046\u3048\u3067\u89e3\u6790\u3092\u9032\u3081\u308b\u3068\u3044\u3046\u3053\u3068\u306b\u306a\u308b\u3002<\/p>\n\n\n\n
Discussion\u3067\u305d\u306e\u9650\u754c\u306b\u89e6\u308c\u308b\u307b\u3046\u304c\u3088\u3044\u3060\u308d\u3046\u3002<\/p>\n\n\n\n
\u5b9f\u969b\u306e\u3068\u3053\u308d\u3001\u5bfe\u6570\u8868\u793a\u306f\u308f\u304b\u308a\u306b\u304f\u3044\u3002<\/p>\n\n\n\n
X\u8ef8\u3092\u5bfe\u6570\u76ee\u76db\u3067\u771f\u6570\u30e9\u30d9\u30eb\u306b\u3059\u308b\u3002<\/p>\n\n\n\n
\u3053\u308c\u3060\u3068\u308f\u304b\u308a\u3084\u3059\u3044\u3002<\/p>\n\n\n\n
funnel<\/span>(<\/span>res.reml,<\/span> atransf=<\/span>exp,<\/span> at=<\/span>log<\/span>(<\/span>c<\/span>(<\/span>0.25<\/span>,<\/span> 0.5<\/span>,<\/span> 1<\/span>,<\/span> 2<\/span>,<\/span> 4<\/span>)))<\/span>\n<\/code><\/pre>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u3068\u30d5\u30a9\u30ec\u30b9\u30c8\u30d7\u30ed\u30c3\u30c8\u306e\u9055\u3044<\/h2>\n\n\n\n
\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u3068\u30d5\u30a9\u30ec\u30b9\u30c8\u30d7\u30ed\u30c3\u30c8\u306f\u3001\u3068\u3082\u306b\u30e1\u30bf\u30a2\u30ca\u30ea\u30b7\u30b9\u3067\u767b\u5834\u3059\u308b\u304c\u3001\u5168\u304f\u306e\u5225\u7269\u3060\u3002<\/p>\n\n\n\n
\u30d5\u30a9\u30ec\u30b9\u30c8\u30d7\u30ed\u30c3\u30c8\u306f\u30e1\u30bf\u30a2\u30ca\u30ea\u30b7\u30b9\u306e\u7d71\u5408\u7d50\u679c\u306e\u4e3b\u8981\u306a\u30b0\u30e9\u30d5\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u3053\u306e\u3069\u3053\u304c\u30d5\u30a9\u30ec\u30b9\u30c8\uff08\u68ee\uff09\u306a\u306e\u304b\u3068\u3044\u3046\u3068\u3001\u5168\u4f53\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\u4e00\u3064\u4e00\u3064\u306e\u7814\u7a76\u3092\u6728\u306b\u898b\u7acb\u3066\u3066\u3001\u5168\u4f53\u304c\u68ee\uff08\u30d5\u30a9\u30ec\u30b9\u30c8\uff09\u3068\u3044\u3046\u3053\u3068\u3060\u3002<\/p>\n\n\n\n
\u75ab\u5b66\u306e\u6388\u696d\u3067\u6700\u3082\u5370\u8c61\u6df1\u3044\u8a00\u8449\u3092\u601d\u3044\u51fa\u3059\u3002<\/p>\n\n\n\n
\u300c\u75ab\u5b66\u306f\u68ee\u3092\u898b\u3066\u6728\u3092\u77e5\u308b\u5b66\u554f\u3067\u3042\u308b\u300d<\/p>\n\n\n\n
\u3053\u306e\u6642\u306e\u68ee\u306f\u4eba\u3005\u306e\u96c6\u56e3\u3001\u6728\u306f\u500b\u3005\u4eba\u30fb\u4e00\u4eba\u3072\u3068\u308a\u306e\u60a3\u8005\u3055\u3093\u306e\u3053\u3068\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\u3053\u3093\u306a\u3068\u3053\u308d\u3067\u3082\u3064\u306a\u304c\u3063\u3066\u3044\u305f\u304b\u3068\u601d\u3046\u3002<\/p>\n\n\n\n
\u3068\u3044\u3046\u3053\u3068\u3067\u3001\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u306f\u51fa\u7248\u30d0\u30a4\u30a2\u30b9\u3092\u30c1\u30a7\u30c3\u30af\u3059\u308b\u30b0\u30e9\u30d5\u3067\u3001\u30d5\u30a9\u30ec\u30b9\u30c8\u30d7\u30ed\u30c3\u30c8\u306f\u30e1\u30bf\u30a2\u30ca\u30ea\u30b7\u30b9\u306e\u500b\u3005\u306e\u8a66\u9a13\u53ca\u3073\u7d71\u5408\u7d50\u679c\u3092\u8868\u793a\u3059\u308b\u30b0\u30e9\u30d5\u3068\u3044\u3046\u5168\u304f\u5225\u306e\u3082\u306e\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\u307e\u3068\u3081<\/h2>\n\n\n\n
\u51fa\u7248\u30d0\u30a4\u30a2\u30b9\u3092\u30c1\u30a7\u30c3\u30af\u3059\u308b\u30d5\u30a1\u30f3\u30cd\u30eb\u30d7\u30ed\u30c3\u30c8\u3092\u89e3\u8aac\u3057\u305f\u3002<\/p>\n\n\n\n
R\u3067\u66f8\u304f\u65b9\u6cd5\u3082\u7d39\u4ecb\u3057\u305f\u3002<\/p>\n\n\n\n
\u30d5\u30a9\u30ec\u30b9\u30c8\u30d7\u30ed\u30c3\u30c8\u3068\u306e\u9055\u3044\u3082\u8aac\u660e\u3057\u305f\u3002<\/p>\n\n\n\n
\u4f55\u3089\u304b\u53c2\u8003\u306b\u306a\u308c\u3070\u5e78\u3044\u3002<\/p>\n\n\n\n
\u53c2\u8003\u66f8\u7c4d<\/h2>\n\n\n\n
\"\u30e1\u30bf\u30fb\u30a2\u30ca\u30ea\u30b7\u30b9\u5165\u9580\u2015\u30a8\u30d3\u30c7\u30f3\u30b9\u306e\u7d71\u5408\u3092\u3081\u3056\u3059\u7d71\u8a08\u624b\u6cd5<\/a><\/figure>\n\n\n\n
<\/div>\n\n\n\n
\n
\u30e1\u30bf\u30fb\u30a2\u30ca\u30ea\u30b7\u30b9\u5165\u9580\u2015\u30a8\u30d3\u30c7\u30f3\u30b9\u306e\u7d71\u5408\u3092\u3081\u3056\u3059\u7d71\u8a08\u624b\u6cd5 (\u533b\u5b66\u7d71\u8a08\u5b66\u30b7\u30ea\u30fc\u30ba)<\/a><\/p>\n