Derivation of a stationary distribution of a Markov chain

The name of the pictureThe name of the pictureThe name of the pictureClash Royale CLAN TAG#URR8PPP











up vote
0
down vote

favorite












Consider the following markov chain in index notation



$$P_j(t+1)=sum_i=1^nM_ijP_i(t)$$



where $M_ij$ is a column transition matrix .



I would like to develop the steady state given by $pi M=pi $ in an index notation for the Column transition matrix.



This is what I have so far



$$P_j(t+1)=sum_i=1^n M_ijP_i(t)=sum_i=1^nleft(M_ijP_i(t)+M_iiP_i(t)right)=$$



$$=sum_i=1^nM_ijP_i(t)+sum_i=1^nM_iiP_i(t)=sum_i=1,ineq j^nM_ijP_i(t)+left(1-sum_j=1,ineq j^nM_jiright)P_j(t)=$$



$$=sum_i=1,ineq j^nM_ijP_i(t)+P_j(t)-sum_j=1,ineq j^nM_jiP_j(t)=$$



$$longrightarrow$$



$$P_j(t+1)-P_j(t)=sum_i=1,ineq j^nM_ijP_i(t)-sum_j=1,ineq j^nM_jiP_j(t)$$
$$=sum_i=1,ineq j^nM_ijP_i(t)-sum_j=1,ineq j^nM_jiP_j(t)+sum_i=1^nM_ii P_i(t)-sum_j=1^nM_jj P_j(t)$$



$$=sum_i=1^nM_ijP_i(t)-sum_j=1^nM_jiP_j(t)$$



At steady state



$$sum_i=1^nM_ijP_i(t)=sum_j=1^nM_jiP_j(t)$$



Which is a trivial solution, I don't understand how to obtain the steady state solution $pi=Mpi$ in matrix form the equation above.



Please advise, note that it is important for me, if possible, to remain in mattrix form for a column transition matrix.







share|cite|improve this question





















  • I'm not really sure what you mean by "develop the steady state...in an index notation." Note that the $n$-dimensional vector whose entries are identically $1$ is always a solution to $pi=Mpi$.
    – Math1000
    Jul 25 at 11:56










  • $mathbf1$ will always work as pointed out above. More generally $pi$ can be any eigenvector of $M$ with eigenvalue 1. Basically you just need to find the nullspace of $M-I$
    – Casey
    Jul 25 at 20:14














up vote
0
down vote

favorite












Consider the following markov chain in index notation



$$P_j(t+1)=sum_i=1^nM_ijP_i(t)$$



where $M_ij$ is a column transition matrix .



I would like to develop the steady state given by $pi M=pi $ in an index notation for the Column transition matrix.



This is what I have so far



$$P_j(t+1)=sum_i=1^n M_ijP_i(t)=sum_i=1^nleft(M_ijP_i(t)+M_iiP_i(t)right)=$$



$$=sum_i=1^nM_ijP_i(t)+sum_i=1^nM_iiP_i(t)=sum_i=1,ineq j^nM_ijP_i(t)+left(1-sum_j=1,ineq j^nM_jiright)P_j(t)=$$



$$=sum_i=1,ineq j^nM_ijP_i(t)+P_j(t)-sum_j=1,ineq j^nM_jiP_j(t)=$$



$$longrightarrow$$



$$P_j(t+1)-P_j(t)=sum_i=1,ineq j^nM_ijP_i(t)-sum_j=1,ineq j^nM_jiP_j(t)$$
$$=sum_i=1,ineq j^nM_ijP_i(t)-sum_j=1,ineq j^nM_jiP_j(t)+sum_i=1^nM_ii P_i(t)-sum_j=1^nM_jj P_j(t)$$



$$=sum_i=1^nM_ijP_i(t)-sum_j=1^nM_jiP_j(t)$$



At steady state



$$sum_i=1^nM_ijP_i(t)=sum_j=1^nM_jiP_j(t)$$



Which is a trivial solution, I don't understand how to obtain the steady state solution $pi=Mpi$ in matrix form the equation above.



Please advise, note that it is important for me, if possible, to remain in mattrix form for a column transition matrix.







share|cite|improve this question





















  • I'm not really sure what you mean by "develop the steady state...in an index notation." Note that the $n$-dimensional vector whose entries are identically $1$ is always a solution to $pi=Mpi$.
    – Math1000
    Jul 25 at 11:56










  • $mathbf1$ will always work as pointed out above. More generally $pi$ can be any eigenvector of $M$ with eigenvalue 1. Basically you just need to find the nullspace of $M-I$
    – Casey
    Jul 25 at 20:14












up vote
0
down vote

favorite









up vote
0
down vote

favorite











Consider the following markov chain in index notation



$$P_j(t+1)=sum_i=1^nM_ijP_i(t)$$



where $M_ij$ is a column transition matrix .



I would like to develop the steady state given by $pi M=pi $ in an index notation for the Column transition matrix.



This is what I have so far



$$P_j(t+1)=sum_i=1^n M_ijP_i(t)=sum_i=1^nleft(M_ijP_i(t)+M_iiP_i(t)right)=$$



$$=sum_i=1^nM_ijP_i(t)+sum_i=1^nM_iiP_i(t)=sum_i=1,ineq j^nM_ijP_i(t)+left(1-sum_j=1,ineq j^nM_jiright)P_j(t)=$$



$$=sum_i=1,ineq j^nM_ijP_i(t)+P_j(t)-sum_j=1,ineq j^nM_jiP_j(t)=$$



