4

これを行うと、2 つの係数 (切片と年) が得られます。

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 式を指定する方法はありますか? 最初のデータ フレームで 2 番目の数式を使用すると、"." が原因で間違った答えが得られます。「年」も含みます。2 番目のデータ フレームでは、列が 1 つしかないためごまかしています。

4

1 に答える 1

11

あなたが探している呪文は次のとおりです。

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 の紹介ドキュメントには、数式インターフェイスの文法に関する非常に有益なセクションが含まれています。

于 2012-01-29T23:55:27.463 に答える