Я оценил смешанную многоуровневую модель с двумя перехватами, используя функцию lme r-пакета nlme.
После этого я проверил автокорреляцию с помощью визуальной проверки с использованием функции plot (ACF).
В первый раз я указал maxlag = 16.
Теперь у меня две проблемы: во-первых, параметр maxlag как-то застрял, то есть все последующие графики отображаются с maxlag = 16, дажекогда maxlag установлен на другие значения.2. График обрезается при y = 0,8, даже если значение задержки 0, очевидно, равно 1.
. Далее я делюсь соответствующим реплексом в надежде получить ответы или информацию о том, как решить эти две проблемы.
Ссылка на набор данных и, если предпочтительнее, также копирование-вставка в следующий код-скрипт:
#read.dataset:
datafclr <-read.csv("datafclr.csv", header = TRUE, sep = ",", dec = ".", fill = TRUE)
#required packages:
library("Matrix")
library("nlme")
#model-estimation:
tim2 <- lme(fixed=EERTmn ~ male + female +
(male:time7c) + (female:time7c) +
(male:IERT_Cp) + (female:IERT_Cp) +
(male:IERT_Cp_Partner) + (female:IERT_Cp_Partner)-1,
control=list(maxIter=100000), data=datafclr,
random=~male + female -1|dyade/female, correlation=corAR1(), na.action=na.omit)
summary(tim2)
#checking for autocorrelation:
plot(ACF(tim2, maxlag = 16), alpha = 0.01)
Результат на следующем графике:
Это приводит к тонкому графику
Когда я изменяю maxlag:
plot(ACF(tim2, maxlag = 10), alpha = 0.01)
Это приводит к тому же графику
Заранее большое спасибо!
Лучший, Патрик