Привет всем, я сейчас работаю с программированием на R, чтобы найти среднее значение цепей Маркова, чтобы я мог определить дисперсию между цепями.
Проблема объясняется как «Выполнить J цепочек длины L из различных начальных значений и отбросить первый D.»Затем мне дали следующую информацию:
Chain 1= x_1^(0),...,x_1^(D),x_1^(D+1),...,x_1^(D+L-1)
....
Chain J= x_J^(0),...,x_J^(D),x_J^(D+1),...,x_J^(D+L-1)
Среднее значение j-й цепочки показано как:
СреднееСредство показано как:
Меня просят найти межцепную дисперсию:
Итак, вот моя попытка его кодирования:
L=5
D=50
chain1=c(1,2,3,4,6)
chain2=c(2,3,4,5,7)
t=D
for(t in 1:(D+L-1)){
xbar1=1/L*sum(chain1)
xbar2=1/L*sum(chain2)
}
xbarj=xbar1+xbar2
J=2
for(j in 1:J){
overallmean=1/J*sum(xbarj)
}
for(j in 1:J){
B=L/(J-1)*sum(xbarj-overallmean)^2
}
Я получаю ответ, но я не уверен, что это правильный способ показать его.Мой код написан на языке R, но вы можете показать его на любом языке программирования, если это более полезно для обучения.