VCOV для LDA в R - PullRequest
Новая
       66

VCOV для LDA в R

1 голос
/ 06 февраля 2020

Я выполнил линейный дискриминантный анализ с помощью команды lda: 

n=1000
library(wakefield)
y<-c("D", "ND")[r_sample_binary(n)]
x_1<-rnorm(n,0,1)
x_2<-rnorm(n,0,1)
x_3<-rnorm(n,3,2)
x_4<-rnorm(n,2,2)
x_5<-rnorm(n,0,1)
df<-data.frame(x_1,x_2,x_3,x_4,x_5,y)
y<-as.factor(y)
train=(x_1<1) # training sample criterion
df.test=df[!train,] # define test sample
y.test=y[!train]
library(MASS) # load the needed library
lda.mod=lda(y~x_1+x_2+x_3+x_4+x_5, data=df, subset=train)

В дополнение к выводу lda.mod, я хотел бы также иметь оценочную ковариационную матрицу дисперсии коэффициентов (I знаю, что знаю в этом случае, но я просто хочу смоделировать его расчет). При использовании vcov я получаю сообщение об ошибке: 

vcov(lda.mod)
Error in UseMethod("vcov") : 
  no applicable method for 'vcov' applied to an object of class "qda"

Есть ли решение этой проблемы?

...