Hello, I am trying to create an arima forecast model using fpp3 package in R. I am trying to use an ARIMA model, it looks like my data has some season component, but hard to tell. Here are the ACF + PACF visuals of the 3 groups - (A, B,C). I am trying to forecast number of clients in each group for the next 1 year and so, I am using the fpp3 package.

if someone could please help me write my code and help me figure out the best pdq and PDQ values based on my ACF and PACF visuals, that would be much appreciated. Thank you

You could also use auto.arima function if you are not sure. There is not enough data to find good arima models here as we dont know if data is stationary, or has constant variance. You can do hegy test for seasonality. If variance changes over time you would need a garch model to model the variance as well...

You need to take one difference for the model to be stationary. But you need hegy test to see if seasonality exists. I recommend studying time series a little.

I did an nsdiffs instead of ndiffs and got this value below, does this mean that my data is seasonal and I need to take 2 difference for the model? I tried to do the hegy method in r using the forecast package but it did not work for me:

> data %>%
+ features(volume, unitroot_nsdiffs)
# A tibble: 3 × 2
provider_group nsdiffs
<chr> <int>
1 A 2
2 B 2
3 C 2

I have also used the following function in provider group B to see if the sesonal component can be determined:

data %>%
filter(provider_group =="B") %>%
gg_tsdisplay(difference(volume, 12) %>% difference(), plot_type = "partial", lag = 24)