更改泊松 GLM R 中的 Y 截距答案

【问题标题】：Change Y intercept in Poisson GLM R更改泊松 GLM R 中的 Y 截距
【发布时间】：2018-09-09 21:04:50
【问题描述】：

背景：我有以下数据运行glm 函数：

location = c("DH", "Bos", "Beth")
count = c(166, 57, 38)

#make into df
df = data.frame(location, count) 

#poisson
summary(glm(count ~ location, family=poisson))

输出：

Coefficients:

            Estimate Std. Error z value Pr(>|z|)    
(Intercept)   3.6376     0.1622  22.424  < 2e-16 ***
locationBos   0.4055     0.2094   1.936   0.0529 .  
locationDH    1.4744     0.1798   8.199 2.43e-16 ***

问题：我想更改(Intercept)，这样我就可以获得与Bos相关的所有值

我看了Change reference group using glm with binomial family 和How to force R to use a specified factor level as reference in a regression?。我尝试了那里的方法，但它不起作用，我不知道为什么。

试过了：

df1 <- within(df, location <- relevel(location, ref = 1))

#poisson
summary(glm(count ~ location, family=poisson, data = df1))

所需的输出：

Coefficients:

            Estimate Std. Error z value Pr(>|z|)    
(Intercept)    ...
locationBeth   ...
locationDH     ...

问题：我该如何解决这个问题？

【问题讨论】：

标签： r output regression intercept coefficients

【解决方案1】：

我认为您的问题是您正在修改数据框，但在您的模型中您没有使用数据框。使用模型中的data 参数来使用数据框中的数据。

location = c("DH", "Bos", "Beth")
count = c(166, 57, 38)
# make into df
df = data.frame(location, count)

请注意，location 本身就是一个 character 向量。默认情况下，data.frame() 在数据框中将其强制为factor。这样转换之后，我们就可以使用relevel来指定参考电平了。

df$location = relevel(df$location, ref = "Bos") # set Bos as reference
summary(glm(count ~ location, family=poisson, data = df))

    # Call:
    # glm(formula = count ~ location, family = poisson, data = df)
    # ...
    # Coefficients:
    #              Estimate Std. Error z value Pr(>|z|)    
    # (Intercept)    4.0431     0.1325  30.524  < 2e-16 ***
    # locationBeth  -0.4055     0.2094  -1.936   0.0529 .  
    # locationDH     1.0689     0.1535   6.963 3.33e-12 ***
    # ...

【讨论】：