这些R glm错误消息的含义是什么：“Error: no valid set of coefficients has been found: please supply starting values”？
以下有两个相关的问题，但它们并不是我的重复问题，因为第一个问题具有特定于数据集的解决方案，第二个问题涉及当start与offset一起提供时，glm的失败。 https://dev59.com/2Izda4cB1Zd3GeqPsOHE 我有以下数据集： library(data.tabl...
如何从glm对象中获取Z-统计量的值？
我该如何从glm对象中获取Z统计量的向量数值？例如，我有以下代码： fit <- glm(y ~ 0 + x,binomial) 我如何以与使用fit$coef获取系数估计值相同的方式访问列Pr(>|z|)？
将R中glm函数计算的标准误差估计值与SAS PROC GENMOD中计算的进行比较
我正在将一个SAS PROC GENMOD的示例转换为R，使用R中的glm。SAS代码如下： proc genmod data=data0 namelen=30; model boxcoxy=boxcoxxy ~ AGEGRP4 + AGEGRP5 + AGEGRP6 + AGEGRP7 +...
修改R中的glm函数以采用用户指定的链接函数
在R中的glm中，Gamma系列的默认链接函数为inverse、identity和log。现在对于我的特定问题，我需要使用响应变量Y的伽玛回归和一个修改过的链接函数log(E(Y)-1))。因此，我考虑修改R中的一些glm相关函数。有几个函数可能与此相关，我正在寻求任何有过这方面经验的人的帮助...
在函数内使用局部协变量更新模型的update()函数
我需要在函数内更新一个回归模型。理想情况下，该函数应适用于任何类型的模型（lm、glm、multinom、clm）。更确切地说，我需要添加一个或多个在函数内定义的协变量。以下为示例：MyUpdate <- function(model){ randData <- data...
在R中决定glm逻辑回归模型的阈值
我有一些带有预测变量和二元目标的数据。例如：df <- data.frame(a=sort(sample(1:100,30)), b= sort(sample(1:100,30)), target=c(rep(0,11),rep(1,4),rep(0,...
使用modelr::add_predictions进行glm模型预测
我正在尝试使用tidyverse和modelr包计算一组数据的逻辑回归预测。显然，我在add_predictions中做错了些什么，因为我没有像在stats中使用'predict'函数那样收到逻辑函数的“响应”。这应该很简单，但我想不出来，多次搜索都没有结果。library(tidyverse...
使用predict.glm函数时，如果type参数设置为"response"，se.fit参数设置为TRUE，会返回哪些标准误差？
我将在提供的数据上拟合模型，该数据来源于这个优秀示例，示例展示了如何在执行逻辑回归后计算响应的95％置信区间： foo <- mtcars[,c("mpg","vs")]; names(foo) <- c("x","y") mod <- glm(y ~ x, data = ...
MATLAB的glmfit与fitglm比较
我正在尝试使用MATLAB进行逻辑回归以进行分类。 在MATLAB的统计工具箱中似乎有两种不同的方法来构建广义线性模型'glmfit'和'fitglm'。 我无法弄清楚两者之间的区别。 是否有一种比另一种更可取？ 以下是函数描述链接。 http://uk.mathworks.com/hel...
如何矢量化方程？
我在观看Andrew Ng教授关于广义线性模型的讲座后，尝试实现Softmax回归算法来解决K分类问题。我认为我理解他所说的一切，直到我开始编写代码来实现Softmax回归的代价函数，如下所示： 我遇到的问题是试图找出一种向量化的方法。虽然我之前已经成功地对线性回归和逻辑回归进行了向量化...