2017-02-08 8 views
0

私はカウントデータを持っており、データには自己相関と過分散の問題があるため、動的負の二項回帰を使用して時系列解析を行う必要があります。 私は使用できるRパッケージをオンラインで検索しましたが、見つけられませんでした。時系列の動的負の二項回帰

私は助けていただきありがとうございます。

私のデータの例:

>St1 
[1] 17 9 28 7 23 16 17 12 11 16 19 29 5 40 13 27 13 11 10 14 13 23 21 24 9 42 14 22 17 9 

>Years 
[1] 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 
[23] 2007 2008 2009 2010 2011 2012 2013 2014 

>library(AER) 
>library(stats) 

>rd <- glm(St1 ~ Years, family = poisson) 
>dispersiontest(rd) 

Overdispersion test 
data: rd 
z = 2.6479, p-value = 0.00405 
alternative hypothesis: true dispersion is greater than 1 
sample estimates: 
dispersion 
    4.305539 

#Autocorrelation 
>Box.test (St1, lag=ceiling(log(length(St1))), type = "Ljung") 

    Box-Ljung test 

data: St1 
X-squared = 13.612, df = 4, p-value = 0.008641 
+1

おそらく、「動的二項回帰」を示すメソドロジ指向のWebページへのリンクを用意する必要があります。作品 "ダイナミック"は、本当に_ _ _ _ _少ない_> _ –

+0

ありがとう、多くの意味を取ることができます。私の理解は、動的な負の二項回帰は、自己相関と過分散を説明することができるということです。こちらをご覧くださいhttp://journals.sagepub.com/doi/pdf/10.1177/1471082X14535530 – Ali

答えて

1

は、これは基本的にパッケージを見つけるための要求である(そして、そのような要求は、オフトピックと考えられています)。だから私はそれをコーディングフレーバを持つ質問に変換できるかどうかを見極めるつもりです。私が私のコメントで言ったように、検索用語として「ダイナミック」を使用しようとすると、しばしば失望してしまいます。なぜなら、誰もが断絶された目的のためにこの言葉を使いたいからです。コンソールからの検索で表示される機能を確認してください。

install.packages("sos") 
sos::findFn(" dynamic negative binomial") 
found 20 matches 
Downloaded 20 links in 13 packages. 

役に立つと思われるものはありません。しかし、あなたの引用を見て、すべてのモデルがautoreggressionコンポーネントを持っていたと思われたので、この検索は....

sos::findFn(" autoregressive negative binomial") 
found 28 matches; retrieving 2 pages 
2 
Downloaded 27 links in 16 packages. 

は、検索:"Fitting Longitudinal Data with Negative Binomial Marginal...""Generalized Linear Autoregressive Moving Average Models with..."。だから、これを "暗黙の質問"への私の答えと考えてください。sos -packageでRコンソールから効果的な検索を行うには?

関連する問題