glm()模型的交叉验证
我想对我之前在R中建立的一些glm模型进行10倍交叉验证。我对boot包中的cv.glm()函数有点困惑，尽管我已经阅读了很多帮助文件。当我提供以下公式时:library(boot) cv.glm(data, glmfit, K=10) 这里的 "data" 参数是指整个数据集还是仅指测试集？...
"predict.glm(, type="terms")" 实际上是做什么的？（涉及IT技术）
我对R中的predict.glm函数的工作方式感到困惑。根据帮助文档， “terms”选项返回一个矩阵，给出模型公式中每个项在线性预测量上的拟合值。 因此，如果我的模型形式为f(y) = X * beta，则命令 predict(model, X, type='terms') 预计会...
理解glm$residuals和resid(glm)
你能告诉我在quasipoisson对象中，glm$residuals和resid(glm)返回什么吗？例如，我如何使用glm$y和glm$linear.predictors来创建它们。 glm$residuals返回的是拟合后残差。 n missing unique M...
在R语言中进行分数响应回归
我正在尝试对数据进行建模，其中响应变量在0和1之间，因此我决定在R中使用分数响应模型。根据我目前的理解，分数响应模型类似于逻辑回归，但它使用准似然方法来确定参数。我不确定自己是否理解正确。 到目前为止，我已经尝试了包“frm”中的“frm”和以下数据上的“glm”，该数据与此“OP”相同：O...
R：VGAM中vglm函数的不兼容维度错误
简述 我正在使用R中的VGAM包进行Tobit回归分析--这里有一个玩具数据集，它一直给我一个错误，我无法诊断出来： library(data.table) library(VGAM) > sessionInfo()$otherPkgs $VGAM Package: VGAM ...
如何保存glm结果，不包含数据或只包含预测系数？
当我使用以下的R代码时，model_glm=glm(V1~. , data=xx,family="binomial"); save(file="modelfile",model_glm); 模型文件的大小将与数据一样大，在我的情况下将达到1 GB。我如何从model_glm的结果中删除数据部分...
使用glm进行逻辑回归时的默认起始值拟合问题
我想知道在glm中如何指定默认的起始值。 这个帖子表明默认值被设置为零。而这个帖子表示有一个算法来处理起始值，然而相关链接已经失效。 我尝试用算法跟踪来拟合简单的逻辑回归模型：set.seed(123) x <- rnorm(100) p <- 1/(1 + exp(-x))...
GLM起始值不被接受的问题（对数链接）
我想使用对数链接和偏移量运行高斯广义线性模型。 以下问题可能会出现： y <- c(1,1,0,0) t <- c(5,3,2,4) 没问题： exp(coef(glm(y~1 + offset(log(t)), family=poisson))) 使用 family...
在glm函数的优化算法的每个步骤中是否有一种方法可以获得系数？
在R中执行逻辑回归时，可以使用coefficients()函数，在优化算法收敛（或未收敛）后获得系数： library(MASS) data(menarche) glm.out = glm(cbind(Menarche, Total-Menarche) ~ Age, ...