{"id":102,"date":"2023-06-14T21:54:12","date_gmt":"2023-06-14T12:54:12","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/poisson-regression-in-ezr\/"},"modified":"2025-01-18T00:01:10","modified_gmt":"2025-01-17T15:01:10","slug":"poisson-regression-in-ezr","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/poisson-regression-in-ezr\/","title":{"rendered":"EZR \u3067\u30dd\u30a2\u30bd\u30f3\u56de\u5e30\u5206\u6790\u3092\u884c\u3046\u65b9\u6cd5"},"content":{"rendered":"\n

\u30dd\u30a2\u30bd\u30f3\u56de\u5e30\u3092 EZR \u3067\u884c\u3046\u65b9\u6cd5\u306e\u89e3\u8aac<\/p>\n\n\n\n\n\n\n\n

\u30dd\u30a2\u30bd\u30f3\u56de\u5e30\u30fb\u30dd\u30a2\u30bd\u30f3\u5206\u5e03\u3068\u306f<\/h2>\n\n\n\n

\u30dd\u30a2\u30bd\u30f3\u56de\u5e30\u3068\u306f\u3001\u307e\u308c\u306b\u3057\u304b\u8d77\u3053\u3089\u306a\u3044\u73fe\u8c61\u3092\u6570\u3048\u305f\u30ab\u30a6\u30f3\u30c8\u30c7\u30fc\u30bf\u3092\u76ee\u7684\u5909\u6570\u306b\u3057\u305f\u56de\u5e30\u5206\u6790\u306e\u3053\u3068<\/p>\n\n\n\n

\u30ab\u30a6\u30f3\u30c8\u30c7\u30fc\u30bf\u304c\u30dd\u30a2\u30bd\u30f3\u5206\u5e03\u306b\u5f93\u3046\u3068\u4eee\u5b9a\u3057\u3066\u3044\u308b<\/p>\n\n\n\n

\u5e73\u5747 $ \\lambda $ \u56de\u8d77\u3053\u308b\u4e8b\u8c61\u304c\u3001$ X $ \u56de\u8d77\u3053\u308b\u78ba\u7387\u5206\u5e03\u304c\u30dd\u30a2\u30bd\u30f3\u5206\u5e03\u3067\u3042\u308b<\/p>\n\n\n\n

\u305d\u306e\u3068\u304d\u306b\u3001\u3042\u308b\u4e8b\u8c61\u304c $ k $ \u56de\u8d77\u3053\u308b\u78ba\u7387 $ P(X = k) $ \u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u66f8\u3051\u308b<\/p>\n\n\n\n

$$ \\displaystyle P(X = k) = \\frac{e^{-\\lambda} \\lambda^k}{k!} $$<\/p>\n\n\n\n\n\n\n\n

\u4f8b\u3048\u3070\u3001$ \\lambda = 3 $, $ k = 0, 1, 2, \\cdots, 10 $ \u306e\u3068\u304d\u306e\u30b0\u30e9\u30d5\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u306a\u308b<\/p>\n\n\n\n

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

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

\u4e0a\u8a18\u30b0\u30e9\u30d5\u306e R \u30b9\u30af\u30ea\u30d7\u30c8\u306f\u4ee5\u4e0b\u306e\u3068\u304a\u308a<\/p>\n\n\n\n

### Graph of Poisson Distribution<\/span>\nlambda <-<\/span> 3<\/span>\nk <-<\/span> 0<\/span>:<\/span>10<\/span>\nP <-<\/span> exp<\/span>(<\/span>-<\/span>lambda)<\/span>*<\/span>lambda^<\/span>k \/<\/span> factorial<\/span>(<\/span>k)<\/span>\nplot<\/span>(<\/span>k,<\/span> P,<\/span> type=<\/span>'h'<\/span>)<\/span>\nlibrary<\/span>(<\/span>ggplot2)<\/span>\nggplot<\/span>(<\/span>data=<\/span>data.frame<\/span>(<\/span>k,<\/span>P),<\/span> aes<\/span>(<\/span>x=<\/span>k,<\/span> y=<\/span>P))<\/span>+<\/span>\ngeom_point<\/span>(<\/span>size=<\/span>3<\/span>)<\/span>\n<\/code><\/pre>\n\n\n\n
\u30b5\u30f3\u30d7\u30eb\u30c7\u30fc\u30bf\u306e\u8aad\u307f\u8fbc\u307f<\/h2>\n\n\n\n
\u4ee5\u4e0b\u306e\u30b5\u30a4\u30c8\u304b\u3089\u3001data3a.csv \u3092\u30c0\u30a6\u30f3\u30ed\u30fc\u30c9\u3059\u308b<\/p>\n\n\n\n
\u751f\u614b\u5b66\u30c7\u30fc\u30bf\u89e3\u6790 – \u672c\/\u30c7\u30fc\u30bf\u89e3\u6790\u306e\u305f\u3081\u306e\u7d71\u8a08\u30e2\u30c7\u30ea\u30f3\u30b0\u5165\u9580<\/a><\/p>\n\n\n\n
CSV \u3092\u8aad\u307f\u8fbc\u3080\u30e1\u30cb\u30e5\u30fc\u3092\u4f7f\u3063\u3066\u3001\u8aad\u307f\u8fbc\u307f<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
data3a \u3068\u3044\u3046\u30c7\u30fc\u30bf\u30bb\u30c3\u30c8\u540d\u306b\u3057\u3066\u3001\u8aad\u307f\u8fbc\u3080<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u30ab\u30a6\u30f3\u30c8\u30c7\u30fc\u30bf\u306e\u30a2\u30a6\u30c8\u30ab\u30e0 y \u3068\u3001\u9023\u7d9a\u30c7\u30fc\u30bf\u306e x\u3001\u30ab\u30c6\u30b4\u30ea\u30ab\u30eb\u30c7\u30fc\u30bf\u306e f \u306e 3 \u3064\u306e\u30c7\u30fc\u30bf\u542b\u307e\u308c\u308b\u3053\u3068\u304c\u308f\u304b\u308b<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
f \u304c\u3001chr\uff08\u6587\u5b57\u578b\uff09\u3067\u8a8d\u8b58\u3055\u308c\u3066\u3044\u308b\u306e\u3067\u3001\u56e0\u5b50\u578b\u306b\u5909\u66f4\u3057\u3066\u304a\u304f<\/p>\n\n\n\n
\u3059\u3079\u3066\u306e\u6587\u5b57\u5909\u6570\u3092\u56e0\u5b50\u306b\u5909\u63db \u30e1\u30cb\u30e5\u30fc\u3092\u4f7f\u3063\u3066\u3001\u5909\u63db\u3059\u308b<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
Factor \u3068\u3044\u3046\u8868\u8a18\u306b\u5909\u308f\u3063\u305f<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u30ab\u30c6\u30b4\u30ea\u30ab\u30eb\u30c7\u30fc\u30bf\u306f\u3001\u3053\u306e\u56e0\u5b50\u578b\u3067\u306a\u3044\u3068\u3001\u7d71\u8a08\u89e3\u6790\u306b\u4f7f\u3048\u306a\u3044\u306e\u3067\u6ce8\u610f<\/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>
\u30c7\u30fc\u30bf\u306e\u56f3\u793a<\/h2>\n\n\n\n
\u30b0\u30e9\u30d5\u3068\u8868\u306e\u6563\u5e03\u56f3\u30e1\u30cb\u30e5\u30fc\u3092\u4f7f\u3063\u3066\u3001f \u306e\u7fa4\u5225\u306e\u6563\u5e03\u56f3\u3092\u66f8\u3044\u3066\u307f\u308b<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u4eca\u56de\u306f\u3001\u7dda\u5f62\u56de\u5e30\u5206\u6790\u3067\u306f\u306a\u3044\u306e\u3067\u3001\u56de\u5e30\u76f4\u7dda\u306b\u76ee\u3092\u596a\u308f\u308c\u306a\u3044\u305f\u3081\u306b\u3001\u6700\u5c0f 2 \u4e57\u76f4\u7dda\uff08\u56de\u5e30\u76f4\u7dda\uff09\u306e\u30aa\u30d7\u30b7\u30e7\u30f3\u306f\u5916\u3057\u3066\u304a\u304f\u307b\u3046\u304c\u3088\u3044<\/p>\n\n\n\n
\u6563\u5e03\u56f3\u306f\u4ee5\u4e0b\u306e\u3068\u304a\u308a\u3060<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u30dd\u30a2\u30bd\u30f3\u56de\u5e30\u306e\u3042\u3066\u306f\u3081\uff08\u9023\u7d9a\u30c7\u30fc\u30bf\uff09<\/h2>\n\n\n\n
\u3053\u306e\u30c7\u30fc\u30bf\u306b\u30dd\u30a2\u30bd\u30f3\u56de\u5e30\u3092\u3042\u3066\u306f\u3081\u3066\u307f\u308b<\/p>\n\n\n\n
\u307e\u305a\u306f\u3001x \u3067 y \u3092\u4e88\u6e2c\uff08\u8aac\u660e\uff09\u3067\u304d\u308b\u304b\u3068\u3044\u3046\u8a08\u7b97<\/p>\n\n\n\n
\u6a19\u6e96\u30e1\u30cb\u30e5\u30fc \u2192 \u7d71\u8a08\u91cf \u2192 \u30e2\u30c7\u30eb\u3078\u306e\u9069\u5408 \u2192 \u4e00\u822c\u5316\u7dda\u578b\u30e2\u30c7\u30eb \u3092\u9078\u629e<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u76ee\u7684\u5909\u6570\u306b y\u3001\u8aac\u660e\u5909\u6570\u306b x\u3001\u30ea\u30f3\u30af\u95a2\u6570\u65cf\u306f\u3001poisson\uff08\u30dd\u30a2\u30bd\u30f3\uff09\u3092\u30c0\u30d6\u30eb\u30af\u30ea\u30c3\u30af\u3057\u3066\u9078\u629e\u3001\u30ea\u30f3\u30af\u95a2\u6570\u306f\u81ea\u52d5\u3067 log \u304c\u9078\u3070\u308c\u305f\u72b6\u614b\u306b\u306a\u308b\u306e\u3067\u305d\u306e\u307e\u307e<\/p>\n\n\n\n
\u3053\u306e\u3088\u3046\u306b\u8a2d\u5b9a\u3059\u308b<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
OK \u3092\u30af\u30ea\u30c3\u30af\u3059\u308b\u3068\u8a08\u7b97\u3055\u308c\u308b<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
x \u304c 1 \u4e0a\u6607\u3059\u308b\u3054\u3068\u306b\u3001\u5bfe\u6570\u30ab\u30a6\u30f3\u30c8\u304c\u30010.07566 \u4e0a\u6607\u3059\u308b\u3068\u3044\u3046\u8a08\u7b97\u7d50\u679c\u3067\u3001\u7d71\u8a08\u5b66\u7684\u6709\u610f\uff08\u691c\u5b9a\u306e\u540d\u524d\u306f Wald \u691c\u5b9a\uff09<\/p>\n\n\n\n
\u5bfe\u6570\u30ab\u30a6\u30f3\u30c8\u4e0a\u6607\u3092\u771f\u6570\u30ab\u30a6\u30f3\u30c8\u306b\u76f4\u3059\u3068\u3001x \u304c 1 \u4e0a\u6607\u3059\u308b\u3054\u3068\u306b\u7d04 1.08 \u500d\u306e\u3079\u304d\u4e57\uff08$ e^{1.08 x} $\uff09\u306b\u306a\u308b\u3068\u3044\u3046\u8a08\u7b97\u306b\u306a\u308b<\/p>\n\n\n\n
\u30dd\u30a2\u30bd\u30f3\u56de\u5e30\u306e\u3042\u3066\u306f\u3081\uff08\u30ab\u30c6\u30b4\u30ea\u30ab\u30eb\u30c7\u30fc\u30bf\uff09<\/h2>\n\n\n\n
\u5909\u6570 f \uff08\u30ab\u30c6\u30b4\u30ea\u30ab\u30eb\u30c7\u30fc\u30bf\uff09\u3092\u8aac\u660e\u5909\u6570\u306b\u3059\u308b\u3068\u3069\u3046\u306a\u308b\u304b<\/p>\n\n\n\n
\u30e1\u30cb\u30e5\u30fc\u5185\u306e\u8a2d\u5b9a\u306f\u4ee5\u4e0b\u306e\u3068\u304a\u308a<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
\u8aac\u660e\u5909\u6570\u306e\u90e8\u5206\u3060\u3051\u3001x \u304b\u3089 f \u306b\u5909\u3048\u305f<\/p>\n\n\n\n
\u7d50\u679c\u306f\u3001\u4ee5\u4e0b\u306e\u3068\u304a\u308a<\/p>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
f \u306e C \u3068\u6bd4\u8f03\u3057\u3066\u3001T \u306e\u4fc2\u6570\u63a8\u5b9a\u5024 0.01277 \u304c\u8a08\u7b97\u3055\u308c\u3066\u3044\u308b<\/p>\n\n\n\n
\u7d71\u8a08\u5b66\u7684\u6709\u610f\u3067\u306f\u306a\u3044<\/p>\n\n\n\n
\u30dd\u30a2\u30bd\u30f3\u56de\u5e30\u306e\u3042\u3066\u306f\u3081\uff08\u9023\u7d9a\u30c7\u30fc\u30bf\uff0b\u30ab\u30c6\u30b4\u30ea\u30ab\u30eb\u30c7\u30fc\u30bf\uff09<\/h2>\n\n\n\n
\u3067\u306f\u3001\u5909\u6570 x \u3068 f \u3092\u3068\u3082\u306b\u8aac\u660e\u5909\u6570\u3068\u3057\u305f\u3089\u3069\u3046\u3060\u308d\u3046\u304b<\/p>\n\n\n\n
f \u306e C \u3068 T \u306e\u5dee\u3092\u898b\u305f\u3044\u3068\u304d\u306b\u3001x \u304c\u4ea4\u7d61\u56e0\u5b50\u3067\u3001\u8abf\u6574\u3057\u306a\u3044\u3068\u3044\u3051\u306a\u3044\u3068\u3044\u3046\u72b6\u6cc1\u3068\u8003\u3048\u3066\u8a08\u7b97\u3057\u3066\u307f\u308b<\/p>\n\n\n\n
\u3061\u306a\u307f\u306b\u3001C \u3068 T \u3067\u306f\u3001x \u306e\u5024\u304c\u7570\u306a\u3063\u3066\u304a\u308a\u3001\u4ea4\u7d61\u56e0\u5b50<\/p>\n\n\n\n