{"id":380,"date":"2021-11-07T12:33:00","date_gmt":"2021-11-07T03:33:00","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/increase-decimal-places-in-conf-int-of-prop-diff-ezr\/"},"modified":"2024-10-03T12:17:58","modified_gmt":"2024-10-03T03:17:58","slug":"increase-decimal-places-in-conf-int-of-prop-diff-ezr","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/increase-decimal-places-in-conf-int-of-prop-diff-ezr\/","title":{"rendered":"EZR \u306e 2 \u7fa4\u306e\u6bd4\u7387\u306e\u5dee\u306e\u4fe1\u983c\u533a\u9593\u8a08\u7b97\u3067\u5c0f\u6570\u70b9\u4ee5\u4e0b\u306e\u6841\u3092\u5897\u3084\u3059\u65b9\u6cd5"},"content":{"rendered":"\n

EZR\u30672\u7fa4\u306e\u6bd4\u7387\u306e\u5dee\u306e\u4fe1\u983c\u533a\u9593\u3092\u6c42\u3081\u308b\u3068\u5c0f\u6570\u70b9\u4ee5\u4e0b\uff13\u6841\u306b\u4e38\u3081\u3089\u308c\u305f\u6570\u5024\u304c\u51fa\u529b\u3055\u308c\u308b\u3002<\/p>\n\n\n\n

\u5c0f\u6570\u70b9\u4ee5\u4e0b\u306e\u6841\u6570\u3092\u5897\u3084\u3059\u65b9\u6cd5\u3002<\/p>\n\n\n\n\n\n\n\n

EZR\u3067\uff12\u7fa4\u306e\u6bd4\u7387\u306e\u5dee\u306e\u4fe1\u983c\u533a\u9593\u3092\u6c42\u3081\u308b<\/h2>\n\n\n\n

\u300c\u7d71\u8a08\u89e3\u6790\u300d\u2192\u300c\u540d\u7fa9\u5909\u6570\u306e\u89e3\u6790\u300d\u2192\u300c2\u7fa4\u306e\u6bd4\u7387\u306e\u5dee\u306e\u4fe1\u983c\u533a\u9593\u306e\u8a08\u7b97\u300d\u30e1\u30cb\u30e5\u30fc\u3092\u4f7f\u3063\u3066\u3001\u8a08\u7b97\u3059\u308b\u3002<\/p>\n\n\n\n

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

\u4f8b\u3048\u3070\u300147000\u5206\u306e460\u306839000\u5206\u306e200\u306e\u5dee\u306e95\uff05\u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97\u3059\u308b\u3068\u3001R\u30b9\u30af\u30ea\u30d7\u30c8\u7a93\u306b\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u66f8\u304b\u308c\u3066\u3001\u51fa\u529b\u7a93\u306b\u7d50\u679c\u304c\u8868\u793a\u3055\u308c\u308b\u3002<\/p>\n\n\n\n

R\u30b9\u30af\u30ea\u30d7\u30c8<\/p>\n\n\n\n

#####2\u7fa4\u306e\u6bd4\u7387\u306e\u5dee\u306e\u4fe1\u983c\u533a\u9593\u306e\u8a08\u7b97#####<\/span>\nprop.diff.conf<\/span>(<\/span>460<\/span>,<\/span> 47000<\/span>,<\/span> 200<\/span>,<\/span> 39000<\/span>,<\/span> 95<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
\u51fa\u529b<\/p>\n\n\n\n
><\/span> prop.diff.conf<\/span>(<\/span>460<\/span>,<\/span> 47000<\/span>,<\/span> 200<\/span>,<\/span> 39000<\/span>,<\/span> 95<\/span>)<\/span>\n[<\/span>1<\/span>]<\/span> \u6bd4\u7387\u306e\u5dee     :<\/span> 0.005<\/span>\n[<\/span>1<\/span>]<\/span> 95<\/span>% \u4fe1\u983c\u533a\u9593              :<\/span> 0.004<\/span> -<\/span> 0.006<\/span>\n<\/code><\/pre>\n\n\n\n
EZR\u3067\uff12\u7fa4\u306e\u6bd4\u7387\u306e\u5dee\u306e\u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97\u3059\u308b\u30b9\u30af\u30ea\u30d7\u30c8\u3092\u78ba\u8a8d\u3057\u3066\u4e00\u90e8\u5909\u66f4\u3059\u308b<\/h2>\n\n\n\n
R\u30b9\u30af\u30ea\u30d7\u30c8\u7a93\u306b\u3001prop.diff.conf \u3068\u66f8\u3044\u3066\u3001\u5b9f\u884c\u3057\u3001prop.diff.conf\u95a2\u6570\u306e\u4e2d\u8eab\u3092\u8868\u793a\u3055\u305b\u308b\u3002<\/p>\n\n\n\n
prop.diff.conf()\u95a2\u6570\u306e\u4e2d\u8eab\u306f\u4ee5\u4e0b\u306e\u901a\u308a\u3002<\/p>\n\n\n\n
><\/span> prop.diff.conf\nfunction<\/span> (<\/span>r1,<\/span> n1,<\/span> r2,<\/span> n2,<\/span> conf)<\/span>\n{<\/span>\nalpha <-<\/span> 1<\/span> -<\/span> conf\/<\/span>100<\/span>\np1 <-<\/span> r1\/<\/span>n1\np2 <-<\/span> r2\/<\/span>n2\nD <-<\/span> p1 -<\/span> p2\nSE <-<\/span> sqrt<\/span>(<\/span>p1 *<\/span> (<\/span>1<\/span> -<\/span> p1)<\/span>\/<\/span>n1 +<\/span> p2 *<\/span> (<\/span>1<\/span> -<\/span> p2)<\/span>\/<\/span>n2)<\/span>\npl <-<\/span> D -<\/span> qnorm<\/span>(<\/span>1<\/span> -<\/span> alpha\/<\/span>2<\/span>)<\/span> *<\/span> SE\npu <-<\/span> D +<\/span> qnorm<\/span>(<\/span>1<\/span> -<\/span> alpha\/<\/span>2<\/span>)<\/span> *<\/span> SE\nprint<\/span>(<\/span>paste<\/span>(<\/span>gettext<\/span>(<\/span>domain =<\/span> \"R-RcmdrPlugin.EZR\"<\/span>,<\/span> \"Difference :\"<\/span>),<\/span>\n\" \"<\/span>,<\/span> round<\/span>(<\/span>D,<\/span> 3<\/span>),<\/span> sep =<\/span> \"\"<\/span>),<\/span> quote =<\/span> F<\/span>)<\/span>\nprint<\/span>(<\/span>paste<\/span>(<\/span>conf,<\/span> gettext<\/span>(<\/span>domain =<\/span> \"R-RcmdrPlugin.EZR\"<\/span>,<\/span> \"% confidence interval :\"<\/span>),<\/span>\n\" \"<\/span>,<\/span> round<\/span>(<\/span>pl,<\/span> 3<\/span>),<\/span> \" - \"<\/span>,<\/span> round<\/span>(<\/span>pu,<\/span> 3<\/span>),<\/span> sep =<\/span> \"\"<\/span>),<\/span> quote =<\/span> F<\/span>)<\/span>\n}<\/span>\n<\/code><\/pre>\n\n\n\n
\u4e0b\u304b\u3089 2 \u884c\u76ee\u3068 4 \u884c\u76ee\u306b\u3001round(, 3)\u3068\u3042\u308b\u3068\u3053\u308d\u8a08 3\u304b\u6240\u3092\u3001\u4f8b\u3048\u3070round(,8)\u3068\u3057\u3066\u3001\u5c0f\u6570\u70b9\u4ee5\u4e0b8\u6841\u307e\u3067\u8868\u793a\u3055\u305b\u308b\u3088\u3046\u306b\u5909\u66f4\u3059\u308b\u3002<\/p>\n\n\n\n
function<\/span> (<\/span>r1,<\/span> n1,<\/span> r2,<\/span> n2,<\/span> conf)<\/span>\n{<\/span>\nalpha <-<\/span> 1<\/span> -<\/span> conf\/<\/span>100<\/span>\np1 <-<\/span> r1\/<\/span>n1\np2 <-<\/span> r2\/<\/span>n2\nD <-<\/span> p1 -<\/span> p2\nSE <-<\/span> sqrt<\/span>(<\/span>p1 *<\/span> (<\/span>1<\/span> -<\/span> p1)<\/span>\/<\/span>n1 +<\/span> p2 *<\/span> (<\/span>1<\/span> -<\/span> p2)<\/span>\/<\/span>n2)<\/span>\npl <-<\/span> D -<\/span> qnorm<\/span>(<\/span>1<\/span> -<\/span> alpha\/<\/span>2<\/span>)<\/span> *<\/span> SE\npu <-<\/span> D +<\/span> qnorm<\/span>(<\/span>1<\/span> -<\/span> alpha\/<\/span>2<\/span>)<\/span> *<\/span> SE\nprint<\/span>(<\/span>paste<\/span>(<\/span>gettext<\/span>(<\/span>domain =<\/span> \"R-RcmdrPlugin.EZR\"<\/span>,<\/span> \"Difference :\"<\/span>),<\/span>\n\" \"<\/span>,<\/span> round<\/span>(<\/span>D,<\/span> 8<\/span>),<\/span> sep =<\/span> \"\"<\/span>),<\/span> quote =<\/span> F<\/span>)<\/span>\nprint<\/span>(<\/span>paste<\/span>(<\/span>conf,<\/span> gettext<\/span>(<\/span>domain =<\/span> \"R-RcmdrPlugin.EZR\"<\/span>,<\/span> \"% confidence interval :\"<\/span>),<\/span>\n\" \"<\/span>,<\/span> round<\/span>(<\/span>pl,<\/span> 8<\/span>),<\/span> \" - \"<\/span>,<\/span> round<\/span>(<\/span>pu,<\/span> 8<\/span>),<\/span> sep =<\/span> \"\"<\/span>),<\/span> quote =<\/span> F<\/span>)<\/span>\n}<\/span>\n<\/code><\/pre>\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>
EZR\u3067\uff12\u7fa4\u306e\u6bd4\u7387\u306e\u5dee\u306e\u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97\u3059\u308b\u4fee\u6b63\u30b9\u30af\u30ea\u30d7\u30c8\u3092\u5b9f\u884c\u3059\u308b<\/h2>\n\n\n\n
\u65b0\u3057\u3044\u95a2\u6570\u540d\u3001\u4f8b\u3048\u3070 prop.diff.conf.rev \u3068\u3044\u3046\u540d\u524d\u306e\u95a2\u6570\u306b\u4e0a\u8a18\u306e\u5185\u5bb9\u3092\u683c\u7d0d\u3059\u308b\u3068\u3044\u3046\u30b9\u30af\u30ea\u30d7\u30c8\u3092R\u30b9\u30af\u30ea\u30d7\u30c8\u7a93\u306b\u8cbc\u308a\u4ed8\u3051\u3066\u5b9f\u884c\u3059\u308b\u3002<\/p>\n\n\n\n
prop.diff.conf.rev <-<\/span>\nfunction<\/span> (<\/span>r1,<\/span> n1,<\/span> r2,<\/span> n2,<\/span> conf)<\/span>\n{<\/span>\nalpha <-<\/span> 1<\/span> -<\/span> conf\/<\/span>100<\/span>\np1 <-<\/span> r1\/<\/span>n1\np2 <-<\/span> r2\/<\/span>n2\nD <-<\/span> p1 -<\/span> p2\nSE <-<\/span> sqrt<\/span>(<\/span>p1 *<\/span> (<\/span>1<\/span> -<\/span> p1)<\/span>\/<\/span>n1 +<\/span> p2 *<\/span> (<\/span>1<\/span> -<\/span> p2)<\/span>\/<\/span>n2)<\/span>\npl <-<\/span> D -<\/span> qnorm<\/span>(<\/span>1<\/span> -<\/span> alpha\/<\/span>2<\/span>)<\/span> *<\/span> SE\npu <-<\/span> D +<\/span> qnorm<\/span>(<\/span>1<\/span> -<\/span> alpha\/<\/span>2<\/span>)<\/span> *<\/span> SE\nprint<\/span>(<\/span>paste<\/span>(<\/span>gettext<\/span>(<\/span>domain =<\/span> \"R-RcmdrPlugin.EZR\"<\/span>,<\/span> \"Difference :\"<\/span>),<\/span>\n\" \"<\/span>,<\/span> round<\/span>(<\/span>D,<\/span> 8<\/span>),<\/span> sep =<\/span> \"\"<\/span>),<\/span> quote =<\/span> F<\/span>)<\/span>\nprint<\/span>(<\/span>paste<\/span>(<\/span>conf,<\/span> gettext<\/span>(<\/span>domain =<\/span> \"R-RcmdrPlugin.EZR\"<\/span>,<\/span> \"% confidence interval :\"<\/span>),<\/span>\n\" \"<\/span>,<\/span> round<\/span>(<\/span>pl,<\/span> 8<\/span>),<\/span> \" - \"<\/span>,<\/span> round<\/span>(<\/span>pu,<\/span> 8<\/span>),<\/span> sep =<\/span> \"\"<\/span>),<\/span> quote =<\/span> F<\/span>)<\/span>\n}<\/span>\n<\/code><\/pre>\n\n\n\n
\u65b0\u3057\u3044\u95a2\u6570 prop.diff.conf.rev() \u3092\u4f7f\u3063\u3066\u3001\u5148\u307b\u3069\u3068\u540c\u3058\u8a08\u7b97\u3092\u3057\u3066\u307f\u308b\u3002<\/p>\n\n\n\n
R\u30b9\u30af\u30ea\u30d7\u30c8\u7a93\u306b\u3001prop.diff.conf.rev(460, 47000, 200, 39000, 95)<\/code> \u3068\u66f8\u3044\u3066\u3001\u5b9f\u884c\u3092\u30af\u30ea\u30c3\u30af\u3059\u308b\u3002<\/p>\n\n\n\n
\u3059\u308b\u3068\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u7d50\u679c\u304c\u8868\u793a\u3055\u308c\u308b\u3002\u5c0f\u6570\u70b9\u4ee5\u4e0b\uff18\u6841\u307e\u3067\u8868\u793a\u3055\u308c\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
><\/span> prop.diff.conf.rev<\/span>(<\/span>460<\/span>,<\/span> 47000<\/span>,<\/span> 200<\/span>,<\/span> 39000<\/span>,<\/span> 95<\/span>)<\/span>\n[<\/span>1<\/span>]<\/span> \u6bd4\u7387\u306e\u5dee     :<\/span> 0.00465903<\/span>\n[<\/span>1<\/span>]<\/span> 95<\/span>% \u4fe1\u983c\u533a\u9593              :<\/span> 0.0035212<\/span> -<\/span> 0.00579685<\/span>\n<\/code><\/pre>\n\n\n\n
\u307e\u3068\u3081<\/h2>\n\n\n\n
EZR\u3067\uff12\u7fa4\u306e\u6bd4\u7387\u306e\u5dee\u306e\u4fe1\u983c\u533a\u9593\u3092\u8a08\u7b97\u3059\u308b\u5834\u5408\u3001\u7d50\u679c\u306e\u5024\u306e\u5c0f\u6570\u70b9\u4ee5\u4e0b\u306e\u6841\u6570\u3092\u5909\u66f4\u3059\u308b\u306b\u306f\u3001\u30e1\u30cb\u30e5\u30fc\u3084\u30c0\u30a4\u30a2\u30ed\u30b0\u30dc\u30c3\u30af\u30b9\u304b\u3089\u306f\u5909\u66f4\u3067\u304d\u306a\u3044\u304c\u3001R\u30b9\u30af\u30ea\u30d7\u30c8\u3092\u5c11\u3057\u4fee\u6b63\u3059\u308b\u3068\u53ef\u80fd\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\u304a\u3059\u3059\u3081\u66f8\u7c4d<\/h2>\n\n\n\n
EZR\u516c\u5f0f\u30de\u30cb\u30e5\u30a2\u30eb<\/p>\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