Displaying 1 result from an estimated 1 matches for "synt_seri".
Did you mean:
synt_series
2007 Jul 16
1
question about ar1 time series
...hite-noise" is a "Gaussian-noise."
#rho1 (or alpha in another text-books ;-)) < 1 (in fact 0 < rho1 < 1) so that
#the system can be stationary.
#Where var_serie is the variance of the serie
cat("\n Hello, this is creat_AR1_synt_ser.R. \n These are the input
parameters: synt_series(ar1_length, rho1, ...), where rho1 is the
correlat. coef.\n")
ar1 <- function(x, rho1, af)
{
return(x*rho1 + runif(1, -af, af))
}
#Spin-up for the AR1 series. For this case is enough with this amount
spinup <- function(x0, rho1, af)
{
xt <- x0
for (i in 1:100) {
xtp...