# How to statistically distinguish between two types of time series (bimodal vs. not bimodal)?

I have two different types of time series. The first group of series is much more bimodal, and the second group is much flatter. For a given time series, I want to test whether it more likely belongs with the first group or the second group. I don't really care about the specific values of the time series. It is the shape of the time series that matters and how much of a difference there is between the higher part of the series and the lower part of the series. Are there any specific tests for this? I think ARIMA modeling requires a stationary timeseries, which these are not.

Ones that look like this:

And ones that look like this:

• They look different, but bimodality usually refers to modes in a frequency distributions, not peaks in a (seasonal?) time series. You might find it interesting and possibly even helpful to plot on logit scale, i.e. log[reflectance / (1 - reflectance)]. – Nick Cox Mar 5 '18 at 15:12
• Hi Nick, thanks for the help. I suppose I was thinking 'bimodal' if you ignore the time component, and just look at the distribution of the reflectance values, but perhaps that is not correct. Unfortunately, if I plot using logit (with reflectance on the x, and log[reflectance / (1 - reflectance)] on the y), they all look pretty linear unfortunately. Would testing for similarities in the variance of timeseries (like an F test) be advisable? – Ana Mar 5 '18 at 15:32
• The point of a logit scale is that they are more similar in that space, suggesting that distinguishing the two is not necessarily a good idea. What is the underlying physical or biological idea? (I am guessing wildly at data from satellite or drone sensing.) – Nick Cox Mar 5 '18 at 15:34
• I am suggesting logit reflectance versus doy (day of year, presumably). But it is true that logit $y$ does not bend much outside (0.1, 0.9). – Nick Cox Mar 5 '18 at 15:36
• It is studying surface reflectance (from satellite data) over time of ice melting to determine ice breakup dates. It is much easier to see when ice has broken up in the first group of timeseries, but much more difficult to be as confident in the breakup detection in the second group. I'm trying to find a way to flag timeseries that look more like the second set, so that we might be more skeptical of ice breakup dates from those time series if that makes sense. – Ana Mar 5 '18 at 15:39

You can expand the time series on some functional basis. You can do this for example with the fda package in R. Now each time-series is represented as a vector of numbers. You could try a classification with these coefficient as input.