library(stats)
set.seed(1)
ts.sim5 <- arima.sim(list(order = c(0,0,2), ma = c(0.5,0.4)), n = 1000)
pacf(ts.sim5, lag.max=10)
acf(ts.sim5, lag.max=10)