Я хочу смоделировать данные временного ряда со средним значением 0, а ковариация cov(ui, uj) = p^2(j-i)^(-a) и a - это константа.
p = 0.5
a = 3
covmat = matrix(NA, nrow = 2, ncol = 2)
for(i in 1: nrow(covmat)){
for(j in 1: ncol(covmat)){
covmat[i,] = p^2*(j-i)^(-a)
}
}
covmat
Почему я получаю Inf в последнем ряду?