我想把下面几行SAS code 在R里实现，请教牛人指点一二。
请问SAS code里ddfm=kr 和 lsmeans 在R 里有类似的option吗？
proc glm data=sample;
class a b c;
model y = a b c*a;
random c*a /test;
lsmeans a /stderr pdiff;
run;
proc mixed data=sample;
class a b c;
model y = a b /ddfm = kr;
random c*a;
lsmeans a /pdiff;
run;
谢谢先！！！
