Доверительные интервалы для прогнозов из линейного дискриминантного анализа

Question

Я хочу нарисовать область прогноза 95% модели LDA.

Я могу нарисовать область прогноза, но без информации о достоверности.

# train model ----
mdl1 <- MASS::lda(y~ x1+ x2,data=traindata)
summary(mdl1)

# partition matrix
klaR::partimat(y~  x1+ x2, data=traindata, method="lda",prec=250,
               image.colors=c("pink","lightgreen","lightblue"),
               na.correct=NA, col.wrong=NA,
               gs=NA,
               col.mean=NA,
               main="Main title")

Я бы хотел получить что-то вроде этой статьи: https://www.sciencedirect.com/science/article/pii/S0048969712016191