Code
bonds <- matrix(c(12, 688, 68, 232), nrow = 2, byrow = TRUE,
dimnames = list(c("Investment", "Speculative"),
c("Default", "No default")))
data.frame(
quantity = c("P(D)", "P(D | S)", "P(D | I)", "P(S | D)"),
value = round(c(sum(bonds[, 1]) / sum(bonds),
bonds["Speculative", "Default"] / sum(bonds["Speculative", ]),
bonds["Investment", "Default"] / sum(bonds["Investment", ]),
bonds["Speculative", "Default"] / sum(bonds[, "Default"])), 4)
) quantity value
1 P(D) 0.0800
2 P(D | S) 0.2267
3 P(D | I) 0.0171
4 P(S | D) 0.8500
