I am working with time series values which are all in the closed interval [0, 1]; these values represent relative frequencies, i.e., empirical probabilities. I would like to create a model such that all forecasted values are within [0, 1], but it would also be fine if the model's output was strictly within the open interval (0, 1).
This answered question tackles the lower bound aspect of my question, but not the upper bound aspect: How to achieve strictly positive forecasts?
I'd like to use the R forecast
package if possible to achieve this, but I am open to other suggestions.