I'm having a bit of trouble with my state space model fit (using KFAS package) when I include a seasonal component in the model. The computations don't seem to end when I include (what I believe is) the frequency of the seasonality in the seasonal component of the model.
So heres what my data Looks like.. The response variable is count data, basically the number of animals at a site. I have 10 of these sites, and each month from 1997-2014 I make an observation at each site and record the number of animals. So 10 observations in total for 10 sites in 1 month. So in total I have 2160 data points.
Now I want to capture the 'seasonality' in the data, because I know that I'll have higher counts during summer months than I will during the winter months. So in 1 year, I've taken 120 observations.. what i've been doing up until this post is setting the frequency of my seasonality as 120. To me it seems correct, but the computation seems to be very intensive so I have a feeling I may be wrong.
Any suggestions?
I also had another idea in mind, but im not sure if its a valid approach. Since my response is count data, it follows a poisson distribution (non-gaussian), if I were to decompose the data into trend,error and seasonal components, could I be able to model just the trend+error and assume its under a gaussian distribution? (thus completely removing the seasonality from my analysis)
Thanks for your help.