R 中線性回歸 model 后的置信區間

Question

我試圖在擬合線性回歸 model 后獲得置信區間。
我的數據包含：

Growthrate  Strains
<dbl>       <fct>

其中growthrate代表生長速率值，Strains代表菌株名稱

我的代碼：

all.lm <- lm(Growthrate~Strains,all) 

confint(all.lm, 'Strains',level=0.95)

Output：

          2.5 % 97.5 %  
Strains    NA     NA

我不明白為什么它的印刷不適用。 在這一點上的任何幫助都非常感謝。

Answer 1

問題是Strains是一個分類變量，因此沒有稱為Strains的系數。 以下是使用 R 中包含的iris數據集的示例：

data(iris)
iris.lm <- lm(Petal.Width~Species, iris)
confint(iris.lm, level=.95)   # All coefficients
#                       2.5 %    97.5 %
# (Intercept)       0.1888041 0.3031959
# Speciesversicolor 0.9991128 1.1608872
# Speciesvirginica  1.6991128 1.8608872

confint(iris.lm, "Species", level=.95)  # Species is not a coefficient
#         2.5 % 97.5 %
# Species    NA     NA

您需要按名稱或編號指定系數：

confint(iris.lm, 2:3, level=.95)
#                       2.5 %   97.5 %
# Speciesversicolor 0.9991128 1.160887
# Speciesvirginica  1.6991128 1.860887

cfnames <- names(coef(iris.lm))[2:3]
cfnames
# [1] "Speciesversicolor" "Speciesvirginica" 
confint(iris.lm, cfnames, level=.95)
#                       2.5 %   97.5 %
# Speciesversicolor 0.9991128 1.160887
# Speciesvirginica  1.6991128 1.860887