Thanks for your comment. Iam going to provide the reprex right down here. Sorry for my mistake.
t <- seq(0.01, 10, 0.01)
Z = rnorm(length(t), 0, 1)
#coefficients
#intercept
a <- -0.005314
#AR-coefficients
phi.1 <- 0.95038
phi.2 <- -0.25543
#MA-coefficients
theta.1 <- 0.04334
theta.2 <- 0.08012
theta.3 <- -0.07681
theta.4 <- -0.01477
#ts construction
Y <- vector()
Y[1:4] = Z[1:4]
for (i in 5:length(t)) {
Y[i] = a + phi.1 * Y[i-1] + phi.2 * Y[i-2] + Z[i] + theta.4 * Z[i-4] + theta.3 * Z[i-3] +theta.2 * Z[i-2] + theta.1 * Z[i-1]
}
acf(Y)

pacf(Y)
The given code reproduces a time series with acf cutting at lag=4 and pacf cutting at lag=2.
And there is no way to get arround using a function resp. a list to summarize the sums? I was hoping to just use a nested loop to be able to reduce the code and also summarize the variables in the form of
phi <- c( 0.95038 , -0.25543 )
theta <- c(0.04334, 0.08012, -0.07681, -0.01477)
and then adress the components of the parameter-vectors one by the other.
Thanks for your help and kind regards
FactOREO