r - 如何产生杠杆统计？

Question

我知道如何使用everagePlot()来生成图，但是我找不到像megastat输出中那样为每个观察值生成统计信息的方法。

Answer 1

我认为您正在寻找帽子值(value)。

使用hatvalues(fit)。经验法则是检查比平均帽子值高2-3倍的任何观测值。我不知道要在一个漂亮的数据框中提供此信息的特定功能或包装，但您自己做起来很简单。这是一个例子:

fit <- lm(hp ~ cyl + mpg, data=mtcars) #a fake model

hatvalues(fit)

hv <- as.data.frame(hatvalues(fit))
mn <-mean(hatvalues(fit))
hv$warn <- ifelse(hv[, 'hatvalues(fit)']>3*mn, 'x3',
   ifelse(hv[, 'hatvalues(fit)']>2*mn, 'x3', '-' ))

hv

对于较大的数据集，您可以使用 subset和/或 order来仅查看hat值的某些值范围:

subset(hv, warn=="x3")
subset(hv, warn%in%c("x2", "x3"))
hv[order(hv['hatvalues(fit)']), ]

实际上，我在 R in Action一书中遇到了一个不错的绘图函数，但由于这是一本受版权保护的书，因此我不会展示Kabacoff的知识产权。但是该图对于中等规模的数据集会更好。

这是一个不错的帽子图，但您可能还需要调查:

plot(hatvalues(fit), type = "h")