GLM起始值不被接受的问题（对数链接）
我想使用对数链接和偏移量运行高斯广义线性模型。 以下问题可能会出现： y <- c(1,1,0,0) t <- c(5,3,2,4) 没问题： exp(coef(glm(y~1 + offset(log(t)), family=poisson))) 使用 family...
在lm()中强制某些参数具有正系数
我想知道如何在lm()中约束某些参数以具有正系数。有一些包或函数（例如display）可以使所有系数和截距都为正。 例如，在这个例子中，我只想强制x1和x2具有正系数。 x1=c(NA,rnorm(99)*10) x2=c(NA,NA,rnorm(98)*10) x...
GLMER - 用于二项数据（cbind计数数据）的预测
我正在尝试预测随时间变化（x轴上的天数）我的二项数据上运行的glmer模型的值。Total Alive和Total Dead是计数数据。这是我的模型以及相应的步骤： full.model.dredge<-glmer(cbind(Total.Alive,Total.Dead)~(CO2....
在R语言中进行分数响应回归
我正在尝试对数据进行建模，其中响应变量在0和1之间，因此我决定在R中使用分数响应模型。根据我目前的理解，分数响应模型类似于逻辑回归，但它使用准似然方法来确定参数。我不确定自己是否理解正确。 到目前为止，我已经尝试了包“frm”中的“frm”和以下数据上的“glm”，该数据与此“OP”相同：O...
如何从非常大的数据集 (尤其是从 rxGlm 输出) 绘制交互效应
我目前正在使用巨大的数据集计算glm模型。无论是glm还是speedglm，都需要几天时间来计算。 我目前拥有约3M个观察值和共400个变量，只有其中一部分用于回归。在回归中，我使用了4个整数自变量(iv1、iv2、iv3、iv4)，1个二元自变量作为因子(iv5)，交互项(x * y，其中...
将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 +...
如何保存glm结果，不包含数据或只包含预测系数？
当我使用以下的R代码时，model_glm=glm(V1~. , data=xx,family="binomial"); save(file="modelfile",model_glm); 模型文件的大小将与数据一样大，在我的情况下将达到1 GB。我如何从model_glm的结果中删除数据部分...
Python中GLM的Anova检验
我正在尝试获取GLM中每个协变量的F统计量和p值。在Python中，我使用stats mode.formula.api来进行GLM分析。 formula = 'PropNo_Pred ~ Geography + log10BMI + Cat_OpCavity + CatLes_neles +...
如何在R中使用anova()函数比较两个glm模型时提取p值
因此，我正在尝试比较两个模型fit1和fit2。 最初，我只是做了anova（fit1，fit2），这产生了我理解的输出（包括p值）。 然而，当我将我的模型从基于lm（）的模型切换到基于glm（）的模型时，anova（fit1，fit2）现在产生了剩余自由度，剩余偏差和Df偏差，我很难解释...
在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,...