汽車保險公司事故分析 (一)
背景：
汽車保險公司發現，每 10 年間，某社團：
前一年曾發生事故的成員有 23% 次年也會發生事故。
前一年為未曾發生事故的成員有 11% 次年發生事故 。
狀態圖 (Transition Diagram)：
機率矩陣圖 (Transition Probability Matrix)：
Si 機率：
GNU-R：
P.L <- c(0.23, 0.11) P <- array(0, dim=c(2,2)) P[1,1] = P.L[1] P[1,2] = 1 - P.L[1] P[2,1] = P.L[2] P[2,2] = 1 - P.L[2] S.0 <- c(0.05, 0.95) S.1 <- S.0 %*% P print(S.1)
穩定狀態向量 (Steady-State Vector)：
P.L <- c(0.23, 0.11) P <- array(0, dim=c(2,2)) P[1,1] = P.L[1] P[1,2] = 1 - P.L[1] P[2,1] = P.L[2] P[2,2] = 1 - P.L[2] S.0 <- c(0.05, 0.95) S.2 <- S.0 %*% P %*% P print(S.2)
GNU-R：
myEq <- array(0, c(2, 2)) myEq[1, 1] = 1 myEq[1, 2] = 1 myEq[2, 1] = 1 myEq[2, 2] = -3 myEqVal <- array(0, c(2, 1)) myEqVal[1, 1] = 77 myEqVal[2, 1] = -89 myAnswer <- solve(myEq, myEqVal) print(myAnswer)
