{"id":283,"date":"2022-07-17T22:09:10","date_gmt":"2022-07-17T13:09:10","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/paired-dichotomous-data-comparing-groups-analyzed-by-glmm-in-ezr\/"},"modified":"2024-09-23T15:19:29","modified_gmt":"2024-09-23T06:19:29","slug":"paired-dichotomous-data-comparing-groups-analyzed-by-glmm-in-ezr","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/paired-dichotomous-data-comparing-groups-analyzed-by-glmm-in-ezr\/","title":{"rendered":"EZR \u3067\u5bfe\u5fdc\u306e\u3042\u308b 2 \u5024\u30ab\u30c6\u30b4\u30ea\u30ab\u30eb\u30c7\u30fc\u30bf\u3092\u7fa4\u9593\u6bd4\u8f03\u3059\u308b\u65b9\u6cd5"},"content":{"rendered":"\n

\u5bfe\u5fdc\u306e\u3042\u308b\uff12\u5024\u30c7\u30fc\u30bf\u3092\u7fa4\u9593\u6bd4\u8f03\u3057\u305f\u3044\u3068\u304d\u3069\u3046\u3059\u308c\u3070\u3088\u3044\u306e\u304b\uff1f<\/p>\n\n\n\n\n\n\n\n

\u5bfe\u5fdc\u306e\u3042\u308b\uff12\u5024\u30c7\u30fc\u30bf\u306e\u7fa4\u9593\u6bd4\u8f03\u3068\u306f\u3069\u3046\u3044\u3046\u72b6\u6cc1\u304b\uff1f<\/h2>\n\n\n\n

\u5bfe\u5fdc\u306e\u3042\u308b\uff12\u5024\u30c7\u30fc\u30bf\u306e\u7fa4\u9593\u6bd4\u8f03\u3068\u306f\u3069\u3046\u3044\u3046\u72b6\u6cc1\u3060\u308d\u3046\u304b\u3002<\/p>\n\n\n\n

\u305d\u308c\u306f\u4ecb\u5165\u524d\u5f8c\u3067\u3001Yes\/No\u3068\u304b\u3001\u3042\u308a\u30fb\u306a\u3057\u3068\u304b\u3001\u305d\u3046\u3044\u3046\uff12\u5024\u306e\u30c7\u30fc\u30bf\u304c\u3042\u3063\u3066\u3001\u3055\u3089\u306b\u4ecb\u5165\u304c\u7570\u306a\u308b\u7fa4\u304c\u8907\u6570\u5b58\u5728\u3059\u308b\u5834\u5408\u3069\u3046\u3059\u308b\u304b\u3068\u3044\u3046\u3053\u3068\u3067\u3042\u308b\u3002<\/p>\n\n\n\n

\u7fa4\u9593\u6bd4\u8f03\u304c\u306a\u3051\u308c\u3070\u30de\u30af\u30cd\u30de\u30fc\u691c\u5b9a\u3092\u5b9f\u65bd\u3057\u305f\u3044\u72b6\u6cc1\u3067\u3042\u308b\u3002<\/p>\n\n\n\n

\u5bfe\u5fdc\u306e\u3042\u308b\uff12\u5024\u30c7\u30fc\u30bf\u3092\u7fa4\u9593\u6bd4\u8f03\u306b\u306f\u4e00\u822c\u5316\u7dda\u5f62\u6df7\u5408\u30e2\u30c7\u30eb\u3092\u4f7f\u3046<\/h2>\n\n\n\n

\u5bfe\u5fdc\u306e\u3042\u308b\uff12\u5024\u30c7\u30fc\u30bf\u3092\u7fa4\u9593\u6bd4\u8f03\u3059\u308b\u4e00\u3064\u306e\u65b9\u6cd5\u306f\u3001\u4e00\u822c\u5316\u7dda\u5f62\u6df7\u5408\u30e2\u30c7\u30eb\u3092\u7528\u3044\u305f\u65b9\u6cd5\u3067\u3042\u308b\u3002<\/p>\n\n\n\n

\u7dda\u5f62\u6df7\u5408\u30e2\u30c7\u30eb\u306f\u3001\u5bfe\u5fdc\u3042\u308b\u9023\u7d9a\u30c7\u30fc\u30bf\u3092\u8907\u6570\u56de\u6e2c\u5b9a\u3057\u305f\u5834\u5408\u306e\u30e2\u30c7\u30eb\u3067\u3001\u305d\u308c\u3092\uff12\u5024\u306e\u7e70\u308a\u8fd4\u3057\u6e2c\u5b9a\u306b\u3082\u5bfe\u5fdc\u3057\u305f\u3082\u306e\u304c\u4e00\u822c\u5316\u7dda\u5f62\u6df7\u5408\u30e2\u30c7\u30eb\u3067\u3042\u308b\u3002<\/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>

\u5bfe\u5fdc\u306e\u3042\u308b\uff12\u5024\u30c7\u30fc\u30bf\u306e\u7fa4\u9593\u6bd4\u8f03\u30b5\u30f3\u30d7\u30eb\u30c7\u30fc\u30bf<\/h2>\n\n\n\n

\u5bfe\u5fdc\u3042\u308b\uff12\u5024\u30c7\u30fc\u30bf\u3067\u7fa4\u9593\u6bd4\u8f03\u3067\u304d\u308b\u30b5\u30f3\u30d7\u30eb\u30c7\u30fc\u30bf\u3092\u751f\u6210\u3057\u3066\u307f\u305f\u3002<\/p>\n\n\n\n

\u4e00\u69d8\u5206\u5e03\u306e\u30e9\u30f3\u30c0\u30e0\u5909\u6570\u3092\u767a\u751f\u3055\u305b\u3066\u300150 \u4f8b 2 \u7fa4\u3067\u3001\u3042\u308a\u306a\u3057\u30fbYesNo\u306a\u3069 2 \u5024\u306e\u8a55\u4fa1\u3092\u3001\u4ecb\u5165\u7b49\u306e\u4e8b\u524d\u4e8b\u5f8c\u306e 2 \u56de\u7e70\u308a\u8fd4\u3057\u305f\u5834\u5408\u306e\u30b5\u30f3\u30d7\u30eb\u30c7\u30fc\u30bf\u3092\u4f5c\u6210\u3057\u305f\u3002<\/p>\n\n\n\n

\u4e0b\u8a18\u306f R \u306e\u30b9\u30af\u30ea\u30d7\u30c8\u3060\u304c\u3001GLMM1.csv\u3068\u3044\u3046\u540d\u524d\u3067\u66f8\u304d\u51fa\u3057
\u3066\u3044\u308b\u3002<\/p>\n\n\n\n

\u66f8\u304d\u51fa\u3059\u30d1\u30b9\u306f\u3001PC\u306e\u74b0\u5883\u306b\u3088\u3063\u3066\u7570\u306a\u308b\u306e\u3067\u3001~~~\u3068\u66f8\u3044\u3066\u3044\u308b\u3002<\/p>\n\n\n\n

group <-<\/span> c<\/span>(<\/span>rep<\/span>(<\/span>1<\/span>,<\/span>50<\/span>),<\/span>rep<\/span>(<\/span>2<\/span>,<\/span>50<\/span>))<\/span>\npre <-<\/span> round<\/span>(<\/span>runif<\/span>(<\/span>100<\/span>))<\/span>\npost <-<\/span> round<\/span>(<\/span>runif<\/span>(<\/span>100<\/span>))<\/span>\ndat <-<\/span> data.frame<\/span>(<\/span>group,<\/span> pre,<\/span> post)<\/span>\nwrite.csv<\/span>(<\/span>dat,<\/span> \"C:<\/span>\\\\<\/span>Users<\/span>\\\\<\/span>~~~<\/span>\\\\<\/span>GLMM1.csv\"<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
\u3053\u308c\u3092\u7e26\u306b\u7a4d\u307f\u91cd\u306d\u305f\u30c7\u30fc\u30bf\u30bb\u30c3\u30c8\u3092\u4f5c\u6210\u3059\u308b\u5fc5\u8981\u304c\u3042\u308b\u3002<\/p>\n\n\n\n
\u4ee5\u4e0b\u306fEZR\u306bCSV\u30d5\u30a1\u30a4\u30eb\u3092\u53d6\u308a\u8fbc\u3093\u3060\u5f8c\u3001\u7e26\u306b\u7a4d\u307f\u91cd\u306d\u305f\u30c7\u30fc\u30bf\u3092\u4f5c\u6210\u3057\u305f\u969b\u306e\u30b9\u30af\u30ea\u30d7\u30c8\u3092\u53c2\u8003\u307e\u3067\u306b\u8cbc\u308a\u4ed8\u3051\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
EZR\u3067\u306f\u3001\u30a2\u30af\u30c6\u30a3\u30d6\u30c7\u30fc\u30bf\u30bb\u30c3\u30c8 \u2192 \u5909\u6570\u306e\u64cd\u4f5c \u2192 \u8907\u6570\u306e\u5909\u6570\u3092\u7e26\u306b\u7a4d\u307f\u91cd\u306d\u305f\u30c7\u30fc\u30bf\u30bb\u30c3\u30c8\u3092\u4f5c\u6210\u3059\u308b \u30e1\u30cb\u30e5\u30fc\u3092\u7528\u3044\u308b\u3068\u3001\u7e26\u306b\u7a4d\u307f\u91cd\u306d\u305f\u30c7\u30fc\u30bf\u30bb\u30c3\u30c8\u3092\u4f5c\u6210\u3067\u304d\u308b\u3002<\/p>\n\n\n\n
#####\u8907\u6570\u306e\u5909\u6570\u3092\u7e26\u306b\u7a4d\u307f\u91cd\u306d\u305f\u30c7\u30fc\u30bf\u30bb\u30c3\u30c8\u3092\u4f5c\u6210\u3059\u308b#####<\/span>\nTempDF <-<\/span>  GLMM1\n#Convert to long format<\/span>\nn <-<\/span> length<\/span>(<\/span>TempDF[,<\/span>1<\/span>])<\/span>\nTempDF2 <-<\/span> reshape<\/span>(<\/span>TempDF,<\/span> varying=<\/span>list<\/span>(<\/span>c<\/span>(<\/span>\"pre\"<\/span>,<\/span> \"post\"<\/span>)),<\/span>\nv.names=<\/span>\"outcome\"<\/span>,<\/span> timevar=<\/span>\"time\"<\/span>,<\/span> direction=<\/span>\"long\"<\/span>)<\/span>\nRepeatNumber <-<\/span> c<\/span>(<\/span>\"pre\"<\/span>,<\/span> \"post\"<\/span>)<\/span>\nTempDF2$<\/span>time <-<\/span> RepeatNumber[<\/span>TempDF2$<\/span>time]<\/span>\nGLMMvert1 <-<\/span> TempDF2\n#####\u9023\u7d9a\u5909\u6570\u3092\u56e0\u5b50\u306b\u5909\u63db\u3059\u308b#####<\/span>\nGLMMvert1$<\/span>group <-<\/span> as.factor<\/span>(<\/span>GLMMvert1$<\/span>group)<\/span>\n<\/code><\/pre>\n\n\n\n
\u5bfe\u5fdc\u306e\u3042\u308b\uff12\u5024\u30c7\u30fc\u30bf\u306e\u7fa4\u9593\u6bd4\u8f03\u3092\u4e00\u822c\u5316\u7dda\u5f62\u6df7\u5408\u30e2\u30c7\u30eb\u3067\u5b9f\u65bd\u3059\u308b<\/h2>\n\n\n\n
EZR\u306e\u6a19\u6e96\u30e1\u30cb\u30e5\u30fc\u304b\u3089\u3001\u7d71\u8a08\u91cf \u2192 \u30e2\u30c7\u30eb\u3078\u306e\u9069\u5408 \u2192 \u4e00\u822c\u5316\u7dda\u5f62\u6df7\u5408\u30e2\u30c7\u30eb\u3068\u9032\u3080\u3068\u89e3\u6790\u3067\u304d\u308b\u3002<\/p>\n\n\n\n
\u6700\u521d\u306blibrary(lme4)\u3067\u3001lme4\u30d1\u30c3\u30b1\u30fc\u30b8\u3092\u547c\u3073\u51fa\u3057\u3066\u304a\u304f\u3002<\/p>\n\n\n\n
\u30e1\u30cb\u30e5\u30fc\u304b\u3089\u89e3\u6790\u3059\u308b\u3068\u4ee5\u4e0b\u306e\u30b9\u30af\u30ea\u30d7\u30c8\u304cR\u30b9\u30af\u30ea\u30d7\u30c8\u7a93\u306b\u51fa\u529b\u3055\u308c\u308b\u3002<\/p>\n\n\n\n
\u30e2\u30c7\u30eb\u306e\u6307\u5b9a\u306f\u3001\u5f93\u5c5e\u5909\u6570\u304coutcome, \u72ec\u7acb\u5909\u6570\u5074\u306e\u67a0\u306b\u306f\u3001group + time + (1|id) \u3068\u66f8\u304d\u5165\u308c\u308b\u3002<\/p>\n\n\n\n
group\u306f\u7fa4\u306e\u5909\u6570\u3001time\u304c\u524d\u304b\u5f8c\u304b\u306e\u6642\u70b9\u306e\u5909\u6570\u3001id\u306f\u540c\u4e00\u4f8b\u3067\u3042\u308b\u3053\u3068\u3092\u793a\u3059\u6307\u6a19\u3067\u3001\u3053\u308c\u304c\u5909\u91cf\u52b9\u679c\u3068\u3057\u3066\u6295\u5165\u3055\u308c\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
library<\/span>(<\/span>lme4)<\/span>\nGLMM.2 <-<\/span> glmer<\/span>(<\/span>outcome ~<\/span> group +<\/span> time +<\/span> (<\/span>1<\/span> |<\/span> id),<\/span> family=<\/span>binomial<\/span>(<\/span>logit),<\/span>\ndata=<\/span>GLMMvert1)<\/span>\nsummary<\/span>(<\/span>GLMM.2)<\/span>\nexp<\/span>(<\/span>coef<\/span>(<\/span>GLMM.2))<\/span>  # Exponentiated coefficients (\"odds ratios\")<\/span>\n<\/code><\/pre>\n\n\n\n
\u7d50\u679c\u306e\u30ad\u30e2\u306e\u90e8\u5206\u306f\u4ee5\u4e0b\u306e\u3068\u304a\u308a\u3067\u3042\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
pre\u3068post\u304c\u9006\u306b\u306a\u3063\u3066\u3057\u307e\u3063\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
post\u306e\u63a8\u5b9a\u5024\u304c\u51fa\u529b\u3055\u308c\u308b\u3088\u3046\u306b\u3001\u30a2\u30af\u30c6\u30a3\u30d6\u30c7\u30fc\u30bf\u30bb\u30c3\u30c8 \u2192 \u5909\u6570\u306e\u64cd\u4f5c \u2192 \u56e0\u5b50\u6c34\u6e96\u3092\u518d\u9806\u5e8f\u5316\u3059\u308b \u30e1\u30cb\u30e5\u30fc\u3067\u5909\u66f4\u3057\u3066\u304a\u304f\u3068\u3088\u304b\u3063\u305f\u3002<\/p>\n\n\n\n
\u3044\u305a\u308c\u306b\u3057\u3066\u3082 group \u306e\u52b9\u679c\u306e\u691c\u5b9a\u3068\u3001time \u306e\u52b9\u679c\u306e\u691c\u5b9a\u304c\u540c\u6642\u306b\u5b9f\u65bd\u3055\u308c\u3066\u3044\u308b\u3002<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u63a8\u5b9a\u5024\u3092\u30aa\u30c3\u30ba\u6bd4\u306b\u5909\u66f4\u3057\u305f\u5024\u3082\u8a08\u7b97\u3055\u308c\u308b\u3002<\/p>\n\n\n\n
95 % \u4fe1\u983c\u533a\u9593\u306f\u3001confint() \u3067\u8a08\u7b97\u3057\u3066\u3001\u771f\u6570\u5909\u63db exp() \u3067\u3001\u7b97\u51fa\u3067\u304d\u308b\uff08\u4ee5\u4e0b\u3001GLMM.1 \u306f\u4e0a\u8a18\u306e GLMM.2 \u3068\u540c\u3058\u3082\u306e\u3067\u3042\u308b\uff09<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u3082\u3057\u3082\u3001\u3053\u306e\u65b9\u6cd5\u3067 95 \uff05 \u4fe1\u983c\u533a\u9593\u304c\u8a08\u7b97\u3067\u304d\u306a\u3044\u5834\u5408\u306f\u3001\u4ee5\u4e0b\u306e\u30b9\u30af\u30ea\u30d7\u30c8\u3092\u4f7f\u3063\u3066\u624b\u52d5\u3067\u8a08\u7b97\u3059\u308b\u3053\u3068\u304c\u3067\u304d\u308b<\/p>\n\n\n\n
se <-<\/span> sqrt<\/span>(<\/span>diag<\/span>(<\/span>vcov<\/span>(<\/span>GLMM.1)))<\/span>\nupperCI <-<\/span> exp<\/span>(<\/span>coef<\/span>(<\/span>GLMM.1)<\/span> +<\/span> 1.96<\/span>*<\/span>se)<\/span>\nlowerCI <-<\/span> exp<\/span>(<\/span>coef<\/span>(<\/span>GLMM.1)<\/span> -<\/span> 1.96<\/span>*<\/span>se)<\/span>\n<\/code><\/pre>\n\n\n\n
vcov \u306f\u5206\u6563\u30fb\u5171\u5206\u6563\u884c\u5217\u3092\u6c42\u3081\u308b\u95a2\u6570\u3001diag \u306f\u884c\u5217\u306e\u5bfe\u89d2\u6210\u5206\u3092\u629c\u304d\u51fa\u3059\u95a2\u6570\u3067\u3001\u3053\u308c\u3067\u5206\u6563\u304c\u53d6\u308a\u51fa\u305b\u3066\u3001\u5e73\u65b9\u6839 sqrt \u3092\u8a08\u7b97\u3059\u308b\u3068\u6a19\u6e96\u8aa4\u5dee\u304c\u5f97\u3089\u308c\u308b<\/p>\n\n\n\n
\u6a19\u6e96\u8aa4\u5dee\u306b 1.96 \u304b\u3051\u305f\u5024\u3092\u504f\u56de\u5e30\u4fc2\u6570\u306e\u70b9\u63a8\u5b9a\u5024\u304b\u3089\u5f15\u304f\u3068 95 \uff05 \u4fe1\u983c\u533a\u9593\u4e0b\u9650\u3001\u8db3\u3059\u3068 95 \uff05 \u4fe1\u983c\u533a\u9593\u4e0a\u9650\u304c\u8a08\u7b97\u3055\u308c\u308b<\/p>\n\n\n\n
\u8a08\u7b97\u7d50\u679c\u306f\u4ee5\u4e0b\u306e\u3068\u304a\u308a<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u4e0a\u8a18\u306e confint \u3092\u4f7f\u3063\u305f\u6642\u3068\u3001\u307b\u3068\u3093\u3069\u540c\u3058\u5024\u306b\u306a\u308b<\/p>\n\n\n\n
confint \u3067\u8a08\u7b97\u3067\u304d\u306a\u3044\u5834\u5408\u306f\u3001\u7c21\u6613\u7684\u306b\u3053\u3061\u3089\u306e\u65b9\u6cd5\u3067\u3088\u3044\u3068\u601d\u308f\u308c\u308b<\/p>\n\n\n\n
\u307e\u3068\u3081<\/h2>\n\n\n\n
\u5bfe\u5fdc\u306e\u3042\u308b\uff12\u5024\u30c7\u30fc\u30bf\u3092\u7fa4\u9593\u6bd4\u8f03\u3057\u305f\u3044\u3068\u304d\u3001\u4e00\u822c\u5316\u7dda\u5f62\u6df7\u5408\u30e2\u30c7\u30eb\u3092\u4f7f\u3046\u65b9\u6cd5\u304c\u3042\u308b\u3002<\/p>\n\n\n\n
\u30c0\u30df\u30fc\u306e\u30b5\u30f3\u30d7\u30eb\u30c7\u30fc\u30bf\u3092\u7528\u3044\u3066\u3001EZR\u3067\u306e\u3084\u308a\u65b9\u3092\u89e3\u8aac\u3057\u305f\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>\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