[R] Gamma regression doesn't converge

2010-05-30 Thread Sherrie Jin
When I ran a Gamma regression in SAS, the algorithm converged. When I ran it in R, it keeps uncoverged even if I used 1 iterations. What was wrong? I used the following code in R: glm(y ~ x1 x2 x3, control=glm.control(maxit=1), data, family=Gamma(link="log")) [[alternative HTML ver

Re: [R] How can I fit a fixed-effect linear model or generalized linear model with method="ml"?

2010-05-30 Thread Sherrie Jin
the standard errors with method "reml"? I'm pretty new to R and haven't write any code by myself. Thanks. On Sun, May 30, 2010 at 8:01 AM, Prof Brian Ripley wrote: > On Sun, 30 May 2010, Achim Zeileis wrote: > > On Sun, 30 May 2010, Sherrie Jin wrote: >> >>

[R] How can I fit a fixed-effect linear model or generalized linear model with method="ml"?

2010-05-30 Thread Sherrie Jin
Hi, I want to fit a linear model (without any random effect) with method "ml". I tried to use "glm" I found that there is no option for "ml" or "reml" and the default one is "reml". THen I tried to use "lme" but it requires a random effect. How can I fix this problem? Of course, it's not necessa

[R] Fixed-effect linear model with method="ml"?

2010-05-29 Thread Sherrie Jin
Hi, When I try to fit a linear fixed-effect model, I was using "glm". I found that there is no option for "ml" or "reml" and the default one is "reml". How can I use "ml" instead? Of course, it's not necessary to be "glm", I am just looking for commands which allow me to fit a linear model (witho