Hi all,
I would like to predict prediction interval for a SAR model (function errorsarlm in R - package spdep). While the function predict.lm allows to set interval='prediction' parameter to predict the upper and lower boundary of predictions, the function predict.sarlm does not. Is it because it is incorrect to estimate the prediction interval of a SAR or simply because it hasn't implemented? Is it possible to calculate those prediction intervals manually? Below a replicable example:
library(spdep)
data(oldcol)
lw