161得票6回答
如何简洁地从数据框中写出具有多个变量的公式?

假设我有一个响应变量和一组包含三个协变量的数据(作为玩具示例):y = c(1,4,6) d = data.frame(x1 = c(4,-1,3), x2 = c(3,9,8), x3 = c(4,-4,-2)) 我想对数据进行线性回归拟合:fit = lm(y ~ d$x1 + d$x2 ...

77得票2回答
逻辑回归预测的置信区间

在 R 中,predict.lm 函数基于线性回归的结果计算预测,并提供计算这些预测的置信区间选项。根据手册,这些区间基于拟合误差方差,而不是系数的误差区间。 另一方面,predict.glm 函数基于逻辑回归和泊松回归(以及其他几种)计算预测,没有置信区间选项。我甚至很难想象如何计算此类置...

66得票4回答
警告:二项式glm的成功次数不是整数!(调查包)

我正在使用 twang 包创建倾向得分,这些得分用作二项式glm中的权重,使用 survey::svyglm 函数。代码大致如下:pscore <- ps(ppci ~ var1+var2+.........., data=dt....) dt$w <- get.weights(...

62得票3回答
如何调试“对具有2个或更多级别的因子才能应用对比”的错误?

这是我正在使用的所有变量:str(ad.train) $ Date : Factor w/ 427 levels "2012-03-24","2012-03-29",..: 4 7 12 14 19 21 24 29 31 ...

53得票5回答
从glm中提取p值

我正在运行多个回归,并且只对一个特定变量的系数和p值感兴趣。因此,在我的脚本中,我希望能够从glm摘要中提取p值(获取系数本身很容易)。我知道查看p值的唯一方法是使用summary(myReg)。还有其他方法吗?fit <- glm(y ~ x1 + x2, myData) x1Coef...

32得票1回答
caret train() 和 predict.glm() 的预测结果非常不同。

我正在尝试使用10折交叉验证来估计逻辑回归。 #import libraries library(car); library(caret); library(e1071); library(verification) #data import and preparation data(Ch...

29得票1回答
如何在R中从glm中排除特定变量?

我有50个变量。这是我在glm中如何使用它们的。var = glm(Stuff ~ ., data=mydata, family=binomial) 但是我想排除其中的2个。那么我如何具体地排除这2个呢?我希望能有像这样的东西:var = glm(Stuff ~ . # notthisstuf...

24得票1回答
多项式模型可以使用广义线性模型进行估计吗?

在分析分类数据时,我们通常使用逻辑回归来估计二项结果和一个或多个协变量之间的关系。我理解这是广义线性模型(GLM)的一种类型。在R中,可以使用glm函数并使用参数family=binomial来实现这一点。另一方面,在分类数据分析中,还有多项式模型。这些不是GLMs吗?它们不能使用glm函数在...

22得票1回答
逻辑回归 - glm中的cbind命令

我正在R中做逻辑回归。有人能否澄清一下运行这两行代码有什么区别?1. glm(Response ~ Temperature, data=temp, family = binomial(link="logit")) 2. glm(cbind(Respon...

22得票3回答
从glm中提取标准误差

我进行了一次glm,我只想提取每个系数的标准误差。我在网上看到了se.coef()函数,但它不起作用,它返回"错误:找不到函数" se.coef ""。