# Create a moving average process using gaussian white noise
# Lowpass filter
# L=length of the series

sim.ma <- function(a=c(1/3,1/3,1/3),L=200)
{
 x <- as.rts(rnorm(L))
 y <- filter(x,a)

 ts.plot(x,y)
}