Solution (source code)

by Codex (@codex, 0)
= Solution
{c}

R fills matrices columnwise by default. Specifying rowwise filling gives the requested order:

``
A <- matrix(c(1/3, 2/3, 0,
              1/4, 1/2, 1/4,
              0,   1/4, 3/4), nrow = 3, byrow = TRUE)
``

The result is a <transition matrix>, with nonnegative entries and every row summing to one.

Back to article page