{"id":1886,"date":"2019-07-04T05:23:27","date_gmt":"2019-07-03T20:23:27","guid":{"rendered":"http:\/\/141.164.34.82\/?p=1886"},"modified":"2019-07-04T05:33:55","modified_gmt":"2019-07-03T20:33:55","slug":"%ec%9d%b8%ea%b3%bc%ea%b4%80%ea%b3%84-%ec%b6%94%ec%a0%95%ec%97%90%ec%84%9c-%eb%aa%a8%ed%98%95-%ec%84%a4%ec%a0%95%ec%9d%98-%ec%a4%91%ec%9a%94%ec%84%b1","status":"publish","type":"post","link":"http:\/\/ds.sumeun.org\/?p=1886","title":{"rendered":"\uc778\uacfc\uad00\uacc4 \ucd94\uc815\uc5d0\uc11c \ubaa8\ud615 \uc124\uc815\uc758 \uc911\uc694\uc131"},"content":{"rendered":"<h1>\uc2ec\uc2a8\uacfc \ubc84\ud06c\uc2a8<\/h1>\n<h2>\uc2ec\uc2a8\uc758 \uc5ed\uc124(Simpson&#39;s paradox)<\/h2>\n<p>R\uc758 <code>datasets::UCBAdmissions<\/code>\uc740 1973\ub144 UC Berkeley\uc758 \ub300\ud559\uc6d0 \uc785\ud559 \uc2dc\ud5d8 \uacb0\uacfc\ub97c \ubcf4\uc5ec \uc900\ub2e4. <\/p>\n<p>1973\ub144 UC \ubc84\ud074\ub9ac(Berkeley) \ub300\ud559\uc740 \uc131\ucc28\ubcc4\ub85c \uace0\uc18c\ub97c \ub2f9\ud55c\ub2e4. \uc9c0\uc6d0\uc790 \uc911 \ub0a8\uc131\uc740 44%\uac00 \ud569\uaca9\ud588\uc9c0\ub9cc, \uc5ec\uc131\uc740 35%\ub9cc\uc774 \ud569\uaca9\uc744 \ud588\uae30 \ub54c\ubb38\uc774\ub2e4. \ub2e4\uc74c\uc740 \ubc84\ud074\ub9ac \ub300\ud559\uc6d0 \uc785\ud559 \uacb0\uacfc\ub97c \uadf8\ub798\ud504\ub85c \ubcf4\uc5ec\uc8fc\uace0 \uc788\ub2e4. <\/p>\n<pre><code class=\"r\"># from help(UCBAdmissions)\r\nrequire(graphics)\r\n## Data aggregated over departments\r\napply(UCBAdmissions, c(1, 2), sum)\r\n<\/code><\/pre>\n<pre>##           Gender\r\n## Admit      Male Female\r\n##   Admitted 1198    557\r\n##   Rejected 1493   1278\r\n<\/pre>\n<pre><code class=\"r\">library(vcd)\r\nmosaic(apply(UCBAdmissions, c(1,2), sum), shade=TRUE, legend=TRUE)\r\n<\/code><\/pre>\n<p><img src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAfgAAAH4CAMAAACR9g9NAAAA9lBMVEUAAAAAADoAAGYAOjoAOmYAOpAAZrY6AAA6ADo6AGY6Ojo6OmY6OpA6ZmY6ZpA6ZrY6kJA6kLY6kNtKb+NmAABmADpmAGZmOgBmOjpmOmZmOpBmZgBmZjpmZmZmZpBmZrZmkNtmtv+QOgCQOjqQOmaQZgCQZjqQZmaQZpCQkDqQkGaQtpCQtraQttuQ29uQ2\/+dqOK2ZgC2Zjq2kGa2tma2tpC2tra2ttu225C229u22\/+2\/7a2\/9u2\/\/\/TP2rbkDrbkGbbtpDbtrbb27bb29vb2\/\/b\/7bb\/9vb\/\/\/i4uLklaX\/tmb\/25D\/27b\/29v\/\/7b\/\/9v\/\/\/9yvoo5AAAACXBIWXMAAAsSAAALEgHS3X78AAASaElEQVR4nO2dC3vbthVAaTeJnMVtOs+1u67ZJqfNunarp27L+liqdu1Sb+Ii8f\/\/mRF8iBQNygJc4oK853ztJ0W6N2R8jAdBEEgyUEkifQIgA+KVgnilIF4piFcK4pWCeKUgXimIVwrilYJ4pSBeKYhXCuKVgnilIF4piFcK4pWCeKUgXimIVwrilYJ4pSBeKYhXCuKVgnilIF4piFcK4pWCeKUgXimIVwrilYJ4pSBeKYhXCuKVgnilIF4piFcK4pWCeKUgXimIVwrilTIJ8ZsXSfLWny1fLI5fBT+ZkTAF8ZurhzebxdH17W8Q38sUxK9q5+ll8uA6S2dvXx79PsuWyYPf5OLLD9fnb82Ov5Y9z6iYgvi6XKezs03+Pp09\/OnKvJytz4sX86F5K3yacTER8evzJJkv85K\/Sua5avOR+VP1kn+4Pn8kfZpxMRHxmZG7TAxnLfHL\/KX8EPEdpiA+d16KL9v62yU+B\/EdpiA+79V\/\/b+rso7PPVfimzbefIj4DlMQX1zHHz29yVaz5Ghel\/i6V19+iPgOkxAP7iBeKYhXCuKVgnilRCg+nZkOeHN\/JX0c5WBrOitGhvq\/j\/O0a2IUf\/LOTbZ+\/3Hs4u84rUhPuyZG8Y8\/vM7SZ\/nPzRSqefETXJ8nsd1kqcWWp7Z5\/mle\/FdFFRD1addEKf5v8+y7zx6\/Wr93bX585ie4OMuWkY3A1OLLU9tcPSraqDT2066JUvyXTzcffVn+XNcXr6qfZf5O+sx2KNr4R1l1apvn15n53\/wxi\/m0a6IU\/89PXv+uKFCLJK8qzU\/wPP8h26bYCFKV+OrUdsTHfNo1UYp\/9cWnZ8XPbV7VmTEWm3RbJZmXlvi4T7smTvGro+uqlUxPrqvGcvXwRvrMdmi18fmptcTHfdo1cYqvrC+T5MHlvOoex1Zltnr1R9ftEh\/3addEKB5CgHilIF4piFcK4pWCeKUgXimIVwrileIiPlFD6H\/sUHb3yXSJ\/ZUSCvHfhgPxkYD4Tqy0kFAgvhMrLSQUiO\/ESgsJxSjE3zW9+y6ZLrHSQkIxDvFmPkA1xc8DxFsYj\/jN1byew11N6n5yefxqlRSflF+kpx+bLxB\/AOMRn55cV3O460nds3kxDdDM6i4nfpnlAizzvxBvYRziizZ+nrXmcJeTultTvLdTAG8\/1YN4C+MQ38z5K+dw15O6i9+Jo6L8m2mAiD+ccYmvSvt2UnfxaV65U+KdGZX4uilvJnWbFt383\/oC8QcxLvHVHO5mUndZ6W979Yg\/mFGIvyeIt4D4Tqy0kFAgvhMrLSQUIxVfteRmNM\/2eUemw1+M+LjFVyDem1GIL0blq3H6cnC+ekyz6tunO2vJNBGI72cc4mfzeh2WanC+fjA7qcU3i7LUw\/ctmYi\/zTjEV6PyudNqcL4SnFf1rWv3ZlGWjkzE32Y84qtx+u3gfDl+04jfjt+XEYjfz3jEN2utlIOzuyV+d\/x+994s4i2MRnw1HF8NzjdtfLEux3FrLZl6+B7xexmP+GqcfrHts2+uTK\/eDNy\/e9FeS2ZBr\/4ARiH+niDeAuI7sdJCQoH4Tqy0kFAgvhMrLSQUiO\/ESgsJxbTE1zdt0ll7DXXEW5iW+GX5OEVxbd9cyiPewqTEp08uC\/Hp6U17pA\/xFsYj\/u7FNTbPP7+ixB\/IiMTvnLON5Vndxq\/Pk9agLeItTEh8XsFX4s1jdqtjqvp9TEj8Mqmfojf3aFq35RFvYULis+3lHCX+bqYnvpp215qLgXgL0xLfk+l1lImD+L6jTBzV4ot5fIZWf0DQRVBUi8+KCV1ZezY24nWILwd228O7gi6Colz85qoo8a3hXUEXQVEu3szYNNvmWo4ycbSL7z\/KxEF831EmjnbxeVX\/8HXraTvE6xCfd+7S05sVnTtt4s0Dtqc3XM6NXbzzRIyyxHM5N3bx1WRL4789Grfnby0u59q\/JIIugjIp8fVkS1N189DkfkYk\/i81d062LGbZHjQDhyHbUYnvidhOtjx06tX27hxt\/JjFtyZbOpb49lGkhYRiQuKbyZaHt\/FU9RMQn20v58w9t0N69VT1UxJfLZF12HU8Vf00xPdkusRKCwmFavHri2\/OmXOnULwtVlpIKLSLX1HiVYq\/tfAt4pWIp1evU3y27OxRjXgt4mnjVYqnjdcqnjZep3jaeKXi14zc6RR\/O1ZaSCgQ34mVFhIK1eJZGEGp+IyFEbSKZ+qVUvEsjKBUPAsjaBV\/K1ZaSCgQ34mVFhIKxJs7dFzO6RO\/u+wt4nWILzaeXszbsdJCQqFafHn9jnh14osx2zni9Yk3LGjjdYo3gzj06lWK34mVFhIKxHdipYWEAvGdWGkhoUB8J1ZaSCgQ34mVFhIKxHdipYWEAvGdWGkhoUB8J1ZaSCgQ34mVFhIKxHdipYWEAvGdWGkhoUB8J1ZaSCgQ34mVFhIKxHdipYWEAvGdWGkhoZiS+HSWJI570iB+AuLNejbpiZlN47AnDeLHL35lVBcT6NiT5m5GJP4\/NX170mT1KlaH71CB+FGJ740pH349eE8axE9E\/Pq8XMOKEn83UxKfzqoZ8rTxdzMh8VvvB+9JYzuKtJBQTEh8uSzt3GlPGsRPQHx\/pkustJBQIL4TKy0kFIjvxEoLCQXiO7HSQkKB+E6stJBQIL4TKy0kFIjvxEoLCQXiO7HSQkKB+E6stJBQIL4TKy0kFIjvxEoLCQXiO7HSQkKB+E6stJBQIL4TKy0kFIjvxEoLCQXiO7HSQkKB+E6stJBQIL4TKy0kFIjvxKoh9D\/WUVpLSRDxEB2IVwrilUIbv5foz91RWktJGPH\/HSeFeOkrhX0gfhgQb8t0iZU26AnibZkusdIGPUG8LdMlVtqgJ4i3ZbrEShv0BPG2TJdYaYOeIN6W6RIrbdATxNsyXWKlDXqCeFumS6y0QU8Qb8t0iZU26AnibZkusdIGPUG8LdMlVtqgJ4i3ZbrEShv0BPG2TJdYaYOeIN6W6RIrbdATxNsyXWKlDXqCeFumS6y0QU8Qb8ts\/2H9a7PGrVnt1h4rbdATxNsym7ebv3548vazZ88ujxEfliHFNyvY9q9lu6qm\/p31HUXaoCeaxTc70ezbk6as6vuPIm3QE83im1Xq+9er33zy0x+eGX57Y\/kbED8cP4t4+0ztZl+K\/h0q1hffnBe5tPGBGbDENzvR7N2T5s2Phn\/3HUXaoCeaxR9S4u8+irRBTzSLP6SNz0nf\/0XOE6r6sAzaqz\/b9ur796TZXCWIF2D46\/g79qTpHbND\/KCID9lmX\/ySzp0A8uIXXM5JIC5+ff6EARwBxMVvPqKNl0Be\/BcP\/vTy5ct\/UOLDIi5+zZCtCOLiN\/96+ZISHx7piRildsQHR1h8VdFT1QdHvsS\/yDt3f798SokPi3gbXw7Zrj9AfFjkxZ8\/\/fHHH15Q1QdGXHy2ZLKlBPLiszd5p76nokf8YEQgfn+stEFPEG\/L7H6w+Z4HKkITgfjv38\/b+FN69WGRFv\/Di9z6249779AhfiDkR+6efLVv+hXiB0JefPLOV4gPj3RVX9b1R5\/3xkob9ATxtszOn7+\/5CZNcGIQbx6TZ159YOIQvydW2qAniLdlusRKG\/QE8bZMl1hpg54g3pbpEitt0BPE2zJdYqUNeoJ4W6ZLrLRBTxBvy3SJlTboCeJtmS6x0gY9Qbwt0yVW2qAniLdlusRKG\/QE8bZMl1hpg54g3pbpEitt0BPE2zJdYqUNeoJ4W6ZLrLRBTxBvy3SJlTboCeJtmS6x0gY9Qbwt0yVW2qAniLdlusRKG\/QE8bZMl1hpg54g3pbpEitt0BPE2zJdYqUNeqJafDpLknnxbs8ixoiXYUDxZmuC9MRsT7BvMyLEyzCg+JVRvTBFfu9GBYgX4R7im3+kfTMiA1uTxMrPIr43xmxNkd2xGRHiRRhI\/CJJHpk+XbmmESU+PoYs8ems7NPTxkfIgOK33vduRnTwUcaFZvHlAnbzOzYjamNZ2RbxAzFs564nc893C9MraNcOgvLuA+Jtmf1fVSvbXlDihyYy8eX137JpFxA\/EJGJNxtPFleCt48yLhBvy\/Q6yrhAvC3T6yjjAvG2zD3f5VX9w9fNIB\/ihyIy8XnnLj29WdG5G5zIxOcXcrl4LueGJzLxZYnncm54IhNfXs61h3cF5d0HxNsyvY4yLhBvy+z\/iiHbUEQlfnt3jjZ+cKISn93esgDxAxGZ+D1HGReIt2Xu+W6VzLPlESN3gxOZ+Go69mM6d0MTmXjzzE05HfvWUcYF4m2Ze74rOvatfUoQPxCxie8\/yrhAvC3T6yjjAvG2zD3fcT8+EJGJ5358KCITz\/34UEQmnvvxoYhMPPfjQxGb+P6jjAvE2zK9jjIuEG\/L7PtiffFNfUe+fpgG8QMRlfiGbb8e8QMRqfjszU33KOMC8bbMPd8xcheIyMQzcheKyMQzcheKyMQzcheKyMQzcheK2MT3H2VcIN6WeWfEG6r6oYlLfDoz8+02V6xzNzhRiS9XwFzRxgcgKvHmMm59nszbsYLy7gPibZl9XxTiW8N2GeIHY1jx5fMRB+9JU4jffWwS8QMxrPhlWW0fuicN4sMxqPj0yeW8FnrIevWsXh2Oe4j\/tqZvT5rN88\/Lqp4dKuLjZxHfE7E8q9p49qSJj4HEmz1p8nK+ocTHyoAlvtyhwixDzp408TFkVb+9nGNPmvgYXrzDnjSID8ew4nsyXWKlDXqCeFumS6y0QU8Qb8t0iZU26AnibZkusdIGPUG8LdMlVtqgJ4i3ZbrEShv0BPG2TJdYaYOeIN6W6RIrbdATxNsyXWKlDXqCeFumS6y0QU8Qb8t0iZU26AnibZkusdIGPUG8LdMlVtqgJ4i3ZbrEShv0BPG2TJdYaYOeIN6W6RIrbdATxNsyXWKlDXqCeFumS6y0QU8Qb8t0iZU26AnibZkusdIGPUG8LdMlVtqgJ4i3ZbrEShv0BPG2TJdYaYOeIN6W6RIrbdATxNsyXWKlDXqCeFumS6y0QU8Qb8t0iZU26AnibZkusdIGPUG8LdMlVtqgJ4i3ZbrEShv0BPG2TJdYaYOeIN6W6RIrbdATxNsyXWKlDXqCeFumS6y0QU8Qb8t0iZU26AnibZkusdIGPUG8LdMlVtqgJ4i3ZbrEShv0RLX45XY54nLLEcRHxJDiF\/VWE2Y50\/Zeci5HkTboiWbxm+f16rXFWrYXLGkaE\/cQf4tOQLH6fLlePSU+OvzF30l6cl2XepY0jY6BxJtly8s3psibX4HVMVV9TAxY4gsK8WxUEB8Dije6Nx+Zck6Jj48hS3x+HX90XSxbvirebWU6\/B2IH4ihq3qrTJdYaYOeIN4m0yVW2qAniLfJdImVNugJ4m0yXWKlDXqCeJtMl1hpg56oEt\/aYWy\/TIe\/E\/EDgfhhmKb49cUfq+tyc99t8\/w6nZnbMbl4477ab+y497cA8RHgJ\/784U05EmduwqSnr9+7rqzX4hf1DoOr4r7dfCcd8RHgKX5e33dbnpn\/smIb0Ua8GZjvbASP+LjwE29uuCzm5c7BP5nfgIWZZNUSb27FtwZpER8d3uKrEr95\/tnpjakBdqr6prRT1UeKZ1X\/qLj3ZlgmZ0Wrnp5cV3X88rho41etqRcdmYiXx7PEf7ityM0tV3Mf7sHl3PjP3717UfTqe2t6xMeAfxvvD+IjAPHDME3x9wTxEXBf8eWYXetd80EryIzlbq7qZh\/xEXBP8abSL7p39bvmg4ZVMXy7mNcdfcRHgLf4dFbMnzXjsovWu\/qD1mD94ugzM3rPkzRR4Sc+r7abq\/Smq1e9My\/bwfril8Rc559+TFUfEV7iV+1r9M3V2e4787I7WF+Iz2sIcy8P8XFw7xK\/Pq+9V++Kl2qwvnqgpizx24cqEB8B92zj65fuB7u35oo79B8gPibu2avv8Z7tDtYXl3MLqvqYuKf4ZXnzLTdbvatedgfrC\/FNPx\/xEcDI3TAg3ibTJXasRH\/uQ9ndJ1PgmBABiFcK4pWCeKUgXimIVwrilYJ4pSBeKYhXCuKVgnilIF4piFcK4pWCeKUgXimIVwrilYJ4pSBeKYhXCuKVgnilIF4piFcK4pWCeKUgXimIVwrilYJ4pSBeKYhXCuKVgnilIF4piFcK4pWCeKUgXimIVwrilYJ4pSBeKYhXCuKVgnilIF4piFcK4pWCeKUgXimIVwrilYJ4pSBeKYhXCuKVgnilIF4piFcK4pWCeKUgXimIVwrilYJ4pSBeKYhXCuKVgnilIF4piFcK4pWCeKUgXimIVwrilfJ\/Z6qcUXhVmesAAAAASUVORK5CYII=\" alt=\"plot of chunk unnamed-chunk-2\"\/><\/p>\n<pre><code class=\"r\">#mosaicplot(apply(UCBAdmissions, c(1, 2), sum),\r\n<\/code><\/pre>\n<p>\ub9cc\uc57d \ud559\uacfc\ubcc4\ub85c \uc790\ub8cc\ub97c \ubcf4\uba74 \ub2e4\uc74c\uacfc \uac19\ub2e4.<\/p>\n<pre><code class=\"r\">round(prop.table(ftable(Admit ~ Dept + Gender, UCBAdmissions), margin=1), 2)\r\n<\/code><\/pre>\n<pre>##             Admit Admitted Rejected\r\n## Dept Gender                        \r\n## A    Male             0.62     0.38\r\n##      Female           0.82     0.18\r\n## B    Male             0.63     0.37\r\n##      Female           0.68     0.32\r\n## C    Male             0.37     0.63\r\n##      Female           0.34     0.66\r\n## D    Male             0.33     0.67\r\n##      Female           0.35     0.65\r\n## E    Male             0.28     0.72\r\n##      Female           0.24     0.76\r\n## F    Male             0.06     0.94\r\n##      Female           0.07     0.93\r\n<\/pre>\n<pre><code class=\"r\">for (idim in 1:dim(UCBAdmissions)[3]) {\r\n  mosaic(UCBAdmissions[,,idim], shade=TRUE, legend=TRUE, main=dimnames(UCBAdmissions)[[3]][idim])  \r\n}\r\n<\/code><\/pre>\n<p><img src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAWgAAAEgCAMAAACuMBnmAAAA81BMVEUAAAAAADoAAGYAOjoAOmYAOpAAZrY6AAA6ADo6AGY6Ojo6OmY6OpA6ZmY6ZpA6ZrY6kJA6kLY6kNtmAABmADpmAGZmOgBmOjpmOmZmOpBmZgBmZjpmZmZmZpBmZrZmkNtmtv+QOgCQOjqQOmaQZgCQZjqQZmaQZpCQkDqQkGaQtpCQtraQttuQ29uQ2\/+dqOK2ZgC2Zjq2kGa2tma2tpC2tra2ttu225C229u22\/+2\/7a2\/9u2\/\/\/bkDrbkGbbtmbbtpDbtrbb27bb29vb2\/\/b\/7bb\/9vb\/\/\/i4uLklaX\/tmb\/25D\/27b\/29v\/\/7b\/\/9v\/\/\/8OxA93AAAACXBIWXMAAAsSAAALEgHS3X78AAALkElEQVR4nO2dDXfbthWGKdeO7NlpWqtWtsXbZHfpR7pqSjev7WZtWbNGK9mQ+v+\/prgASEIfkAiCvRCh9zknx45EXNOPwUsAAsBkCVhIQp\/AsQDRTEA0ExDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNRP9E56PBNPQ5tKB\/otMkuQ59Di3on+jZye9PHkOfhDu9E52Pzv6TTEKfhTu9E50m1\/noSeizcKdvoos7kTdmPbwd9k10Nnwia3Xo83Cmb6LnlJ9Fnl6EPhFXeiZaZg7R8ujf7bBnorNhoujd7bBnoufac9K7pnS\/ROcjLXjeu9thv0SnZcrIhn27HfZLdHUTLO761pTulegqc1Du6NntsFei+wxEMwHRTEA0ExDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTAQVXbxMkg\/+suWNWe+Weu8lpOji7mxRbN1NBqI7JS0dZ+PkdLrMhk\/Hgz\/RothT2gdMvZiPPhiefB\/wJLsipOiy3mbD60J8nw3P3t3Rl2taiqxfNFYl95vAovNRkkzmomanyUSopZfof\/qLeLGXO4JtI3iNFjLVbifXhui5+KJehOgOSGlTiFTV6KXMIGs1ekkbGkC0N6LV8f3PdypnCK9adJ2j6UWI7gJqRw+eLZbpMBlMyhpdtjrUixAN3IBoJiCaCYhmAqKZCCVa7pprDB5lF4172mrbTPsGgw6hGAkm+vzDxTJ\/ftFG9J5DIdoku3gxXWa3wglV0Im0k48abTlailSHF\/dfiuqtnrDgHKoKueci6YBwov82Wf771cVj\/smU1JCdmegLNuielKLV4cXdE5mHsjahVkJS8V+PcKK\/fVZ8+q1ylt88ak\/iu\/1Fh3JbY314cT9d0j+tyS1UfTaPNCIwKa8EfWlcjk8eU7Wzr3oju\/pz0nLr6nCi\/\/X5j3+Uv+CMfhWyMxICG2w5qmu0PnxFtGuolZDZ+VRfCeWlMZzIPxddG+KN9GxBowRpu+1nw4l+fP3ltXQy0dd70yqYVZcBfTFEu4eqQsocPVkaV4K6NIwLhb6lV1reawOKTgdTnVFFXdKJtUl1MXK0ONwQ7R5qPWR1JZSXhvwbDGT9pp\/ST9Ha8jxJTscT3VRocrkbVgZTs0a7h9oIqWtzdWnIV8WfrL81+rCo5OlUXF8adF3QP+MNiG5PJU9fCfWloZJI1eqA6IMHopmA6O3oBEGdmG2vuwPRO4HobpGdbd39Vn1u3UbUt8RsZfCrPsIBiCaos112v1Wfu+z1JKXoesSq7JW7AdFE2dkWDnWfWwsVqcNo0tUjVu50LTrhw\/Lj2py1ThSy+131uVUzuhZddcvVEa5i2pzXrnj\/50KJ\/us67UXXA1GqD7hao1e75e5DeBBNqBore9m6z13naDmIcmIMfpW9ckcxbc5rV7z+itbd71nVpijuqNVB\/fGPbszBr9khtDr6KZoBiGYCopmAaCYgmgmIZgKimYBoJiCaCYiuT72T8Slr9C6DLfst+le9KiC6PnWIPhTR6tMsIh8ljuN3XqLlWDlRT\/qOWbT8NOuchvzpM1vHT7N8a\/SMpskbPzNm0alcdkNVmj4T4B0mVR9LGB9O9Fl0Hdre6lAfGGZXC9cFAp6iiztZo+t8FYfoPb+v\/ICFVzR9CkHrHOp4UYvOR2pBEX+N3owXs2g5\/YPgz9Gb8SIWXXmWGYS51SFSx9mPxkUUs2i1+1C5ooCzHS3\/tCJfpUdzM2yPf\/NOiI6ueecnZXt0v+KqRh9P8649nTTvjHQF0bboXceDaEt0v+KRdsH9pGyP7lO4Gr1Djt4f3a\/4xqJriLZF9yuO1NE4uk9hpA6H6H7FkToaR+86HkRbovsUzm\/+OYroM8PDFb0tHkRbonuWT1GjG0b3K77xeQ5E26L7FUero3F0z\/Lzte0PIdoW3bP8HDm6YXS\/4sjRjaP7FUeObhzdszxydNPorgXy35rjdTl6hk2jux1efPPi\/Ont7e3YtjszRNuiOx6f6oFR26bWEG2L7lpApQ79fVwT0Q9JdPH5u89uiT\/oof6YJqIfkuj1gVF8lNU4unOJ92+J\/6n\/HNtE9PZ0MlPpWCai13NIOWaTZs9\/I7g8wuZdvRaLY1WWqMLHKrqe588x43\/Lcwrm0TXvtq\/KqleusKxhef2xcTOUPRhzW8M4RG8\/oF6LxbIqa2Y07+TmkXKJ4+bZxieat0bno8uqw6LSyNGI5s3RxadGfPn89KMRXa\/F4liVVbw+\/erh4eEfVStydjQ5WreemVZlbYxA0983tlZH6xg7orsWKP778LBSo9fiQbQlutvhWjNEu0d3O3zLCPRaPIi2RHc7XNTol+Jm+PfxM9Rox+iuBfQI9O8g2jG6a4F89Ozt2x9eInW4RncuMceHs62iuxd5Lxod1sY6RNuidx0Poi3R2xQq3mACjXN09yJvnoscfYVWh2N0x+N\/eCksP7WPEEK0Lbrb4aJnePndrofIQ7Qtutvh1AX\/8DuIbhHdtYDMHYOvrfEg2hK9RZk3YwwquUdvVar4xjqvgw\/Lj2uv4gBFRwlEMwHRTByw6L2ft\/SKAxa93DLjPwzGT26v4pBFB5zxbyJF\/ySJVHTAGf\/HJTrgjP8jE70RD6It0buOB9GW6J7lw209f1yiA249H0B0sGdlBd16nl90yGdlBdx6nl90yGdlBdx6vnvRP5Xs6GSGelbWZrwoRFuPCfasrKi64LtEz2S3LNSzsoLuH80rWhLwWVkBdwnjFx30WVkb8SIWHfJZWctUtODngzh6hg1uhu3xTR2qtXMR\/c3QG+8OC6WtSLrghyxaNTyMjwwh2ha963gQbYnedTyItkT3LB\/TePQhi45qPPqQRUc1Hn3IoqMajz5k0ZGOR3tK2Rq963gQbYnedbyYRadG38xxlLSjh5LVk5ViFk1uy8HR1HUGbUc1ump3xCya0BV5NnjFWaMN3uv7YZ9Fb5sMvE413M+aOpZR9QwbkA2roXdm0TH1DHcy07eh6iNZZtEx9QybUe5bGaRGx9Ez3MvKbI4AOTqWnuF+5nILZ6WYW\/RGvJhF+9DNeb0\/itThhd95ZUPqIIXZBPaYRKsJO+nR5GgfPMc65JydifESRFvwF706qRKiLfiLXp3mGJVo1ybcLjoXHQbjJ3vpWOVwRMeyu0F+84V+UgFNMS\/up2r9lRBNrvX0Ub9f8lBTGi\/56GyRSpFCspBN45HKcil6Vk6IVo8onewJuAlEE\/loIhUv5WOh1ZOhabysEk23\/B2b0DUAognZdJpNaCw0u3pHxmeUDw3RlCQHTotW1oBogkTrGl3cv7paUA1fSR11bUbq8CEfPakmec+Ta5mVs\/OpzhnzE5mjU8fFFKtANJF\/8qJKDHIR8jxJTsdytYr47iPVAfbKHBAtcVwz2AaIJiCai3oOkt4qQg0A18vc1GKdGv3GvHFvDaKJeg6S3ipC3gPPFvV2EfPVhkb5xqxx8wOiS9TAht4qQi72vnmsliJnl+PJ0uiH6zd0k7AJEF1SLzkW1VnX6HJxfXH\/NaWOqh9errqX\/ZhmlRqiFcYcJLlVhMrB5QSD+TXlCqMfrt+gpmDDWg3RJWXLQ24VQQbTk0ddccUXKVr1w1VHvZ7i0SxPQ3SFEqa2MNA1Vqdi\/dQqY1TJ3EcCoptTz0HSW0XoGm0+zlfl6LRs7ck36H8rTzu2A9GSag6S3iqifK78ajva6IfX7eiGPXOIZgKimYBoJiCaCYhmAqKZgGgmIJoJiGYCopmAaCYgmgmIZgKimYBoJiCaCYhmAqKZgGgmIJoJiGYCopmAaCYgmgmIZuIXnxWvuWMpLq0AAAAASUVORK5CYII=\" alt=\"plot of chunk unnamed-chunk-3\"\/><img src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAWgAAAEgCAMAAACuMBnmAAAA7VBMVEUAAAAAADoAAGYAOjoAOmYAOpAAZrY6AAA6ADo6AGY6Ojo6OmY6OpA6ZmY6ZpA6ZrY6kJA6kLY6kNtmAABmADpmAGZmOgBmOjpmOmZmOpBmZgBmZjpmZmZmZpBmZrZmkLZmkNtmtv+QOgCQOjqQOmaQZgCQZjqQZmaQZpCQkDqQkGaQtpCQtraQttuQ29uQ2\/+2ZgC2Zjq2kGa2tma2tpC2tra2ttu225C229u22\/+2\/7a2\/9u2\/\/\/bkDrbkGbbtpDbtrbb27bb29vb2\/\/b\/7bb\/9vb\/\/\/i4uL\/tmb\/25D\/27b\/29v\/\/7b\/\/9v\/\/\/+OcnvkAAAACXBIWXMAAAsSAAALEgHS3X78AAALQUlEQVR4nO2dDXfbthWGKceO4tmJW6lR1sX7oNymH9kqZpvXdrOWNWs0Mhb1\/3\/OcAFQBGXBIgj6koTf5xwfyTRxTT8GLwEIBKMNYCHq+gAeCxDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNBEQzAdFMDEz0MpI8+eOq6yNxZZiio+hp10fiyuBEx\/TyaapeB8QwRW\/SaNLxkbgC0UwMU\/Sn2WjR9aE4MjjRitHQUvRQRUcXXR+JK4MTLavy7Tvk6IeluBjm86Objg\/FkYGK3iQQ\/aCgRjOhRefJ4PrggxNdtO\/Qjn5QtOjRi5+7PhJXBiZ6uEA0ExDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNRKei8zdR9OQve34wuDu9D9Ol6Hx+ssqTfffAQnSrpIXjbBYdLzbZ+Pls9Ce6O\/b4d0K02riePhkfDe422T10Kbqot9l4kov32fjk45xeJuupfKGN9LbDQ2yPjkWvp1EUL0XNTqNYqKVN9J1+ERvX06EtY2Ch8xotZKo76SeG6KV4URshugVSWnwjVTV6IzPITo0WQHQLiFbHz5\/mKmcIr1p0maNpI0S3AbWjRxerTTqmhey06KLVoTZCNHADopmAaCYgmgmIZqJb0dmY2hTlEFJ2Wru\/nY1lh8b+8\/qhWOhY9LMXq8361WkT0Qd2hWiT7PT1YpNdCidUQWNpZz2N6owjFSLV7vnV96J6p7KKO4c6fHq0QNei\/xZv\/v329Gb9xYLUkJ1E9AhrdFIK0Wr3fP5U5qGsSSgdjAo+HF2L\/vEi\/+pH5Wz98kZ7Eu8OFx3Lh4To3fOrxYa+tCy3UFp0Po+Lc0CfFGezoxtxltAW9YPs\/Ouo4fMyuhb9r29\/\/YP8MxP6g8jOtN7Ko7pG690rol1D6WDZs4U+B4qTYhzLfxSdFeIH6cmKRgnES5M\/tWvRN+++n0gnsT7f61XBUrTa3RDtHqrI0fHGOAfUSWGcIvSWtjS8ynYuOh0tdEYVNUon1jqVxsjRYndDtHuoypVVnQPFSSH\/ByNZvyn+kEVry8soOp7FuqlQ53Q33IwWZo12D2UE07V5e1LIreKfNfQa3Re28nQqLk8KOiPoy\/gBRDdnK0+fA+VJoZLIttUB0b0HopmA6P3oBEGdmH3b3YHoe4HodpGdbd39Vn1u3TrUl8SsMuxV7uEARBPU2S6636rPXfR3okJ0OVZV9MrdgGii6GwLh7rPrYWK1GE06cqxKnfaFh11RPm7mxy1ThSy+73tc6tmdCl62y1Xe7iKaXJc98X7XydE29\/dXHQ5BKX6gNUaXe2Wuw\/hQTShaqzsZes+d5mj5fDJkTHsVfTKHcU0Oa774g1XtO5+J9s2RT6nVgf1xz97aQ57JX1odQxTNAMQzQREMwHRTEA0ExDNBETXQbT8dMO5nObh2JKG6BrQYKkaRpKfIz81NtQX0\/IxBSmaOuNl\/0T3ztFhaXzoe0eqiOx8ZczMo0HUyoZ60RsflyXegEVXoxnQyEbhNRvTLb0nEP0QoqsVmCbSMNdoOYRLlLOQgxS9k5KTuIMcndDsbeMCHKTofF5MtNY5o9xQP3rj45Ko0XJjzDxI0brZrOcI0lgqezua\/rWiRpe\/NEzRLeAbM59HNPG+jAfRluhtx4NoS\/S240G0JbpneZE6Tn41WpQQbYvuV1xcDEXjPcXF8HB0v+I0d+d8FXzzrgVaqdFo3tWI7lleNu+MtjtE26K3HQ+iLdH9ij+SLngLtDN6hxx9OLpf8Tt3AUO0LbpfcaSO2tF9CiN1OET3K47UUTt62\/Eg2hLdp\/D65T+nj+EzwzZAjd5z6H0UnaJG14zuV\/zO5AaItkX3K45WR+3onuWXO6vyQbQtumf5JXJ0zeh+xZGja0f3K44cXTu6Z3nk6LrR\/Yqv0TOsG921wPq3OwOjO\/Eg2hLdbff8r6+fPb+8vJzZ1mWGaFt0x\/1TPQItU\/NjmYjeBg1Tx5ZHMRG9DVxTx7cfv7kkfq8+U8FHWbWju+2+OwKNiei1ozuXuP1A\/Fd\/h4nodaO3HQ+iLdGdS2SvfiM4e1TNu\/LeIP2O4WYhkSt2RS9Db96V99jrdxw33e\/2CdPqgyGCFF3ev6nfsdzQ+e7z8mIo1zRMzAVpwxBtvelev2O5RTkpm3eqdgcoeucn5T32+h3HTffr6VnZYZGP9Q5fdCc1Ov9qJzUlyNG1orseTv7u+M\/X19f\/MBo3+Tz4Vsdk2+rQz0B7+Jvu7wz12442INHGTfeM7ej\/XF\/v1Oj9RxuS6BZwHb1TmiHaPbrb7nuG+nfiQbQlutvuoka\/ERfDv88uUKMdo7sW0EP9X0K0Y3TXAuvpxYcPv7xB6nCN7lxiWX44uy8eRFuiuxe5FY0OaxsSom3R244H0ZboTQrl7zGBxjm6e5H3r0SOPkerwzG64\/6\/vBGWn9sHriDaFt1td9EzPPvpvsfHQ7Qtutvu1AV\/8RNEN4juWkDmjtEP1ngQbYneoMz72T2DSh1R\/u7mKnonmqZJ2ybQDJg+ig4SiGaix6L3zfjvPEk3V9Ff0ZvezPgvmx2Biu7LjP\/gRfdlxn\/wovsy4z980XfiQbQletvxINoS3bN8T5aeD150X5aeD150X5aeD150X5aeD150X5aeD1\/0nXgQbYnuVxxd8NrRfQr3Z\/3owEX3Z5Ww4EXfjQfRluie5dMo3ixH6Bkeju5XXN\/meIqL4cHofsXpNn91424RD6It0T3Ly4aHMckDom3R244H0ZbobceDaEt0z\/IYj64b3a84xqNrR\/crjvHo2tH9imM8unZ0z\/IYj64bve14EG2J3nY8iLZE9ym8XRK2nKwUpuhywZlU9YOzsfWeB1v0xsdVPZKi3RGk6HLhRho+E+\/W09iccFgreuPjqnKrf2uQoquLgol3LAsMVgmpZ7h3cjtR9SpqNL\/okHqGdsyFG7PxaKFSx4hTdEg9w\/0k4kJfrcD0TlwMP7\/ir9Fh9Azt7CzcqNYIvef24X2gZ1jrjywWbtRJhCr1Q6\/keCheiKKNlRyXarnsNHrwlRz3cht26mgDv+NSHaQeLAIbuGiVu9Lgc3QbeI51yOVne7BQ92MQXe0gQbQFf9HV5iREW4Doe3BcXf5e2pkfHcTqBnfpjehgWL\/8Tj8SgoY18quFaLeKazwNiJ7eFEvOO4707wDRxHp6skqlSCFZyKaBX2W5EJ1sO+HytIkPBLwLRBM07Jmr0bjlRD+CmwYmt6KpbeU4irQDRBOyjZrEakj0IxlP6MJjiKarkdsA9A4QTZBoXaPzq7fnK6rhldRR1makDh\/W06fb2fTLaCKzcvZsoXPG8kjm6NRxvK4KRBPrL15vE4MQLJd9P57Felz0MzXS4JU5IFri+EFrEyCagGguyqlI+bzouqhuobodSs2haXgdlEA0YTxDNonVVS\/VPcGlMqtfml8RIZqoPrCQSEZv5ffZ2SwuXzxyDEQT5iM4vzZTR371A6UO\/bJx\/ejbAKKJcipSNo6ldi16OZE5Wr\/4XDQhmrjzUNmNEi2+I8P6xSNDQ7TCyNFfmqL146q2T61KbE+uOgxEE8YzZJNK6qg074rhkEZAtKTyUFnVrtsj2qdfA9FMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNBEQzAdFMQDQTEM0ERDMB0UxANBMQzQREMwHRTEA0ExDNBEQzAdFM\/B8rdGsILf0NBgAAAABJRU5ErkJggg==\" alt=\"plot of chunk unnamed-chunk-3\"\/><img src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAWgAAAEgCAMAAACuMBnmAAAA8FBMVEUAAAAAADoAAGYAOjoAOmYAOpAAZrY6AAA6ADo6AGY6OgA6Ojo6OmY6OpA6ZmY6ZpA6ZrY6kJA6kLY6kNtmAABmADpmAGZmOgBmOjpmOmZmOpBmZgBmZjpmZmZmZpBmZrZmkNtmtv+QOgCQOjqQOmaQZgCQZjqQZmaQZpCQkDqQkGaQtpCQtraQttuQ29uQ2\/+2ZgC2Zjq2kGa2tma2tpC2tra2ttu225C229u22\/+2\/7a2\/9u2\/\/\/bkDrbkGbbtmbbtpDbtrbb27bb29vb2\/\/b\/7bb\/9vb\/\/\/i4uL\/tmb\/25D\/27b\/29v\/\/7b\/\/9v\/\/\/8i0buQAAAACXBIWXMAAAsSAAALEgHS3X78AAALo0lEQVR4nO2dDXvbthWFade27M5O0lq1sy7eRjn1+i2mm9d0k9a0WeORNcn\/\/2+GC4AiKQsWQMBXInXOkzyiaeIYegVe4kIgGJUQi6JNV2BXBNBMAmgmATSTAJpJAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkkAzSSAZhJAMwmgmQTQTAJoJgE0kwCaSQDNJIBmEkAzCaCZBNBMAmgmATST+gf63cso2nv+06ar4aq+gS4mkdTedNM1cVTfQCfRwduyvH8T7c82XRU39Qx0Gh3eyY0kOt9wVRzVM9BJFKuN7PS7zdbEVf0CXUz6FjEW6hfofKwjR\/8E0EwCaCb1C3QjRv\/SM+L9Al33OvJxzy6LPQOdjdCP5pHIDEUH+v4GmeETC2MdXCp+\/hijd5BZAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkkAzSSAZhJAMwmgmQTQTAJoJgE0kwCaSQDNJIBmEkAzCaCZBNBMAmgmATSTAJpJAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkkAzSSAZhJAMwmgmQTQTAJoJm0UdHETRR+tWkkm6dt9seu1SdDF5PCuSFbdAwvQQZVWjLPL6GBaZqNnl3t\/Lct5dPAnAVrtzMcfjfZ7d5vsCm0SdNVus9F5Ibaz0eGHCb2c0xoRemfvloswacOg83EUxXPRstMoFmhpF\/2kX8TOfHy0wRoG1MZbtIA5l+sVnDdAz8WL2gnQAZTScjKpatGljCBLLbqklWYA2lui1\/HT7xMVMwRXDbqO0bQToEOI+tF7L+7KdBTtxVWLrnodaidAQ24CaCYBNJMAmkkAzaSnA52NqL9QDw9lJ9a5dDaSyYr59\/ZWW6MnBH38\/K7MX550Ab3mUIBuKjt5NS2zK8GEGmgs6eRjqxUYK5Dq8OL6a9G8U9nEna0s67rmFAqgpwT997j8zzcns\/yzKaEhOonI9iwSkAq0OryYHMk4lHWxsq2r+INk\/nR6StA\/vihe\/6iY5RczzUlsrS9KDeyo1IcX19OS\/msQbla2dZ3RiEBcnSf6xDm93J+JM4n2qF9kZ19E1Xq\/jnpK0P\/+8re\/yLeQUGWJzthuVVHdovXhLdCuVrZ1JbrHU32eVCfOKJYfJp054hfp4R2NEqTd1gp\/StCzN1+fSyaxPt9tm2C2OA3opQHa3cqyrjJGx2XjPFEnTuM0ok3a0\/FK\/KSg072pjqiitejAatMgGjFaHN4A7W5lXVf5ujhPqhNHfgZ7sn1THbYVtKY8j6KDy1h3FWxO98b73ps2W7S7lXVd1R\/UrXlx4si94gPd5hbdJy3g6VBcnzh01tD\/xi8AursW8PR5Up84Kogseh0AvfUCaCYB9GrpAEFJzKr97gLoRwXQYSWTbZ1+q5xb9yD1JTFrDY3VRzgIoEmUbFfpt8q5q5woqkDX41lVVu4mL9AykSL1fX5clWwLhjrn1kBF6Gh06erxLHf5tuiEBnEbn2\/kpg5FbCxdpQOFTL8XObfqRtegF2m5OsJRnqBVztoY4Yn+56LIvYiVpask6PpNqByw3aLbabn7MIsn6GIiW3T9Z\/sLWmfZOueuY7QcYtlvDI1VWbmbfEOHfHZV48rQY9A6\/U4WfQrx3kSvg\/LxTy6aQ2PJNvQ6+gmaQQDNpACh4\/C3Rn8HoA3yvxhmZ3dp3y+GDPLv3gnQve\/eMShIi+59945BQbp3jU7lMEGLnp9+j3M14CDetmNyiF6HhWiwVA8zJHH14pizIAW3fJcqP6Fv49tv2FZhRu8GEaPN41PiMqRn0sh3HMupYbyh48FH22fQRg+KEgo0TRujiTSjWNJ3ce9cL6khhQ4z6LpFSyXx0g4r9871KgcXOowedYyWSuL8c17QuxI6aDBY9TooiBSvZ9Tr4A0dD\/2GCFr3o\/U8QroKut9v4Bc6Lv6lgsf+sGN0CKFFr6j6NoJO0aIt3f2KP7j0ArTJ3a\/4bvQ6QsjXc750cx5Am9w9y88Roy3d\/YojRlu7+xVHjLZ29yyPGG3r7lc8R2Zo6x7aD6AN7qH9ANrg7lN4xUR0gDa5e5Z\/MBEdoA3ufsV346usEMJE9BrFNoMeyET0papvI+gHfgBtcA\/tB9AG9wAec3TvLNx9DdL2GlIAbXL3Ki1vbUya96UDtMndp7DqPwO0lbtXafl0D4C2cvd2SNox2v3GbcciNpadtO2gKWnp++oGUtsPeiDaMtD5H5eGkQajrQJd\/PDq+NnV1dWlKVg4RtPA8XlAMTrVb8e0qLVDFyJyO9zBtSOKrQKtQ4fZD6AN7m6HF19++NsV6c+G6e4AbXJ3O3zF1PMlP4A2uDuXuH9P+q\/JD6AN7qH9ANrg7lwie\/kHoVOEDkd31wLFJALoLu6uBdbkhABtcncu8eZTXAy7uDuXSNC96+TuWiAfnyJh6eLuWoDuhH7MD6AN7q4FijcH397e3v4TLdrR3bXAgzn+S34AbXB3LVD8cnuLFt3B3e1wjRmg3d3dDm\/P8feb8Q\/Qj0i06BtxMfzH5Qvdoj1m\/AP0o9Jz\/D+\/a\/7UacY\/QD+qfPzi\/ftfb3Sw8Jnx3yPQ9ZKZeqveYe3uXJ9568tZjxn\/\/QFdL5mptxpraFq7u1foXnQ6jJ\/mIEHXy7HpraX12azcO9fL4DcI0EtTROr1BPUW1wKDxbvFBBqPpee3GPTSb+olM\/VWvcPe3b1C716KD1uvruez9Hx\/QG+iRf96Iyg\/W4Qnn6Xn+wN6AzFadGtO3zb6zT5Lz\/cHdL1kpt6qd9i7u1WGku7nbxst2GPp+f6AbiyZydePlrFj7\/v1tR0S6ADq4vnucjGMhBTc2r1TqeKH09bo3cBjdAh5evqsEgbQXn4AbXD3LJ9GcTlv3P8G0CZ3v+I6XzrBxXCtu19x9YDyoafgIeTrKTsejbkHAG1yD+0H0Ab30H4AbXD3LL8T49Eh5H0x3IXx6BDyzwx3YDw6hIK06KGPR4dQgBi9A+PRAYRex4qqA\/Sugl7cGV5PVgJok3sYm0W\/A6BN7oF87vX1EKBN7p7lkRnauvsVR2Zo7e5XHJmhtbtfcWSG1u6e5ZEZ2rqH9gNog3sQl3uEjrXuXqWzEX1f2FwEFqBN7j6F1ezVFDHaxt2nMHXr8nHUcaFugLaWBN2+wwCgTe4+hSXo9jRHgDa5+xQGaAd3n8JY3cDBPbTfIEE3b1mR0w1Fv7bVB7Bw71wvg98QQbdu\/Z4LwtQDyI6nTu6d62XwGyLo5m2F2ellXKZEPXFq0gBdV934uIDGjbLF9fcTBZjjXvDH\/PoL2qzGrd\/zczUlXK9UYi+AXqMkio5a94Ir0PnYjTNA26iO0dWqMNnIrc8B0FZq3fpNLdqdM0Bbqb4XXPf1ZMNGr2MbBdBMAmgmATSTAJpJAM0kgH5EjutTPargoN0e8ehyuL1rKG0z6H4qv\/hKPwuahjWK66ka2BegibVetMrvCcYATcrHh3epBCkgC9g041tRrkAnVRKeuieFUgBNysexRFzSOCj9K9WM5AVoGrvze6AxQJPkGGgSqyHRD0Q8oW+cG6Dpa+jmg+adBdAkAq1bdHH9zdkdtfBW6KhbM0KHj\/Lx0WIZnTkNNxPe46mOGfN9GaNTx7Ub2wJoUv7Zq0VgkN9uz6Po4DIm3mLrEzXF0CtyALSU4xetXQTQJIDmUj0VKVUT3NQUe\/lVinjVO4tJ9\/gB0KR6KpIMy0eyidO9ZmqSjN5JP3W+IgI0qb3COfXqaILBxUx3+fROr4wFoEntNfvrFi3zlLjamZ19gdDhp+ZTKLIRwdTfex\/rREbupEkG2RlCh4faLZoeSyEIp3q4TsXpbs+qqAXQpKWnUCRxs4nrWaNJTE9iA2gv1VORNGHdoumn4vWswp4gdPiqnoo0V6N0qXqZt158Rv8BmkkAzSSAZhJAMwmgmQTQTAJoJgE0kwCaSQDNJIBmEkAzCaCZBNBMAmgmATSTAJpJAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkn\/B8li9L2VoCTnAAAAAElFTkSuQmCC\" alt=\"plot of chunk unnamed-chunk-3\"\/><img src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAWgAAAEgCAMAAACuMBnmAAAA7VBMVEUAAAAAADoAAGYAOjoAOmYAOpAAZrY6AAA6ADo6AGY6Ojo6OmY6OpA6ZmY6ZpA6ZrY6kJA6kLY6kNtmAABmADpmAGZmOgBmOjpmOmZmOpBmZgBmZjpmZmZmZpBmZrZmkLZmkNtmtv+QOgCQOjqQOmaQZgCQZjqQZmaQZpCQkDqQkGaQtpCQtraQttuQ29uQ2\/+2ZgC2Zjq2kGa2tma2tpC2tra2ttu225C229u22\/+2\/7a2\/9u2\/\/\/bkDrbkGbbtpDbtrbb27bb29vb2\/\/b\/7bb\/9vb\/\/\/i4uL\/tmb\/25D\/27b\/29v\/\/7b\/\/9v\/\/\/+OcnvkAAAACXBIWXMAAAsSAAALEgHS3X78AAALZ0lEQVR4nO2dDXvbthWFKce27NmpO6t21sbdRiXLumablLbup7UsWaOSDcn\/\/3OGC4AESQsWQcBXonTO0z6iaOEIegNe4oIgGBUQi6JNV2BfBNBMAmgmATSTAJpJAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkkAzSSAZhJAMwmgmQTQTAJoJgE0kwCaSQDNJIBmEkAzCaCZBNBMAmgmATSTBgZ6EUk9uVhuuiauGiboKBrNNl0VRw0OdEwvv0+jo4G16WGCLvKp3hiMBgq6SMfHm62Jq4YKOp8e3G22Ko4CaCYBNJMAmklDBZ1NBta\/Gypo9DoeV+hHM0mDfofM8JFVjXUM7FQ4VNBP\/jKw9jw40MMVQDMJoJkE0EwCaCYBNJMAmkkAzSSAZhJAMwmgmQTQTAJoJgE0kwCaSQDNJIBmEkAzCaCZBNBMAmgmATSTAJpJAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkkAzSSAZhJAMwmgmQTQTAJoJgE0kwCaSRsFnb+Koif\/WvGH+eDu11yrTYLOp0fLfL5qvS+ADqqkZJxeR4ezIh0\/vR79le6OPfxCgFY7s8mT8cFPG6xkKG0SdNlu0\/FlLrbT8dGHKb1cZhP5Qjtpc4NVDKcNg84mURQvRMtOoligpV30Tr+IndlkYOty2LTxFi1gqjvpL2ugF+JF7QToAEpo8Y1EtehCRpBWiy5oSR+A9pbodfz0+1TFDMFVgzYxmnYCdAhRP3p0sSyScTSKyxZd9jrUToCG3ATQTAJoJgE0kwCaSQyg1aKLZpwoPe2cVKdjmbXY\/97datPiAH3yybLInp32Ab3mowBdV3r6fFakN4IJNdBY0skmnZYILEGqj+cvvhbNO5FN3Nnq4e9Zc+QEEAvob+PiP69P77LPZoSG6MxF2tchEylBq4\/n02MZh9I+Vuu\/hzwfTyygf7jIX\/6gmGVXd5qT2FpflFracaE\/nr+YFfS\/JuJmtaaKdzQiEJeHhz5ezq4P7hK1PKf6Q3r+t6jn6rMsoH\/56tc\/y98yp1oTnUm3R6noFq0\/3gDtarX+e9KTmT48yuNlHMt\/QzpgxB+SoyWNEiT9lp9lAX335utLySTWx3vXJphWhwG91EC7Wz38PTJGx0Xt8FDHS+3ooU3a0\/MEzAM6Gc10RBXNRgfWLi2jFqPFx2ug3a06fU91eJTHi\/w3GMn2TV+99aA15UUUHV7HuqvQ5XCvARjN6i3a3arb9+jWXB0vcq\/4dxxEix6AKng6FJvjhQ4W+r\/2B4DurwqePjzM8aKCSNXrAOitF0AzCaBXSwcISmJW7XcXQD8ogA4rmWzr9Fvl3LrjqE+JaWNEzHzCQV6gZf9+kA\/raImS7TL9Vjl3mQpFJWgzjFVm5W7ybdFzGlv0Hj3btMpkWzDUObcGKkJHrUtnhrHc5QlapVK18YbIV94efX6GDhQy\/a5ybtWNNqCrtFx9wlGeoPOpbNFmrCH6zU+Rr0d\/0Ka1qByw2aKbabn76Ipv6MinEY0YG7\/BgtZZts65TYyWIysHtRGxMit3U+hex4BB6\/R7XvUpRCMSvQ7Kxz+9qo+Izdl7Hav8BgmaQQFCx9GvtdMwQFvkfzJMz5fJ0E+GDPLv3gnQ9e4dQK9WkBY9+O4dg4J072p9HYC2CL2OLhI9P92Y9DQPNcLkouAp+C6CpsFSNZ4jLzHriR8nTml4mNG7HY\/R1JBMfkLJoZwf69Skg7Tout9wQdvHp8T5vjYzT49VOg7iIXQYFFYPyhNKrulYjdyp4TQH9971KnYtdNhBN1u03MomjpN8ETpWVL3t0YzRFJzlJRk39971svjtImgKEyoy6yDiztkzdFz9rILHwW7HaN2P1tMHR3LmX+TYkUaLXlH1x+gi+nome9GiA8j3ZNjuTAK0zd2v+H70OkLI13PR6k4CtM3ds\/wCMbqju19xxOjO7n7FEaM7u3uWR4zu6u5XPNuPzDCAkBmuqDpA7yvoFRPRAdrm7lm+PREdoG3ufsX341JWCGEiukGxzaB3ZCJ6q+rbCPqeH0Bb3EP7AbTFPYDHAt27Du6+BklzRSOAtrl7lZZ33DXmoAG0zd2nsOo\/A3Qnd6\/S8qETTdC+8vboj2KLQZPm\/mtHbYe2HTQlLUNf3UBq+0HviACaSVsGOvtTa7yu5RdAfjb9UWwT6Pyb5ydPb25urm1R2bd795tvD29HQOtJjfZFrQHa5u5aQIUOux9AW9zdPp5\/9eHvN6QvLQuDALTN3e3jK+b4t\/wA2uLuXOLje9L\/bH4AbXEP7QfQFnfnEumzPwidIXQ4ursWyKcRQPdxdy2w5uEFAG1zdy7x5o84GfZxdy4xR\/eul7trgWxyhoSlj7trgfwlYnQvd9cC+ZvDf97e3n6PFu3o7lrg3s0ULb+dBG0Wr9JbZkdnd9fq5P+9vd23Fm0Wr9JbZkd3d7fKaMwa9CPM+N9K0GZhFL3VWimlk7tbZe6hDT3jfztAt66KmaV+9FZz7Z9u7m6VES36lTgZfnd9oUNW6Bn\/2wG69RezeJXeqq9m1dXdtToa7ecKdPAZ\/1sJegMtmoLHxfv3716VoSP0jP+tBL2BGF2U6xns1cVZs3iV3jI7uru7V+ij6HRY+5A7Cbq2eBVfP7ql0EvPbyfoAOrlmb8tJ9AEX3oeoI3ePhMx+rzq3oVdeh6gtd69EpSfVifc4EvPA7SSOAmc\/VhPUEIvPQ\/QSpSCf\/LjA9cNAdrm7lpAxo7Rv\/W7\/UjBQ6iP59trPaj0COtHA3RD+TdnjUfD1\/wA2uIe2g+gLe6e5ZMoLhYjZIbr3f2K68HDU5wM17r7FVfPzUYK3sHds7zseNQuiQO0zT20H0Bb3EP7AbTF3bM8xqO7uvsVx3h0Z3e\/4hiP7uzuVxzj0Z3dPctjPLqre2g\/gLa4h\/YDaIu7T+HqhmUzWQmgbe5hbKp+B0Db3AP5fNTnQ4C2uXuWR2bY1d2vODLDzu5+xZEZdnb3K47MsLO7Z3lkhl3dQ\/sBtMU9iMtHhI617l6l0zFdL6wvAgvQNnefwupWjgQxuou7T2Hq1mWTKOSK6AC9ShJ083Y7gLa5+xSWoJvTHAHa5u5TGKAd3H0K78vqBiGEhGVF1QF6U6DNjbJJefg63goO0Kuq3vYwC84QXrmVWBc7srn3rpfFbxdBNxczoK356DVadO+qW5\/L0Fyew7RtJ\/fe9bL4DRe0XfUFZ9KxupEEoAP\/oHkUHTdb9L3bSbqBCVyvnQNNai04I5\/CBtCBfxDJLDhjgghAB\/5BUmYFmkX5QFKADvyDQgmgmQTQTAJoJgE0kwCaSQD9gFy7cA8pOOgA8rMJ+GO2GfQwlV39o8xDzpdF\/mKWjqMoJtDEWi+V6feoaIAmZZOjZSJBCsgCNk2tV5RL0PMqCZeHTbzG8L4AmpRNYolYaHFJ\/xVq6ncFmsY31jy8ao0AmiTHieaxGhL9QMTndG2wBpqu94+cVuZuCaBJBFq36PzF6\/MltfBG6DCtGaHDR9nkuFqvaBFdyqicnsx0zFgcyBidOK4Y3RRAk7LPnleBQQCWy74fXsd6XPRTNZfTK3IAtJTjgxH6CKBJAM0lMxWJnhUhYrK6l6F8p5PE6l0PATSp9uwreeVVNnG6qW+uuxdqYtLcvbNRCaBJ5jK37uTJCQZXd\/qdnphUvuslgCaZiRsyM4nLFq3fFUU53tGnB60E0KTaLIITlbro694nZSKje9b9WzVAk1rPvprHxDTR573mfJm+cRqgSfemItVn2wF0ODWnIuUv73SL1u8KBdq86yGAlro3FSlRL82JSQuPETyAZhJAMwmgmQTQTAJoJgE0kwCaSQDNJIBmEkAzCaCZBNBMAmgmATSTAJpJAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkkAzSSAZtL\/AQembkBUTIqeAAAAAElFTkSuQmCC\" alt=\"plot of chunk unnamed-chunk-3\"\/><img src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAWgAAAEgCAMAAACuMBnmAAAA6lBMVEUAAAAAADoAAGYAOjoAOmYAOpAAZrY6AAA6ADo6AGY6Ojo6OmY6OpA6ZmY6ZpA6ZrY6kJA6kLY6kNtmAABmADpmAGZmOgBmOjpmOmZmOpBmZgBmZjpmZmZmZpBmZrZmkNtmtv+QOgCQOjqQOmaQZgCQZjqQZmaQZpCQkDqQkGaQtpCQtraQttuQ29uQ2\/+2ZgC2Zjq2kGa2tma2tpC2tra2ttu225C229u22\/+2\/7a2\/9u2\/\/\/bkDrbkGbbtpDbtrbb27bb29vb2\/\/b\/7bb\/9vb\/\/\/i4uL\/tmb\/25D\/27b\/29v\/\/7b\/\/9v\/\/\/+wcmXmAAAACXBIWXMAAAsSAAALEgHS3X78AAALP0lEQVR4nO2dDXvbthWFKdeO4tlpWqtR2sVrabdZ12yV0tX92KymTRuNXET9\/79TXAAkQdkIAQK6FJlznuSRRBJH0GvwEoDAq2QLsSjpuwLviwCaSQDNJIBmEkAzCaCZBNBMAmgmATSTAJpJAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkkAzSSAZhJAMwmgmQTQTAJoJgE0kwCaSQDNJIBm0sBAr5JSD\/quiqcAmkmDA532XYWOAmgmATSTAJpJgwOtNVn0XRVPATSTBgcaoYNFAM0kgGYSQDMJoJlUz3Uc3fZdFz8BNJMGBnq4AmgmATSTAJpJAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkkAzSSAZhJAMwmgmQTQTAJoJgE0kwCaSQDNJIBmEkAzCaCZBNBMAmgmATSTAJpJAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkkAzSSAZhJAM6lX0MXzJPngX\/fsWA7tLsJ29Qm6uDpZF8v7Mm8AdFRlJeN8nhwvtvn00XzyBd0de\/xXAVpt3Mw+mB793GMlY6lP0GW7zacXhXieT0\/eXNHDxWYmH2gjPe2xivHUM+jNLEnSlWjZWZIKtLSJXukHsXEzG1omQYt6b9ECprqT\/sIAvRIPaiNAR1BGyTcy1aK3MoLstGghgI4g0ev4+f9XKmYIrhp0HaNpI0DHEPWjJ4\/X22yaTNKyRZe9DrURoCE\/ATSTAJpJAM0kgGbS3kDnU+ou1LND+anzUDqfyrGKfb+71eFof6Affrjebp6edgHdcihAm8pPny22+aVgQg00lXQ2M6d0dSVIdXhx\/Y1o3pls4t5WrpVtOYciaI+g\/51uf3lxerv5ZEFoiM5SDPYcxh8laHV4cfVAxqG8i5VzZcU7kvv+tEfQPz4uvvxRMds8udWcxLP2olP58x\/68OJ6saX\/moOflXNlb2lGIC1PFH3mnM2PbjOVNFLtyM+\/SjpmR90j6P9+\/cfn8hMsqa5EZ+aWyly3aH14A7SvlXNlie7DhT5RyjNnmsq\/Jp06Ykd2sqZZAvHQ5S32CPr25TcXkkmqz3fXJphXpwE9GKD9rVwrK2N0ujVOFHXmGOcRPaUtHS\/F+wSdTRY6oorGogOrS3swYrQ43ADtb+VeWflYnSjlmSP\/BhPZvqkSBwpaU14lyfE81V0Fl9Pd+NiThdmi\/a3cK6veUbfm6syRW8Vf9IBb9KBUwdOhuD5z6LSh\/8aOHkDLM22IyZx3VcHTJ0p95qggUvU6+mvRS+rlR+zRjlaBoFVQi9wHGKUCQRdXskXHu\/4finSAoEHMfdv9FRo6iqsB\/oijuw4H9DgkB9t6+K3G3LoLqS+JeWNurD7CQwBNosF2OfxWY+5yUJSUoOsJrXJU7qcIoePkj71Oe3GoHGwLhnrMrYGK0GF06eoJLX+FXwzz87UxGE685F3AzdNbOlDI4Xc15lbd6Bp0NSxXR3gqvHsnQBvdu+R\/Hkp8Czh6ekuCrj+EGgM2W3RzWO4\/zxKlRRvdu+GC1qNsPeauY7ScYzky5sbKUbmfonTvjHcdMGg9\/F5WfQrx2USvg8bjHz0x58aWh9DrGCZoBsUeggO0RXFm74YeoxkUpUWbfgB9vxA6mITQwSSEDqePOdNtSd3CRDN6vtN46N45iCZL62kkOVjJfL++CwsdT\/6jgsfRuGM0nbf114pyIc\/kBVp056pb56fy83W9Mk83be7QkY2nRSdWDwoWJejykRn0nbnZUYI2W3Q5ncQNekS9DjtoM0Yv9Spq7tCx2lm9PUrQ9F2\/Ds20OE2KHfT7EKN1P9pYXYoYvR\/QEYQYfU\/VDxD0+xGjYyi0RY9oZHjQoO\/6AbTFPbYfQFvcQwrfsxAdoG3ugeV3F6IDtM09rPiYvso6aNB3FqIDtM09sPzuQnSAtrnH9gNoi3tsP98ltn4F3Dw7Vv3gQa8Gf5+h1IGDzmJmGehVhwxarn1fpu0HDkEHDFr1nwHayT2otEz\/3ADtc8mKfh0c98Vw2YjRjr21xONYP40XNA1avCeVADrYD6At7rH9ANriHtsPoC3usf0A2uIe2w+gLe6x\/QDa4h7bD6At7rH9ANri7ltg8+k701UBtM3d7\/Diu2cPH11eXs5tM9AAbXP3PD7TUze2pNYAbXP3LaBCh90PoC3ufocXX7\/5+yXpb5bEIABtc\/c7\/J47C3f8ANri7l3i7WvS7zY\/gLa4x\/YDaIu7d4n86V+EzhA6PN19CxRXCUB3cfct0JLCGKBt7t4lXn6Mi2EXd+8Sy7p7133FP0C3ajM7MwcsHVf8A3Srii\/NGN11xT9At6p4efzPm5ubH1SL7rriH6BbtXMLZ8cV\/wDdquLXm5u6Rd\/1A2iLu9\/hGjNA+7v7HX6nQ7ebeh6gbe5+h4sW\/VxcDL+fP64uhjup5wHa4u5bQHfoPluXrzqlngfoVm1mj1+\/\/u25Dh1dU88DdLtWjS9nO6aeB2gHvRWdDutPCYwTdJUEtnxWb3B271wv9b7vxRC8TgKrnzWzwrq5d6lS8WrenL0beYyuEwzqZ42ssI7u\/hV69VSgPa96HbbaDhn0zg1edcpM\/ayRFdbR3bM2vz0XFXhk\/2OOA\/TOnjoJrH5mZoV1dferjAgWZz81fu0oEdFqMvKRYQ8tmqLyhz\/VoPXbn477YthLjJaxY\/KtfqV+wXLsQ\/A6Cax+ZmSFdXbvUKNX83pSSXY8jPVhowRtJoHl7UcX32Fdh697bD+AtrgHlsd8tKt7WHHMRzu7hxXHfLSze1hxzEc7uweWx3y0q3tsP4C2uMf2A2iLe0jh6taherESQNvc49hU\/Q6AtrlH8nmrr4cAbXMPLI+Roat7WHGMDJ3dw4pjZOjsHlYcI0Nn98DyGBm6usf2A2iLexSXtwgdre5BpfMpfV+IJLAu7iGF1ZfBGWK0i3tIYerWbWZJavoBtMU9pLAE3VywA9A295DCEnRzmSNA29xDCgO0h3tIYWQ38HCP7QfQFvfYfgBtcY\/tB9AW99h+AG1xj+0H0Bb32H4AbXGP7QfQFvfYfqMEbS7wl3eTqGlLL\/fO9bL4jRF040bZVZIK8KmZSMrJvXO9LH5jBG3ehJWfzdMtww2drX4DBm39aUSDa3H97RVAh4G2y7hRdnVBcUSGDr+f2gXoFi2T5EHjztnyYvjxNUDHVh2jjawwLWmH74CJXKdRgm7cKCtDh2jeHPk63uU3RtDGnbO6H50lzBlo7vqNEnQEATSTAJpJAM0kgGYSQDMJoN8hz2w+71R00I7yOdZPET\/MIYMepjZP\/pGoOSKa1iiuF\/k0SVICTax1ih\/Pmf4dATRpMztZZxKkgCxg0w19inIJelkOwtVPlKYthncF0CSa9iyu9Two\/duqG84q0DS54TmLtCOAJsk50GWqpkTfEPElLSg0QNMqQ78J6B0BNIlA6xZdXL84X1MLb4SOujUjdIRoM3tQZUlcJRcyKucPFzpmrI5kjM485+uaAmjS5pNnVWAQgOUE\/\/E8Jd7i2UfqDpKgyAHQUp5ftHYRQJMAmkv1UqSsvH9Bxw155dO76WXXYQtAk+qlSBIvPctKpOIaWO5e+nc2KgE0qZkPWo4EJy\/U63KoIjYWfgsMmgJoUnPlUd229Qu9W45aujZqgCaZOfvzqZ5eUinQaaveTR2\/zq0aoEnNFm3m05ejFHN31zgN0KSdnP0Spnotf7zY3A3QIaqXItVBRJJVkULvpn3N3zb2EEBL1UuRVuUsnQStA0bdj+48DgdoJgE0kwCaSQDNJIBmEkAzCaCZBNBMAmgmATSTAJpJAM0kgGYSQDMJoJkE0EwCaCYBNJMAmkkAzSSAZhJAMwmgmQTQTAJoJv0JdzTpfcI7VLEAAAAASUVORK5CYII=\" alt=\"plot of chunk unnamed-chunk-3\"\/><img src=\"data:image\/png;base64,iVBORw0KGgoAAAANSUhEUgAAAWgAAAEgCAMAAACuMBnmAAAA6lBMVEUAAAAAADoAAGYAOjoAOmYAOpAAZrY6AAA6ADo6AGY6Ojo6OmY6OpA6ZmY6ZpA6ZrY6kJA6kLY6kNtmAABmADpmAGZmOgBmOjpmOmZmOpBmZgBmZjpmZmZmZpBmZrZmkNtmtv+QOgCQOjqQOmaQZgCQZjqQZmaQZpCQkDqQkGaQtpCQtraQttuQ29uQ2\/+2ZgC2Zjq2kGa2tma2tpC2tra2ttu225C229u22\/+2\/7a2\/9u2\/\/\/bkDrbkGbbtpDbtrbb27bb29vb2\/\/b\/7bb\/9vb\/\/\/i4uL\/tmb\/25D\/27b\/29v\/\/7b\/\/9v\/\/\/+wcmXmAAAACXBIWXMAAAsSAAALEgHS3X78AAALYUlEQVR4nO2dgX\/auBXHTS4pzZJeb+FKtzXbnHTd7boV97bs7raw7rors1fM\/\/\/vTE+SLRsQWPbLMzi\/36cNxlg\/xBf5IcniEa0gEUV9V+CxCKCFBNBCAmghAbSQAFpIAC0kgBYSQAsJoIUE0EICaCEBtJAAWkgALSSAFhJACwmghQTQQgJoIQG0kABaSAAtJIAWEkALCaCFBNBCAmghHRvoeWR1ct93VcIE0EI6PtBx31VoJ4AWEkALCaCFdHygjc4WfdckUAAtpOMDjdAhIoAWEkALCaCFBNBCAmghATS0WwAtJIAWEkALCaCFBNBCAmghAbSQAFpIAC0kgBYSQAsJoIUE0EICaCEBtJAAWkgALSSAFhJACwmghQTQQgJoIQG0kABaSAAtJIAWEkALCaCFBNBCAmghAbSQAFpIAC0kgBYSQAsJoIUE0EICaCH1Cjp\/E0Vf\/GXLA8mxpabarz5B5zdnizwZzTYfAWhWpQXjbBqdzlbZ+Nl09Hv6cuzpbxRos3M5+WJ88mOPleRSn6CLdpuNr3K1nY3PPt3QzdVyom9oJ232WEU+9Qx6OYmieK5adhrFCi3tonv2Ru1cTp70WENG9d6iFUyTJ+mqAnqubsxOgGZQSpkKUtOiVzqCrLVoJYBmkOp1\/Pi\/GxMzFFcL2sVo2gnQHKJ+9Oj5YpWOo1FctOii12F2AjQUJoAWEkALCaCF1Ak0DTeOMj15D+raopMr9We+o2eQjelBN0uUnTd+U7KxHrP4H29u1b86gl6+uC\/\/blf29MvFavnyvA3oPYc+JtD5jW7RO5IMZ+evZqvsWjGhBhprOirkNAk2BUhzeH77rWreqW7iwVb7nmnPucOgrqEjv1E13DWmyM7\/Gq\/+9fb8fvn1jNAQHRVvdkUbV9QgNIfnN090HMraWDV5JnJ9OD14ryM7\/\/55\/vp7w0zFGMtpV7Qpi471m2gPz29nK\/pveYRZ7a3kPbWZuDhB7BlzMT25T81HvXkgu\/xD1DKVpADof37z8+\/0K0mozkSHOivbLqxsFDXvjjm8BjrUqskzZU9n9gQpzphxrN9FOmXUA+nZgmYJ0nbJ2BlCx9nPu845VeH3315pJrE935s2waw8DfQzOdDhVvueScfoeFU5QcwZUzl\/aJP2tPwI7v5hmF0udr3Jql7paGYjqmo0NrA2aReVGK0Or4AOt2r4TOUJUpwx+j0Y6fZNT94faHrfLxc7u3caMdVuHkWn09h2FZqc7pWXP5pVW3S4VdNnsq25PGP0XvVOHkiL3tW9OwqV8GwodmcMnS70v\/JAL6BN9+7YOTt49gRxZ4wJImWvoz\/QUEM9+BAcMuKZvTv62LEhGyBoELNtf7hYWvRwdWCgjx63Hmzb4bcZc9uuo\/1IzGpzYu6IACF0kGiwXQy\/zZi7GAxFBWg3kVWMysOE0EEqBtvq5dgxtwWqQkelS+cmssLF3b2Luqhr+cIkWDZQ6OF3OeY23WgHuhyWmyNCwbSpV6Hli3+Y4OHm3qP\/tlfUsXxpEiwN2p2cZgxYb9H1YXn4\/Ap7iz5a0HaUbcfcLkbruZWTypxYMSoPBNOmXhWlg2nRxfA7KfsU+Q31Omg8\/tWL6pxYIt7r2HL55zhBC4i71wHQHnWt13zt0jFAe9QZ9CBitIDYY3THLnCn8oXJIQojQyFxx+hhSvX8yo6zntBLgxd2dm3R6yPDQYrYltNIczuHFzivdFBzHWWQPbAYTQGyGJ9kF1MzRS07YNn069pr6Nr3aP+C\/O9ZdrkoPvbz23f2WoBki96yEP2YQXs9aGbDgp5fmYsu2ThwAo97IfogQbsWrbaKq1uB09Lcl7IGCdrFaPvFab03CVpWyr0QfZCg6VWWZy216Eosaezeul7F00a1heiDBG370aZR69AxD14sjF7Hlqo\/RBcRoLdU\/VBBz4feveNQZ8+0Hq0A2ufeqbS+8F7r5gC0z71LYdN\/BuhG7p1K69xTAN3IvbNDghjdyJ3BI79Br2O\/O7cfQHvcuf0A2uPO7QfQHnduP4D2uHP7AbTHndsPoD3u3H4A7XHn9gNojzu3H0B73Ln9ANrjzu0H0B53bj+A9rhz+wG0x53bD6A97tx+AO1x5\/YDaI87tx9Ae9y5\/QDa487tB9Ae99ACy1\/tzDkD0D73sMPz7149fXZ9fT3Vl2Mfy4p\/DoV6phZt8a23R7Hin0MtQ0d571Gs+OdQaOj45tMfr0m\/NWv8H8eKfw4Fem4k93kcK\/4ZFO75+SPpPz4\/gPa4c\/sBtMc9uET28hdKFy501FLPA7TPPbSAIlsBvZF6HqB97qEF6mPCjdTzAO1zDy7x\/peVD8ON1PMA7XMPLpFsdu8q2fYA2uceWmA5uagMWDb9ANrjHlogf40heCv30AL5+9M\/393d\/Z1a9Jb80QDtcw8tUM+ixJrJ8XBBu+RVdquazaqhe2h18n\/f3RUtemdtBwTaJa+yW7VsVg3dwypjMTvQaRTTD9Jvqe2AQLvEKHarms2qqXtYZdavqdgEOOcD+zD0Jq+yW5VsVo3dwyqjWvQb9WH4t+nzYj46XtXygw8D9NojLuGM3ZLIQGM7dL+2aHUTr+QXHCToHlo0kX3+8eNPb3y5GwcJuocYvaonytpV2wGBdsmr7FYtm1VD9\/AKfVadjjImP475aJe8Sq4fXRfmoxu7tymUf5i6kSHmo5u5hxf58FLF6MtyuQHmo5u5Bx7\/0xtF+Zn7wMV8dFP3sMPVh8DFD7sSzgO0zz3scBqffPkDQLdwDy2gY8fond4u1y25xUoA7XNvUebDdCOrf9nGAdrn3qpU\/t3FWvD4bD8PAdrn3rH84xgZMggjQ4fikEFjZNjYvVtxjAwbu3csj5FhU3duP4D2uLO4fEbo2OveqXQ2ppELksA2ce9S2FzRSRGjm7h3KUzduuUkiqt+AO1x71JYg65fdQdon3uXwhp0fcoUoH3uXQoDdIB7l8LIbhDgzu0H0B53bj+A9rhz+wG0x53bD6A97tx+AO1x5\/YDaI87tx9Ae9y5\/QDa487tB9Aed24\/gPa4c\/sBtMed2w+gPe7cfgDtcef2A2iPO7cfQHvcuf0A2uPO7QfQHnduP4D2uHP7DRJ09YuyJqFD6FfBAXpb1dc9agln5nodS7r+5ZK97q3r5fEbIuhqMoPsYqpAJ6O3aNGtq76h4pFKeo789h1CR0fQflUSzsyvEKMfBnQSRU9qGWgA+mFAk1yMdllhAJr5BZFqCWfQoh8OdCUDDUA\/KGgGAbSQAFpIAC0kgBYSQAsJoHcotAu3S+ygu6ujC+OLOWTQx6nliz9FJts4TWvkt7NsHEUxgSbWNlVm4AT0mgCatJycLVINUkFWsCmnjqFcgE6KQbj5idI4+CkAmrScxBrxiuZB6d\/K5HwpQdPc3Y4kdA0E0CQ9B5rEZkr0ExFP6Dt9FdD0Rb9RUGbuNQE0iUDbFp3fvr1cUAuvhQ7XmhE6umg5eVL+UME8utJROXs6szFjfqJjdBqYMbougCYtv35VBgYFWE\/wn05j4q22vjJJHDpFDoDWCvxhhDYCaBJAS8ktRaKrgiompyaTgL1nr6uYYUw7ATSpshQpKS9U0Y6k4Errk6jd6wjeRgBNcpe5bSePpHaU9\/T6pJTeiqRlkwZoklu4oUcmhqVq0cW9cn1S+2gO0CS3FIlCg27H2Xg0K++V65No4UE7ATRp7bevTHgodiRxuT5pOWnLGaC11n77ysbh8qZYn5SN2\/Y5ANrILUWiIJK\/vq\/8nJ75NWNq0V04A7SRW4o0N7N09RvbAWw5n6QF0EICaCEBtJAAWkgALSSAFhJACwmghQTQQgJoIQG0kABaSAAtJIAWEkALCaCFBNBCAmghAbSQAFpIAC0kgBYSQAsJoIUE0EL6Pzr95Jm37U0KAAAAAElFTkSuQmCC\" alt=\"plot of chunk unnamed-chunk-3\"\/><\/p>\n<p>\uc804\uccb4 \uc9c0\uc6d0\uc790\ub97c \ubcf4\uba74 \ub0a8\uc131\uc758 \ud569\uaca9\ub960\uc774 \uc5ec\uc131\uc758 \ud569\uaca9\ub960\ubcf4\ub2e4 \ub192\ub2e4. \ud558\uc9c0\ub9cc \ud559\uacfc\ubcc4\ub85c \ub098\ub220 \ubcf4\uba74 A \ud559\uacfc\uc758 \uacbd\uc6b0 \uc5ec\uc131\uc758 \ud569\uaca9\ub960\uc774 \ub192\uace0, \ub098\uba38\uc9c0 \ud559\uacfc\ub294 \uac70\uc758 \ube44\uc2b7\ud558\ub2e4. \uc5b4\ub5bb\uac8c \uc774\ub7f0 \uc77c\uc774 \uc77c\uc5b4\ub0a0 \uc218 \uc788\ub294\uac00?<\/p>\n<p>\uacb0\ub860\ubd80\ud130 \ub9d0\ud558\uc790\uba74, <strong>\uc131\uc740 \ud559\uacfc \uc120\ud0dd\uc5d0 \uc601\ud5a5\uc744 \ubbf8\uce58\uace0 \ub2e4\uc2dc \ud559\uacfc\uc5d0 \uc758\ud574 \ud569\uaca9\ub960\uc774 \uc601\ud5a5\uc744 \ubc1b\ub294\ub2e4. \uadf8\ub9ac\uace0 \uc5ec\uc131\uc740 \ub0a8\uc131\uc5d0 \ube44\ud574 \ud569\uaca9\ub960\uc774 \ub0ae\uc740 \ud559\uacfc\uc5d0 \ub300\ud55c \uc9c0\uc6d0\uc728\uc774 \ub192\uc558\ub2e4.<\/strong><\/p>\n<p>\uc0ac\uc2e4 \uc5ec\ub7ec \ubc88 \uc598\uae30\ud574\ub3c4 \uc774\ud574\ud558\uae30 \uc5b4\ub824\uc6b4 \uac1c\ub150\uc774\ub2e4. \ub2e4\uc2dc \ud55c\ubc88 \uc124\uba85\ud574\ubcf4\uc790.<\/p>\n<p>\uac01 \ud559\uacfc\ubcc4\ub85c \uc0b4\ud3b4\ubcf4\uba74, \uc131\ubcc4\uc5d0 \ub530\ub978 \ud569\uaca9\ub960\uc758 \ucc28\uc774\uac00 \uc5c6\uc5c8\ub2e4. \ud558\uc9c0\ub9cc \uc5ec\uc131\uc758 \uacbd\uc6b0 \ud569\uaca9\ub960\uc774 \ub0ae\uc740 \ud559\uacfc(\uc778\ubb38\ub300, \uc0ac\ud68c\ub300\uc758 \uacbd\uc6b0 \uc815\uc6d0\uc5d0 \ube44\ud574 \uc9c0\uc6d0\uc790\uac00 \ub9ce\uc544 \ud569\uaca9\ub960\uc774 \ub0ae\ub2e4)\uc5d0 \ub300\ud55c \uc9c0\uc6d0\uc774 \ub9ce\uc544\uace0, \ub0a8\uc131\uc758 \uacbd\uc6b0 \ud569\uaca9\ub960\uc774 \ub192\uc740 \ud559\uacfc(\uc790\uc5f0\ub300, \uacf5\ub300 \ub4f1)\uc5d0 \ub300\ud55c \uc9c0\uc6d0\uc774 \ub9ce\ub2e4. \ub530\ub77c\uc11c \ud559\uacfc\ubcc4\ub85c \uc131\uc5d0 \ub530\ub978 \ud569\uaca9\ub960\uc774 \uc5c6\uc5c8\uc9c0\ub9cc, \uc5ec\uc131\uc758 \uacbd\uc6b0 \ud569\uaca9\ub960\uc774 \ub0ae\uc740 \ud559\uacfc\uc5d0 \ub300\ud55c \uc9c0\uc6d0\uc774 \ub9ce\uc544\uc11c \uc804\uccb4 \uc5ec\uc131\uc5d0 \ub300\ud55c \ud569\uaca9\ub960\uc744 \uacc4\uc0b0\ud574\ubcf4\uba74 \ub0a8\uc131\uc5d0 \ube44\ud574 \ub0ae\uc558\ub358 \uac83\uc774\ub2e4.<\/p>\n<p>\uc774\ub294 \ub2e4\uc74c\uacfc \uac19\uc774 \uadf8\ub798\ud504\ub85c \uadf8\ub824\ubcfc \uc218\ub3c4 \uc788\ub2e4. \ud559\uacfc\ub97c \uc0dd\uac01\ud558\uc9c0 \uc54a\uace0 \ub0a8\/\ub140\uc640 \ud569\uaca9\ub960\ub9cc \ubcf8\ub2e4\uba74 \ub2e4\uc74c\uacfc \uac19\uc740 \uc778\uacfc\uad00\uacc4 \ubaa8\ud615\uc744 \uadf8\ub824 \ubcfc \uc218 \uc788\ub2e4. <\/p>\n<p><img src=\"http:\/\/141.164.34.82\/wp-content\/uploads\/2019\/07\/SimsonDia01.png\" alt=\"\uc798\ubabb\ub41c \ubc84\ud074\ub9ac\ub300 \uc785\ud559 \uc778\uacfc\uad00\uacc4 \ubaa8\ud615\"\/><\/p>\n<p>\ud558\uc9c0\ub9cc \uc774\ub294 \ud559\uacfc\ub97c \uace0\ub824\ud558\uc9c0 \uc54a\uc558\uae30 \ub54c\ubb38\uc774\uba70, \ud559\uacfc\ub97c \ubaa8\ud615\uc5d0 \ucd94\uac00\ud558\uc5ec \uc2e4\uc81c \uc778\uacfc\uad00\uacc4 \ubaa8\ud615\uc744 \uadf8\ub824\ubcf4\uba74 \ub2e4\uc74c\uacfc \uac19\ub2e4. \ud559\uacfc\ub294 \ud569\uaca9\uc728\uc5d0 \uc601\ud5a5\uc744 \uc8fc\uace0, \uc131\ubcc4\uc774 \ud559\uacfc \uc120\ud0dd\uc5d0 \uc601\ud5a5\uc744 \uc8fc\uc9c0\uba74, \uc131\ubcc4\uc774 \ud569\uaca9\uc728\uc5d0 \uc601\ud5a5\uc744 \uc8fc\uc9c0\ub294 \uc54a\ub294\ub2e4.<\/p>\n<p><img src=\"http:\/\/141.164.34.82\/wp-content\/uploads\/2019\/07\/Simpson_Dia02.png\" alt=\"\ubc84\ud074\ub9ac\ub300 \uc785\ud559 \uc6d0\uc778\uad00\uacc4 \ubaa8\ud615\"\/><\/p>\n<p>\uc774\ub294 <strong>\ub0b4\uc0dd\uc131<\/strong>\uc5d0\uc11c\ub3c4 \uc0b4\ud3b4\ubd24\ub4ef\uc774 \ud1b5\uc81c\ubcc0\uc218\ub97c \uc801\uc808\ud788 \ud3ec\ud568\ud558\uc9c0 \uc54a\uc544\uc11c \uc624\ucc28\uc640 \uc124\uba85\ubcc0\uc218\uc5d0 \uc0c1\uad00\uc774 \uc0dd\uacbc\uae30 \ub54c\ubb38\uc774\ub2e4. (\ucc38 \uc778\uacfc\uad00\uacc4 \ubaa8\ud615\uc5d0\uc11c \\(\ud569\uaca9\ub960 = f(0 \\cdot \ub0a8\ub140 + e)\\) \ub77c\uace0 \ud560 \ub54c, \ub0a8\uc790\uc77c \ub54c \\(e\\) \ub294 \uc591\uc218\uc774\uace0, \uc5ec\uc790\uc77c \ub54c \\(e\\) \ub294 \uc74c\uc218\uc774\ub2e4.)<\/p>\n<p>\ub2e4\uc74c\uc740 R\uc758 <code>glm<\/code>\uc73c\ub85c \ubd84\uc11d\ud55c \uacb0\uacfc\uc774\ub2e4. \ud1b5\uc81c\ubcc0\uc218\ub85c \ubb34\uc5c7\uc744 \ud3ec\ud568\uc2dc\ud0a4\ub290\ub0d0\uc5d0 \ub530\ub77c <code>GenderFemale<\/code>\uc758 \uacc4\uc218\uac00 \ub2ec\ub77c\uc9d0\uc744 \ud655\uc778\ud574\ubcf4\uc790. <\/p>\n<pre><code class=\"r\">require(tidyr)\r\nUCBAdf &lt;- as.data.frame(UCBAdmissions)\r\ndf &lt;- spread(UCBAdf, key=&#39;Admit&#39;, value=&#39;Freq&#39;)\r\nsummary(glm(cbind(Admitted, Rejected) ~ Gender, family=binomial, data=df))\r\n<\/code><\/pre>\n<pre>## \r\n## Call:\r\n## glm(formula = cbind(Admitted, Rejected) ~ Gender, family = binomial, \r\n##     data = df)\r\n## \r\n## Deviance Residuals: \r\n##      Min        1Q    Median        3Q       Max  \r\n## -16.7915   -4.7613   -0.4365    5.1025   11.2022  \r\n## \r\n## Coefficients:\r\n##              Estimate Std. Error z value Pr(&gt;|z|)    \r\n## (Intercept)  -0.22013    0.03879  -5.675 1.38e-08 ***\r\n## GenderFemale -0.61035    0.06389  -9.553  &lt; 2e-16 ***\r\n## ---\r\n## Signif. codes:  0 &#39;***&#39; 0.001 &#39;**&#39; 0.01 &#39;*&#39; 0.05 &#39;.&#39; 0.1 &#39; &#39; 1\r\n## \r\n## (Dispersion parameter for binomial family taken to be 1)\r\n## \r\n##     Null deviance: 877.06  on 11  degrees of freedom\r\n## Residual deviance: 783.61  on 10  degrees of freedom\r\n## AIC: 856.55\r\n## \r\n## Number of Fisher Scoring iterations: 4\r\n<\/pre>\n<pre><code class=\"r\">summary(glm(cbind(Admitted, Rejected) ~ Gender + Dept, family=binomial, data=df))\r\n<\/code><\/pre>\n<pre>## \r\n## Call:\r\n## glm(formula = cbind(Admitted, Rejected) ~ Gender + Dept, family = binomial, \r\n##     data = df)\r\n## \r\n## Deviance Residuals: \r\n##       1        2        3        4        5        6        7        8  \r\n## -1.2487  -0.0560   1.2533   0.0826   1.2205  -0.2076   3.7189   0.2706  \r\n##       9       10       11       12  \r\n## -0.9243  -0.0858  -0.8509   0.2052  \r\n## \r\n## Coefficients:\r\n##              Estimate Std. Error z value Pr(&gt;|z|)    \r\n## (Intercept)   0.58205    0.06899   8.436   &lt;2e-16 ***\r\n## GenderFemale  0.09987    0.08085   1.235    0.217    \r\n## DeptB        -0.04340    0.10984  -0.395    0.693    \r\n## DeptC        -1.26260    0.10663 -11.841   &lt;2e-16 ***\r\n## DeptD        -1.29461    0.10582 -12.234   &lt;2e-16 ***\r\n## DeptE        -1.73931    0.12611 -13.792   &lt;2e-16 ***\r\n## DeptF        -3.30648    0.16998 -19.452   &lt;2e-16 ***\r\n## ---\r\n## Signif. codes:  0 &#39;***&#39; 0.001 &#39;**&#39; 0.01 &#39;*&#39; 0.05 &#39;.&#39; 0.1 &#39; &#39; 1\r\n## \r\n## (Dispersion parameter for binomial family taken to be 1)\r\n## \r\n##     Null deviance: 877.056  on 11  degrees of freedom\r\n## Residual deviance:  20.204  on  5  degrees of freedom\r\n## AIC: 103.14\r\n## \r\n## Number of Fisher Scoring iterations: 4\r\n<\/pre>\n<pre><code class=\"r\">#summary(glm(cbind(Admitted, Rejected) ~ Gender * Dept, family=binomial, data=df))\r\n<\/code><\/pre>\n<p>\uc774\ub54c \ub2e4\uc74c\uacfc \uac19\uc774 \uc624\ud574\ud558\uae30 \uc27d\ub2e4.<\/p>\n<p><strong>\uadf8\ub807\ub2e4\uba74 \uc778\uacfc \uad00\uacc4\ub97c \ucd94\uc815\ud558\uae30 \uc704\ud574\uc11c\ub294 \ub420 \uc218 \uc788\uc73c\uba74 \ub9ce\uc740 \ud1b5\uc81c \ubcc0\uc218\ub97c \ud3ec\ud568\ud558\uc5ec \ubaa8\ud615\uc744 \uad6c\uc131\ud558\uba74 \ub418\ub294\uac70 \uc544\ub2cc\uac00?<\/strong><\/p>\n<p>\uc74c.. \ub531\ud788 \uadf8\ub807\uc9c0 \uc54a\ub2e4. \ub2e4\uc74c\uc758 \uc608\ub97c \ubcf4\uc790.<\/p>\n<h2>\ubc84\ud06c\uc2a8\uc758 \uc5ed\uc124(Berkson&#39;s paradox)<\/h2>\n<p>\uc2ec\uc2a8\uc758 \uc5ed\uc124\uc5d0 \uc775\uc219\ud55c \uc0ac\ub78c\ub3c4 \ubc84\ud06c\uc2a8\uc758 \uc5ed\uc124\uc740 \uc0dd\uc18c\ud560 \uc218 \uc788\ub2e4. \ubc84\ud06c\uc2a8\uc758 \uc5ed\uc124\uc740 <strong>\ud1b5\uc81c \ubcc0\uc218<\/strong>\uc640 \uc778\uacfc\uad00\uacc4 \ucd94\uc815\uc5d0\uc11c \uc2ec\uc2a8 \uc5ed\uc124\uc758 \uc815\ubc18\ub300 \uc9c0\uc810\uc5d0 \ub193\uc5ec\uc788\ub2e4.<\/p>\n<p>\uc2ec\uc2a8 \uc5ed\uc124\uc740 <strong>\ud1b5\uc81c \ubcc0\uc218\ub97c \ud3ec\ud568\ud558\uc9c0 \uc54a\uc544 \uc0dd\uae30\ub294<\/strong> \uac00\uc9dc \uc778\uacfc \uad00\uacc4\ub97c \ubcf4\uc5ec\uc900\ub2e4\uba74, \ubc84\ud06c\uc2a8\uc758 \uc5ed\uc124\uc740 <strong>\ud1b5\uc81c \ubcc0\uc218\ub97c \ud3ec\ud568\ud574\uc11c \uc0dd\uae30\ub294<\/strong> \uac00\uc9dc \uc778\uacfc\uad00\uacc4\ub97c \ubcf4\uc5ec\uc900\ub2e4. <\/p>\n<p>\ub2e4\uc74c\uc758 \uc608\ub97c \ubcf4\uc790.<\/p>\n<h3>\uc2dd\ub2f9\uc758 \ud3c9\uc810 = \ub9db + \uc778\ud14c\ub9ac\uc5b4<\/h3>\n<p>\uc2dd\ub2f9\uc758 \ud3c9\uc810\uc740 \ud06c\uac8c \ub450 \uac00\uc9c0 \uc694\uc18c(\ub9db\uacfc \uc778\ud14c\ub9ac\uc5b4)\uc5d0 \ub300\ud55c \ud3c9\uc810\uc758 \ud569\uc73c\ub85c \uc774\ub8e8\uc5b4\uc9c4\ub2e4. \ub530\ub77c\uc11c \uc778\uacfc\uad00\uacc4 \ubaa8\ud615\uc740 \ub2e4\uc74c\uacfc \uac19\ub2e4.<\/p>\n<p>\\[\ud3c9\uc810 = \ub9db + \uc778\ud14c\ub9ac\uc5b4 + e\\]<\/p>\n<p>\uc774\ub97c \uadf8\ub798\ud504\ub85c \ud45c\ud604\ud558\uba74 \ub2e4\uc74c\uacfc \uac19\ub2e4.<\/p>\n<p><img src=\"http:\/\/141.164.34.82\/wp-content\/uploads\/2019\/07\/Berkson.png\" alt=\"\ud3c9\uc810 \uc778\uacfc\uad00\uacc4 \ubaa8\ud615\"\/><\/p>\n<p>\uc774\ub54c \ub9db\uacfc \uc778\ud14c\ub9ac\uc5b4\ub294 \uc11c\ub85c \uc778\uacfc \uad00\uacc4\uac00 \uc5c6\uae30 \ub54c\ubb38\uc5d0 \ub2e4\uc74c\uacfc \uac19\uc774 \uc778\uacfc\uad00\uacc4 \uacc4\uc218\ub294 \uc720\uc758\ubbf8\ud558\uc9c0 \uc54a\ub2e4. <\/p>\n<pre><code class=\"r\">N &lt;- 1000\r\ntaste &lt;- rnorm(N)\r\ninterior &lt;- rnorm(N)\r\nscore &lt;- taste + interior + (e=rnorm(N))\r\ndat &lt;- data.frame(taste, interior, score)\r\nsummary(lm(taste ~ interior, data=dat))\r\n<\/code><\/pre>\n<pre>## \r\n## Call:\r\n## lm(formula = taste ~ interior, data = dat)\r\n## \r\n## Residuals:\r\n##      Min       1Q   Median       3Q      Max \r\n## -2.58375 -0.66466  0.00948  0.66956  2.45230 \r\n## \r\n## Coefficients:\r\n##             Estimate Std. Error t value Pr(&gt;|t|)\r\n## (Intercept) -0.02091    0.03014  -0.694    0.488\r\n## interior    -0.01450    0.03093  -0.469    0.639\r\n## \r\n## Residual standard error: 0.9525 on 998 degrees of freedom\r\n## Multiple R-squared:  0.0002203,  Adjusted R-squared:  -0.0007815 \r\n## F-statistic: 0.2199 on 1 and 998 DF,  p-value: 0.6392\r\n<\/pre>\n<p>\ud558\uc9c0\ub9cc <code>score<\/code>(\ud3c9\uc810)\uc744 \ud1b5\uc81c\ubcc0\uc218\ub97c \ub123\uc5b4 \ud68c\uadc0 \ubd84\uc11d\uc744 \ud558\uba74, \ub2e4\uc74c\uacfc \uac19\uc774 <code>interior<\/code>(\uc778\ud14c\ub9ac\uc5b4)\ub294 <code>taste<\/code>(\ub9db)\uc744 \uac10\uc18c\uc2dc\ud0a4\ub294 \uac83\uc73c\ub85c \ub098\ud0c0\ub09c\ub2e4.<\/p>\n<pre><code class=\"r\">summary(lm(taste ~ interior + score, data=dat))\r\n<\/code><\/pre>\n<pre>## \r\n## Call:\r\n## lm(formula = taste ~ interior + score, data = dat)\r\n## \r\n## Residuals:\r\n##     Min      1Q  Median      3Q     Max \r\n## -2.0488 -0.4223 -0.0136  0.4515  2.2760 \r\n## \r\n## Coefficients:\r\n##             Estimate Std. Error t value Pr(&gt;|t|)    \r\n## (Intercept) -0.02000    0.02147  -0.932    0.352    \r\n## interior    -0.46642    0.02637 -17.687   &lt;2e-16 ***\r\n## score        0.47923    0.01538  31.162   &lt;2e-16 ***\r\n## ---\r\n## Signif. codes:  0 &#39;***&#39; 0.001 &#39;**&#39; 0.01 &#39;*&#39; 0.05 &#39;.&#39; 0.1 &#39; &#39; 1\r\n## \r\n## Residual standard error: 0.6783 on 997 degrees of freedom\r\n## Multiple R-squared:  0.4935, Adjusted R-squared:  0.4925 \r\n## F-statistic: 485.8 on 2 and 997 DF,  p-value: &lt; 2.2e-16\r\n<\/pre>\n<p><img src=\"http:\/\/141.164.34.82\/wp-content\/uploads\/2019\/07\/Berkson0.png\" alt=\"\ub9db \uc608\uce21 \ubaa8\ud615, \uc798\ubabb\ub41c \ub9db \uc778\uacfc\uad00\uacc4 \ubaa8\ud615\"\/><\/p>\n<p>\uc774\uc5d0 \ub300\ud55c \uac04\ub2e8\ud55c \uc124\uba85\uc740 \ub2e4\uc74c\uacfc \uac19\ub2e4. \ud3c9\uc810\uc774 \uac19\uc740 \uacbd\uc6b0, \uc778\ud14c\ub9ac\uc5b4\uac00 \ub192\ub2e4\uba74, \ub9db\uc740 \ub0ae\uc544\uc9c0\uac8c \ub9c8\ub828\uc774\ub2e4. \uc774\ub294 \uc608\uce21\ubaa8\ud615\uc73c\ub85c\ub294 \ud6cc\ub96d\ud558\ub2e4. \ud3c9\uc810\uc774 \uc778\ud14c\ub9ac\uc5b4\uc640 \ub9db \uc810\uc218\uc758 \ud569\uc774\uae30 \ub54c\ubb38\uc5d0 \uc911\uac04 \ud3c9\uc810\uc758 \uc2dd\ub2f9\uc740 (\uc911 \uc778\ud14c\ub9ac\uc5b4, \uc911 \ub9db), (\uace0 \uc778\ud14c\ub9ac\uc5b4, \uc800 \ub9db), (\uc800 \uc778\ud14c\ub9ac\uc5b4, \uace0 \ub9db) \uc911\uc758 \ud558\ub098\uc77c \uac83\uc774\ub2e4. \uc774\ub7f0 \uc0c1\ud669(\ud3c9\uc810\uc744 \ud1b5\uc81c\ud55c \uc0c1\ud0dc)\uc5d0\uc11c \uc778\ud14c\ub9ac\uc5b4\uac00 \ub0ae\uc744 \uc218\ub85d \ub9db\uc774 \uc88b\ub2e4\uace0 \uae30\ub300\ud560 \uc218 \uc788\ub2e4.<\/p>\n<p>\uc5ec\uce5c : &ldquo;\uc790\uae30\uc57c, \uc774 \uc2dd\ub2f9\uc740 \ud3c9\uc810\uc740 \uc5c4\uccad \ub192\uc740\ub370 \uc778\ud14c\ub9ac\uc5b4\ub294 \ubcc4\ub85c\ub124. \ub2e4\ub978 \uc2dd\ub2f9 \uac08\uae4c?&rdquo;<\/p>\n<p>\uc774\uae00\uc744 \uc77d\uc740 \ub0a8\uce5c : &ldquo;\ubc84\ud06c\uc2a8\uc758 \uc5ed\uc124\uc774\ub77c\uace0 \uc788\uc5b4. \uc5ec\uae34 \uc544\ub9c8 \ub9db\uc774 \ub05d\ub0b4\uc904\uaebc\uc57c!&rdquo;<\/p>\n<h2>\ub9e4\uac1c \ud6a8\uacfc \ubaa8\ud615<\/h2>\n<p>\uc774\ub7f0 \ubb38\uc81c\ub294 \ub9e4\uac1c \ud6a8\uacfc \ubaa8\ud615\uc5d0\uc11c\ub3c4 \ub9c8\ucc2c\uac00\uc9c0\uc774\ub2e4. \ub9e4\uac1c\ubcc0\uc218\ub9cc\uc744 \ub123\uc5b4 \ud68c\uadc0\ubd84\uc11d\uc744 \ud558\uba74 \ub9e4\uac1c \ubcc0\uc218\uc758 \uc778\uacfc\uad00\uacc4 \uacc4\uc218\ub97c \uacfc\ub300 \ub610\ub294 \uacfc\uc18c \ucd94\uc815\ud558\uac8c \ub418\uace0, \ub3c5\ub9bd\ubcc0\uc218\uc5d0 \ub9e4\uac1c\ubcc0\uc218\ub97c \ud1b5\uc81c\ubcc0\uc218\ub85c \ub193\uc5b4 \ud68c\uadc0\ubd84\uc11d\uc744 \ud558\uba74 \ucd1d\ud6a8\uacfc\uac00 \uc544\ub2c8\ub77c \uc9c1\uc811\ud6a8\uacfc(\ucd1d\ud6a8\uacfc-\uac04\uc811\ud6a8\uacfc)\ub9cc\uc744 \ucd94\uc815\ud558\uac8c \ub41c\ub2e4. <\/p>\n<p>\uc774\ub7f0 \uc5f0\uc720\ub85c \ub9ce\uc740 \uacbd\uc81c\ud559 \ub17c\ubb38\ub4e4\uc774 \uc5ec\ub7ec \ubcc0\uc218\ub97c \ud1b5\uc81c\ubcc0\uc218\ub85c \ub123\uc5b4\uc11c \ubd84\uc11d\ud574\ubcf4\uae30\ub3c4 \ud558\uace0, \ube7c\ub193\uace0 \ud574\ubcf4\uae30\ub3c4 \ud558\uba74\uc11c <strong>\uac00\uc815<\/strong>\uc5d0 \ub530\ub77c \uacb0\ub860\uc774 \uc5b4\ub5bb\uac8c \ub2ec\ub77c\uc9c0\ub294 \uc0b4\ud3b4\ubcf8\ub2e4. \ud558\uc9c0\ub9cc \uac1c\uc778\uc801\uc73c\ub85c \uadf8\ub807\uac8c \ud1b5\uc81c \ubcc0\uc218\ub97c \uc120\ud0dd\ud558\ub294 \uac83\ub9cc\uc73c\ub85c \uc778\uacfc \uad00\uacc4\ub97c \uc815\ud655\ud558\uac8c \ucd94\uc815\ud560 \uc218 \uc788\uc744\uc9c0\ub294 \uc758\ubb38\uc774\ub2e4. \uc77c\ub2e8 \ud1b5\uc81c \ubcc0\uc218\ub85c \ud3ec\ud568\uc2dc\ud0a4\uace0 \uc548 \ud558\uace0\ub294 \uc774\ub860\uc801 \ubc30\uacbd\uc5d0 \uc758\ud574 \ub2ec\ub77c\uc838\uc57c \ud558\ub294\ub370 \uadf8\ub7f0 \uac83\uc740 \ubcc4\ub85c \uc5c6\uace0, \ud1b5\uc81c\ub418\uc5b4\uc57c \ud558\uc9c0\ub9cc \ud1b5\uc81c\ub418\uc9c0 \ubabb\ud55c(\ud639\uc740 \uad00\ucc30\ub418\uc9c0 \ubabb\ud55c) \ubcc0\uc218\uac00 \uc788\ub2e4\uba74 \uacb0\ub860\uc740 \ub2e4\uc2dc \ubc14\ub014 \uc218 \uc788\uc73c\uba70, \ud1b5\uc81c\ubcc0\uc218\uc640 \uacb0\uacfc \ubcc0\uc218\uc758 \uad00\uacc4\uac00 \uc120\ud615\uc801\uc774\ub77c\ub294 \ubcf4\uc7a5\ub3c4 \uc5c6\uc73c\uba70&hellip;<\/p>\n<h2>\uacb0\ub860<\/h2>\n<p><strong>\ud1b5\uc81c \ubcc0\uc218\ub97c \ub9ce\uc774 \ub123\ub294\ub2e4\uace0 \uc88b\uc740 \uac83\uc740 \uc544\ub2c8\ub2e4.<\/strong><\/p>\n<p>\uc608\uce21 \ubaa8\ud615\uc774\ub77c\uba74 \uacfc\uc801\ud569\uc774 \uc77c\uc5b4\ub098\uc9c0 \uc54a\ub294 \ud55c\ub3c4\uc5d0\uc11c \uc124\uba85 \ubcc0\uc218\ub97c \ub9ce\uc774 \ub123\uc5b4\uc11c \ub098\uc060 \uac83\uc740 \uc5c6\uc9c0\ub9cc,<\/p>\n<p>\uc778\uacfc\uad00\uacc4 \ubaa8\ud615\uc774\ub77c\uba74 \ud1b5\uc81c \ubcc0\uc218\uc5d0 \ub530\ub77c \uacb0\uacfc\uac00 \ubc14\ub014 \uc218\ub3c4 \uc788\ub2e4.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>\uc2ec\uc2a8\uacfc \ubc84\ud06c\uc2a8 \uc2ec\uc2a8\uc758 \uc5ed\uc124(Simpson&#39;s paradox) R\uc758 datasets::UCBAdmissions\uc740 1973\ub144 UC Berkeley\uc758 \ub300\ud559\uc6d0 \uc785\ud559 \uc2dc\ud5d8 \uacb0\uacfc\ub97c \ubcf4\uc5ec \uc900\ub2e4. 1973\ub144 UC \ubc84\ud074\ub9ac(Berkeley) \ub300\ud559\uc740 \uc131\ucc28\ubcc4\ub85c \uace0\uc18c\ub97c \ub2f9\ud55c\ub2e4. \uc9c0\uc6d0\uc790 \uc911 \ub0a8\uc131\uc740 44%\uac00 \ud569\uaca9\ud588\uc9c0\ub9cc, \uc5ec\uc131\uc740 35%\ub9cc\uc774 \ud569\uaca9\uc744 \ud588\uae30 \ub54c\ubb38\uc774\ub2e4. \ub2e4\uc74c\uc740 \ubc84\ud074\ub9ac \ub300\ud559\uc6d0 \uc785\ud559 \uacb0\uacfc\ub97c \uadf8\ub798\ud504\ub85c \ubcf4\uc5ec\uc8fc\uace0 \uc788\ub2e4. # from help(UCBAdmissions) require(graphics) ## Data aggregated over departments apply(UCBAdmissions, c(1, 2), sum) ## [&hellip;]<\/p>\n","protected":false},"author":1,"featured_media":1884,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":[],"categories":[146,319],"tags":[354,351,353,352,337],"jetpack_featured_media_url":"http:\/\/ds.sumeun.org\/wp-content\/uploads\/2019\/07\/Simpson_Dia02.png","_links":{"self":[{"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/posts\/1886"}],"collection":[{"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=1886"}],"version-history":[{"count":2,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/posts\/1886\/revisions"}],"predecessor-version":[{"id":1888,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/posts\/1886\/revisions\/1888"}],"wp:featuredmedia":[{"embeddable":true,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=\/wp\/v2\/media\/1884"}],"wp:attachment":[{"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=1886"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=1886"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/ds.sumeun.org\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=1886"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}