{"id":519,"date":"2018-08-11T17:11:00","date_gmt":"2018-08-11T08:11:00","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/dunnett-test-in-r\/"},"modified":"2024-10-13T19:46:30","modified_gmt":"2024-10-13T10:46:30","slug":"dunnett-test-in-r","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/dunnett-test-in-r\/","title":{"rendered":"R \u3067\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u3092\u884c\u3046\u65b9\u6cd5"},"content":{"rendered":"\n<p>\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u306f\u3001\u6bd4\u8f03\u5bfe\u7167\u7fa4\u3068\u3044\u304f\u3064\u304b\u306e\u5b9f\u9a13\u7fa4\u3092\u591a\u91cd\u6bd4\u8f03\u3059\u308b\u65b9\u6cd5\u3002<\/p>\n\n\n\n<p>R\u3067\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u3092\u3059\u308b\u306b\u306f\u3069\u3046\u3057\u305f\u3089\u3088\u3044\u304b\uff1f<\/p>\n\n\n\n<!--more-->\n\n\n\n<h2 class=\"wp-block-heading\" id=\"R\u3067\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u3092\u3059\u308b\u306b\u306f\">R\u3067\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u3092\u3059\u308b\u306b\u306f\uff1f<\/h2>\n\n\n\n<p>\u307e\u305amultcomp\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3002<\/p>\n\n\n\n<p>\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u306f\u4e00\u56de\u3060\u3051\u3067OK\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synIdentifier\">install.packages<\/span><span class=\"synSpecial\">(<\/span><span class=\"synConstant\">\"multcomp\"<\/span><span class=\"synSpecial\">)<\/span>\n<\/code><\/pre>\n\n\n\n<p>\u6b21\u306bmultcomp\u3092\u547c\u3073\u51fa\u3059\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synPreProc\">library<\/span><span class=\"synSpecial\">(<\/span>multcomp<span class=\"synSpecial\">)<\/span>\n<\/code><\/pre>\n\n\n\n<p>\u4eca\u56de\u3001\u4f8b\u3068\u3057\u3066warpbreaks\u3068\u3044\u3046\u30c7\u30fc\u30bf\u3092\u4f7f\u3046\u3002<\/p>\n\n\n\n<p>\u6a5f\u7e54\u308a\u306b\u304a\u3044\u3066\u3001tension\uff08\u7dca\u5f35\u3001\u5f35\u308a\u3001\u30c6\u30f3\u30b7\u30e7\u30f3\uff09\u306e\u9055\u3044\u306b\u3088\u3063\u3066\u3001warpbreaks\u304c\u3044\u304f\u3064\u8d77\u304d\u308b\u304b\u3068\u3044\u3046\u30c7\u30fc\u30bf\u3092\u5206\u6790\u3059\u308b\u3002<\/p>\n\n\n\n<p>\u3061\u306a\u307f\u306bWarp\u306f\u7e26\u7cf8\u3002\u6a2a\u7cf8\u306fWeft\u3068\u3044\u3046\u3002<\/p>\n\n\n\n<p>\u624b\u9806\u3068\u3057\u3066\u306f\u3001\u307e\u305a\u3001\u4e00\u5143\u914d\u7f6e\u5206\u6563\u5206\u6790ANOVA\u3092\u884c\u3046\u3002<\/p>\n\n\n\n<p>\u6b21\u306b\u3001ANOVA\u306e\u30aa\u30d6\u30b8\u30a7\u30af\u30c8\u3092\u4f7f\u3063\u3066\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u3092\u3059\u308b\u3002<\/p>\n\n\n\n<p>\u3053\u306e\u6d41\u308c\u306b\u306a\u308b\u3002<\/p>\n\n\n\n<p>amod\u304cANOVA\u306e\u30aa\u30d6\u30b8\u30a7\u30af\u30c8\u3060\u3002<\/p>\n\n\n\n<p>General Linear Hypotheses\u3092\u5206\u6790\u3059\u308bglht()\u3067\u3001\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u3092\u5b9f\u884c\u3059\u308b\u3002<\/p>\n\n\n\n<p>linfct=\u306e\u3042\u3068\u306bmcp()\u306e\u4e2d\u306b\u3001\u591a\u91cd\u6bd4\u8f03\u3059\u308b\u5909\u6570\u3068\u65b9\u6cd5\u3092\u66f8\u304f\u3002<\/p>\n\n\n\n<p>\u4eca\u56de\u306ftension\u3092Dunnett\u3067\u6bd4\u8f03\u3059\u308b\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synIdentifier\">boxplot<\/span><span class=\"synSpecial\">(<\/span>breaks <span class=\"synStatement\">~<\/span> tension<span class=\"synSpecial\">,<\/span> data <span class=\"synStatement\">=<\/span> warpbreaks<span class=\"synSpecial\">)<\/span>\namod <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">aov<\/span><span class=\"synSpecial\">(<\/span>breaks <span class=\"synStatement\">~<\/span> tension<span class=\"synSpecial\">,<\/span> data <span class=\"synStatement\">=<\/span> warpbreaks<span class=\"synSpecial\">)<\/span>\nglht.res <span class=\"synStatement\">&lt;-<\/span> <span class=\"synIdentifier\">glht<\/span><span class=\"synSpecial\">(<\/span>amod<span class=\"synSpecial\">,<\/span> linfct <span class=\"synStatement\">=<\/span> <span class=\"synIdentifier\">mcp<\/span><span class=\"synSpecial\">(<\/span>tension <span class=\"synStatement\">=<\/span> <span class=\"synConstant\">\"Dunnett\"<\/span><span class=\"synSpecial\">))<\/span>\n<span class=\"synIdentifier\">summary<\/span><span class=\"synSpecial\">(<\/span>glht.res<span class=\"synSpecial\">)<\/span>\n<\/code><\/pre>\n\n\n\n<p>tension\uff08\u7dca\u5f35\u3001\u5f35\u308a\uff09\u306fL\uff08Low\uff09\u3001M\uff08Middle\uff09\u3001H\uff08High\uff09\u306e\u4e09\u7a2e\u985e\u3002<\/p>\n\n\n\n<p>\u7bb1\u3072\u3052\u56f3\u306b\u3057\u3066\u307f\u308b\u3068\u3001L\u304c\u591a\u304f\u3066\u3001M\u3001H\u3078tension\u304c\u9ad8\u304f\u306a\u308b\u3068\u6570\u304c\u6e1b\u308b\u50be\u5411\u304c\u3042\u308b\u3002<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img decoding=\"async\" width=\"672\" height=\"672\" src=\"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2018\/08\/boxplot_example_warpbreak-1.png\" alt=\"\" class=\"wp-image-2831\" srcset=\"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2018\/08\/boxplot_example_warpbreak-1.png 672w, https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2018\/08\/boxplot_example_warpbreak-1-300x300.png 300w, https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2018\/08\/boxplot_example_warpbreak-1-150x150.png 150w\" sizes=\"(max-width: 672px) 100vw, 672px\" \/><\/figure>\n\n\n\n\n\n\n\n<p>\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u306e\u7d50\u679c\u306f\u3001M\u3082H\u3082L\u306b\u6bd4\u3079\u3066\u7d71\u8a08\u5b66\u7684\u6709\u610f\u3067\u3001breaks\u306e\u6570\u306f\u3069\u3061\u3089\u3082L\u306b\u6bd4\u3079\u3066\u5c11\u306a\u3044\u3002<\/p>\n\n\n\n<p>M\u306fL\u306b\u6bd4\u3079\u3066\u5e73\u574710\u500b\u5c11\u306a\u3044\u3002<\/p>\n\n\n\n<p>H\u306fL\u306b\u6bd4\u3079\u3066\u5e73\u5747\u7d0415\u500b\u5c11\u306a\u3044\u3002<\/p>\n\n\n\n<p>warpbreaks\u3092\u4e88\u9632\u3059\u308b\u306a\u3089\u3001L\u306f\u3084\u3081\u3066\u304a\u3044\u305f\u307b\u3046\u304c\u3088\u3044\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synStatement\">><\/span> <span class=\"synIdentifier\">summary<\/span><span class=\"synSpecial\">(<\/span>glht.res<span class=\"synSpecial\">)<\/span>\nSimultaneous Tests <span class=\"synStatement\">for<\/span> General Linear Hypotheses\nMultiple Comparisons of Means<span class=\"synSpecial\">:<\/span> Dunnett Contrasts\nFit<span class=\"synSpecial\">:<\/span> <span class=\"synIdentifier\">aov<\/span><span class=\"synSpecial\">(<\/span>formula <span class=\"synStatement\">=<\/span> breaks <span class=\"synStatement\">~<\/span> tension<span class=\"synSpecial\">,<\/span> data <span class=\"synStatement\">=<\/span> warpbreaks<span class=\"synSpecial\">)<\/span>\nLinear Hypotheses<span class=\"synSpecial\">:<\/span>\n           Estimate Std. Error t value <span class=\"synIdentifier\">Pr<\/span><span class=\"synSpecial\">(<\/span><span class=\"synStatement\">>|<\/span>t<span class=\"synStatement\">|<\/span><span class=\"synSpecial\">)<\/span>\nM <span class=\"synStatement\">-<\/span> L <span class=\"synStatement\">==<\/span> <span class=\"synConstant\">0<\/span>   <span class=\"synStatement\">-<\/span><span class=\"synConstant\">10.00<\/span>       <span class=\"synConstant\">3.96<\/span>  <span class=\"synStatement\">-<\/span><span class=\"synConstant\">2.525<\/span> <span class=\"synConstant\">0.027536<\/span> <span class=\"synStatement\">*<\/span>\nH <span class=\"synStatement\">-<\/span> L <span class=\"synStatement\">==<\/span> <span class=\"synConstant\">0<\/span>   <span class=\"synStatement\">-<\/span><span class=\"synConstant\">14.72<\/span>       <span class=\"synConstant\">3.96<\/span>  <span class=\"synStatement\">-<\/span><span class=\"synConstant\">3.718<\/span> <span class=\"synConstant\">0.000978<\/span> <span class=\"synError\">***<\/span>\n<span class=\"synStatement\">---<\/span>\nSignif. codes<span class=\"synSpecial\">:<\/span>  <span class=\"synConstant\">0<\/span> \u2018<span class=\"synError\">***<\/span>\u2019 <span class=\"synConstant\">0.001<\/span> \u2018<span class=\"synStatement\">**<\/span>\u2019 <span class=\"synConstant\">0.01<\/span> \u2018<span class=\"synStatement\">*<\/span>\u2019 <span class=\"synConstant\">0.05<\/span> \u2018.\u2019 <span class=\"synConstant\">0.1<\/span> \u2018 \u2019 <span class=\"synConstant\">1<\/span>\n<span class=\"synSpecial\">(<\/span>Adjusted p values reported <span class=\"synStatement\">--<\/span> single<span class=\"synStatement\">-<\/span>step method<span class=\"synSpecial\">)<\/span>\n<\/code><\/pre>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"R\u3067\u30c0\u30cd\u30c3\u30c8\u6bd4\u8f03\u306e\u5dee\u306e\u4fe1\u983c\u533a\u9593\u3092\u6c42\u3081\u308b\u306b\u306f\">R\u3067\u30c0\u30cd\u30c3\u30c8\u6bd4\u8f03\u306e\u5dee\u306e\u4fe1\u983c\u533a\u9593\u3092\u6c42\u3081\u308b\u306b\u306f\uff1f<\/h2>\n\n\n\n<p>confint()\u3067\u6c42\u3081\u3089\u308c\u308b\u3002<\/p>\n\n\n\n<p>M\u3068L\u306e\u5dee\u306f\u3001\u5e73\u5747\u304c\uff7010\u3067\u300195%\u4fe1\u983c\u533a\u9593\u4e0b\u9650\u304c\uff7019.0094\u3001\u4e0a\u9650\u304c-0.9906\u3002<\/p>\n\n\n\n<p>H\u3068L\u306e\u5dee\u306f\u3001\u5e73\u5747\u304c\uff7014.7222\u3067\u3001\u4fe1\u983c\u533a\u9593\u4e0b\u9650\u304c\uff7023.7316\u3001\u4e0a\u9650\u304c\uff705.712\u3002<\/p>\n\n\n\n<p>\u3068\u3082\u306b\u4fe1\u983c\u533a\u9593\u306e\u4e0a\u9650\u304c0\u3092\u4e0b\u56de\u3063\u3066\u3044\u3066\u3001\u7d71\u8a08\u5b66\u7684\u6709\u610f\u306b\u7570\u306a\u308b\u3053\u3068\u3068\u4e00\u81f4\u3057\u3066\u3044\u308b\u3002<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code><span class=\"synStatement\">&gt;<\/span> <span class=\"synIdentifier\">confint<\/span><span class=\"synSpecial\">(<\/span>glht.res<span class=\"synSpecial\">)<\/span>\nSimultaneous Confidence Intervals\nMultiple Comparisons of Means<span class=\"synSpecial\">:<\/span> Dunnett Contrasts\nFit<span class=\"synSpecial\">:<\/span> <span class=\"synIdentifier\">aov<\/span><span class=\"synSpecial\">(<\/span>formula <span class=\"synStatement\">=<\/span> breaks <span class=\"synStatement\">~<\/span> tension<span class=\"synSpecial\">,<\/span> data <span class=\"synStatement\">=<\/span> warpbreaks<span class=\"synSpecial\">)<\/span>\nQuantile <span class=\"synStatement\">=<\/span> <span class=\"synConstant\">2.275<\/span>\n<span class=\"synConstant\">95<\/span>% family<span class=\"synStatement\">-<\/span>wise confidence level\nLinear Hypotheses<span class=\"synSpecial\">:<\/span>\nEstimate lwr      upr\nM <span class=\"synStatement\">-<\/span> L <span class=\"synStatement\">==<\/span> <span class=\"synConstant\">0<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">10.0000<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">19.0094<\/span>  <span class=\"synStatement\">-<\/span><span class=\"synConstant\">0.9906<\/span>\nH <span class=\"synStatement\">-<\/span> L <span class=\"synStatement\">==<\/span> <span class=\"synConstant\">0<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">14.7222<\/span> <span class=\"synStatement\">-<\/span><span class=\"synConstant\">23.7316<\/span>  <span class=\"synStatement\">-<\/span><span class=\"synConstant\">5.7129<\/span>\n<\/code><\/pre>\n\n\n\n<div id=\"biost-398251163\" 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=\"\u307e\u3068\u3081\">\u307e\u3068\u3081<\/h2>\n\n\n\n<p>R\u3067\u306f\u3001multcomp \u30d1\u30c3\u30b1\u30fc\u30b8\u306e glht() \u3092\u4f7f\u3048\u3070\u3001\u7c21\u5358\u306b\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u304c\u3067\u304d\u308b\u3002<\/p>\n\n\n\n<p>confint()\u3067\u7fa4\u9593\u5dee\u306e\u4fe1\u983c\u533a\u9593\u3082\u6c42\u3081\u3089\u308c\u308b\u3002<\/p>\n\n\n\n<h2 class=\"wp-block-heading\" id=\"EZR\u3067\u591a\u91cd\u6bd4\u8f03\u7121\u6599\u7d71\u8a08\u30bd\u30d5\u30c8EZR\u3067\u7c21\u5358\u7d71\u8a08\u52d5\u753b\">EZR\u3067\u591a\u91cd\u6bd4\u8f03\u3010\u7121\u6599\u7d71\u8a08\u30bd\u30d5\u30c8EZR\u3067\u7c21\u5358\u7d71\u8a08\u3011\u3010\u52d5\u753b\u3011<\/h2>\n\n\n\n<p><iframe width=\"560\" height=\"315\" src=\"https:\/\/www.youtube.com\/embed\/ixy3hM1sMsY?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen=\"allowfullscreen\" title=\"EZR\u3067\u591a\u91cd\u6bd4\u8f03\u3010\u7121\u6599\u7d71\u8a08\u30bd\u30d5\u30c8EZR\u3067\u7c21\u5358\u7d71\u8a08\u3011\"><\/iframe><cite class=\"hatena-citation\"><a href=\"https:\/\/youtu.be\/ixy3hM1sMsY\">youtu.be<\/a><\/cite><\/p>\n","protected":false},"excerpt":{"rendered":"<p>\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u306f\u3001\u6bd4\u8f03\u5bfe\u7167\u7fa4\u3068\u3044\u304f\u3064\u304b\u306e\u5b9f\u9a13\u7fa4\u3092\u591a\u91cd\u6bd4\u8f03\u3059\u308b\u65b9\u6cd5\u3002 R\u3067\u30c0\u30cd\u30c3\u30c8\u691c\u5b9a\u3092\u3059\u308b\u306b\u306f\u3069\u3046\u3057\u305f\u3089\u3088\u3044\u304b\uff1f<\/p>\n","protected":false},"author":2,"featured_media":2831,"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":[4,5,49],"tags":[],"class_list":["post-519","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-ezr","category-r","category-49"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2018\/08\/boxplot_example_warpbreak-1.png","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/519","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=519"}],"version-history":[{"count":2,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/519\/revisions"}],"predecessor-version":[{"id":2833,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/519\/revisions\/2833"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/media\/2831"}],"wp:attachment":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/media?parent=519"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/categories?post=519"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/tags?post=519"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}