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 (without any random effect) with "ml" method. Any suggestions? Thanks. Sherrie [[alternative HTML version deleted]]