Я пытаюсь использовать LDA в модели с использованием определенных вручную приоров, чтобы увидеть изменение частоты ошибок при изменении приоров.Я попытался написать цикл внутри функции, чтобы создать последовательность приоров, но получаю ошибку.
Library(MASS)
Classif.f <- function(liver) {
p1 <- seq(from = 0.3, to= 0.7, by = .05)
p2 <- 1-p1
lda.obj <- lda(group ~ ., data = liver, prior=c(p2,p1))
lda.obj
}
Error in lda.default(x, grouping, ...) : invalid 'prior'
из данных, я знаю, что p1= .42
и p2= .58
, поэтому я пытался исследовать значения выше и ниже этих априорных значений.
Пример данных:
R1 R2 R3 R4 R5 bev group
85 92 45 27 31 0 Severe
85 64 59 32 23 0 Not Severe
86 54 33 16 54 0 Not Severe
91 78 34 24 36 0 Not Severe
87 70 12 28 10 0 Not Severe
98 55 13 17 17 0 Not Severe