library(monobin)
library(LGDtoolkit)
data(lgd.ds.c)
#build dummy model
rf <- c("rf_01", "rf_02", "rf_16", "rf_03", "rf_09")
for (i in 1:length(rf)) {
rf_l <- rf[i]
lgd.ds.c[, rf_l] <- sts.bin(x = lgd.ds.c[, rf_l],
y = lgd.ds.c[, "lgd"])[[2]]
}
str(lgd.ds.c)
frm <- paste0("lgd ~ ", paste(rf, collapse = " + "))
model <- lm(formula = as.formula(frm), data = lgd.ds.c)
summary(model)$coefficients
#create lgd pools
lgd.ds.c$pred <- unname(predict(model))
lgd.ds.c$pool <- sts.bin(x = lgd.ds.c$pred,
y = lgd.ds.c$lgd)[[2]]
#test homogeneity on development sample
#(the same procedure can be applied on application portfolio)
homogeneity(app.port = lgd.ds.c,
loss = "lgd",
pools = "pool",
segment = "rf_03",
segment.num = 3,
method = "t.test",
alpha = 0.05)
Run the code above in your browser using DataLab