我试图将我的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
答案 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)