Dear all,
  I want to fit a random effect model with only one random factor. I do not
want to
include the intercept term either.
  The model I using now is

  lmer(values ~ (1|lot), data=tmp)

The results are as below. How do I take out the intercept term? Or
if this is not possible for the lmer function, is it possible using lme
function in the "nlme" package?
Thank you very much in advance.
    Hanna


Linear mixed model fit by REML
Formula: values ~ (1 | lot)
   Data: resamp
    AIC    BIC logLik deviance REMLdev
 -14.21 -9.459   10.1   -23.88  -20.21
Random effects:
 Groups   Name        Variance Std.Dev.
 lot      (Intercept) 0.036077 0.18994
 Residual             0.017278 0.13144
Number of obs: 36, groups: lot, 10
Fixed effects:
            Estimate Std. Error t value
(Intercept) 99.78693    0.06421    1554

	[[alternative HTML version deleted]]

