Forecasting the Appearance Frequency of Rafflesia arnoldii in Bengkulu, Indonesia, Using Discrete-valued Time Series Modeling
DOI:
https://doi.org/10.47654/v30y2026i2p39-67Keywords:
Rafflesia arnoldii, Discrete-valued, INARMA-INGARCH models, Bootstrap AnalysisAbstract
Purpose: This study aims to develop a probabilistic forecasting model to estimate the frequency of Rafflesia arnoldii appearances, a rare flower species endemic to Bengkulu Province, Indonesia. Motivated by the need for effective conservation and ecotourism planning, the study addresses a gap in the literature by applying discrete-valued time series models to ecological count data.
Design/methodology/approach: This study analyzes monthly data on the frequency of Rafflesia arnoldii appearances from January 2016 to December 2023. We applied and compared multiple discrete-valued time series models, including Poisson INAR, Negative Binomial INAR, Zero-Inflated INAR, INARMA, NB-INARMA, and INGARCH. The models were fitted using Maximum Likelihood Estimation and the EM algorithm. Model performance was evaluated using AIC, BIC, RMSE, and MAAPE, with post-estimation diagnostics conducted via Ljung-Box, Jarque-Bera, and Lagrange Multiplier tests.
Findings: The analysis found that the optimal models were INAR(1)-INGARCH(1,1) for South Bengkulu, and INAR(1)-INGARCH(1,1) for Kepahiang. These models effectively captured the discrete, overdispersed, and zero-inflated characteristics of the data, producing reliable short-term forecasts. The variation in model performance across regencies highlights underlying ecological differences and reinforces the need for region-specific forecasting strategies.
Research limitations/implications: Limitations include a focus on four regencies and the exclusion of exogenous environmental variables such as rainfall and temperature. Implications suggest that future studies should expand the data scope, incorporate environmental covariates, and explore spatio-temporal extensions to improve predictive accuracy. In addition, future extensions may explore integer-valued time series frameworks designed to accommodate non-stationary count data, including signed-thinning or generalized INAR formulations, to further enhance model flexibility and predictive performance.
Practical implications: The forecasting framework can inform the scheduling of conservation monitoring and guide ecotourism planning by anticipating peak flowering periods. This study’s primary relevance to Decision Sciences lies in providing a quantitative, data-driven framework for conservation resource allocation and ecotourism planning under uncertainty, transforming ecological data into actionable management strategies.
Originality/value: This study pioneers the application of discrete-valued time series models in forecasting the occurrences of a rare and sporadically flowering plant species. It contributes to the Decision Sciences literature by demonstrating how probabilistic modeling can be applied to biodiversity management under uncertainty.
References
Akaike, H. (1974). A New Look at The Statistical Model Identification. IEEE Transactions on Automatic Control, 19(6), 716 723. https://doi.org/10.1109/TAC.1974.1100705
Al-Osh, M. A., & Alzaid, A. A. (1987). First-Order Integer-Valued Autoregressive (INAR(1)) process. Journal of Time Series Analysis, 8(3), 261 275. https://doi.org/10.1111/j.1467-9892.1987.tb00438.x
Aly, E.-E. A. A., & Nadjib, B. (1994). On Some Integer-Valued Autoregressive Moving Average Models. Journal of Multivariate Analysis, 50(1), 132 151. https://doi.org/10.1006/jmva.1994.1038
Box, G. E. P., Jenkins, G. M., Reinsel, G. C., & Ljung, G. M. (2016). Time Series Analysis Forecasting and Control (5th ed.). John Wiley & Sons, Inc. https://doi.org/10.1111/jtsa.12194
Cameron, A. C., & Trivedi, P. K. (1990). Regression-Based Tests for Overdispersion in The Poisson Model. Journal of Econometrics, 46(3), 347 364. https://doi.org/10.1016/0304-4076(90)90014-K
Chen, C., Davis, R. A., & Brockwell, P. J. (1996). Order Determination for Multivariate Autoregressive Processes Using Resampling Methods. Journal of Multivariate Analysis, 57(2), 175 190. https://doi.org/10.1006/jmva.1996.0028
Collins, N. M., Sayer, J. A., & Whitmore, T. C. (Eds.). (1991). The conservation atlas of tropical forests: Asia and the Pacifics. Palgrave Macmillan London. https://doi.org/10.1007/978-1-349-12030-7
Dungey, M., Martin, V. L., Tang, C., & Tremayne, A. (2019). A Threshold Mixed Count Time Series Model: Estimation and Application. Studies in Nonlinear Dynamics and Econometrics, 24(2), 1 18. https://doi.org/10.1515/snde-2018-0029
Engle, R. F. (1982). Autoregressive Conditional Heteroscedasticity with Estimates of the Variance of United Kingdom Inflation. Econometrica: Journal of the Econometric Society, 987-1007. https://doi.org/10.2307/1912773
Ferland, R., Latour, A., & Oraichi, D. (2006). Integer-Valued GARCH Process. Journal of Time Series Analysis, 27(6), 923-942. https://doi.org/10.1111/j.1467-9892.2006.00496.x
Gujarati, D. N., & Porter, D. C. (2009). Basic Econometrics. McGraw-Hill/Irwin.
Hardin, J. W., & Hilbe, J. M. (2018). Generalized Linear Models and Extensions. Stata press.
Harvey, A. C., & Fernandes, C. (1989). Time Series Models for Count or Qualitative Observations. Journal of Business & Economic Statistics, 7(4), 407-417. https://doi.org/10.2307/1391639
Kim, H. Y., & Park, Y. (2008). A non-stationary integer-valued autoregressive model. Statistical papers, 49(3), 485-502. https://doi.org/10.1007/s00362-006-0028-1
Kim, J. M., & Jun, S. (2017). Integer-Valued GARCH Processes for Apple Technology Analysis. Industrial Management & Data Systems, 117(10), 2381-2399. http://dx.doi.org/10.1108/IMDS-01-2017-0023
Kutoyants, Y. A. (2023). Introduction to the Statistics of Poisson Processes and Applications. Springer Nature. https://doi.org/10.1007/978-3-031-37054-0
Mohammadipour, M., & Boylan, J. E. (2012). Forecast horizon aggregation in integer autoregressive moving average (INARMA) models. Omega, 40(6), 703-712. https://doi.org/10.1016/j.omega.2011.08.008
Nie, H., Liu, G., Liu, X., & Wang, Y. (2012). Hybrid of ARIMA and SVMs for Short-Term Load Forecasting. Energy Procedia, 16, 1455 1460. https://doi.org/10.1016/j.egypro.2012.01.229
Olofsson, P., & Andersson, M. (2012). Probability, Statistics, and Stochastic Processes. John Wiley & Sons. https://doi.org/10.1002/9781118231296
Priatna, D. R., Zuhud, E. A., & Alikodra, H. S. (1989). Kajian Ekologis Rafflesia Patma Blume di Cagar Alam Leuweung Sancang Jawa Barat. Media Konservasi, 2(2). https://dx.doi.org/10.29243/medkon.2.2.%p
Reboredo, J. C., Queiruga, J. R. B., Ferreiro, J. O., & Santias, F. R. (2023). Forecasting Emergency Department Arrivals Using INGARCH Models. Health Economics Review, 13(51). https://doi.org/10.1186/s13561-023-00456-5
Sammaknejad, N., Zhao, Y., & Huang, B. (2019). A review of the expectation maximization algorithm in data-driven process identification. Journal of Process Control, 73, 123 136. https://doi.org/10.1016/j.jprocont.2018.12.010
Schwarz, G. (1978). Estimating The Dimension of a Model. The Annals of Statistics, 461-464. https://doi.org/10.1214/aos/1176344136
Smith, J. Q. (1979). A Generalization of the Bayesian Steady Forecasting Model. Journal of the Royal Statistical Society: Series B (Methodological), 41(3), 375-387. https://doi.org/10.1111/j.2517-6161.1979.tb01092.x
Susatya. (2011). Raflesia Pesona Bunga Terbesar di Dunia. Direktorat Kawasan Konservasi dan Bina Hutan Lindung Dengan Pendanaan dari DIPA TA 2011.
Tijms, H. C. (2003). A First Course in Stochastic Models. John Wiley Sons Ltd. https://doi.org/10.1002/047001363X
Tolod, J. R., Galindon, J. M. M., Atienza, R. R., Duya, M. V., Fernando, E. S., & Ong, P. S. (2021). Flower and fruit development and life history of Rafflesia consueloae (Rafflesiaceae). Philippine Journal of Science, 150(S1), 321 334. https://doi.org/10.56899/150.S1.23
Tsay, R. S. (2010). Analysis of Financial Time Series. John Wiley & Sons. https://doi.org/10.1002/9780470644560
Wei, W. W. S. (2006). Time Series Analysis: Univariate and Multivariate Methods (2nd ed.). Pearson Education Inc.
Weiss, C. H. (2018). An Introduction to Discrete-Valued Time Series. John Wiley & Sons. https://doi.org/10.1002/9781119097013
Weiß, C. H. (2021). Stationary count time series models. Wiley Interdisciplinary Reviews: Computational Statistics, 13(1), e1502. https://doi.org/10.1002/wics.1502
Weiß, C. H., Feld, M. H.-J. M., Mamode Khan, N., & Sunecher, Y. (2019). INARMA Modeling of Count Time Series. Stats, 2(2), 284 320. https://doi.org/10.3390/stats2020022
Weiß, C. H., Homburg, A., & Puig, P. (2016). Testing for Zero Inflation and Overdispersion in INAR(1) Models. Springer. https://doi.org/10.1007/s00362-016-0851-y
Yang, Z., Hardin, J. W., & Addy, C. L. (2009). A Score Test for Overdispersion in Poisson Regression Based on the Generalized Poisson-2 Model. Journal of Statistical Planning and Inference, 139(4), 1514-1521. http://dx.doi.org/10.1016/j.jspi.2008.08.018
Zeger, S. L., & Qaqish, B. (1988). Markov Regression Models for Time Series: A Quasi-Likelihood Approach. Biometrics, 1019-1031. https://doi.org/10.2307/2531732
Zhang, H., Liu, Y., & Li, B. (2014). Notes on Discrete Compound Poisson Model with Applications to Risk Theory. Insurance: Mathematics and Economics, 59, 325-336. https://doi.org/10.1016/j.insmatheco.2014.09.012.
Published
Issue
Section
License
Copyright (c) 2026 Jose Rizal (Corresponding Author); Department of Mathematics, The University of Bengkulu, Bengkulu, Indonesia, Gusman Juliadi, Indah Wahyuliani, Cinta Rizki Oktarina (Author)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
Scientific and Business World