{"id":456,"date":"2019-01-04T12:50:39","date_gmt":"2019-01-04T03:50:39","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/meta-analysis-of-correlation-coefficients-in-r\/"},"modified":"2024-10-12T17:20:15","modified_gmt":"2024-10-12T08:20:15","slug":"meta-analysis-of-correlation-coefficients-in-r","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/meta-analysis-of-correlation-coefficients-in-r\/","title":{"rendered":"R \u3067\u76f8\u95a2\u4fc2\u6570\u306e\u30e1\u30bf\u30a2\u30ca\u30ea\u30b7\u30b9\u3092\u884c\u3046\u65b9\u6cd5"},"content":{"rendered":"\n<p>\u76f8\u95a2\u4fc2\u6570\u3092\u7d71\u5408\u3057\u305f\u3044\u5834\u5408\u306f\u3069\u3046\u3084\u308b\u304b\uff1f<\/p>\n\n\n\n<p>R \u3067\u306e\u65b9\u6cd5\u3002<\/p>\n\n\n\n<!--more-->\n\n\n\n<h2 class=\"wp-block-heading\" id=\"\u500b\u3005\u306e\u7814\u7a76\u306e\u76f8\u95a2\u4fc2\u6570\u306895\u4fe1\u983c\u533a\u9593\u306e\u6e96\u5099\">\u500b\u3005\u306e\u7814\u7a76\u306e\u76f8\u95a2\u4fc2\u6570\u306895%\u4fe1\u983c\u533a\u9593\u306e\u6e96\u5099<\/h2>\n\n\n\n<p>\u4f7f\u3046\u30c7\u30fc\u30bf\u306f\u4ee5\u4e0b\u306e\u901a\u308a\u3002<\/p>\n\n\n\n<p>r \u304c\u76f8\u95a2\u4fc2\u6570\u3002<\/p>\n\n\n\n<p>n \u304c\u30b5\u30f3\u30d7\u30eb\u30b5\u30a4\u30ba\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code>r <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">0.307<\/span><span class=\"synSpecial\">,<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">0.01<\/span><span class=\"synSpecial\">,<\/span><span class=\"synConstant\">0.300<\/span><span class=\"synSpecial\">,<\/span><span class=\"synConstant\">0.119<\/span><span class=\"synSpecial\">,<\/span><span class=\"synConstant\">0.194<\/span><span class=\"synSpecial\">,<\/span><span class=\"synConstant\">0.248<\/span><span class=\"synSpecial\">)<\/span>\nn <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">107<\/span><span class=\"synSpecial\">,<\/span><span class=\"synConstant\">1524<\/span><span class=\"synSpecial\">,<\/span><span class=\"synConstant\">154<\/span><span class=\"synSpecial\">,<\/span><span class=\"synConstant\">6165<\/span><span class=\"synSpecial\">,<\/span><span class=\"synConstant\">4138<\/span><span class=\"synSpecial\">,<\/span><span class=\"synConstant\">1559<\/span><span class=\"synSpecial\">)<\/span>\n<\/code><\/pre>\n\n\n\n<p>z \u304cFisher\u306e\u5206\u6563\u5b89\u5b9a\u5316\u5909\u63db\u5f8c\u306e\u76f8\u95a2\u4fc2\u6570\u3002<\/p>\n\n\n\n<p>Hyperbolic arc-tangent \u95a2\u6570 atanh() \u3092\u4f7f\u3046\u3068\u7c21\u5358\u306b\u5909\u63db\u3067\u304d\u308b\u3002<\/p>\n\n\n\n<p>z \u306b\u5909\u63db\u3057\u305f\u5f8c\u306e\u6a19\u6e96\u8aa4\u5dee\u306f $ \\sqrt{\\frac{1}{n &#8211; 3}} $ \u306b\u306a\u308b\u3002<\/p>\n\n\n\n<p>zl, zu \u306f z \u306e95%\u4fe1\u983c\u533a\u9593\u4e0b\u9650\u3068\u4e0a\u9650\u3002<\/p>\n\n\n\n<p>rl, ru \u306f zl, zu\u3092\u9006\u5909\u63db\u3057\u305f\uff08\u5143\u306e\u4e16\u754c\u306b\u623b\u3057\u305f\uff09\u5024\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code>z <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">atanh<\/span><span class=\"synSpecial\">(<\/span>r<span class=\"synSpecial\">)<\/span>\nse <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sqrt<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">1<\/span><span class=\"synStatement\">\/<\/span><span class=\"synSpecial\">(<\/span>n<span class=\"synStatement\">-<\/span><span class=\"synConstant\">3<\/span><span class=\"synSpecial\">))<\/span>\nZa <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">qnorm<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">1<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">0.05<\/span><span class=\"synStatement\">\/<\/span><span class=\"synConstant\">2<\/span><span class=\"synSpecial\">)<\/span>\nzl <span class=\"synStatement\">&lt;-<\/span> z<span class=\"synStatement\">-<\/span>Za<span class=\"synStatement\">*<\/span>se\nzu <span class=\"synStatement\">&lt;-<\/span> z<span class=\"synStatement\">+<\/span>Za<span class=\"synStatement\">*<\/span>se\nrl <span class=\"synStatement\">&lt;-<\/span> <span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">exp<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>zl<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">)<\/span><span class=\"synStatement\">\/<\/span><span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">exp<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>zl<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">+<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">)<\/span>\nru <span class=\"synStatement\">&lt;-<\/span> <span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">exp<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>zu<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">)<\/span><span class=\"synStatement\">\/<\/span><span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">exp<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>zu<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">+<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">)<\/span>\n<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"Fisher\u306e\u5206\u6563\u5b89\u5b9a\u5316\u5909\u63db\u306e\u9006\u5909\u63db\u306ftanh\">Fisher\u306e\u5206\u6563\u5b89\u5b9a\u5316\u5909\u63db\u306e\u9006\u5909\u63db\u306ftanh()<\/h2>\n\n\n\n<p>Hyperbolic arc-tangent\u306e\u9006\u306f\u3001hyperbolic tangent\u3002<\/p>\n\n\n\n<p>\u306a\u306e\u3067\u3001rl, ru \u306e\u8a08\u7b97\u5f0f\u306ftanh()\u3067\u7f6e\u304d\u63db\u3048\u3089\u308c\u308b\u3002<\/p>\n\n\n\n<p>\u4ee5\u964d\u306ftanh()\u3092\u4f7f\u3046\u3002<\/p>\n\n\n\n<p>\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u3001\u30b7\u30f3\u30d7\u30eb\u306a\u30b9\u30af\u30ea\u30d7\u30c8\u3067\u66f8\u3051\u308b\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synStatement\">&gt;<\/span> <span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">exp<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>zl<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">)<\/span><span class=\"synStatement\">\/<\/span><span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">exp<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>zl<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">+<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synSpecial\">&#91;<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">]<\/span>  <span class=\"synConstant\">0.12439197<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">0.06018300<\/span>  <span class=\"synConstant\">0.14890459<\/span>  <span class=\"synConstant\">0.09431717<\/span>  <span class=\"synConstant\">0.16450217<\/span>  <span class=\"synConstant\">0.20082642<\/span>\n<span class=\"synStatement\">&gt;<\/span>  <span class=\"synIdentifier\">tanh<\/span><span class=\"synSpecial\">(<\/span>zl<span class=\"synSpecial\">)<\/span>\n<span class=\"synSpecial\">&#91;<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">]<\/span>  <span class=\"synConstant\">0.12439197<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">0.06018300<\/span>  <span class=\"synConstant\">0.14890459<\/span>  <span class=\"synConstant\">0.09431717<\/span>  <span class=\"synConstant\">0.16450217<\/span>  <span class=\"synConstant\">0.20082642<\/span>\n<\/code><\/pre>\n\n\n\n<div id=\"biost-1417229356\" class=\"biost- biost-entity-placement\"><p style=\"text-align: center;\"><span style=\"font-size: 20px;\"><strong><a href=\"https:\/\/best-biostatistics.com\/kmhl\">\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 href=\"https:\/\/best-biostatistics.com\/kmhl\"><img class=\"aligncenter wp-image-2794 size-full\" src=\"https:\/\/best-biostatistics.com\/wp\/wp-content\/uploads\/2023\/11\/bn_r_03.png\" alt=\"\" width=\"500\" height=\"327\" \/><\/a>\r\n<p style=\"text-align: center;\"><span style=\"color: #ff0000; font-size: 20px;\"><strong><span class=\"marker2\">\u21911\u4e07\u4eba\u4ee5\u4e0a\u306e\u533b\u7642\u5f93\u4e8b\u8005\u304c\u8cfc\u8aad\u4e2d<\/span><\/strong><\/span><\/p><\/div><h2 class=\"wp-block-heading\" id=\"\u500b\u3005\u306e\u7814\u7a76\u306e\u76f8\u95a2\u4fc2\u6570\u306895\u4fe1\u983c\u533a\u9593\u3092\u30b0\u30e9\u30d5\u306b\u3059\u308b\">\u500b\u3005\u306e\u7814\u7a76\u306e\u76f8\u95a2\u4fc2\u6570\u306895%\u4fe1\u983c\u533a\u9593\u3092\u30b0\u30e9\u30d5\u306b\u3059\u308b<\/h2>\n\n\n\n<p>\u500b\u3005\u306e\u7814\u7a76\u3092\u30b0\u30e9\u30d5\u306b\u3059\u308b\u306b\u306f\u4ee5\u4e0b\u306e\u901a\u308a\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synComment\"># ----- individual graph -----<\/span>\nk <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">length<\/span><span class=\"synSpecial\">(<\/span>n<span class=\"synSpecial\">)<\/span>\nid <span class=\"synStatement\">&lt;-<\/span> k<span class=\"synSpecial\">:<\/span><span class=\"synConstant\">1<\/span>\n<span class=\"synIdentifier\">plot<\/span><span class=\"synSpecial\">(<\/span>r<span class=\"synSpecial\">,<\/span> id<span class=\"synSpecial\">,<\/span> ylim<span class=\"synStatement\">=<\/span><span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">4<\/span><span class=\"synSpecial\">,<\/span> k<span class=\"synStatement\">+<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">),<\/span> pch<span class=\"synStatement\">=<\/span><span class=\"synConstant\">\" \"<\/span><span class=\"synSpecial\">,<\/span> xlim<span class=\"synStatement\">=<\/span><span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">,<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">),<\/span> yaxt<span class=\"synStatement\">=<\/span><span class=\"synConstant\">\"n\"<\/span><span class=\"synSpecial\">,<\/span>\nylab<span class=\"synStatement\">=<\/span><span class=\"synConstant\">\"Citation\"<\/span><span class=\"synSpecial\">,<\/span> xlab<span class=\"synStatement\">=<\/span><span class=\"synConstant\">\"Correlation coefficients\"<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">title<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">\" Correlation coefficient model \"<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">symbols<\/span><span class=\"synSpecial\">(<\/span>r<span class=\"synSpecial\">,<\/span> id<span class=\"synSpecial\">,<\/span> squares<span class=\"synStatement\">=<\/span><span class=\"synIdentifier\">sqrt<\/span><span class=\"synSpecial\">(<\/span>n<span class=\"synStatement\">-<\/span><span class=\"synConstant\">3<\/span><span class=\"synSpecial\">),<\/span> add<span class=\"synStatement\">=<\/span><span class=\"synConstant\">T<\/span><span class=\"synSpecial\">,<\/span> inches<span class=\"synStatement\">=<\/span><span class=\"synConstant\">0.25<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">abline<\/span><span class=\"synSpecial\">(<\/span>v<span class=\"synStatement\">=<\/span><span class=\"synConstant\">0<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synStatement\">for<\/span> <span class=\"synSpecial\">(<\/span>i <span class=\"synStatement\">in<\/span> <span class=\"synConstant\">1<\/span><span class=\"synSpecial\">:<\/span>k<span class=\"synSpecial\">){<\/span>\nj <span class=\"synStatement\">&lt;-<\/span> k<span class=\"synStatement\">-<\/span>i<span class=\"synStatement\">+<\/span><span class=\"synConstant\">1<\/span>\nx <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span>rl<span class=\"synSpecial\">&#91;<\/span>i<span class=\"synSpecial\">],<\/span> ru<span class=\"synSpecial\">&#91;<\/span>i<span class=\"synSpecial\">])<\/span>\ny <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span>j<span class=\"synSpecial\">,<\/span> j<span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">lines<\/span><span class=\"synSpecial\">(<\/span>x<span class=\"synSpecial\">,<\/span> y<span class=\"synSpecial\">,<\/span> type<span class=\"synStatement\">=<\/span><span class=\"synConstant\">\"l\"<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">text<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">0.5<\/span><span class=\"synSpecial\">,<\/span> i<span class=\"synSpecial\">,<\/span> j<span class=\"synSpecial\">)<\/span>\n<span class=\"synSpecial\">}<\/span>\n<\/code><\/pre>\n\n\n\n<p>\u30b0\u30e9\u30d5\u306f\u3053\u3093\u306a\u98a8\u306b\u51fa\u529b\u3055\u308c\u308b\u3002<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img decoding=\"async\" width=\"553\" height=\"552\" src=\"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104121705.png\" alt=\"\" class=\"wp-image-2687\" title=\"\" srcset=\"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104121705.png 553w, https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104121705-300x300.png 300w, https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104121705-150x150.png 150w\" sizes=\"(max-width: 553px) 100vw, 553px\" \/><\/figure>\n\n\n\n<p><span itemscope=\"\" itemtype=\"http:\/\/schema.org\/Photograph\"><\/span><\/p>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"\u56fa\u5b9a\u52b9\u679c\u3067\u76f8\u95a2\u4fc2\u6570\u3092\u7d71\u5408\u3059\u308b\">\u56fa\u5b9a\u52b9\u679c\u3067\u76f8\u95a2\u4fc2\u6570\u3092\u7d71\u5408\u3059\u308b<\/h2>\n\n\n\n<p>$ w_i = n_i &#8211; 3 $ \u3067\u91cd\u307f\u3065\u3051\u5e73\u5747\u3092\u53d6\u308b\u3068\u7d71\u5408\u76f8\u95a2\u4fc2\u6570\u306e\u63a8\u5b9a\u5024\uff08Fisher\u306e\u5206\u6563\u5b89\u5b9a\u5316\u5909\u63db\u5f8c\uff09\u3002<\/p>\n\n\n\n<p>tanh()\u3067\u623b\u3059\u3068\u6c42\u3081\u308b\u70b9\u63a8\u5b9a\u5024\u306b\u306a\u308b\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synComment\"># ----- fixed effects -----<\/span>\nw <span class=\"synStatement\">&lt;-<\/span> n<span class=\"synStatement\">-<\/span><span class=\"synConstant\">3<\/span>\nsw <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>w<span class=\"synSpecial\">)<\/span>\nvarz <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>w<span class=\"synStatement\">*<\/span>z<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">\/<\/span>sw\nvarr <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">tanh<\/span><span class=\"synSpecial\">(<\/span>varz<span class=\"synSpecial\">)<\/span>\nvarrl <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">tanh<\/span><span class=\"synSpecial\">(<\/span>varz<span class=\"synStatement\">-<\/span>Za<span class=\"synStatement\">*<\/span><span class=\"synIdentifier\">sqrt<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">1<\/span><span class=\"synStatement\">\/<\/span>sw<span class=\"synSpecial\">))<\/span>\nvarru <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">tanh<\/span><span class=\"synSpecial\">(<\/span>varz<span class=\"synStatement\">+<\/span>Za<span class=\"synStatement\">*<\/span><span class=\"synIdentifier\">sqrt<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">1<\/span><span class=\"synStatement\">\/<\/span>sw<span class=\"synSpecial\">))<\/span>\nq1 <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>w<span class=\"synStatement\">*<\/span><span class=\"synSpecial\">(<\/span>z<span class=\"synStatement\">-<\/span>varz<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">^<\/span><span class=\"synConstant\">2<\/span><span class=\"synSpecial\">)<\/span>\ndf1 <span class=\"synStatement\">&lt;-<\/span> k<span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span>\npval1 <span class=\"synStatement\">&lt;-<\/span> <span class=\"synConstant\">1<\/span><span class=\"synStatement\">-<\/span><span class=\"synIdentifier\">pchisq<\/span><span class=\"synSpecial\">(<\/span>q1<span class=\"synSpecial\">,<\/span> df1<span class=\"synSpecial\">)<\/span>\nq2 <span class=\"synStatement\">&lt;-<\/span> varz<span class=\"synStatement\">^<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>sw\ndf2 <span class=\"synStatement\">&lt;-<\/span> <span class=\"synConstant\">1<\/span>\npval2 <span class=\"synStatement\">&lt;-<\/span> <span class=\"synConstant\">1<\/span><span class=\"synStatement\">-<\/span><span class=\"synIdentifier\">pchisq<\/span><span class=\"synSpecial\">(<\/span>q2<span class=\"synSpecial\">,<\/span> df2<span class=\"synSpecial\">)<\/span>\n<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"\u5909\u91cf\u52b9\u679c\u3067\u76f8\u95a2\u4fc2\u6570\u3092\u7d71\u5408\u3059\u308b\">\u5909\u91cf\u52b9\u679c\u3067\u76f8\u95a2\u4fc2\u6570\u3092\u7d71\u5408\u3059\u308b<\/h2>\n\n\n\n<p>Dersimonian-Laird\u306e\u65b9\u6cd5\u306e\u65b9\u6cd5\u3067\u7d71\u5408\u3059\u308b\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synComment\"># ----- random effects -----<\/span>\ntau2 <span class=\"synStatement\">&lt;-<\/span> <span class=\"synSpecial\">(<\/span>q1<span class=\"synStatement\">-<\/span><span class=\"synSpecial\">(<\/span>k<span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">))<\/span><span class=\"synStatement\">\/<\/span><span class=\"synSpecial\">(<\/span>sw<span class=\"synStatement\">-<\/span><span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>w<span class=\"synStatement\">*<\/span>w<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">\/<\/span>sw<span class=\"synSpecial\">)<\/span>\ntau2 <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">max<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">0<\/span><span class=\"synSpecial\">,<\/span> tau2<span class=\"synSpecial\">)<\/span>\nwx <span class=\"synStatement\">&lt;-<\/span> <span class=\"synConstant\">1<\/span><span class=\"synStatement\">\/<\/span><span class=\"synSpecial\">(<\/span>tau2<span class=\"synStatement\">+<\/span>se<span class=\"synStatement\">*<\/span>se<span class=\"synSpecial\">)<\/span>\nswx <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>wx<span class=\"synSpecial\">)<\/span>\nvarzd <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>wx<span class=\"synStatement\">*<\/span>z<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">\/<\/span>swx\nvarrd <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">tanh<\/span><span class=\"synSpecial\">(<\/span>varzd<span class=\"synSpecial\">)<\/span>\nvarrdl <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">tanh<\/span><span class=\"synSpecial\">(<\/span>varzd<span class=\"synStatement\">-<\/span>Za<span class=\"synStatement\">*<\/span><span class=\"synIdentifier\">sqrt<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">1<\/span><span class=\"synStatement\">\/<\/span>swx<span class=\"synSpecial\">))<\/span>\nvarrdu <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">tanh<\/span><span class=\"synSpecial\">(<\/span>varzd<span class=\"synStatement\">+<\/span>Za<span class=\"synStatement\">*<\/span><span class=\"synIdentifier\">sqrt<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">1<\/span><span class=\"synStatement\">\/<\/span>swx<span class=\"synSpecial\">))<\/span>\nqx2 <span class=\"synStatement\">&lt;-<\/span> varzd<span class=\"synStatement\">^<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>swx\npvalx2 <span class=\"synStatement\">&lt;-<\/span> <span class=\"synConstant\">1<\/span><span class=\"synStatement\">-<\/span><span class=\"synIdentifier\">pchisq<\/span><span class=\"synSpecial\">(<\/span>qx2<span class=\"synSpecial\">,<\/span> df2<span class=\"synSpecial\">)<\/span>\n<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"\u56fa\u5b9a\u52b9\u679c\u3068\u5909\u91cf\u52b9\u679c\u306e\u7d71\u5408\u5024\u3092\u30b0\u30e9\u30d5\u306b\u66f8\u304d\u5165\u308c\u308b\">\u56fa\u5b9a\u52b9\u679c\u3068\u5909\u91cf\u52b9\u679c\u306e\u7d71\u5408\u5024\u3092\u30b0\u30e9\u30d5\u306b\u66f8\u304d\u5165\u308c\u308b<\/h2>\n\n\n\n<p>\u56fa\u5b9a\u52b9\u679c\uff08\u30b0\u30e9\u30d5\u4e2d\u306fCombined: fixed\uff09\u3068\u5909\u91cf\u52b9\u679c\uff08\u30b0\u30e9\u30d5\u4e2d\u306fCombined: random\uff09\u3092\u30b0\u30e9\u30d5\u306b\u66f8\u304d\u5165\u308c\u308b\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synComment\"># ----- graph -----<\/span>\nx <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span>varrl<span class=\"synSpecial\">,<\/span> varru<span class=\"synSpecial\">)<\/span>\ny <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">,<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">lines<\/span><span class=\"synSpecial\">(<\/span>x<span class=\"synSpecial\">,<\/span> y<span class=\"synSpecial\">,<\/span> type<span class=\"synStatement\">=<\/span><span class=\"synConstant\">\"b\"<\/span><span class=\"synSpecial\">)<\/span>\nx <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span>varrdl<span class=\"synSpecial\">,<\/span> varrdu<span class=\"synSpecial\">)<\/span>\ny <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">2<\/span><span class=\"synSpecial\">,<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">2<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">lines<\/span><span class=\"synSpecial\">(<\/span>x<span class=\"synSpecial\">,<\/span> y<span class=\"synSpecial\">,<\/span> type<span class=\"synStatement\">=<\/span><span class=\"synConstant\">\"b\"<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">abline<\/span><span class=\"synSpecial\">(<\/span>v<span class=\"synStatement\">=<\/span><span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span>varr<span class=\"synSpecial\">,<\/span> varrd<span class=\"synSpecial\">),<\/span> lty<span class=\"synStatement\">=<\/span><span class=\"synConstant\">2<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">text<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">0.5<\/span><span class=\"synSpecial\">,<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">,<\/span> <span class=\"synConstant\">\"Combined: fixed\"<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">text<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">0.5<\/span><span class=\"synSpecial\">,<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">2<\/span><span class=\"synSpecial\">,<\/span> <span class=\"synConstant\">\"Combined: random\"<\/span><span class=\"synSpecial\">)<\/span>\n<\/code><\/pre>\n\n\n\n<p>\u30b0\u30e9\u30d5\u306f\u3053\u3093\u306a\u611f\u3058\u306b\u306a\u308b\u3002<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img decoding=\"async\" width=\"553\" height=\"552\" src=\"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104122328.png\" alt=\"\" class=\"wp-image-2688\" title=\"\" srcset=\"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104122328.png 553w, https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104122328-300x300.png 300w, https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104122328-150x150.png 150w\" sizes=\"(max-width: 553px) 100vw, 553px\" \/><\/figure>\n\n\n\n<p><span itemscope=\"\" itemtype=\"http:\/\/schema.org\/Photograph\"><\/span><\/p>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"\u5236\u9650\u4ed8\u304d\u6700\u5c24\u63a8\u5b9a\u91cf\u3092\u7528\u3044\u305f\u65b9\u6cd5REML\u3067\u76f8\u95a2\u4fc2\u6570\u3092\u7d71\u5408\u3059\u308b\">\u5236\u9650\u4ed8\u304d\u6700\u5c24\u63a8\u5b9a\u91cf\u3092\u7528\u3044\u305f\u65b9\u6cd5\uff08REML\uff09\u3067\u76f8\u95a2\u4fc2\u6570\u3092\u7d71\u5408\u3059\u308b<\/h2>\n\n\n\n<p>Restricted Maximum Likelihood estimator\uff08REML\uff09\u3092\u9069\u7528\u3057\u305f\u65b9\u6cd5\u3002<\/p>\n\n\n\n<p>\u4ee5\u4e0b\u306e\u30b9\u30af\u30ea\u30d7\u30c8\u3067\u7e70\u308a\u8fd4\u3057\u53ce\u675f\u8a08\u7b97\u3092\u884c\u3044\u89e3\u3092\u6c42\u3081\u308b\u3002<\/p>\n\n\n\n<p>$ \\hat{\\tau}^2 $ \u306e\u521d\u671f\u5024\u306b\u306fDerSimonian-Laird\u306e\u65b9\u6cd5\u3067\u6c42\u3081\u305f\u5024\u3092\u7528\u3044\u308b\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synComment\"># ----- REML method -----<\/span>\nintau <span class=\"synStatement\">&lt;-<\/span> tau2\ntau <span class=\"synStatement\">&lt;-<\/span> intau\n<span class=\"synComment\">#<\/span>\nnrep <span class=\"synStatement\">&lt;-<\/span> <span class=\"synConstant\">10<\/span>\nnewt <span class=\"synStatement\">&lt;-<\/span> <span class=\"synConstant\">1<\/span><span class=\"synSpecial\">:<\/span>nrep\n<span class=\"synStatement\">for<\/span> <span class=\"synSpecial\">(<\/span>i <span class=\"synStatement\">in<\/span> <span class=\"synConstant\">1<\/span><span class=\"synSpecial\">:<\/span>nrep<span class=\"synSpecial\">){<\/span>\nwb <span class=\"synStatement\">&lt;-<\/span> <span class=\"synConstant\">1<\/span><span class=\"synStatement\">\/<\/span><span class=\"synSpecial\">(<\/span>tau<span class=\"synStatement\">+<\/span>se<span class=\"synStatement\">*<\/span>se<span class=\"synSpecial\">)<\/span>\nzmb <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>wb<span class=\"synStatement\">*<\/span>z<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">\/<\/span><span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>wb<span class=\"synSpecial\">)<\/span>\nqf <span class=\"synStatement\">&lt;-<\/span> k<span class=\"synStatement\">\/<\/span><span class=\"synSpecial\">(<\/span>k<span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">)<\/span><span class=\"synStatement\">*<\/span><span class=\"synSpecial\">(<\/span>z<span class=\"synStatement\">-<\/span>zmb<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">^<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">-<\/span>se<span class=\"synStatement\">*<\/span>se\ndkx <span class=\"synStatement\">&lt;-<\/span> <span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synStatement\">*<\/span><span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>z<span class=\"synStatement\">*<\/span>wb<span class=\"synStatement\">*<\/span>wb<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">+<\/span>zmb<span class=\"synStatement\">*<\/span><span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>wb<span class=\"synStatement\">*<\/span>wb<span class=\"synSpecial\">))<\/span><span class=\"synStatement\">\/<\/span><span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>wb<span class=\"synSpecial\">)<\/span>\nqf2 <span class=\"synStatement\">&lt;-<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>k<span class=\"synStatement\">\/<\/span><span class=\"synSpecial\">(<\/span>k<span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">)<\/span><span class=\"synStatement\">*<\/span><span class=\"synSpecial\">(<\/span>z<span class=\"synStatement\">-<\/span>zmb<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">*<\/span>dkx\nh <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>wb<span class=\"synStatement\">*<\/span>wb<span class=\"synStatement\">*<\/span><span class=\"synSpecial\">(<\/span>qf<span class=\"synStatement\">-<\/span>tau<span class=\"synSpecial\">))<\/span>\ndh <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>wb<span class=\"synStatement\">*<\/span>wb<span class=\"synStatement\">*<\/span>wb<span class=\"synStatement\">*<\/span><span class=\"synSpecial\">(<\/span>qf<span class=\"synStatement\">-<\/span>tau<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">+<\/span>wb<span class=\"synStatement\">*<\/span>wb<span class=\"synStatement\">*<\/span><span class=\"synSpecial\">(<\/span>qf2<span class=\"synStatement\">-<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">))<\/span>\nnewt<span class=\"synSpecial\">&#91;<\/span>i<span class=\"synSpecial\">]<\/span> <span class=\"synStatement\">&lt;-<\/span> tau <span class=\"synStatement\">-<\/span> h<span class=\"synStatement\">\/<\/span>dh\nrel <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">abs<\/span><span class=\"synSpecial\">((<\/span>newt<span class=\"synSpecial\">&#91;<\/span>i<span class=\"synSpecial\">]<\/span><span class=\"synStatement\">-<\/span>tau<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">\/<\/span>tau<span class=\"synSpecial\">)<\/span>\ntau <span class=\"synStatement\">&lt;-<\/span> newt<span class=\"synSpecial\">&#91;<\/span>i<span class=\"synSpecial\">]<\/span>\n<span class=\"synSpecial\">}<\/span>\nwg <span class=\"synStatement\">&lt;-<\/span> <span class=\"synConstant\">1<\/span><span class=\"synStatement\">\/<\/span><span class=\"synSpecial\">(<\/span>tau<span class=\"synStatement\">+<\/span>se<span class=\"synStatement\">*<\/span>se<span class=\"synSpecial\">)<\/span>\nswg <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>wg<span class=\"synSpecial\">)<\/span>\nzRM <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">sum<\/span><span class=\"synSpecial\">(<\/span>z<span class=\"synStatement\">*<\/span>wg<span class=\"synSpecial\">)<\/span><span class=\"synStatement\">\/<\/span>swg\nrRM <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">tanh<\/span><span class=\"synSpecial\">(<\/span>zRM<span class=\"synSpecial\">)<\/span>\nrRMl <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">tanh<\/span><span class=\"synSpecial\">(<\/span>zRM<span class=\"synStatement\">-<\/span>Za<span class=\"synStatement\">*<\/span><span class=\"synIdentifier\">sqrt<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">1<\/span><span class=\"synStatement\">\/<\/span>swg<span class=\"synSpecial\">))<\/span>\nrRMu <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">tanh<\/span><span class=\"synSpecial\">(<\/span>zRM<span class=\"synStatement\">+<\/span>Za<span class=\"synStatement\">*<\/span><span class=\"synIdentifier\">sqrt<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">1<\/span><span class=\"synStatement\">\/<\/span>swg<span class=\"synSpecial\">))<\/span>\nqx2RM <span class=\"synStatement\">&lt;-<\/span> zRM<span class=\"synStatement\">^<\/span><span class=\"synConstant\">2<\/span><span class=\"synStatement\">*<\/span>swg\npvalx2RM <span class=\"synStatement\">&lt;-<\/span> <span class=\"synConstant\">1<\/span><span class=\"synStatement\">-<\/span><span class=\"synIdentifier\">pchisq<\/span><span class=\"synSpecial\">(<\/span>qx2RM<span class=\"synSpecial\">,<\/span> df2<span class=\"synSpecial\">)<\/span>\ntau2 <span class=\"synStatement\">&lt;-<\/span> tau\n<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"REML\u306e\u7d50\u679c\u3092\u30b0\u30e9\u30d5\u306b\u66f8\u304d\u5165\u308c\u308b\">REML\u306e\u7d50\u679c\u3092\u30b0\u30e9\u30d5\u306b\u66f8\u304d\u5165\u308c\u308b<\/h2>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synComment\"># ----- REML graph -----<\/span>\nx <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span>rRMl<span class=\"synSpecial\">,<\/span> rRMu<span class=\"synSpecial\">)<\/span>\ny <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">3<\/span><span class=\"synSpecial\">,<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">3<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">lines<\/span><span class=\"synSpecial\">(<\/span>x<span class=\"synSpecial\">,<\/span> y<span class=\"synSpecial\">,<\/span> type<span class=\"synStatement\">=<\/span><span class=\"synConstant\">\"b\"<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">abline<\/span><span class=\"synSpecial\">(<\/span>v<span class=\"synStatement\">=<\/span>rRM<span class=\"synSpecial\">,<\/span> lty<span class=\"synStatement\">=<\/span><span class=\"synConstant\">2<\/span><span class=\"synSpecial\">)<\/span>\n<span class=\"synIdentifier\">text<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">-<\/span><span class=\"synConstant\">0.5<\/span><span class=\"synSpecial\">,<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">3<\/span><span class=\"synSpecial\">,<\/span> <span class=\"synConstant\">\"Combined: REML\"<\/span><span class=\"synSpecial\">)<\/span>\n<\/code><\/pre>\n\n\n\n<p>\u30b0\u30e9\u30d5\u306f\u3053\u3061\u3089\u300295%\u4fe1\u983c\u533a\u9593\u306f\u82e5\u5e72\u5e83\u304f\u306a\u308b\u304c\u3001\u7d71\u8a08\u5b66\u7684\u6709\u610f\u3067\u306f\u3042\u3063\u305f\u3002<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img decoding=\"async\" width=\"553\" height=\"552\" src=\"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104122855.png\" alt=\"\" class=\"wp-image-2689\" title=\"\" srcset=\"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104122855.png 553w, https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104122855-300x300.png 300w, https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104122855-150x150.png 150w\" sizes=\"(max-width: 553px) 100vw, 553px\" \/><\/figure>\n\n\n\n<p><span itemscope=\"\" itemtype=\"http:\/\/schema.org\/Photograph\"><\/span><\/p>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"\u7d50\u679c\u306e\u6570\u5024\u3092\u4e26\u3079\u3066\u6bd4\u3079\u3066\u307f\u308b\">\u7d50\u679c\u306e\u6570\u5024\u3092\u4e26\u3079\u3066\u6bd4\u3079\u3066\u307f\u308b\u3002<\/h2>\n\n\n\n<p>\u70b9\u63a8\u5b9a\u5024\u306f\u56fa\u5b9a\u52b9\u679c\u2192DerSimonian-Laird\u306e\u65b9\u6cd5\u2192REML\u3068\u3060\u3093\u3060\u3093\u306b\u5927\u304d\u304f\u306a\u308a\u300195%\u4fe1\u983c\u533a\u9593\u3082\u3060\u3093\u3060\u3093\u5e83\u304f\u306a\u3063\u3066\u3044\u3063\u305f\u3002<\/p>\n\n\n\n<p>\u6700\u7d42\u7684\u306b\u3001\u7d71\u8a08\u5b66\u7684\u6709\u610f\uff08$ P_{REML} = 0.0003 $\uff09<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synStatement\">&gt;<\/span> <span class=\"synComment\">## Fixed effects<\/span>\n<span class=\"synStatement\">&gt;<\/span> <span class=\"synType\">list<\/span><span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">round<\/span><span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span>rFE<span class=\"synStatement\">=<\/span>varr<span class=\"synSpecial\">,<\/span> LL<span class=\"synStatement\">=<\/span>varrl<span class=\"synSpecial\">,<\/span> UL<span class=\"synStatement\">=<\/span>varru<span class=\"synSpecial\">,<\/span> Q1<span class=\"synStatement\">=<\/span>q1<span class=\"synSpecial\">,<\/span> df1<span class=\"synStatement\">=<\/span>df1<span class=\"synSpecial\">,<\/span> p1<span class=\"synStatement\">=<\/span>pval1<span class=\"synSpecial\">,<\/span> Q2<span class=\"synStatement\">=<\/span>q2<span class=\"synSpecial\">,<\/span> df2<span class=\"synStatement\">=<\/span>df2<span class=\"synSpecial\">,<\/span> p2<span class=\"synStatement\">=<\/span>pval2<span class=\"synSpecial\">),<\/span><span class=\"synConstant\">4<\/span><span class=\"synSpecial\">))<\/span>\n<span class=\"synSpecial\">&#91;&#91;<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">]]<\/span>\nrFE       LL       UL       Q1      df1       p1       Q2      df2       p2\n<span class=\"synConstant\">0.1463<\/span>   <span class=\"synConstant\">0.1298<\/span>   <span class=\"synConstant\">0.1627<\/span>  <span class=\"synConstant\">76.8331<\/span>   <span class=\"synConstant\">5.0000<\/span>   <span class=\"synConstant\">0.0000<\/span> <span class=\"synConstant\">295.8090<\/span>   <span class=\"synConstant\">1.0000<\/span>   <span class=\"synConstant\">0.0000<\/span>\n<span class=\"synStatement\">&gt;<\/span>\n<span class=\"synStatement\">&gt;<\/span> <span class=\"synComment\">## Random effects<\/span>\n<span class=\"synStatement\">&gt;<\/span> <span class=\"synType\">list<\/span><span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">round<\/span><span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span>rDL<span class=\"synStatement\">=<\/span>varrd<span class=\"synSpecial\">,<\/span> LL<span class=\"synStatement\">=<\/span>varrdl<span class=\"synSpecial\">,<\/span> UL<span class=\"synStatement\">=<\/span>varrdu<span class=\"synSpecial\">,<\/span> Q2DL<span class=\"synStatement\">=<\/span>qx2<span class=\"synSpecial\">,<\/span> df2<span class=\"synStatement\">=<\/span>df2<span class=\"synSpecial\">,<\/span> p2DL<span class=\"synStatement\">=<\/span>pvalx2<span class=\"synSpecial\">,<\/span> <span class=\"synConstant\">\"tau^2\"<\/span><span class=\"synStatement\">=<\/span>intau<span class=\"synSpecial\">),<\/span><span class=\"synConstant\">4<\/span><span class=\"synSpecial\">))<\/span>\n<span class=\"synSpecial\">&#91;&#91;<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">]]<\/span>\nrDL      LL      UL    Q2DL     df2    p2DL   tau<span class=\"synStatement\">^<\/span><span class=\"synConstant\">2<\/span>\n<span class=\"synConstant\">0.1741<\/span>  <span class=\"synConstant\">0.0966<\/span>  <span class=\"synConstant\">0.2496<\/span> <span class=\"synConstant\">19.0392<\/span>  <span class=\"synConstant\">1.0000<\/span>  <span class=\"synConstant\">0.0000<\/span>  <span class=\"synConstant\">0.0078<\/span>\n<span class=\"synStatement\">&gt;<\/span>\n<span class=\"synStatement\">&gt;<\/span> <span class=\"synComment\">## REML<\/span>\n<span class=\"synStatement\">&gt;<\/span> <span class=\"synType\">list<\/span><span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">round<\/span><span class=\"synSpecial\">(<\/span><span class=\"synIdentifier\">c<\/span><span class=\"synSpecial\">(<\/span>rRM<span class=\"synStatement\">=<\/span>rRM<span class=\"synSpecial\">,<\/span> LL<span class=\"synStatement\">=<\/span>rRMl<span class=\"synSpecial\">,<\/span> UL<span class=\"synStatement\">=<\/span>rRMu<span class=\"synSpecial\">,<\/span> Q2RM<span class=\"synStatement\">=<\/span>qx2RM<span class=\"synSpecial\">,<\/span> df2<span class=\"synStatement\">=<\/span>df2<span class=\"synSpecial\">,<\/span> p2RM<span class=\"synStatement\">=<\/span>pvalx2RM<span class=\"synSpecial\">,<\/span> <span class=\"synConstant\">\"tau^2\"<\/span><span class=\"synStatement\">=<\/span>tau2<span class=\"synSpecial\">),<\/span><span class=\"synConstant\">4<\/span><span class=\"synSpecial\">))<\/span>\n<span class=\"synSpecial\">&#91;&#91;<\/span><span class=\"synConstant\">1<\/span><span class=\"synSpecial\">]]<\/span>\nrRM      LL      UL    Q2RM     df2    p2RM   tau<span class=\"synStatement\">^<\/span><span class=\"synConstant\">2<\/span>\n<span class=\"synConstant\">0.1800<\/span>  <span class=\"synConstant\">0.0831<\/span>  <span class=\"synConstant\">0.2735<\/span> <span class=\"synConstant\">13.0604<\/span>  <span class=\"synConstant\">1.0000<\/span>  <span class=\"synConstant\">0.0003<\/span>  <span class=\"synConstant\">0.0130<\/span>\n<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"\u307e\u3068\u3081\">\u307e\u3068\u3081<\/h2>\n\n\n\n<p>\u76f8\u95a2\u4fc2\u6570\u306e\u7d71\u5408\u3092 R \u3067\u3084\u3063\u3066\u307f\u305f\u3002<\/p>\n\n\n\n<p>\u76f8\u95a2\u4fc2\u6570\u306bFisher\u306e\u5206\u6563\u5b89\u5b9a\u5316\u5909\u63db\u3092\u307b\u3069\u3053\u3057\u3066\u3001n &#8211; 3 \u3092\u91cd\u307f\u306b\u3057\u3066\u91cd\u307f\u3065\u3051\u5e73\u5747\u3092\u53d6\u308b\u306e\u304c\u4e2d\u6838\u3060\u3002<\/p>\n\n\n\n<p>\u76f8\u95a2\u4fc2\u6570\u304c\u95a2\u5fc3\u4e8b\u306e\u7814\u7a76\u3092\u7d71\u5408\u3059\u308b\u3068\u304d\u306b\u306f\u3053\u306e\u65b9\u6cd5\u304c\u4f7f\u3048\u308b\u3002<\/p>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"\u95a2\u9023\u8a18\u4e8b\">\u95a2\u9023\u8a18\u4e8b<\/h2>\n\n\n\n<p>\u6f38\u8fd1\u5206\u6563\u6cd5\u306b\u3088\u308b\u30aa\u30c3\u30ba\u6bd4\u306e\u7d71\u5408<\/p>\n\n\n<div class=\"swell-block-postLink\">\t\t\t<div class=\"p-blogCard -internal\" data-type=\"type1\" data-onclick=\"clickLink\">\n\t\t\t\t<div class=\"p-blogCard__inner\">\n\t\t\t\t\t<span class=\"p-blogCard__caption\">\u3042\u308f\u305b\u3066\u8aad\u307f\u305f\u3044<\/span>\n\t\t\t\t\t<div class=\"p-blogCard__thumb c-postThumb\"><figure class=\"c-postThumb__figure\"><img decoding=\"async\" src=\"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2018\/07\/20200920210208-300x300.png\" alt=\"\" class=\"c-postThumb__img u-obf-cover\" width=\"320\" height=\"180\"><\/figure><\/div>\t\t\t\t\t<div class=\"p-blogCard__body\">\n\t\t\t\t\t\t<a class=\"p-blogCard__title\" href=\"https:\/\/best-biostatistics.com\/toukei-er\/entry\/meta-analysis-odds-ratio-asymptotic-variance-based-method\/\">R \u3067\u6f38\u8fd1\u5206\u6563\u6cd5\u3067\u30aa\u30c3\u30ba\u6bd4\u3092\u7d71\u5408\u3059\u308b\u30e1\u30bf\u30a2\u30ca\u30ea\u30b7\u30b9<\/a>\n\t\t\t\t\t\t<span class=\"p-blogCard__excerpt\">\u30e1\u30bf\u30a2\u30ca\u30ea\u30b7\u30b9\u306e\u305f\u3081\u306b\u30aa\u30c3\u30ba\u6bd4\u3092\u7d71\u5408\u3057\u305f\u3044\u5834\u5408\u3067\u3001\u305d\u308c\u305e\u308c\u306e\u7814\u7a76\u306e\u30b5\u30f3\u30d7\u30eb\u30b5\u30a4\u30ba\u304c\u5927\u304d\u3044\u5834\u5408\u3001\u30aa\u30c3\u30ba\u6bd4\u306e\u5bfe\u6570\u304c\u6f38\u8fd1\u7684\u306b\u6b63\u898f\u8fd1\u4f3c\u3067\u304d\u308b\u3002 \u500b\u3005\u306e\u7814\u7a76\u306e\u30b5\u30f3\u30d7\u30eb\u30b5&#8230;<\/span>\t\t\t\t\t<\/div>\n\t\t\t\t<\/div>\n\t\t\t<\/div>\n\t\t<\/div>\n\n\n\n\n\n\n<h2 class=\"wp-block-heading\" id=\"\u5f15\u7528\u66f8\u7c4d\">\u5f15\u7528\u66f8\u7c4d<\/h2>\n\n\n\n<p>\u30cd\u30bf\u672c\u306f\u3053\u3061\u3089<\/p>\n\n\n\n<p>4.2 \u76f8\u95a2\u4fc2\u6570<\/p>\n\n\n\n<figure class=\"wp-block-image\"><a class=\"hatena-asin-detail-image-link\" href=\"https:\/\/www.amazon.co.jp\/dp\/4254127545?tag=ttyotani-22&amp;linkCode=ogi&amp;th=1&amp;psc=1\" target=\"_blank\" rel=\"noopener\"><img decoding=\"async\" src=\"https:\/\/m.media-amazon.com\/images\/I\/41PBSEteD0L._SL500_.jpg\" alt=\"\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)\" title=\"\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><\/figure>\n\n\n\n<div class=\"hatena-asin-detail\"><\/div>\n\n\n\n<div class=\"hatena-asin-detail-info\">\n<p class=\"hatena-asin-detail-title\"><a href=\"https:\/\/www.amazon.co.jp\/dp\/4254127545?tag=ttyotani-22&amp;linkCode=ogi&amp;th=1&amp;psc=1\" target=\"_blank\" rel=\"noopener\">\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<ul class=\"hatena-asin-detail-meta\">\n<li><span class=\"hatena-asin-detail-label\">\u4f5c\u8005:<\/span><a href=\"https:\/\/d.hatena.ne.jp\/keyword\/%BD%D3%CF%BA%2C%20%C3%B0%B8%E5\" class=\"keyword\">\u4fca\u90ce, \u4e39\u5f8c<\/a><\/li>\n<li>\u671d\u5009\u66f8\u5e97<\/li>\n<\/ul>\n<p><a href=\"https:\/\/www.amazon.co.jp\/dp\/4254127545?tag=ttyotani-22&amp;linkCode=ogi&amp;th=1&amp;psc=1\" class=\"asin-detail-buy\" target=\"_blank\" rel=\"noopener\">Amazon<\/a><\/p>\n<\/div>\n\n\n\n<p>\u65b0\u7248\u306f\u3053\u3061\u3089<\/p>\n\n\n\n<figure class=\"wp-block-image\"><a class=\"hatena-asin-detail-image-link\" href=\"https:\/\/www.amazon.co.jp\/dp\/425412760X?tag=ttyotani-22&amp;linkCode=ogi&amp;th=1&amp;psc=1\" target=\"_blank\" rel=\"noopener\"><img decoding=\"async\" src=\"https:\/\/m.media-amazon.com\/images\/I\/41P-7dAdGgL._SL500_.jpg\" alt=\"\u65b0\u7248 \u30e1\u30bf\u30fb\u30a2\u30ca\u30ea\u30b7\u30b9\u5165\u9580 \u2500\u30a8\u30d3\u30c7\u30f3\u30b9\u306e\u7d71\u5408\u3092\u3081\u3056\u3059\u7d71\u8a08\u624b\u6cd5\u2500 (\u533b\u5b66\u7d71\u8a08\u5b66\u30b7\u30ea\u30fc\u30ba)\" title=\"\u65b0\u7248 \u30e1\u30bf\u30fb\u30a2\u30ca\u30ea\u30b7\u30b9\u5165\u9580 \u2500\u30a8\u30d3\u30c7\u30f3\u30b9\u306e\u7d71\u5408\u3092\u3081\u3056\u3059\u7d71\u8a08\u624b\u6cd5\u2500 (\u533b\u5b66\u7d71\u8a08\u5b66\u30b7\u30ea\u30fc\u30ba)\"\/><\/a><\/figure>\n\n\n\n<div class=\"hatena-asin-detail\"><\/div>\n\n\n\n<div class=\"hatena-asin-detail-info\">\n<p class=\"hatena-asin-detail-title\"><a href=\"https:\/\/www.amazon.co.jp\/dp\/425412760X?tag=ttyotani-22&amp;linkCode=ogi&amp;th=1&amp;psc=1\" target=\"_blank\" rel=\"noopener\">\u65b0\u7248 \u30e1\u30bf\u30fb\u30a2\u30ca\u30ea\u30b7\u30b9\u5165\u9580 \u2500\u30a8\u30d3\u30c7\u30f3\u30b9\u306e\u7d71\u5408\u3092\u3081\u3056\u3059\u7d71\u8a08\u624b\u6cd5\u2500 (\u533b\u5b66\u7d71\u8a08\u5b66\u30b7\u30ea\u30fc\u30ba)<\/a><\/p>\n<ul class=\"hatena-asin-detail-meta\">\n<li><span class=\"hatena-asin-detail-label\">\u4f5c\u8005:<\/span><a href=\"https:\/\/d.hatena.ne.jp\/keyword\/%C3%B0%B8%E5%20%BD%D3%CF%BA\" class=\"keyword\">\u4e39\u5f8c \u4fca\u90ce<\/a><\/li>\n<li>\u671d\u5009\u66f8\u5e97<\/li>\n<\/ul>\n<p><a href=\"https:\/\/www.amazon.co.jp\/dp\/425412760X?tag=ttyotani-22&amp;linkCode=ogi&amp;th=1&amp;psc=1\" class=\"asin-detail-buy\" target=\"_blank\" rel=\"noopener\">Amazon<\/a><\/p>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>\u76f8\u95a2\u4fc2\u6570\u3092\u7d71\u5408\u3057\u305f\u3044\u5834\u5408\u306f\u3069\u3046\u3084\u308b\u304b\uff1f R \u3067\u306e\u65b9\u6cd5\u3002<\/p>\n","protected":false},"author":2,"featured_media":2689,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"swell_btn_cv_data":"","_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_publicize_message":"","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":false,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2}},"categories":[5,46,36],"tags":[],"class_list":["post-456","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-r","category-46","category-36"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2019\/01\/20190104122855.png","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/456","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/comments?post=456"}],"version-history":[{"count":4,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/456\/revisions"}],"predecessor-version":[{"id":2692,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/456\/revisions\/2692"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/media\/2689"}],"wp:attachment":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/media?parent=456"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/categories?post=456"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/tags?post=456"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}