是否可以拟合仅具有响应变量的线性模型？

Question

如果这样做，我会得到两个系数（截距和年份），

data <- data.frame(accidents=c(3,1,5,0,2,3,4), year=1:7)
glm(accidents ~ year, family=poisson(link = log), data)

Coefficients:
(Intercept)         year  
     0.7155       0.0557

但正确答案是 0.944

data <-data.frame(accidents=c(3,1,5,0,2,3,4))
glm(accidents ~ ., family=poisson(link=log), data)

Coefficients:
(Intercept)  
  0.944 

有没有办法只为响应变量指定 glm 公式？如果我将第二个公式与第一个数据框一起使用，我会得到错误的答案，因为“。”还包括“年”。在第二个数据框中，我作弊了，因为只有一列。

Answer 1

这是您正在寻找的咒语：

glm(accidents ~ 1, family=poisson(link = log), data)

将其与原始数据框架一起使用：

data <- data.frame(accidents=c(3,1,5,0,2,3,4), year=1:7)
coef(glm(accidents ~ 1, family=poisson(link = log), data))
(Intercept) 
  0.9444616 

此外，正如 Ben Bolker 提到的，R 附带的 R 介绍文档包括 有关公式界面语法的内容丰富的部分。