$$longrightarrow$$



$$P_j(t+1)-P_j(t)=sum_i=1,ineq j^nM_ijP_i(t)-sum_j=1,ineq j^nM_jiP_j(t)$$
$$=sum_i=1,ineq j^nM_ijP_i(t)-sum_j=1,ineq j^nM_jiP_j(t)+sum_i=1^nM_ii P_i(t)-sum_j=1^nM_jj P_j(t)$$



$$=sum_i=1^nM_ijP_i(t)-sum_j=1^nM_jiP_j(t)$$



At steady state



$$sum_i=1^nM_ijP_i(t)=sum_j=1^nM_jiP_j(t)$$



Which is a trivial solution, I don't understand how to obtain the steady state solution $pi=Mpi$ in matrix form the equation above.



Please advise, note that it is important for me, if possible, to remain in mattrix form for a column transition matrix.







share|cite|improve this question













Consider the following markov chain in index notation



$$P_j(t+1)=sum_i=1^nM_ijP_i(t)$$



where $M_ij$ is a column transition matrix .



I would like to develop the steady state given by $pi M=pi $ in an index notation for the Column transition matrix.



This is what I have so far



$$P_j(t+1)=sum_i=1^n M_ijP_i(t)=sum_i=1^nleft(M_ijP_i(t)+M_iiP_i(t)right)=$$



$$=sum_i=1^nM_ijP_i(t)+sum_i=1^nM_iiP_i(t)=sum_i=1,ineq j^nM_ijP_i(t)+left(1-sum_j=1,ineq j^nM_jiright)P_j(t)=$$



$$=sum_i=1,ineq j^nM_ijP_i(t)+P_j(t)-sum_j=1,ineq j^nM_jiP_j(t)=$$



$$longrightarrow$$



$$P_j(t+1)-P_j(t)=sum_i=1,ineq j^nM_ijP_i(t)-sum_j=1,ineq j^nM_jiP_j(t)$$
$$=sum_i=1,ineq j^nM_ijP_i(t)-sum_j=1,ineq j^nM_jiP_j(t)+sum_i=1^nM_ii P_i(t)-sum_j=1^nM_jj P_j(t)$$



$$=sum_i=1^nM_ijP_i(t)-sum_j=1^nM_jiP_j(t)$$



At steady state



$$sum_i=1^nM_ijP_i(t)=sum_j=1^nM_jiP_j(t)$$



Which is a trivial solution, I don't understand how to obtain the steady state solution $pi=Mpi$ in matrix form the equation above.



Please advise, note that it is important for me, if possible, to remain in mattrix form for a column transition matrix.









share|cite|improve this question












share|cite|improve this question




share|cite|improve this question








edited Jul 25 at 10:31
























asked Jul 25 at 10:05









jarhead

1148




1148











  • I'm not really sure what you mean by "develop the steady state...in an index notation." Note that the $n$-dimensional vector whose entries are identically $1$ is always a solution to $pi=Mpi$.
    – Math1000
    Jul 25 at 11:56










  • $mathbf1$ will always work as pointed out above. More generally $pi$ can be any eigenvector of $M$ with eigenvalue 1. Basically you just need to find the nullspace of $M-I$
    – Casey
    Jul 25 at 20:14
















  • I'm not really sure what you mean by "develop the steady state...in an index notation." Note that the $n$-dimensional vector whose entries are identically $1$ is always a solution to $pi=Mpi$.
    – Math1000
    Jul 25 at 11:56










  • $mathbf1$ will always work as pointed out above. More generally $pi$ can be any eigenvector of $M$ with eigenvalue 1. Basically you just need to find the nullspace of $M-I$
    – Casey
    Jul 25 at 20:14















I'm not really sure what you mean by "develop the steady state...in an index notation." Note that the $n$-dimensional vector whose entries are identically $1$ is always a solution to $pi=Mpi$.
– Math1000
Jul 25 at 11:56




I'm not really sure what you mean by "develop the steady state...in an index notation." Note that the $n$-dimensional vector whose entries are identically $1$ is always a solution to $pi=Mpi$.
– Math1000
Jul 25 at 11:56












$mathbf1$ will always work as pointed out above. More generally $pi$ can be any eigenvector of $M$ with eigenvalue 1. Basically you just need to find the nullspace of $M-I$
– Casey
Jul 25 at 20:14




$mathbf1$ will always work as pointed out above. More generally $pi$ can be any eigenvector of $M$ with eigenvalue 1. Basically you just need to find the nullspace of $M-I$
– Casey
Jul 25 at 20:14















active

oldest

votes











Your Answer




StackExchange.ifUsing("editor", function ()
return StackExchange.using("mathjaxEditing", function ()
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix)
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
);
);
, "mathjax-editing");

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "69"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
convertImagesToLinks: true,
noModals: false,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);








 

draft saved


draft discarded


















StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2862256%2fderivation-of-a-stationary-distribution-of-a-markov-chain%23new-answer', 'question_page');

);

Post as a guest



































active

oldest

votes













active

oldest

votes









active

oldest

votes






active

oldest

votes










 

draft saved


draft discarded


























 


draft saved


draft discarded














StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2862256%2fderivation-of-a-stationary-distribution-of-a-markov-chain%23new-answer', 'question_page');

);

Post as a guest













































































Comments

Popular posts from this blog

What is the equation of a 3D cone with generalised tilt?

Color the edges and diagonals of a regular polygon

Relationship between determinant of matrix and determinant of adjoint?