Я использую приведенный ниже код, который является реализацией прямого алгоритма. Но это приводит к ошибке:
Ошибка в альфа [t - 1,]: индекс за пределами границ
forward = function(v, a, b, initial_distribution){
T = length(v)
M = nrow(a)
alpha = matrix(0, T, M)
alpha[1, ] = initial_distribution*b[, v[1]]
for(t in 2:T){
tmp = alpha[t-1, ] %*% a
alpha[t, ] = tmp * b[, v[t]]
}
return(alpha)
}
Как я могу ее решить?