如何输出模型系数的估计值?

时间:2012-11-26 17:02:13

标签: r curve-fitting

很抱歉,如果这是一个新手问题,但我不明白。我使用glm()为我的数据拟合sigmoid曲线。这是有效的,我可以绘制输出,我看到一个很好的S形曲线。

但是,如何让R返回它适合的最终值?根据我的理解,R将数据拟合到logit(y) = b0 + b1x,但当我> summary(glm.out)时,我只得到

Call:
glm(formula = e$V2 ~ e$V1, family = binomial(logit), data = e)

Deviance Residuals: 
       1         2         3         4         5         6         7  
-0.00001  -0.06612  -0.15118  -0.34237   0.20874   0.08724  -0.19557  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)
(Intercept)  -24.784     20.509  -1.208    0.227
e$V1           2.073      1.725   1.202    0.229

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 4.60338  on 6  degrees of freedom
Residual deviance: 0.23388  on 5  degrees of freedom
AIC: 5.8525

Number of Fisher Scoring iterations: 8

我如何获得b0和b1?

样本数据集:

Z列可以忽略。

X   Y   Z
0.0 0.0 6
6.5 0.0 3
8.8 0.333333333333  3
10.5    0.2 10
11.1    0.0 3
11.25   0.166666666667  6
12.0    0.2 5
12.75   0.5 6
13.4    0.333333333333  3
13.5    0.2 5
14.25   0.5 6
15.0    0.333333333333  6
15.7    0.666666666667  3
15.75   0.666666666667  6
16.5    0.833333333333  6
17.25   0.555555555556  9
18.0    1.0 3

1 个答案:

答案 0 :(得分:3)

您可以通过fitted()方法获取拟合值,即fitted(glm.out)。但是,您希望估计的系数不是拟合值,并且您希望使用coef()方法,如coef(glm.out)中所示。