Displaying 1 result from an estimated 1 matches for "out2_glm".
Did you mean:
out1_glm
1998 Jan 20
1
R-beta: questions and comments
...times.
=============================================================
In glm the offset does not seem to work.
When I generate data like
n <- 100
beta <- 2
x <- 5*(0:n)/n
y <- x+beta*x+rnorm(n)
dat <- data.frame(y=y,x=x)
out1_glm.fit(x,y,family=gaussian(),offset=x,intercept=FALSE)
out2_glm(y~x-1,family=gaussian,offset=x,data=dat)
out_c(out1$coef,out2$coef)
then out contains two numbers differing by 1, indicating
that glm.fit seems to work and glm does not
honor the offset.
I was trying to implement the Weibull estimation routines
as described in Chapter 13 of Gen. Lin. Models 2nd...