Классификационные функции в линейном дискриминантном анализе в R

после завершения линейного дискриминантного анализа в R с помощью lda(), есть ли удобный способ для извлечения классификация функций для каждой группы?

по ссылке,

их не следует путать с дискриминантными функциями. Классификационные функции могут использоваться для определения, к какой группе, скорее всего, принадлежит каждый случай. Существует столько же классификационных функций, сколько и групп. Каждая функция позволяет вычислите классификационные баллы для каждого случая для каждой группы, применив формулу:

Si = ci + wi1*x1 + wi2*x2 + ... + wim*xm

в этой формуле индекс I обозначает соответствующую группу; индексы 1, 2, ..., m обозначает M переменных; ci-константа для i-й группы, wij-вес для j-й переменной при вычислении классификационного балла для i-й группы; xj-наблюдаемое значение для соответствующего случая для j-й переменной. Si является результирующим классификационный балл.

мы можем использовать классификационные функции для прямого вычисления классификационных баллов для некоторых новых наблюдений.

Я могу построить их с нуля, используя формулы учебника, но для этого требуется перестроить ряд промежуточных шагов из анализа lda. Есть ли способ получить их после факта из объекта lda?

добавлено:

если я все еще не понимаю что-то в Брендоне ответ (извините за путаницу!), похоже, ответ-нет. Предположительно, большинство пользователей могут получить необходимую им информацию от predict(), который предоставляет классификации на основе lda().

9
classification r

3 ответов

предположим, что x-ваш объект LDA:

x$terms

вы можете иметь пик на объекте, посмотрев на его структуру: 

str(x)

обновление: 

Iris <- data.frame(rbind(iris3[,,1], iris3[,,2], iris3[,,3]),Sp = rep(c("s","c","v"), rep(50,3)))
train <- sample(1:150, 75)
table(Iris$Sp[train])
z <- lda(Sp ~ ., Iris, prior = c(1,1,1)/3, subset = train)
predict(z, Iris[-train, ])$class
str(z)
List of 10
 $ prior  : Named num [1:3] 0.333 0.333 0.333
  ..- attr(*, "names")= chr [1:3] "c" "s" "v"
 $ counts : Named int [1:3] 30 25 20
  ..- attr(*, "names")= chr [1:3] "c" "s" "v"
 $ means  : num [1:3, 1:4] 6.03 5.02 6.72 2.81 3.43 ...
  ..- attr(*, "dimnames")=List of 2
  .. ..$ : chr [1:3] "c" "s" "v"
  .. ..$ : chr [1:4] "Sepal.L." "Sepal.W." "Petal.L." "Petal.W."
 $ scaling: num [1:4, 1:2] 0.545 1.655 -1.609 -3.682 -0.443 ...
  ..- attr(*, "dimnames")=List of 2
  .. ..$ : chr [1:4] "Sepal.L." "Sepal.W." "Petal.L." "Petal.W."
  .. ..$ : chr [1:2] "LD1" "LD2"
 $ lev    : chr [1:3] "c" "s" "v"
 $ svd    : num [1:2] 33.66 2.93
 $ N      : int 75
 $ call   : language lda(formula = Sp ~ ., data = Iris, prior = c(1, 1, 1)/3, subset = train)
 $ terms  :Classes 'terms', 'formula' length 3 Sp ~ Sepal.L. + Sepal.W. + Petal.L. + Petal.W.
  .. ..- attr(*, "variables")= language list(Sp, Sepal.L., Sepal.W., Petal.L., Petal.W.)
  .. ..- attr(*, "factors")= int [1:5, 1:4] 0 1 0 0 0 0 0 1 0 0 ...
  .. .. ..- attr(*, "dimnames")=List of 2
  .. .. .. ..$ : chr [1:5] "Sp" "Sepal.L." "Sepal.W." "Petal.L." ...
  .. .. .. ..$ : chr [1:4] "Sepal.L." "Sepal.W." "Petal.L." "Petal.W."
  .. ..- attr(*, "term.labels")= chr [1:4] "Sepal.L." "Sepal.W." "Petal.L." "Petal.W."
  .. ..- attr(*, "order")= int [1:4] 1 1 1 1
  .. ..- attr(*, "intercept")= int 1
  .. ..- attr(*, "response")= int 1
  .. ..- attr(*, ".Environment")=<environment: R_GlobalEnv
  .. ..- attr(*, "predvars")= language list(Sp, Sepal.L., Sepal.W., Petal.L., Petal.W.)
  .. ..- attr(*, "dataClasses")= Named chr [1:5] "factor" "numeric" "numeric" "numeric" ...
  .. .. ..- attr(*, "names")= chr [1:5] "Sp" "Sepal.L." "Sepal.W." "Petal.L." ...
 $ xlevels: Named list()
 - attr(*, "class")= chr "lda"
1

Я думаю, что ваш вопрос был ущербным ... Хорошо, может быть, не дефектный, но несколько вводящий в заблуждение, по крайней мере. Дискриминантная функция(ы) относится к расстояниям между группами, поэтому нет функции, связанной с одной группой, а скорее функция, описывающая расстояния между любыми двумя центроидами группы. Я просто ответил на более поздний вопрос и разместил пример вычисления функции оценки с использованием набора данных iris и использования его для обозначения случаев в 2d-графике предсказатели. В случае анализа 2 групп функция будет больше нуля для одной группы и меньше нуля для другой группы.

0

нет встроенного способа получить необходимую мне информацию, поэтому я написал функцию для этого:

ty.lda <- function(x, groups){
  x.lda <- lda(groups ~ ., as.data.frame(x))

  gr <- length(unique(groups))   ## groups might be factors or numeric
  v <- ncol(x) ## variables
  m <- x.lda$means ## group means

  w <- array(NA, dim = c(v, v, gr))

  for(i in 1:gr){
    tmp <- scale(subset(x, groups == unique(groups)[i]), scale = FALSE)
    w[,,i] <- t(tmp) %*% tmp
  }

  W <- w[,,1]
  for(i in 2:gr)
    W <- W + w[,,i]

  V <- W/(nrow(x) - gr)
  iV <- solve(V)

  class.funs <- matrix(NA, nrow = v + 1, ncol = gr)
  colnames(class.funs) <- paste("group", 1:gr, sep=".")
  rownames(class.funs) <- c("constant", paste("var", 1:v, sep = "."))

  for(i in 1:gr) {
    class.funs[1, i] <- -0.5 * t(m[i,]) %*% iV %*% (m[i,])
    class.funs[2:(v+1) ,i] <- iV %*% (m[i,])
  }

  x.lda$class.funs <- class.funs

  return(x.lda)
}

этот код следует формулам в численной экологии Лежандра и Лежандра (1998), стр. 625, и соответствует результатам работающего примера, начиная со страницы 626.

0