0
Я использую многомерную смешанную модель в R с помощью пакета nlme. Предположим, что x и y - переменные отклика для продольных данных, которые предполагали, что ошибка внутри группы коррелирована. Матрица остаточной ошибки представлена как: 
корреляция в многомерной смешанной модели в r
Итак, мой вопрос заключается в том, как связать корреляцию с функцией lme? Я пробовал команды corr = corComSymm(from =~ 1 | x) или corr = corAR1(from =~ 1 | x), но не работал!
здесь ан например:
# visiting time by months
time = rep(c(0,3,6,9),time = 4, 200)
# subjects
subject = rep(1:50, each = 4)
# first response variable "identity"
x = c(rep(0, 100), rep(1,100))
# second response variable "identity"
y = c(rep(1, 100), rep(0,100))
# values of both reponses variables (x_1, x_2)
value = c(rnorm(100,20,1),rnorm(100,48,1))
# variables refer to reponses variables (x_1, x_2)
variable = factor(c(rep(0,150),rep(1,50)), label=c("X","Y"))
df = data.frame(subject , time, x,y,value, variable)
library(nlme)
# fit the model that each response variable has intercept and slope (time) for each random and fixed effects
# as well as fixed effects slopes for sex and lesion, and each response has different variance
f= lme(value ~ -1 + x + y + x:time + y:time , random = ~ -1 + (x + y) + time:(x + y)|subject ,
weights = varIdent(form=~1| x),corr = corAR1(from = ~ 1|x), control=lmeControl(opt="optim"), data =df)
Error in corAR1(from = ~1 | x) : unused argument (from = ~1 | x)
Любые предложения?
Попробуйте добавить [минимальный воспроизводимый пример] (http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example), обычно включающий в себя набор данных и некоторый код чтобы показать проблему. – aosmith
Я добавил воспроизводимый пример. –
Это четвертое или пятое сообщение, которое я редактировал для уточнения тегов. Вы используете пакет 'nlme'. Вы не используете пакет 'lme4'. Пожалуйста, отметьте ваши вопросы соответствующим образом. –