使用具有高斯族分布的广义线性模型来预测离散量是否可以（合法方法），例如通过将高斯 GLM 的输出四舍五入到最接近的整数来预测离散量？

背景：使用glmulti进行多模型推理

glmulti是一个 R 函数/包，用于通用线性模型的自动模型选择，它在给定因变量和一组预测变量的情况下构建所有可能的通用线性模型，通过经典glm函数拟合它们，然后允许进行多模型推理（例如，使用从 AICc、BIC 派生的模型权重）。glmulti在理论上也可以与任何其他函数一起使用，该函数以与glm相同的格式返回系数、模型的对数似然和自由参数的数量（可能还有其他信息？）。

我的目标：具有鲁棒错误的多模型推理

我想使用glmulti对定量因变量的误差进行稳健建模，以防止异常值的影响。

例如，我可以假设线性模型中的误差分布为t 分布而不是正态分布。使用其峰度参数，t 分布可以有重尾，因此对异常值更稳健（与正态分布相比）。

但是，我并不致力于使用 t 分布方法。我对任何能够返回对数似然的方法感到满意，因此可以与glmulti中的多模型方法一起使用。但这意味着，不幸的是，我不能使用 R 中众所周知的稳健线性模型（例如，来自robustbase的lmRob或来自 robustbase 的lmrob ），因为它们不在对数似然框架下运行，因此无法与glmulti一起使用。

问题：我找不到适用于glmulti的稳健回归函数

发现在对数似然框架下运行的唯一稳健的 RI 线性回归函数是heavyLm（来自重包）；它在分布时对错误进行建模。不幸的是，heavyLm不适用于glmulti（至少不是开箱即用的），因为它没有用于loglik的 S3 方法（可能还有其他东西）。

为了显示：

使用数据集stackloss

正则高斯线性模型：

使用glm的默认高斯链接函数使用glmulti进行多模型推理

具有 t 分布误差的线性模型（默认为 df=4）

glmulti调用heavyLm作为拟合函数的多模型推理

给出以下错误：

如果我定义以下函数，

glmulti 可以得到对数似然，但随后发生下一个错误：

问题：鲁棒线性回归的哪个函数/包适用于 glmulti（即，表现得像 glm）？

可能有一种方法可以定义进一步的函数来让heavyLm与glmulti一起工作，但在开始这段旅程之前，我想问一下是否有人

• 知道一个健壮的线性回归函数，它 (a) 在对数似然框架下运行，并且 (b) 行为类似于glm（因此可以开箱即用地使用glmulti）。
• 得到了heavyLm已经在使用glmulti。

很感谢任何形式的帮助！

我正在使用 step() 函数在 R 中进行反向消除。现在，我正在尝试查看每个自变量如何与它们的 AIC、F 和 P 值一起排名。

接下来，我想

i) 根据 p 值从最重要（顶部）到不太重要（底部）对变量进行排序

ii) 获取每个自变量的 R 平方并显示在最后一列

如何从glm对象中获取 Z - 统计数据的值作为向量？例如，我有

如何Pr(>|z|)以与获取系数估计值相同的方式访问该列fit$coef？

是否可以使用 caret 包中的 rfe 函数在二项式 GLM 上进行特征选择？我想传递预测器和预测的数据框，让它返回最好的模型和特征。如果没有，我应该使用另一个包吗？

我有一个二元结果的回归模型。我用 glmnet 拟合了模型并得到了选定的变量及其系数。

由于 glmnet 不计算变量的重要性，我想将确切的输出（选定的变量及其系数）提供给 glm 以获取信息（标准错误等）。

我搜索了 r 文档，看来我可以使用 glm 中的“方法”选项来指定用户定义的函数。但是我没有这样做，有人可以帮助我吗？

我正在使用该twang包来创建倾向得分，这些得分在二项式 glm 中用作权重，使用survey::svyglm. 代码看起来像这样：

这会产生以下警告：

有谁知道我可能做错了什么？

我不确定这条消息在 stats.SE 上是否会更好，但总的来说，我想我会先在这里尝试。

这是我试图在 R 中运行但出现错误的代码

这是我在 R 的第一天。请帮助我。Conv当我使用view()命令时，我也可以在数据中看到。Conv是包含 1/0 的结果变量。同样在 sas 或 spss 中，我们可以选择在二进制 logit 模型中建模 1 或 0。我们如何在 R 中使用它，或者这个错误与此有关。

我有一个包含多个组的数据集“关键字”。我想将 glm 单独应用于每个组，以创建一个 glm 列表适合每个组。

我可以用for循环来做到这一点，但这不是 R 精神。相反，我尝试使用by函数来做到这一点：

问题是虽然glm通常返回一个 glm 类对象，我可以从中提取各种好东西，但返回一个列表

除拟合参数外，所有信息均已丢失。有没有办法强制保留有关列表的所有信息？

ps，我尝试使用plyr工具箱，但它卡住了，因为我的关键字中有空格。

pps，这篇文章应该有标签“by”，但我无法创建该标签（stackoverflow 的新标签），有人可以重新标记它吗？

假设我们想使用收入、年轻人、城市和地区作为回归变量来模拟美国公立学校的支出（教育）。欲了解更多信息：?Anscombe 模型：教育〜（收入+年轻+城市）*地区

具有 3 个变量和 1 个交互作用（教育~收入+年轻+城市+RegionWest:young）的模型在 BIC 方面似乎是最好的。

问题是，如何在不手动编写公式的情况下从该模型中获取 ML 对象？

在发布之前，我发现包 HH 具有一些有趣的 regsubsets 对象功能，例如summaryHH和lm.regsubsets。

lm.regsubsets做我想要的，但我认为解析交互有一些问题，还有其他选择吗？