吸收移轉矩陣範例二 (Absorbing Transition Matrix)

• • 無肇事紀錄為 S₁ ：

    • 保持在 S₁ 狀態有 5 年。

    • 保持在 S₂ 狀態有 2.5 年。

    • 經過 7.5 年後成為 多次肇事 者。

    • 成為 多次肇事 者機率為 1。

  • 一次肇事紀錄為 S₂ ：

    • 保持在 S₂ 狀態有 3.3 年。

    • 經過 3.3 年後成為 多次肇事 者。

    • 成為 多次肇事 者機率為 1。

GNU-R Script:

myEq <- array(0, c(2, 2)) myEq[1, 1] = 0.3 myEq[1, 2] = 0 myEq[2, 1] = -0.15 myEq[2, 2] = 0.2 N = solve(myEq) myEqVal <- array(0, c(2, 1)) myEqVal[1, 1] = 1 myEqVal[2, 1] = 1 T = N %*% myEqVal myEqVal[1, 1] = 0.3 myEqVal[2, 1] = 0.05 B = N %*% myEqVal print(N) print(T) print(B)