{"id":310,"date":"2022-04-17T00:00:00","date_gmt":"2022-04-16T15:00:00","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/how-to-draw-regression-line-with-conf-int-using-ggplot2-in-ezr\/"},"modified":"2024-09-29T22:05:59","modified_gmt":"2024-09-29T13:05:59","slug":"how-to-draw-regression-line-with-conf-int-using-ggplot2-in-ezr","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/how-to-draw-regression-line-with-conf-int-using-ggplot2-in-ezr\/","title":{"rendered":"R \u306e ggplot \u3067\u66f8\u3044\u305f\u56de\u5e30\u76f4\u7dda\u306b geom_smooth \u3067 95 % \u4fe1\u983c\u533a\u9593\u3092\u4ed8\u3051\u308b\u65b9\u6cd5"},"content":{"rendered":"\n

EZR\u3067 ggplot2 \u3092\u4f7f\u3063\u306695%\u4fe1\u983c\u533a\u9593\u3068\u56de\u5e30\u5f0f\u4ed8\u304d\u56de\u5e30\u76f4\u7dda\u3092\u63cf\u304f\u65b9\u6cd5\u3002<\/p>\n\n\n\n

EZR\u3067\u56de\u5e30\u76f4\u7dda\u306b95%\u4fe1\u983c\u533a\u9593\u3092\u63cf\u304d\u305f\u3044\u5834\u5408\u3069\u3046\u3059\u308b\u304b\uff1f<\/p>\n\n\n\n

ggplot2 \u3092\u4f7f\u3046\u3068\u7c21\u5358\u30fb\u304d\u308c\u3044\u306b\u63cf\u3051\u308b\u3002<\/p>\n\n\n\n\n\n\n\n

ggplot \u3067\u56de\u5e30\u76f4\u7dda\u3092\u66f8\u304f\u65b9\u6cd5 ggplot2 \u3068\u306f\uff1f<\/h2>\n\n\n\n

ggplot2 \u3068\u306f\u3001\u7c21\u5358\u306b\u304d\u308c\u3044\u306a\u30b0\u30e9\u30d5\u304c\u63cf\u3051\u308bR\u306e\u9ad8\u6027\u80fd\u306a\u30d1\u30c3\u30b1\u30fc\u30b8\u3067\u3042\u308b\u3002<\/p>\n\n\n\n

ggplot2\u306fEZR\u306b\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3055\u308c\u3066\u3044\u308b\u306e\u3067\u3001library()\u3067\u3059\u3050\u306b\u547c\u3073\u51fa\u305b\u308b\u3002<\/p>\n\n\n\n

\u3082\u3057\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3055\u308c\u3066\u3044\u306a\u304b\u3063\u305f\u3089\u3001install.packages(\"ggplot2\")<\/code>\u3067\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3059\u308b\u3002<\/p>\n\n\n\n

ggplot \u3067\u56de\u5e30\u76f4\u7dda\u3092\u63cf\u3044\u3066\u307f\u308b<\/h2>\n\n\n\n

\u4ee5\u4e0b\u306e\u30b9\u30af\u30ea\u30d7\u30c8\u3092R\u30b9\u30af\u30ea\u30d7\u30c8\u7a93\u306b\u66f8\u304f\u3002<\/p>\n\n\n\n

\uff13\u3064\u306e\u30d1\u30fc\u30c8\u304b\u3089\u69cb\u6210\u3055\u308c\u3066\u3044\u308b\u3002<\/p>\n\n\n\n

    \n
  1. ggplot()<\/li>\n\n\n\n
  2. geom_point()<\/li>\n\n\n\n
  3. geom_smooth()<\/li>\n<\/ol>\n\n\n\n

    \u3053\u308c\u3089\u3092+\uff08\u30d7\u30e9\u30b9\uff09\u3067\u3064\u306a\u3044\u3067\u3044\u304f\u3001\u3061\u3087\u3063\u3068\u5909\u308f\u3063\u305f\u8a18\u6cd5\u3067\u3042\u308b\u3002<\/p>\n\n\n\n

    library<\/span>(<\/span>ggplot2)<\/span>\nggplot<\/span>(<\/span>data=<\/span>GCSF_CD34,<\/span> aes<\/span>(<\/span>x=<\/span>WBC.entry,<\/span> y=<\/span>LogCD34perBW))<\/span>+<\/span>\ngeom_point<\/span>()<\/span>+<\/span>\ngeom_smooth<\/span>(<\/span>method=<\/span>lm,<\/span> se=<\/span>TRUE<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
    ggplot()\u306e\u30ab\u30c3\u30b3\u5185\u306f\u3001\u30c7\u30fc\u30bf\u306e\u6307\u5b9a\u3068X\u8ef8\u3001Y\u8ef8\u306e\u6307\u5b9a\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
    geom_point()\u306f\u3001\u6563\u5e03\u56f3\u3092\u63cf\u304b\u305b\u308b\u6307\u4ee4\u3002<\/p>\n\n\n\n
    geom_smooth()\u306f\u3001\u5e73\u6ed1\u7dda\u3092\u63cf\u304b\u305b\u308b\u6307\u4ee4\u3067\u3001method=lm\u306f\u76f4\u7dda\u56de\u5e30\u3001se=TRUE\u306f\u4fe1\u983c\u533a\u9593\u3092\u8868\u793a\u3059\u308b\u6307\u5b9a\u3002<\/p>\n\n\n\n
    \u4fe1\u983c\u533a\u9593\u306f\u3001\u7279\u306b\u6307\u5b9a\u3057\u306a\u3051\u308c\u307095\uff05\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
    \u4ee5\u4e0b\u304c\u7d50\u679c\u3002<\/p>\n\n\n\n
    \u30c7\u30b6\u30a4\u30f3\u6027\u306b\u512a\u308c\u305f\u30b0\u30e9\u30d5\u304c\u3042\u3063\u3068\u3044\u3046\u9593\u306b\u63cf\u3051\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>
    ggplot \u3067\u66f8\u3044\u305f\u56de\u5e30\u76f4\u7dda\u306e\u898b\u305f\u76ee\u3092\u5909\u66f4\u3059\u308b<\/h2>\n\n\n\n
    theme_bw()\u3092\u52a0\u3048\u308b\u3068\u3001\u767d\u9ed2\u306e\u30b7\u30f3\u30d7\u30eb\u306a\u898b\u305f\u76ee\u306b\u306a\u308b\u3002<\/p>\n\n\n\n
    ggplot<\/span>(<\/span>data=<\/span>GCSF_CD34,<\/span> aes<\/span>(<\/span>x=<\/span>WBC.entry,<\/span> y=<\/span>LogCD34perBW))<\/span>+<\/span>\ngeom_point<\/span>()<\/span>+<\/span>\ngeom_smooth<\/span>(<\/span>method=<\/span>lm,<\/span> se=<\/span>TRUE<\/span>)<\/span>+<\/span>\ntheme_bw<\/span>()<\/span>\n<\/code><\/pre>\n\n\n\n
    \"\"<\/figure>\n\n\n\n
    <\/span><\/p>\n\n\n\n
    theme_classic()\u3060\u3068\u3001X\u8ef8\u3001Y\u8ef8\u3060\u3051\u306b\u306a\u308b\u3002<\/p>\n\n\n\n
    ggplot<\/span>(<\/span>data=<\/span>GCSF_CD34,<\/span> aes<\/span>(<\/span>x=<\/span>WBC.entry,<\/span> y=<\/span>LogCD34perBW))<\/span>+<\/span>\ngeom_point<\/span>()<\/span>+<\/span>\ngeom_smooth<\/span>(<\/span>method=<\/span>lm,<\/span> se=<\/span>TRUE<\/span>)<\/span>+<\/span>\ntheme_classic<\/span>()<\/span>\n<\/code><\/pre>\n\n\n\n
    \"\"<\/figure>\n\n\n\n
    <\/span><\/p>\n\n\n\n
    geom_smooth()\u5185\u306bcolor=”black”\u3067\u6307\u5b9a\u3059\u308b\u3068\u9ed2\u306b\u5909\u3048\u3089\u308c\u308b\u3002<\/p>\n\n\n\n
    \u307b\u304b\u306e\u8272\u3082\u82f1\u8a9e\u3067\u6307\u5b9a\u3059\u308c\u3070\u305d\u306e\u8272\u306b\u5909\u308f\u308b\u3002<\/p>\n\n\n\n
    ggplot<\/span>(<\/span>data=<\/span>GCSF_CD34,<\/span> aes<\/span>(<\/span>x=<\/span>WBC.entry,<\/span> y=<\/span>LogCD34perBW))<\/span>+<\/span>\ngeom_point<\/span>()<\/span>+<\/span>\ngeom_smooth<\/span>(<\/span>method=<\/span>lm,<\/span> se=<\/span>TRUE<\/span>,<\/span> color=<\/span>\"black\"<\/span>)<\/span>+<\/span>\ntheme_classic<\/span>()<\/span>\n<\/code><\/pre>\n\n\n\n
    \"\"<\/figure>\n\n\n\n
    <\/span><\/p>\n\n\n\n
    \u4fe1\u983c\u533a\u9593\u306e\u30b0\u30ec\u30fc\u3092\u5c11\u3057\u8584\u304f\u3059\u308b\u305f\u3081\u306b\u306f\u3001geom_smooth()\u5185\u306baes()\u3067fill=’confidence’\u3068\u3044\u3046\u3075\u3046\u306bfill\u3092\u6307\u5b9a\u3059\u308b\u3002<\/p>\n\n\n\n
    \u305d\u3057\u3066scale_fill_manual()\u3067\u3001lightgray\u306b\u5909\u66f4\u3059\u308b\u3002<\/p>\n\n\n\n
    \u305d\u3046\u3059\u308b\u3068\u5c11\u3057\u8584\u3044\u30b0\u30ec\u30fc\u306b\u5909\u308f\u308b\u3002<\/p>\n\n\n\n
    ggplot<\/span>(<\/span>data=<\/span>GCSF_CD34,<\/span> aes<\/span>(<\/span>x=<\/span>WBC.entry,<\/span> y=<\/span>LogCD34perBW))<\/span>+<\/span>\ngeom_point<\/span>()<\/span>+<\/span>\ngeom_smooth<\/span>(<\/span>method=<\/span>lm,<\/span> color=<\/span>\"black\"<\/span>,<\/span>\naes<\/span>(<\/span>fill=<\/span>'confidence'<\/span>),<\/span> show.legend=<\/span>FALSE<\/span>)<\/span>+<\/span>\nscale_fill_manual<\/span>(<\/span>'Interval'<\/span>,<\/span> values=<\/span>'lightgray'<\/span>)<\/span>+<\/span>\ntheme_classic<\/span>()<\/span>\n<\/code><\/pre>\n\n\n\n
    \"\"<\/figure>\n\n\n\n
    <\/span><\/p>\n\n\n\n
    \u4fe1\u983c\u533a\u9593\u3092\u5857\u308a\u3064\u3076\u3055\u305a\u3001\u9650\u754c\u3092\u7dda\u3067\u8868\u793a\u3059\u308b\u65b9\u6cd5\u306f\u308f\u304b\u3089\u306a\u304b\u3063\u305f\u3002<\/p>\n\n\n\n
    ggplot \u3067\u66f8\u3044\u305f\u56de\u5e30\u76f4\u7dda\u306e\u5f0f\u3092\u30b0\u30e9\u30d5\u5185\u306b\u66f8\u304d\u5165\u308c\u308b<\/h2>\n\n\n\n
    \u4f8b\u3048\u3070\u3001\u30b0\u30e9\u30d5\u5185\u306b\u56de\u5e30\u5f0f\u3092\u66f8\u304d\u5165\u308c\u305f\u3044\u5834\u5408\u3001\u3069\u3046\u3059\u308b\u304b\uff1f<\/p>\n\n\n\n
    \u307e\u305a\u5358\u56de\u5e30\u3092\u884c\u3063\u3066\u304a\u304f\u3002<\/p>\n\n\n\n
    #####\u7dda\u5f62\u56de\u5e30(\u5358\u56de\u5e30\u3001\u91cd\u56de\u5e30)#####<\/span>\nlibrary<\/span>(<\/span>aod,<\/span> pos=<\/span>18<\/span>)<\/span>\nRegModel.1 <-<\/span> lm<\/span>(<\/span>LogCD34perBW~<\/span>WBC.entry,<\/span> data=<\/span>GCSF_CD34)<\/span>\n<\/code><\/pre>\n\n\n\n
    \u6b21\u306b\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u56de\u5e30\u5f0f\u3092equation\u3068\u3057\u3066\u5b9a\u7fa9\u3059\u308b\u3002<\/p>\n\n\n\n
    equation <-<\/span> sprintf<\/span>(<\/span>\n\"italic(Y) == %.3g + %.3g * italic(X)\"<\/span>,<\/span>\ncoef<\/span>(<\/span>RegModel.1)[<\/span>1<\/span>],<\/span> coef<\/span>(<\/span>RegModel.1)[<\/span>2<\/span>]<\/span>\n)<\/span>\n<\/code><\/pre>\n\n\n\n
    %.3g \u306e\u90e8\u5206\u306f\u3001\u6b21\u306ecoef(RegModel.1)[1])\uff08\u5207\u7247\u306e\u63a8\u5b9a\u5024\uff09\u3068coef(RegModel.1[2])\uff08\u50be\u304d\u306e\u63a8\u5b9a\u5024\uff09\u3092\u9806\u756a\u306b\u4ee3\u5165\u3059\u308b\u610f\u5473\u3067\u30013g\u306f\u6709\u52b9\u6570\u5b573\u6841\u306b\u3059\u308b\u3068\u3044\u3046\u610f\u5473\u5408\u3044\u3002<\/p>\n\n\n\n
    \u6700\u5f8c\u306b\u3001\u5148\u307b\u3069\u307e\u3067\u306e\u30b9\u30af\u30ea\u30d7\u30c8\u306e\u6700\u5f8c\u306b\u4e00\u884c\u3092\u66f8\u304d\u8db3\u3059\u3002<\/p>\n\n\n\n
    \uff0b\u3067\u3064\u306a\u3050\u306e\u3092\u5fd8\u308c\u305a\u306b\u3002<\/p>\n\n\n\n
    ggplot<\/span>(<\/span>data=<\/span>GCSF_CD34,<\/span> aes<\/span>(<\/span>x=<\/span>WBC.entry,<\/span> y=<\/span>LogCD34perBW))<\/span>+<\/span>\ngeom_point<\/span>()<\/span>+<\/span>\ngeom_smooth<\/span>(<\/span>method=<\/span>lm,<\/span> color=<\/span>\"black\"<\/span>,<\/span>\naes<\/span>(<\/span>fill=<\/span>'confidence'<\/span>),<\/span> show.legend=<\/span>FALSE<\/span>)<\/span>+<\/span>\nscale_fill_manual<\/span>(<\/span>'Interval'<\/span>,<\/span> values=<\/span>'lightgray'<\/span>)<\/span>+<\/span>\ntheme_classic<\/span>()<\/span>+<\/span>\nannotate<\/span>(<\/span>\"text\"<\/span>,<\/span> x=<\/span>5000<\/span>,<\/span> y=<\/span>1.0<\/span>,<\/span> label=<\/span>equation,<\/span> parse=<\/span>TRUE<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
    X\u8ef8\u304c5000\u3001Y\u8ef8\u304c1.0\u306e\u5ea7\u6a19\u306b\u3001\u56de\u5e30\u5f0f\u3092\u66f8\u304d\u5165\u308c\u308b\u305f\u3081\u306e\u30b9\u30af\u30ea\u30d7\u30c8\u3002<\/p>\n\n\n\n
    parse=TRUE\u306f\u3001\u6570\u5b66\u8868\u73fe\u3067\u3042\u308b\u3053\u3068\u306e\u6307\u793a\u3002<\/p>\n\n\n\n
    \u30b0\u30e9\u30d5\u3092\u8868\u793a\u3055\u305b\u308b\u3068\u3053\u306e\u3088\u3046\u306b\u306a\u308b\u3002<\/p>\n\n\n\n
    \"\"<\/figure>\n\n\n\n
    <\/span><\/p>\n\n\n\n
    \u307e\u3068\u3081<\/h2>\n\n\n\n
    EZR\u306795\uff05\u4fe1\u983c\u533a\u9593\u4ed8\u304d\u306e\u56de\u5e30\u76f4\u7dda\u3092\u63cf\u304f\u306b\u306f\u3001ggplot2\u3092\u4f7f\u3046\u3068\u3088\u3044\u3002<\/p>\n\n\n\n
    \u7c21\u5358\u306b\u304d\u308c\u3044\u306a\u56f3\u304c\u63cf\u3051\u308b\u3002<\/p>\n\n\n\n
    \u53c2\u8003\u30b5\u30a4\u30c8<\/h2>\n\n\n\n
    RPubs – Confidence and Prediction intervals for Linear Regression<\/a><\/p>\n\n\n\n
    ggplot2: \u304d\u308c\u3044\u306a\u30b0\u30e9\u30d5\u3092\u7c21\u5358\u306b\u5408\u7406\u7684\u306b – Heavy Watal<\/a><\/p>\n\n\n\n
    ggplot2\u3092\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3057\u7f8e\u3057\u3044\u30b0\u30e9\u30d5\u3092\u4f5c\u308b\u307e\u3067\u306e\u57fa\u790e\u77e5\u8b58\u3010\u5165\u9580\u7de8\u3011<\/a><\/p>\n\n\n\n
    statistics – creating confidence area for normally distributed scatterplot in ggplot2 and R – Stack Overflow<\/a><\/p>\n\n\n\n
    5.9 Adding Annotations with Model Coefficients | R Graphics Cookbook, 2nd edition<\/a><\/p>\n\n\n\n
    \u304a\u3059\u3059\u3081\u66f8\u7c4d<\/h2>\n\n\n\n
    \"EZR\u3067\u3084\u3055\u3057\u304f\u5b66\u3076\u7d71\u8a08\u5b66<\/a><\/figure>\n\n\n\n
    <\/div>\n\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