...
hi, there:
First of all, i would check if they are using the same method; as i know, for SAS, it is REML as default; what about R? it may be ML;if that is true, you would expect some difference in the results.\\jingju

多谢jingju的回复，但是就这个例子来说，SAS和R是一样的，都是REML，不知道问题在哪里？

you are right. the default for R is reml too. so, what do you mean the resuts are different? jingju

运行上述程序后R的结果是：
Linear mixed-effects model fit by REML
Data: a
       AIC      BIC    logLik
  244.3839 254.9511 -116.1919

Random effects:
Formula: ~1 + NAP | fBeach
Structure: General positive-definite, Log-Cholesky parametrization
            StdDev   Corr  
(Intercept) 3.549100 (Intr)
NAP         1.715015 -0.99
Residual    2.702785      

Fixed effects: Richness ~ NAP
                Value Std.Error DF   t-value p-value
(Intercept)  6.588729 1.2647708 35  5.209425   0e+00
NAP         -2.830029 0.7229514 35 -3.914549   4e-04
而SAS的结果是：
The Mixed Procedure

                                                     Fit Statistics

                                          -2 Res Log Likelihood           238.0
                                          AIC (smaller is better)         244.0
                                          AICC (smaller is better)        244.6
                                          BIC (smaller is better)         244.6


                                              Type 3 Tests of Fixed Effects

                                                    Num     Den
                                      Effect         DF      DF    F Value    Pr > F

                                      Intercept       1       8      29.53    0.0006
                                      NAP             1       8      14.51    0.0052
可见两种结果在AIC,BIC, 分母自由度，及最终的P值上都有不同！问题到底出在哪里呢？

in sas, df =3 were used but it is 6 in r. so, i doubut the group effect in r can be identical to subject effect in sas.//jingju