library(astsa)
## Warning: package 'astsa' was built under R version 4.3.3
set.seed(123)
sim_data <- sarima.sim(
n = 200,
ar = numeric(0),
ma = c(0.5, 0.3),
sar = c(0.6, -0.2),
sma = c(-0.4, 0.3),
d = 1,
D = 1,
S = 12
)
ts.plot(sim_data, col = "red", main = "Time Series (With Differencing)", ylab = "Value")

acf(sim_data, main = "ACF - With Differencing")

pacf(sim_data, main = "PACF - With Differencing")

data_nodiff <- cumsum(cumsum(sim_data))
ts.plot(data_nodiff, col = "cyan", main = "Time Series (Without Differencing)", ylab = "Value")

acf(data_nodiff, main = "ACF - Without Differencing")

pacf(data_nodiff, main = "PACF - Without Differencing")
