We use turbidity to estimate suspended-sediment concentration (SSC)- our data was serially correlated. We ran an ARMA process and ended up with a AR (2) model. Our equation in log form is:
estlogi(SSC) = 1.35149 + 0.785099logi(Turbidity) - 0.00095 + 0.7492ei-1+ 0.2052ei-2
How can we input the ei-1 & ei-2 values in real time, that is, how do we enter residuals for current or future values when our turbidity/SSC data pairs are in the past (2010 to 2011). Someone did suggest an R script that will estimate K steps ahead for a specified ARMA model and past data- does anyone know this script?