在R语言中进行分数响应回归
我正在尝试对数据进行建模，其中响应变量在0和1之间，因此我决定在R中使用分数响应模型。根据我目前的理解，分数响应模型类似于逻辑回归，但它使用准似然方法来确定参数。我不确定自己是否理解正确。 到目前为止，我已经尝试了包“frm”中的“frm”和以下数据上的“glm”，该数据与此“OP”相同：O...
使用modelr::add_predictions进行glm模型预测
我正在尝试使用tidyverse和modelr包计算一组数据的逻辑回归预测。显然，我在add_predictions中做错了些什么，因为我没有像在stats中使用'predict'函数那样收到逻辑函数的“响应”。这应该很简单，但我想不出来，多次搜索都没有结果。library(tidyverse...
使用glm模型可以运行，但使用bigglm模型则无法运行。
我试图对320,000行数据（6个变量）运行逻辑回归。在数据的一个样本（10000）上进行逐步模型选择，得到了一个包含5个交互项的相当复杂的模型：Y~X1+ X2*X3+ X2*X4+ X2*X5+ X3*X6+ X4*X5。使用glm()函数可以拟合该模型中的10000行数据，但无法处理整个...
如何矢量化方程？
我在观看Andrew Ng教授关于广义线性模型的讲座后，尝试实现Softmax回归算法来解决K分类问题。我认为我理解他所说的一切，直到我开始编写代码来实现Softmax回归的代价函数，如下所示： 我遇到的问题是试图找出一种向量化的方法。虽然我之前已经成功地对线性回归和逻辑回归进行了向量化...
如何在R中拟合自回归泊松混合模型（计数时间序列）？
我的任务是评估各种环境变量如何影响年度人口波动。为此，我需要为时间序列计数拟合泊松自回归模型： 其中Ni,j是观察到的第i个站点在第j年的个体计数，xi,j是第i个站点在第j年的环境变量 - 这些是输入数据，其余是参数：\mui,j是第i个站点在第j年预期的个体数量，\gammaj是每年...
修改R中的glm函数以采用用户指定的链接函数
在R中的glm中，Gamma系列的默认链接函数为inverse、identity和log。现在对于我的特定问题，我需要使用响应变量Y的伽玛回归和一个修改过的链接函数log(E(Y)-1))。因此，我考虑修改R中的一些glm相关函数。有几个函数可能与此相关，我正在寻求任何有过这方面经验的人的帮助...
如何简洁地从数据框中写出具有多个变量的公式？
假设我有一个响应变量和一组包含三个协变量的数据（作为玩具示例）：y = c(1,4,6) d = data.frame(x1 = c(4,-1,3), x2 = c(3,9,8), x3 = c(4,-4,-2)) 我想对数据进行线性回归拟合：fit = lm(y ~ d$x1 + d$x2 ...
为什么从glmnet模型获取回归系数的统计摘要信息是不可取的？
我有一个二元结果的回归模型。我使用glmnet拟合了该模型，并获得了所选择的变量及其系数。 由于glmnet不计算变量重要性，我想将精确的输出（所选变量及其系数）馈送给glm以获取信息（标准误差等）。 我搜索了r文档，似乎可以使用glm中的“method”选项指定用户定义函数。但我没有成功...