Мне сложно определить стационарное распределение марковской модели. Я начинаю понимать теорию и связи: Учитывая стохастическую матрицу, для dermine стационарного распределения мы должны найти собственный вектор для наибольшего собственного значения (то есть 1)Вычисление собственных значений/собственных векторов стохастической матрицы
Я начал с созданием стохастической матрицы
set.seed(6534)
stoma <- matrix(abs(rnorm(25)), nrow=5, ncol=5)
stoma <- (stoma)/rowSums(stoma) # that should make it a stochastic matrix rowSums(stoma) == 1
Затем я использую R eigen
функцию
ew <- eigen(stoma)
Но я не понимаю, результаты
> ew
$values
[1] 1.000000e+00+0.000000e+00i -6.038961e-02+0.000000e+00i -3.991160e-17+0.000000e+00i
[4] -1.900754e-17+1.345763e-17i -1.900754e-17-1.345763e-17i
$vectors
[,1] [,2] [,3] [,4] [,5]
[1,] -0.4472136+0i 0.81018968+0i 0.3647755+0i -0.0112889+0.1658253i -0.0112889-0.1658253i
[2,] -0.4472136+0i 0.45927081+0i -0.7687393+0i 0.5314923-0.1790588i 0.5314923+0.1790588i
[3,] -0.4472136+0i 0.16233945+0i 0.2128250+0i -0.7093859+0.0000000i -0.7093859+0.0000000i
[4,] -0.4472136+0i -0.09217315+0i 0.4214660+0i -0.1305497-0.1261247i -0.1305497+0.1261247i
[5,] -0.4472136+0i -0.31275073+0i -0.2303272+0i 0.3197321+0.1393583i 0.3197321-0.1393583i
Вектор для наибольшего значения (1) имеет все те же значения компонента «-0.4472136». Даже если я изменю семя, чтобы нарисовать различное число, я снова получаю те же значения. Что мне не хватает? Почему компоненты собственного вектора все эквалы? Почему они не суммируются до 1 - так как это должно быть стационарное распределение?
Благодарим за помощь!
если размер корреляционной матрицы меньше, то количество периодов, то матрица вырождена – Qbik