将SAS转换为(随机效应方差分析)

时间:2018-02-08 06:50:16

标签: r sas

我试图将我的SAS代码转换为随机效果ANOVA 这是我的代码:

 proc glm data=A;
 class group;
 model common = group;
 random group;
 run;

'组'是团体会员,常见的是IV。 请将此代码翻译成R代码。

(编辑)的

我的数据如下:

   id  common    group
    1       4        A
    2       2        A
    3       3        A
    4       2        B
    5       2        B
    6       3        C
    7       4        C
    8       3        C

1 个答案:

答案 0 :(得分:0)

我认为您正在寻找lme,代码可以写成:

library(nlme)

#let's say A (a dataframe) has the sample data
A$group <- as.factor(A$group)

#model
fit <- lme(common ~ group, random = ~ 1 | group, data = A, na.action=T)