汽車保險公司事故分析 (二)
背景：
汽車保險公司定義「安全駕駛」為：一年之內無事故。
有 98% 安全駕駛次年也不會發生事故。
有 78% 非 安全駕駛次年次年不會發生事故 。
狀態圖 (Transition Diagram)：
機率矩陣圖 (Transition Probability Matrix)：
Si 機率：
GNU-R：
P.L <- c(0.98, 0.78) 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.1, 0.9) S.1 <- S.0 %*% P print(S.1)
穩定狀態向量 (Steady-State Vector)：
P.L <- c(0.98, 0.78) 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.1, 0.9) S.2 <- S.0 %*% P %*% P print(S.2)
GNU-R：
myEq <- array(0, c(2, 2)) myEq[1, 1] = -0.02 myEq[1, 2] = 0.02 myEq[2, 1] = 1 myEq[2, 2] = 1 myEqVal <- array(0, c(2, 1)) myEqVal[1, 1] = 0 myEqVal[2, 1] = 1 myAnswer <- solve(myEq, myEqVal) print(myAnswer)
