{"id":1012,"date":"2023-07-28T00:10:59","date_gmt":"2023-07-28T00:10:59","guid":{"rendered":"https:\/\/statorials.org\/cn\/%e5%9b%be%e5%9f%ba%e6%b5%8b%e8%af%95\/"},"modified":"2023-07-28T00:10:59","modified_gmt":"2023-07-28T00:10:59","slug":"%e5%9b%be%e5%9f%ba%e6%b5%8b%e8%af%95","status":"publish","type":"post","link":"https:\/\/statorials.org\/cn\/%e5%9b%be%e5%9f%ba%e6%b5%8b%e8%af%95\/","title":{"rendered":"\u5982\u4f55\u5728 r \u4e2d\u6267\u884c tukey \u6d4b\u8bd5"},"content":{"rendered":"

<\/p>\n

\n

\u5355\u5411\u65b9\u5dee\u5206\u6790<\/a>\u7528\u4e8e\u786e\u5b9a\u4e09\u4e2a\u6216\u66f4\u591a\u72ec\u7acb\u7ec4\u7684\u5e73\u5747\u503c\u4e4b\u95f4\u662f\u5426\u5b58\u5728\u7edf\u8ba1\u663e\u7740\u5dee\u5f02\u3002<\/span><\/p>\n

\u5982\u679c\u65b9\u5dee\u5206\u6790\u8868\u7684\u603b\u4f53p \u503c<\/a>\u4f4e\u4e8e\u4e00\u5b9a\u7684\u663e\u7740\u6027\u6c34\u5e73\uff0c\u5219\u6211\u4eec\u6709\u8db3\u591f\u7684\u8bc1\u636e\u8868\u660e\u81f3\u5c11\u5176\u4e2d\u4e00\u4e2a\u7ec4\u5747\u503c\u4e0e\u5176\u4ed6\u5747\u503c\u4e0d\u540c\u3002<\/span><\/p>\n

\u7136\u800c\uff0c\u8fd9\u5e76\u6ca1\u6709\u544a\u8bc9\u6211\u4eec\u54ea\u4e9b<\/em>\u7fa4\u4f53\u5f7c\u6b64\u4e0d\u540c\u3002\u8fd9\u53ea\u662f\u544a\u8bc9\u6211\u4eec\u5e76\u975e\u6240\u6709\u7ec4\u7684\u5e73\u5747\u503c\u90fd\u662f\u76f8\u7b49\u7684\u3002\u4e3a\u4e86\u51c6\u786e\u5730\u4e86\u89e3\u54ea\u4e9b\u7ec4\u5f7c\u6b64\u4e0d\u540c\uff0c\u6211\u4eec\u9700\u8981\u6267\u884c\u4e8b\u540e\u6d4b\u8bd5<\/a>\u3002<\/span><\/p>\n

\u6700\u5e38\u7528\u7684\u4e8b\u540e\u68c0\u9a8c\u4e4b\u4e00\u662fTukey \u68c0\u9a8c<\/strong>\uff0c\u5b83\u5141\u8bb8\u6211\u4eec\u5728\u6bcf\u7ec4\u5747\u503c\u4e4b\u95f4\u8fdb\u884c\u6210\u5bf9\u6bd4\u8f83\uff0c\u540c\u65f6\u63a7\u5236\u65cf\u9519\u8bef\u7387<\/a>\u3002<\/span><\/p>\n

\u672c\u6559\u7a0b\u4ecb\u7ecd\u5982\u4f55\u5728 R \u4e2d\u6267\u884c Tukey \u6d4b\u8bd5\u3002<\/span><\/p>\n

\u6ce8\u610f\uff1a<\/strong>\u5982\u679c\u7814\u7a76\u4e2d\u7684\u4efb\u4f55\u7ec4\u88ab\u89c6\u4e3a\u5bf9\u7167\u7ec4\uff0c\u5219\u5e94\u4f7f\u7528Dunnett \u68c0\u9a8c<\/a>\u4f5c\u4e3a\u4e8b\u540e\u68c0\u9a8c\u3002<\/em><\/span><\/p>\n

\u793a\u4f8b\uff1aR \u4e2d\u7684 Tukey \u6d4b\u8bd5<\/strong><\/span><\/h3>\n

\u6b65\u9aa4 1\uff1a\u62df\u5408\u65b9\u5dee\u5206\u6790\u6a21\u578b\u3002<\/strong><\/span><\/p>\n

\u4ee5\u4e0b\u4ee3\u7801\u5c55\u793a\u4e86\u5982\u4f55\u521b\u5efa\u4e00\u4e2a\u5305\u542b\u4e09\u7ec4\uff08A\u3001B \u548c C\uff09\u7684\u5047\u6570\u636e\u96c6\uff0c\u5e76\u5bf9\u6570\u636e\u62df\u5408\u5355\u5411\u65b9\u5dee\u5206\u6790\u6a21\u578b\u4ee5\u786e\u5b9a\u6bcf\u7ec4\u7684\u5e73\u5747\u503c\u662f\u5426\u76f8\u7b49\uff1a<\/span><\/p>\n

 #make this example reproducible<\/span>\nset.seed(0)\n\n#create data<\/span>\ndata <- data.frame(group = rep<\/span> (c(\"A\", \"B\", \"C\"), each<\/span> = 30),\n                   values = c(runif(30, 0, 3),\n                                   runif(30, 0, 5),\n                                   runif(30, 1, 7)))\n\n#view first six rows of data<\/span>\nhead(data)\n\n  group values\n1 A 2.6900916\n2 A 0.7965260\n3 A 1.1163717\n4 A 1.7185601\n5 A 2.7246234\n6 A 0.6050458\n\n#fit one-way ANOVA model<\/span>\nmodel <- aov<\/span> (values~group, data=data)\n\n#view the model output\n<\/span>summary(model)\n\n            Df Sum Sq Mean Sq F value Pr(>F)    \ngroup 2 98.93 49.46 30.83 7.55e-11 ***\nResiduals 87 139.57 1.60                     \n---\nSignificant. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1\n<\/strong><\/pre>\n
\u6211\u4eec\u53ef\u4ee5\u770b\u5230\u65b9\u5dee\u5206\u6790\u8868\u4e2d\u7684\u603b\u4f53 p \u503c\u4e3a7.55e-11<\/strong> \u3002\u7531\u4e8e\u8fd9\u4e2a\u6570\u5b57\u5c0f\u4e8e0.05\uff0c\u6211\u4eec\u6709\u8db3\u591f\u7684\u8bc1\u636e\u8868\u660e\u6bcf\u7ec4\u4e2d\u7684\u5e73\u5747\u503c\u5e76\u4e0d\u76f8\u7b49\u3002\u56e0\u6b64\uff0c\u6211\u4eec\u53ef\u4ee5\u8fdb\u884c Tukey \u68c0\u9a8c\u6765\u51c6\u786e\u786e\u5b9a\u54ea\u4e9b\u7ec4\u5747\u503c\u4e0d\u540c\u3002<\/span><\/p>\n
\u6b65\u9aa4 2\uff1a\u6267\u884c Tukey \u6d4b\u8bd5\u3002<\/strong><\/span><\/p>\n
\u4ee5\u4e0b\u4ee3\u7801\u663e\u793a\u5982\u4f55\u4f7f\u7528TukeyHSD()<\/strong>\u51fd\u6570\u6267\u884c Tukey \u6d4b\u8bd5\uff1a<\/span><\/p>\n
 #perform Tukey's Test\n<\/span>TukeyHSD(model, conf.level= .95<\/span> ) \n\n  Tukey multiple comparisons of means\n    95% family-wise confidence level\n\nFit: aov(formula = values ~ group, data = data)\n\n$group\n         diff lwr upr p adj\nBA 0.9777414 0.1979466 1.757536 0.0100545\nCA 2.5454024 1.7656076 3.325197 0.0000000\nCB 1.5676610 0.7878662 2.347456 0.0000199\n<\/strong><\/pre>\n
 p \u503c\u8868\u793a\u6bcf\u4e2a\u7a0b\u5e8f\u4e4b\u95f4\u662f\u5426\u5b58\u5728\u7edf\u8ba1\u663e\u7740\u5dee\u5f02\u3002\u7ed3\u679c\u663e\u793a\uff0c\u5404\u65b9\u6848\u7684\u5e73\u5747\u51cf\u80a5\u6548\u679c\u57280.05\u663e\u7740\u6027\u6c34\u5e73\u4e0a\u5b58\u5728\u7edf\u8ba1\u5b66\u4e0a\u7684\u663e\u7740\u5dee\u5f02\u3002<\/span><\/p>\n
\u5c24\u5176\uff1a<\/span><\/p>\n