假设我有一个响应变量和一组包含三个协变量的数据(作为玩具示例):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 ...
在 R 中,predict.lm 函数基于线性回归的结果计算预测,并提供计算这些预测的置信区间选项。根据手册,这些区间基于拟合误差方差,而不是系数的误差区间。 另一方面,predict.glm 函数基于逻辑回归和泊松回归(以及其他几种)计算预测,没有置信区间选项。我甚至很难想象如何计算此类置...
我正在使用 twang 包创建倾向得分,这些得分用作二项式glm中的权重,使用 survey::svyglm 函数。代码大致如下:pscore <- ps(ppci ~ var1+var2+.........., data=dt....) dt$w <- get.weights(...
这是我正在使用的所有变量:str(ad.train) $ Date : Factor w/ 427 levels "2012-03-24","2012-03-29",..: 4 7 12 14 19 21 24 29 31 ...
我正在尝试使用10折交叉验证来估计逻辑回归。 #import libraries library(car); library(caret); library(e1071); library(verification) #data import and preparation data(Ch...
我有50个变量。这是我在glm中如何使用它们的。var = glm(Stuff ~ ., data=mydata, family=binomial) 但是我想排除其中的2个。那么我如何具体地排除这2个呢?我希望能有像这样的东西:var = glm(Stuff ~ . # notthisstuf...
在分析分类数据时,我们通常使用逻辑回归来估计二项结果和一个或多个协变量之间的关系。我理解这是广义线性模型(GLM)的一种类型。在R中,可以使用glm函数并使用参数family=binomial来实现这一点。另一方面,在分类数据分析中,还有多项式模型。这些不是GLMs吗?它们不能使用glm函数在...
我正在R中做逻辑回归。有人能否澄清一下运行这两行代码有什么区别?1. glm(Response ~ Temperature, data=temp, family = binomial(link="logit")) 2. glm(cbind(Respon...
我进行了一次glm,我只想提取每个系数的标准误差。我在网上看到了se.coef()函数,但它不起作用,它返回"错误:找不到函数" se.coef ""。