是否可以只使用響應變量擬合線性模型？

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簡介文檔包含一個關於公式界面語法的信息豐富的部分 。