私は、ゼロインフレーションポアソン(zip)とゼロインフレーション負の二項(zinb)回帰を使用して、カウントデータの時間的傾向を検出しています6つの病院で報告されています)、ゼロと過分散の可能性があります。 私はpscl packageを使用していくつかのコードを書いており、私の目標は病院間の傾向を比較することです。トレンドのゼロインフレーション回帰とゼロインフレーション負の二項回帰の使用
ジップとzinbショー負の動向に対しCounts<- read.csv("data.csv", header = T)
Years= Counts$X
Ho1= Counts$Ho1
Ho2= Counts$Ho2
Ho3= Counts$Ho3
... .........
... ..........
require(pscl)
zip1 <- zeroinfl(Ho1 ~ Years, dist = "poisson")
zinb4 <- zeroinfl(Ho4 ~ Years, dist = "negbin")
しかし、私はいくつかのデータをプロットするとき、それはやや増加傾向を示して
ここでは一例です:
ジップ結果:
zip1
Call:
zeroinfl(formula = Ho1 ~ Years, dist = "poisson")
Count model coefficients (poisson with log link):
(Intercept) Years
-4.836815 0.002837
Zero-inflation model coefficients (binomial with logit link):
(Intercept) Years
467.2323 -0.2353
このモデルの傾向(勾配)は-0.235であり、普通の最小二乗(OLS)を使用したときの傾向= 0.043です。 私の理解では、zipとOLSの両方が少し異なるはずです。
私のコードが正しくないかもしれないと思っていたのですが、何か不足しています。
は、私はあなたがカウントを増やしますYears
が増加すると任意の考えや提案
私が心から援助を受けることに興味があれば、再現可能な例を提供することはできません... –