在R语言中，glmnet()和cv.glmnet()有什么区别？
我正在开发一个项目，旨在展示一组事件对结果的潜在影响。我正在使用glmnet()包，具体使用泊松特征。这是我的代码：# de <- data imported from sql connection x <- model.matrix(~.,data = de[,2...
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中提取标准误差
我进行了一次glm，我只想提取每个系数的标准误差。我在网上看到了se.coef()函数，但它不起作用，它返回"错误：找不到函数" se.coef ""。
如何预测随机效应和固定效应模型？
最近我从STATA转到R，对于STATA命令“xtlogit,fe或re”和“predict”的R等效实现有些困难。请问是否能提供一些帮助，调整以下情形： data <- read.table("http://people.stern.nyu.edu/wgreene/Econome...
Python中GLM的Anova检验
我正在尝试获取GLM中每个协变量的F统计量和p值。在Python中，我使用stats mode.formula.api来进行GLM分析。 formula = 'PropNo_Pred ~ Geography + log10BMI + Cat_OpCavity + CatLes_neles +...
逻辑回归 - glm中的cbind命令
我正在R中做逻辑回归。有人能否澄清一下运行这两行代码有什么区别？1. glm(Response ~ Temperature, data=temp, family = binomial(link="logit")) 2. glm(cbind(Respon...
如何从非常大的数据集 (尤其是从 rxGlm 输出) 绘制交互效应
我目前正在使用巨大的数据集计算glm模型。无论是glm还是speedglm，都需要几天时间来计算。 我目前拥有约3M个观察值和共400个变量，只有其中一部分用于回归。在回归中，我使用了4个整数自变量(iv1、iv2、iv3、iv4)，1个二元自变量作为因子(iv5)，交互项(x * y，其中...
如何在R中拟合自回归泊松混合模型（计数时间序列）？
我的任务是评估各种环境变量如何影响年度人口波动。为此，我需要为时间序列计数拟合泊松自回归模型： 其中Ni,j是观察到的第i个站点在第j年的个体计数，xi,j是第i个站点在第j年的环境变量 - 这些是输入数据，其余是参数：\mui,j是第i个站点在第j年预期的个体数量，\gammaj是每年...
如何调试“对具有2个或更多级别的因子才能应用对比”的错误？
这是我正在使用的所有变量：str(ad.train) $ Date : Factor w/ 427 levels "2012-03-24","2012-03-29",..: 4 7 12 14 19 21 24 29 31 ...
如何在R中从glm中排除特定变量？
我有50个变量。这是我在glm中如何使用它们的。var = glm(Stuff ~ ., data=mydata, family=binomial) 但是我想排除其中的2个。那么我如何具体地排除这2个呢？我希望能有像这样的东西：var = glm(Stuff ~ . # notthisstuf...