我在最新版本的R中计算glm OR置信区间时遇到了问题,但以前从未遇到过这个问题。对于任何family="binomial"的glm,无论模型多么简单,它都可以轻松地让我提取摘要和exp(coef(model)),但是当我尝试提取confint()或exp(confint(model))时,显...
当我只使用 glm 进行训练时,一切正常,并且我甚至没有接近耗尽内存。但是当我运行 train(..., method='glm') 时,我会用完内存。 这是因为 train 在交叉验证的每次迭代(或者其他 trControl 过程)中存储了大量数据吗?我查看了 trainControl,但...
我的任务是评估各种环境变量如何影响年度人口波动。为此,我需要为时间序列计数拟合泊松自回归模型: 其中Ni,j是观察到的第i个站点在第j年的个体计数,xi,j是第i个站点在第j年的环境变量 - 这些是输入数据,其余是参数:\mui,j是第i个站点在第j年预期的个体数量,\gammaj是每年...
我正在使用 twang 包创建倾向得分,这些得分用作二项式glm中的权重,使用 survey::svyglm 函数。代码大致如下:pscore <- ps(ppci ~ var1+var2+.........., data=dt....) dt$w <- get.weights(...
简述 我正在使用R中的VGAM包进行Tobit回归分析--这里有一个玩具数据集,它一直给我一个错误,我无法诊断出来: library(data.table) library(VGAM) > sessionInfo()$otherPkgs $VGAM Package: VGAM ...
最近我从STATA转到R,对于STATA命令“xtlogit,fe或re”和“predict”的R等效实现有些困难。请问是否能提供一些帮助,调整以下情形: data <- read.table("http://people.stern.nyu.edu/wgreene/Econome...
我想对我之前在R中建立的一些glm模型进行10倍交叉验证。我对boot包中的cv.glm()函数有点困惑,尽管我已经阅读了很多帮助文件。当我提供以下公式时:library(boot) cv.glm(data, glmfit, K=10) 这里的 "data" 参数是指整个数据集还是仅指测试集?...
我进行了一次glm,我只想提取每个系数的标准误差。我在网上看到了se.coef()函数,但它不起作用,它返回"错误:找不到函数" se.coef ""。
我有一个目标变量Y和一列可能影响Y的100个维度(假设为X1...X100)。 运行我的glm并查看模型摘要后,我可以看到那些具有统计显著性的变量。我想能够选择这些变量并运行另一个模型来比较性能。是否有一种方法可以解析模型摘要并仅选择那些显著的变量?
我需要在函数内更新一个回归模型。理想情况下,该函数应适用于任何类型的模型(lm、glm、multinom、clm)。更确切地说,我需要添加一个或多个在函数内定义的协变量。以下为示例:MyUpdate <- function(model){ randData <- data...