{"id":559,"date":"2018-06-29T21:56:46","date_gmt":"2018-06-29T12:56:46","guid":{"rendered":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/negative-binomial-model-by-r\/"},"modified":"2024-10-14T11:46:34","modified_gmt":"2024-10-14T02:46:34","slug":"negative-binomial-model-by-r","status":"publish","type":"post","link":"https:\/\/best-biostatistics.com\/toukei-er\/entry\/negative-binomial-model-by-r\/","title":{"rendered":"R \u3067\u8ca0\u306e\u4e8c\u9805\u56de\u5e30\u30e2\u30c7\u30eb\u3092\u5b9f\u884c\u3059\u308b\u65b9\u6cd5"},"content":{"rendered":"\n

\u8ca0\u306e\u4e8c\u9805\u56de\u5e30\u30e2\u30c7\u30eb\u3092 R \u3067\u884c\u3046\u65b9\u6cd5<\/p>\n\n\n\n\n\n\n\n

\u8ca0\u306e\u4e8c\u9805\u30e2\u30c7\u30eb\u3092 R \u3067\u884c\u3046\u65b9\u6cd5<\/h2>\n\n\n\n

\u307e\u305a\u3001\u30c7\u30fc\u30bf\u3092\u8aad\u307f\u8fbc\u3080<\/p>\n\n\n\n

X <-<\/span> c<\/span>(<\/span>rep<\/span>(<\/span>1<\/span>,<\/span> 913<\/span>),<\/span> rep<\/span>(<\/span>2<\/span>,<\/span> 1345<\/span>),<\/span> rep<\/span>(<\/span>3<\/span>,<\/span> 519<\/span>),<\/span> rep<\/span>(<\/span>9<\/span>,<\/span> 863<\/span>))<\/span>\nY <-<\/span> c<\/span>(<\/span>rep<\/span>(<\/span>0<\/span>,<\/span> 910<\/span>),<\/span> rep<\/span>(<\/span>1<\/span>,<\/span> 2<\/span>),<\/span> rep<\/span>(<\/span>2<\/span>,<\/span> 1<\/span>),<\/span> rep<\/span>(<\/span>4<\/span>,<\/span> 0<\/span>),<\/span>\nrep<\/span>(<\/span>0<\/span>,<\/span> 1333<\/span>),<\/span> rep<\/span>(<\/span>1<\/span>,<\/span> 11<\/span>),<\/span> rep<\/span>(<\/span>2<\/span>,<\/span> 1<\/span>),<\/span> rep<\/span>(<\/span>4<\/span>,<\/span> 0<\/span>),<\/span>\nrep<\/span>(<\/span>0<\/span>,<\/span> 509<\/span>),<\/span> rep<\/span>(<\/span>1<\/span>,<\/span> 9<\/span>),<\/span> rep<\/span>(<\/span>2<\/span>,<\/span> 1<\/span>),<\/span> rep<\/span>(<\/span>4<\/span>,<\/span> 0<\/span>),<\/span>\nrep<\/span>(<\/span>0<\/span>,<\/span> 849<\/span>),<\/span> rep<\/span>(<\/span>1<\/span>,<\/span> 13<\/span>),<\/span> rep<\/span>(<\/span>2<\/span>,<\/span> 0<\/span>),<\/span> rep<\/span>(<\/span>4<\/span>,<\/span> 1<\/span>))<\/span>\ndat <-<\/span> data.frame<\/span>(<\/span>X,<\/span> Y)<\/span>\n<\/code><\/pre>\n\n\n\n
Y\u306e\u5206\u5e03\u3092\u78ba\u8a8d<\/p>\n\n\n\n
><\/span> table<\/span>(<\/span>dat$<\/span>Y)<\/span>\n   0<\/span>    1<\/span>    2<\/span>    4<\/span>\n3601<\/span>   35<\/span>    3<\/span>    1<\/span>\n<\/code><\/pre>\n\n\n\n
Y\u306e\u5206\u5e03\u56f3\u306f\u3059\u3054\u3044\u3053\u3068\u306b\u3001\u306a\u3063\u3066\u3044\u308b<\/p>\n\n\n\n
plot<\/span>(<\/span>table<\/span>(<\/span>dat$<\/span>Y))<\/span>\n<\/code><\/pre>\n\n\n\n
\"\"<\/figure>\n\n\n\n
<\/span><\/p>\n\n\n\n
X\u3068Y\u306e\u95a2\u4fc2\u3092\u898b\u3066\u3044\u3053\u3046<\/p>\n\n\n\n
\u307e\u305a\u30af\u30ed\u30b9\u96c6\u8a08\u8868<\/p>\n\n\n\n
><\/span> table<\/span>(<\/span>dat$<\/span>X,<\/span> dat$<\/span>Y)<\/span>\n <\/span>    0<\/span>    1<\/span>    2<\/span>    4<\/span>\n1<\/span>  910<\/span>    2<\/span>    1<\/span>    0<\/span>\n2<\/span> 1333<\/span>   11<\/span>    1<\/span>    0<\/span>\n3<\/span>  509<\/span>    9<\/span>    1<\/span>    0<\/span>\n9<\/span>  849<\/span>   13<\/span>    0<\/span>    1<\/span>\n<\/code><\/pre>\n\n\n\n
Y \u304c\u30bc\u30ed\u3070\u304b\u308a\u306a\u306e\u304c\u76ee\u7acb\u3064<\/p>\n\n\n\n
\u8ca0\u306e\u4e8c\u9805\u56de\u5e30\u30e2\u30c7\u30eb\u306f\u3001MASS \u30d1\u30c3\u30b1\u30fc\u30b8\u306e glm.nb() \u3092\u4f7f\u3046<\/p>\n\n\n\n
library<\/span>(<\/span>MASS)<\/span>\nnb.res <-<\/span> glm.nb<\/span>(<\/span>Y ~<\/span> factor<\/span>(<\/span>X),<\/span> data=<\/span>dat)<\/span>\nsummary<\/span>(<\/span>nb.res)<\/span>\n<\/code><\/pre>\n\n\n\n
\u7d50\u679c\u306f\u3001X\u304c\uff11\u306b\u6bd4\u3079\u3066\u3001\uff13\u3068\u304b\u3001\uff19\u3068\u304b\u306fY\u306b\u95a2\u9023\u304c\u3042\u308b<\/p>\n\n\n\n
><\/span> summary<\/span>(<\/span>nb.res)<\/span>\nCall:<\/span>\nglm.nb<\/span>(<\/span>formula =<\/span> Y ~<\/span> factor<\/span>(<\/span>X),<\/span> data =<\/span> dat,<\/span> init.theta =<\/span> 0.05009536409<\/span>,<\/span>\nlink =<\/span> log)<\/span>\nDeviance Residuals:<\/span>\nMin        1Q    Median        3Q       Max\n-<\/span>0.18802<\/span>  -<\/span>0.18228<\/span>  -<\/span>0.13295<\/span>  -<\/span>0.09165<\/span>   3.10151<\/span>\nCoefficients:<\/span>\nEstimate Std. Error z value Pr<\/span>(<\/span>>|<\/span>z|<\/span>)<\/span>\n(<\/span>Intercept)<\/span>  -<\/span>5.4304<\/span>     0.5214<\/span> -<\/span>10.415<\/span>   <<\/span>2e-16<\/span> ***<\/span>\nfactor<\/span>(<\/span>X)<\/span>2<\/span>    0.7912<\/span>     0.6030<\/span>   1.312<\/span>   0.1895<\/span>\nfactor<\/span>(<\/span>X)<\/span>3<\/span>    1.5764<\/span>     0.6334<\/span>   2.489<\/span>   0.0128<\/span> *<\/span>\nfactor<\/span>(<\/span>X)<\/span>9<\/span>    1.5032<\/span>     0.5948<\/span>   2.527<\/span>   0.0115<\/span> *<\/span>\n---<\/span>\nSignif. codes:<\/span>  0<\/span> \u2018***<\/span>\u2019 0.001<\/span> \u2018**<\/span>\u2019 0.01<\/span> \u2018*<\/span>\u2019 0.05<\/span> \u2018.\u2019 0.1<\/span> \u2018 \u2019 1<\/span>\n(<\/span>Dispersion parameter for<\/span> Negative Binomial<\/span>(<\/span>0.0501<\/span>)<\/span> family taken to be 1<\/span>)<\/span>\nNull deviance:<\/span> 210.18<\/span>  on 3639<\/span>  degrees of freedom\nResidual deviance:<\/span> 199.65<\/span>  on 3636<\/span>  degrees of freedom\nAIC:<\/span> 465.67<\/span>\nNumber of Fisher Scoring iterations:<\/span> 1<\/span>\nTheta:<\/span>  0.0501<\/span>\nStd. Err.:<\/span>  0.0236<\/span>\n2<\/span> x log-<\/span>likelihood:<\/span>  -<\/span>455.6660<\/span>\n<\/code><\/pre>\n\n\n\n
\u307e\u3068\u3081<\/h2>\n\n\n\n
\u8ca0\u306e\u4e8c\u9805\u56de\u5e30\u30e2\u30c7\u30eb\u3092 R \u3067\u5b9f\u884c\u3059\u308b\u65b9\u6cd5\u3092\u7c21\u5358\u306b\u7d39\u4ecb\u3057\u305f<\/p>\n\n\n\n
MASS \u30d1\u30c3\u30b1\u30fc\u30b8\u306e glm.nb() \u95a2\u6570\u3092\u4f7f\u3046\u3068\u5b9f\u884c\u3067\u304d\u308b<\/p>\n","protected":false},"excerpt":{"rendered":"
\u8ca0\u306e\u4e8c\u9805\u56de\u5e30\u30e2\u30c7\u30eb\u3092 R \u3067\u884c\u3046\u65b9\u6cd5<\/p>\n","protected":false},"author":2,"featured_media":3061,"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":[5,87],"tags":[],"class_list":["post-559","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-r","category-87"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"https:\/\/best-biostatistics.com\/toukei-er\/wp-content\/uploads\/2018\/06\/20230107163259.png","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/559","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=559"}],"version-history":[{"count":2,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/559\/revisions"}],"predecessor-version":[{"id":3063,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/posts\/559\/revisions\/3063"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/media\/3061"}],"wp:attachment":[{"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/media?parent=559"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/categories?post=559"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/best-biostatistics.com\/toukei-er\/wp-json\/wp\/v2\/tags?post=559"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}