{"id":294,"date":"2022-06-24T20:40:08","date_gmt":"2022-06-24T11:40:08","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/line-graph-of-mean-and-error-bars-of-repeated-measure-with-ggplot2\/"},"modified":"2025-01-14T22:55:12","modified_gmt":"2025-01-14T13:55:12","slug":"line-graph-of-mean-and-error-bars-of-repeated-measure-with-ggplot2","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/line-graph-of-mean-and-error-bars-of-repeated-measure-with-ggplot2\/","title":{"rendered":"EZR \u3068 R \u306e ggplot \u3067\u30a8\u30e9\u30fc\u30d0\u30fc\u4ed8\u304d\u306e\u5e73\u5747\u5024\u6298\u308c\u7dda\u30b0\u30e9\u30d5\u3092\u66f8\u304f\u65b9\u6cd5"},"content":{"rendered":"\n

ggplot2 \u3068\u3044\u3046\u30b0\u30e9\u30d5\u3092\u304d\u308c\u3044\u306b\u63cf\u753b\u3059\u308b\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u4f7f\u3063\u3066\u3001\u30a8\u30e9\u30fc\u30d0\u30fc\u4ed8\u304d\u306e\u5e73\u5747\u5024\u306e\u6298\u308c\u7dda\u30b0\u30e9\u30d5\u3092\u63cf\u304f\u65b9\u6cd5\u3092\u7d39\u4ecb\u3059\u308b\u3002<\/p>\n\n\n\n

\u53cd\u5fa9\u6e2c\u5b9a\u3057\u305f\u30a8\u30f3\u30c9\u30dd\u30a4\u30f3\u30c8\u306e\u5e73\u5747\u5024\u3068\u3001SD\u3068\u304b95\uff05\u4fe1\u983c\u533a\u9593\u3068\u304b\u3092\u63cf\u753b\u3059\u308b\u65b9\u6cd5\u3002<\/p>\n\n\n\n\n\n\n\n

ggplot \u3067\u30a8\u30e9\u30fc\u30d0\u30fc\u4ed8\u304d\u6298\u308c\u7dda\u30b0\u30e9\u30d5\u3092\u66f8\u304f\u571f\u53f0\u3092\u4f5c\u308b<\/h2>\n\n\n\n

\u307e\u305aEZR\u306e\u30e1\u30cb\u30e5\u30fc\u304b\u3089\u3001\u53cd\u5fa9\u6e2c\u5b9a\u5206\u6563\u5206\u6790\u3092\u5b9f\u884c\u3059\u308b\u3002<\/p>\n\n\n\n

\u4ee5\u4e0b\u306e\u30c7\u30fc\u30bf\u3092\u4fee\u6b63\u3057\u3066\u3001\u30b0\u30e9\u30d5\u3092\u63cf\u304f\u30b9\u30af\u30ea\u30d7\u30c8\u306e\u524d\u307e\u3067\u3067\u4f5c\u6210\u3055\u308c\u308b\u3001TempDF3\u3068\u3044\u3046\u30c7\u30fc\u30bf\u30bb\u30c3\u30c8\u3092\u6d41\u7528\u3059\u308b\u3002<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

<\/span><\/p>\n\n\n\n

TempDF3 \u306e\u5185\u5bb9\u3092\u898b\u3066\u307f\u308b\u3068\u3001\u30ef\u30a4\u30c9\u5f62\u5f0f\u3060\u3063\u305f\u30c7\u30fc\u30bf\u3092\u30ed\u30f3\u30b0\u5f62\u5f0f\u306b\u5909\u3048\u305f\u30c7\u30fc\u30bf\u3067\u3042\u308b\u3068\u8a00\u3048\u308b\u3002<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

<\/span><\/p>\n\n\n\n

\u6700\u521d\u304b\u3089\u3001\u30ed\u30f3\u30b0\u5f62\u5f0f\u3001\u3064\u307e\u308a\u30011 \u75c7\u4f8b\u8907\u6570\u884c\u306e\u30c7\u30fc\u30bf\u306b\u306a\u3063\u3066\u3044\u308b\u5834\u5408\u306f\u3001\u6b21\u306e\u30b9\u30c6\u30c3\u30d7\u304b\u3089\u59cb\u3081\u3066\u3082\u5927\u4e08\u592b\u3068\u601d\u3046\u3002<\/p>\n\n\n\n

\u30b0\u30eb\u30fc\u30d7\u5909\u6570\uff08\u3053\u306e\u4f8b\u306f\u3001FCZ\uff09\u3068\u6642\u70b9\u5909\u6570\uff08\u3053\u306e\u4f8b\u3067\u306f\u3001time2\uff09\u3092\u3054\u81ea\u5206\u306e\u30c7\u30fc\u30bf\u306e\u5909\u6570\u540d\u306b\u7f6e\u304d\u63db\u3048\u308b\u3053\u3068\u306b\u306a\u308b\u3002<\/p>\n\n\n\n

ggplot \u3067\u30a8\u30e9\u30fc\u30d0\u30fc\u4ed8\u304d\u6298\u308c\u7dda\u30b0\u30e9\u30d5\u3092\u66f8\u304f\u524d\u6e96\u5099\u3092\u3059\u308b<\/h2>\n\n\n\n

\u30b0\u30e9\u30d5\u63cf\u753b\u306e\u524d\u306b\u3001\u30c7\u30fc\u30bf\u51e6\u7406\u3092\u4fbf\u5229\u306b\u3057\u3066\u304f\u308c\u308b dplyr \u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb install.packages(“dplyr”) \u3057\u3066\u3001library() \u3067\u547c\u3073\u51fa\u3057\u3066\u304a\u304f\u3002<\/p>\n\n\n\n

library<\/span>(<\/span>dplyr)<\/span>\n<\/code><\/pre>\n\n\n\n
\u6b21\u306b\u3001\u30b0\u30eb\u30fc\u30d7 FCZ \u306e\u30ab\u30c6\u30b4\u30ea\u3054\u3068\u3001\u6642\u70b9 time2 \u3054\u3068\u306b\u5e73\u5747\u3001\u6a19\u6e96\u504f\u5dee\u3001\uff4e\u3001\u6a19\u6e96\u8aa4\u5dee\u3092\u96c6\u8a08\u3059\u308b\u3002<\/p>\n\n\n\n
group_time_mean_sd <- TempDF3 %>%\n  group_by(FCZ, time2) %>%\n  summarize(mean=mean(data, na.rm=T), sd=sd(data, na.rm=T), n=n(), se=sd\/sqrt(n))\n<\/code><\/pre>\n\n\n\n
\u3061\u306a\u307f\u306b\u3001%>%\u3068\u304b\u3001group_by()\u3068\u304b\u304c\u3001dplyr \u30d1\u30c3\u30b1\u30fc\u30b8\u306e\u6a5f\u80fd\u30fb\u95a2\u6570\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\u96c6\u8a08\u7d50\u679c\u306f\u4ee5\u4e0b\u306e\u901a\u308a\u306b\u306a\u308b\u3002<\/p>\n\n\n\n
FCZ \u30b0\u30eb\u30fc\u30d7\u3054\u3068\u3001time2 \u6642\u70b9\u3054\u3068\u306b\u3001\u5e73\u5747\u3001\u6a19\u6e96\u504f\u5dee\u3001\u30b5\u30f3\u30d7\u30eb\u30b5\u30a4\u30ba\u3001\u6a19\u6e96\u8aa4\u5dee\u304c\u96c6\u8a08\u30fb\u8a08\u7b97\u3055\u308c\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\uff1e\uff1e\u3082\u3046\u7d71\u8a08\u3067\u60a9\u3080\u306e\u306f\u7d42\u308f\u308a\u306b\u3057\u307e\u305b\u3093\u304b\uff1f\u00a0<\/a><\/strong><\/span><\/p>\r\n\"\"<\/a>\r\n
\u21911\u4e07\u4eba\u4ee5\u4e0a\u306e\u533b\u7642\u5f93\u4e8b\u8005\u304c\u8cfc\u8aad\u4e2d<\/span><\/strong><\/span><\/p><\/div>
ggplot2 \u3092\u547c\u3073\u51fa\u3057\u3066\u30a8\u30e9\u30fc\u30d0\u30fc\u4ed8\u304d\u6298\u308c\u7dda\u30b0\u30e9\u30d5\u3092\u63cf\u753b\u3059\u308b<\/h2>\n\n\n\n
\u6b21\u306b\u3001ggplot2 \u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u547c\u3073\u51fa\u3059\u3002<\/p>\n\n\n\n
library<\/span>(<\/span>ggplot2)<\/span>\n<\/code><\/pre>\n\n\n\n
ggplot(), geom_line(), geom_errorbar(), geom_point() \u3067\u30b0\u30e9\u30d5\u3092\u63cf\u304f\u3002<\/p>\n\n\n\n
ggplot2\u306f\u3001\uff0b\uff08\u30d7\u30e9\u30b9\uff09\u3067\u95a2\u6570\u3092\u3064\u306a\u3052\u3066\u3044\u304f\u8a18\u6cd5\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\u307e\u305a\u306f\u3001\u6a19\u6e96\u504f\u5dee\u3067\u30a8\u30e9\u30fc\u30d0\u30fc\u3092\u63cf\u304f\u5834\u5408\u3002<\/p>\n\n\n\n
ggplot(data=group_time_mean_sd, aes(x = time2, y = mean, group = FCZ, color = FCZ)) +\n  geom_line() +\n  labs(x = \"Weeks\", y = \"CD ratio\") +\n  geom_errorbar(aes(ymax = mean + sd, ymin = mean - sd), width = 0.2) +\n  geom_point(aes(color = FCZ, shape = FCZ), size = 4)\n<\/code><\/pre>\n\n\n\n
\u5b9f\u884c\u3059\u308b\u3068\u4ee5\u4e0b\u306e\u3088\u3046\u306a\u30b0\u30e9\u30d5\u304c\u63cf\u3051\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u30a8\u30e9\u30fc\u30d0\u30fc\u309295\uff05\u4fe1\u983c\u533a\u9593\u306b\u3059\u308b\u5834\u5408\u306f\u4ee5\u4e0b\u306e\u30b9\u30af\u30ea\u30d7\u30c8\u3002<\/p>\n\n\n\n
mean + sd \u306e\u3068\u3053\u308d\u304c\u3001mean + 1.96*se \u3068\u5909\u308f\u3063\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
ggplot(data=group_time_mean_sd, aes(x = time2, y = mean, group = FCZ, color = FCZ)) +\n  geom_line() +\n  labs(x = \"Weeks\", y = \"CD ratio\") +\n  geom_errorbar(aes(ymax = mean + 1.96*se, ymin = mean - 1.96*se), width = 0.2) +\n  geom_point(aes(color = FCZ, shape = FCZ), size = 4)\n<\/code><\/pre>\n\n\n\n
\u30b0\u30e9\u30d5\u306f\u3053\u3093\u306a\u611f\u3058\u306b\u306a\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u70b9\u3068\u30a8\u30e9\u30fc\u30d0\u30fc\u3092\u5c11\u3057\u3060\u3051\u305a\u3089\u3057\u305f\u3044\u5834\u5408\u306f\u3001position = position_dodge() \u3092\u4f7f\u3046\u3002<\/p>\n\n\n\n
ggplot(data=group_time_mean_sd, aes(x = time2, y = mean, group = FCZ, color = FCZ)) +\n  geom_line(position = position_dodge(0.2)) +\n  labs(x = \"Weeks\", y = \"CD ratio\") +\n  geom_errorbar(aes(ymax = mean + 1.96*se, ymin = mean - 1.96*se), width = 0.2, position = position_dodge(0.2)) +\n  geom_point(aes(color = FCZ, shape = FCZ), size = 4, position = position_dodge(0.2))\n<\/code><\/pre>\n\n\n\n
\u3053\u3093\u306a\u611f\u3058\u3067\u5c11\u3057\u3060\u3051\u305a\u308c\u3066\u3001\u898b\u3084\u3059\u304f\u306a\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u5b66\u8853\u96d1\u8a8c\u7b49\u3078\u6295\u7a3f\u3059\u308b\u5834\u5408\u306b\u3001\u767d\u9ed2\u3067\u8868\u73fe\u3057\u305f\u3051\u308c\u3070\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u66f8\u304f\u3002<\/p>\n\n\n\n
ggplot()\u306e\u4e2d\u306e\u3001aes()\u5185\u306elinetype=\u304c\u7dda\u3092\u5b9f\u7dda\u3001\u7834\u7dda\u3001\u70b9\u7dda\u306a\u3069\u3092\u6307\u5b9a\u3057\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
theme_classic()\u304c\u767d\u9ed2\u3067\u30b7\u30f3\u30d7\u30eb\u306a\u898b\u305f\u76ee\u306b\u3059\u308b\u95a2\u6570\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
#for academic presentation and article\nggplot(data=group_time_mean_sd, aes(x = time2, y = mean, group = FCZ, linetype = FCZ)) +\n  geom_line() +\n  labs(x = \"Weeks\", y = \"CD ratio\") +\n  geom_errorbar(aes(ymax = mean + sd, ymin = mean - sd), width = 0.2) +\n  geom_point(aes(shape = FCZ), size = 4) +\n  theme_classic()\n<\/code><\/pre>\n\n\n\n
\u30b0\u30e9\u30d5\u306f\u3053\u3093\u306a\u611f\u3058\u306b\u306a\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
95\uff05\u4fe1\u983c\u533a\u9593\u3067\u63cf\u753b\u3059\u308b\u3068\u3053\u3093\u306a\u611f\u3058\u3002<\/p>\n\n\n\n
ggplot(data=group_time_mean_sd, aes(x = time2, y = mean, group = FCZ, linetype = FCZ)) +\n  geom_line() +\n  labs(x = \"Weeks\", y = \"CD ratio\") +\n  geom_errorbar(aes(ymax = mean + 1.96*se, ymin = mean - 1.96*se), width = 0.2) +\n  geom_point(aes(shape = FCZ), size = 4) +\n  theme_classic()\n<\/code><\/pre>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u307e\u3068\u3081<\/h2>\n\n\n\n
\u4ee5\u4e0a\u3001\u7d71\u8a08\u30bd\u30d5\u30c8\uff32\u306e ggplot2 \u30d1\u30c3\u30b1\u30fc\u30b8\u3067\u3001\u304d\u308c\u3044\u306a\u30a8\u30e9\u30fc\u30d0\u30fc\u4ed8\u304d\u5e73\u5747\u5024\u306e\u6298\u308c\u7dda\u30b0\u30e9\u30d5\u3092\u63cf\u304f\u65b9\u6cd5\u3092\u7d39\u4ecb\u3057\u305f\u3002<\/p>\n\n\n\n
EZR\u306e\u30e1\u30cb\u30e5\u30fc\u3092\u6d41\u7528\u3057\u3066\u3001\u9014\u4e2d\u304b\u3089\u30b9\u30af\u30ea\u30d7\u30c8\u3092\u66f8\u304f\u3068\u3044\u3046\u65b9\u6cd5\u3067\u3001\u7c21\u5358\u306b\u304d\u308c\u3044\u306a\u30b0\u30e9\u30d5\u304c\u63cf\u3051\u308b\u3002<\/p>\n\n\n\n
\u53c2\u8003\u306b\u306a\u308c\u3070\u3002<\/p>\n\n\n\n
\u53c2\u8003\u30b5\u30a4\u30c8<\/h2>\n\n\n\n
ggplot2\u3092\u4f7f\u3063\u305f\u30b0\u30e9\u30d5\u4f5c\u6210\uff08\u6298\u308c\u7dda\u3001\u6642\u7cfb\u5217\u3001\u30a8\u30e9\u30fc\u30d0\u30fc\uff09 : \u5e72\u304b\u3089\u3073\u305f\u30a6\u30a7\u30c3\u30c8\u6559\u6388\u306e\u72ec\u7fd2\u3067\u30a2\u30fc\u30eb R<\/a><\/p>\n\n\n\n
\u3010R\u3011ggplot2\u3067\u30a8\u30e9\u30fc\u30d0\u30fc\u3092\u4ed8\u3051\u308b\uff5c\u30c9\u30af\u30bf\u30fc\u30d5\u30f3\u30c8(Dr. Hund)<\/a><\/p>\n\n\n\n
15.18 Summarizing Data with Standard Errors and Confidence Intervals | R Graphics Cookbook, 2nd edition<\/a><\/p>\n\n\n\n
\u53c2\u8003\u66f8\u7c4d<\/h2>\n\n\n\n
EZR\u3067\u3084\u3055\u3057\u304f\u5b66\u3076\u7d71\u8a08\u5b66 \u6539\u8a023\u7248 \u301cEBM\u306e\u5b9f\u8df5\u304b\u3089\u81e8\u5e8a\u7814\u7a76\u307e\u3067\u301c<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"
ggplot2 \u3068\u3044\u3046\u30b0\u30e9\u30d5\u3092\u304d\u308c\u3044\u306b\u63cf\u753b\u3059\u308b\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u4f7f\u3063\u3066\u3001\u30a8\u30e9\u30fc\u30d0\u30fc\u4ed8\u304d\u306e\u5e73\u5747\u5024\u306e\u6298\u308c\u7dda\u30b0\u30e9\u30d5\u3092\u63cf\u304f\u65b9\u6cd5\u3092\u7d39\u4ecb\u3059\u308b\u3002 \u53cd\u5fa9\u6e2c\u5b9a\u3057\u305f\u30a8\u30f3\u30c9\u30dd\u30a4\u30f3\u30c8\u306e\u5e73\u5747\u5024\u3068\u3001SD\u3068\u304b95\uff05\u4fe1\u983c\u533a\u9593\u3068\u304b\u3092\u63cf\u753b\u3059\u308b\u65b9\u6cd5\u3002<\/p>\n","protected":false},"author":2,"featured_media":1861,"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,43,5,37,32],"tags":[],"class_list":["post-294","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-ezr","category-ggplot2","category-r","category-37","category-32"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2022\/06\/20230205131645.png","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/294","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=294"}],"version-history":[{"count":2,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/294\/revisions"}],"predecessor-version":[{"id":3421,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/294\/revisions\/3421"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/media\/1861"}],"wp:attachment":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/media?parent=294"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/categories?post=294"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/tags?post=294"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}