[Rd] [R] length 1 offset in glm (& lm)

[Kenneth Knoblauch]

Hi

Quoting Prof Brian Ripley <ripley at stats.ox.ac.uk>:

> On Wed, 25 Feb 2009, Heather Turner wrote:
>
>> This post about length 1 offsets on R help seems to have been ignored
>> (sorry deleted original email - is there a way to continue thread in
>> this case?)
>>
>> https://stat.ethz.ch/pipermail/r-help/2009-February/189352.html
>
> So let's be clear: this was the 'offset' argument' and not the offset()
> function as you refer to below.

It occurs for both

c1 <- structure(list(Contr = c(0.028, 0.043, 0.064, 0.097, 0.146, 0.219
), Correct = c(34L, 57L, 94L, 152L, 160L, 160L), Incorrect = c(126L,
103L, 66L, 8L, 0L, 0L)), .Names = c("Contr", "Correct", "Incorrect"
), row.names = c(NA, 6L), class = "data.frame")

glm(cbind(Correct, Incorrect) ~ Contr + offset(qlogis(0.25)) - 1, binomial,
	data = c1)
Error in model.frame.default(formula = cbind(Correct, Incorrect) ~ Contr +  :
   variable lengths differ (found for 'offset(qlogis(0.25))')
+

Thanks for looking into this when time permits.

Ken

--- remaining deleted



>
> -- 
> Brian D. Ripley,                  ripley at stats.ox.ac.uk
> Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
> University of Oxford,             Tel:  +44 1865 272861 (self)
> 1 South Parks Road,                     +44 1865 272866 (PA)
> Oxford OX1 3TG, UK                Fax:  +44 1865 272595



-- 
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
69500 Bron
France
tel: +33 (0)4 72 91 34 77
fax: +33 (0)4 72 91 34 61
portable: +33 (0)6 84 10 64 10
http://www.sbri.fr/members/kenneth-knoblauch.html