Dear all,
    I am starting to use R for mixed models.
    For example, for the date below,
      I want to fit the model values=sample+ind(sample).
Here "sample" is a fixed effect and "ind"
should be a random effect nested in "sample".

 In Sas, it would be
    proc mixed;
    class ind sample;
    model y=sample;
    random ind(sample);

    In R, should the syntax for this model be lme(values ~ sample, random=
~1|sample/ind, data=y)?
I tried this and it seems that the results are different? Can anyone give
some hint on this.
    Thanks you so much.
      Hannah

   Here is the data structure:


  values ind sample
1  0.03325   1      1
2  0.03305   1      1
3  0.03185   1      1
4  0.03515   1      1
5  0.03375   1      1
6  0.01180   1      2
7  0.01850   1      3
8  0.02915   1      4
9  0.06200   1      5
10 0.03230   2      1
11 0.03345   2      1
12 0.03385   2      1
13 0.03605   2      1
14 0.03225   2      1
15 0.01145   2      2
16 0.01805   2      3
17 0.02950   2      4
18 0.05995   2      5
.......

	[[alternative HTML version deleted]]

