![](/img/trans.png)
[英]Estimating the Variance Inflation Factors (VIF) of a SEM model (lavaan)
[英]Variance inflation VIF for glm caret model in R
我想计算 R 中插入符号 glm model 的方差膨胀因子 (VIF)。 这是我的代码,数据集来自 UCI:
library(caret)
library(tidyverse)
url <- paste0("https://archive.ics.uci.edu/ml/machine-learning-databases/",
"00267/data_banknote_authentication.txt")
dataset <- read_csv(url, col_names = c("varWav","skeWav","curtWav","entropy","class"))
dataset$class <- as.factor(ifelse(dataset$class == 0,"Authentic","Forgery"))
idx <- createDataPartition(dataset$class, p = 0.8, list = FALSE)
train_set <- dataset[idx,]
test_set <- dataset[-idx,]
notes_model <- train(class ~.,
data = train_set,
method = "glm")
但是当我尝试这段代码时,它会返回一个错误:
car::vif(notes_model)
UseMethod("vcov") 中的错误:没有适用于 'vcov' 的方法应用于 class "c('train', 'train.formula')" 的 object
也许我的代码错了? 请,任何帮助将不胜感激。
您可以提取最终训练的 model ,然后计算vif
:
car::vif(notes_model$finalModel)
varWav skeWav curtWav entropy
63.978111 184.323806 356.526156 1.935005
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.