以前一直用的SAS, 现在也想学学R, 但在做混合模型的时候，发现 R 和SAS 的结果总是有差别，例如下面的例子，数据见附件。
R 程序：
a =read.table("RIKZ.txt",header=T)
a$fBeach <- factor(a$Beach)
library (nlme)
mixed <- lme(Richness ~ NAP, random = ~1 +NAP| fBeach,data=a)
summary (mixed)

SAS程序
proc mixed data=rikz;
class beach;
model richness=nap /intercept;
random intercept nap/ subject=beach;
run;
但是R 和SAS的结果却不同，这是怎么回事呢？求高人指点！
该数据和R程序见于Mixed Effects Models and Extensions in Ecology with R， P109
SAS程序是我自己写的，是不是有哪里写的不对啊？