Re: [R] GLM output problem

2016-09-01 Thread Anderson Eduardo
Embarrassing but that's true. I wrote 'binamial' instead of 'binomial'. I tried now with the correct spelling and everything is ok, in fact. > summary(GLM) Call: glm(formula = model, family = binomial(link = logit)) Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept)

Re: [R] GLM output problem

2016-09-01 Thread peter dalgaard
>> And use the parameters returned by GLM to contruct an equation for the >> regression model: >> >> model.eq = -0.446078 + 0.267673*x - 0.014577*I(x^2) > > ## Not what I got with your data. I got: > > Coefficients: > (Intercept)x I(x^2) > -18.5750 5.0403 -0.2845

Re: [R] GLM output problem

2016-08-31 Thread Bert Gunter
Inline. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Aug 31, 2016 at 10:03 AM, Anderson Eduardo wrote: > Hello > > I

[R] GLM output problem

2016-08-31 Thread Anderson Eduardo
Hello I have started to work with GLM and I am facing the following problem: If I take: y = c(0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0) x = 1:18 model = y ~x + I(x^2) GLM = glm(model, family=binamial(link = logit)) And use the parameters returned by GLM to contruct an equation for