We implement the function fregre.gls (and predict.fregre.gls) that estimates (and predicts) the functional regression model with correlated errors, see S1 Appendix. The fregre.gls function calls the gls function of nlme package. Therefore, the correlation structures allowed are those programmed by the original authors of the package .