Re: [R] logit GLM without intercept

2008-04-22 Thread gorgasal
Hi Robert, you can exclude the intercept by including -1 in the formula: value - as.numeric(runif(20).4) ppm - rnorm(20) glm(value~ppm-1,family=binomial) HTH Stephan Dear Statisticians, I would like to analyse my data with a GLM with binomial error distribution and logit link function.

Re: [R] logit GLM without intercept

2008-04-22 Thread Prof Brian Ripley
On Tue, 22 Apr 2008, [EMAIL PROTECTED] wrote: Hi Robert, you can exclude the intercept by including -1 in the formula: value - as.numeric(runif(20).4) ppm - rnorm(20) glm(value~ppm-1,family=binomial) Yes, but that is the same as 0 + ppm, which I tend to think is more intuitive. In so

Re: [R] logit GLM without intercept

2008-04-22 Thread Gavin Simpson
On Tue, 2008-04-22 at 17:04 +0200, [EMAIL PROTECTED] wrote: Hi Robert, you can exclude the intercept by including -1 in the formula: value - as.numeric(runif(20).4) ppm - rnorm(20) glm(value~ppm-1,family=binomial) Note that 0+ppm is the same thing, try it: glm(value ~ 0 + ppm,

[R] logit GLM without intercept

2008-04-21 Thread Robert Junker
Dear Statisticians, I would like to analyse my data with a GLM with binomial error distribution and logit link function. The point is that I want a model fitted without intercept, i.e. the fitted curve should start at y=0.5 for x=0. I tried it with the following code: glm(value~0+ppm,