Диадическая кривая роста в R (lme)

Question

Я установил диадическую кривую роста (модель с двумя перехватами) с моделью стабильности и влияния (Р. Гарсия) и обнаружил запаздывающую синхронность между вариабельностью сердечного ритма врачей и пациентов (ВСР).

У меня естьЯ задавался вопросом, можно ли будет адаптировать синтаксис, чтобы увидеть, когда врачи прогнозируют ВСР пациентов.Другими словами, если пациенты следуют за своими врачами, когда ВСР врачей выше или ниже.Кроме того, пациенты вложены в докторов (18 врачей и 103 пациента, один из которых имеет много дизайнов), поэтому я создал DocID для каждого врача.Как я могу включить это в синтаксис?Я был бы очень признателен за любую помощь или совет.

Мои переменные следующие:

• log_RMSSDms.2: HRV пациентов
• Doc_Pat.1: строковая переменнаяс ролевой функцией DOC или PAT
• Clog_RMSSDms.1_lag: HRV центрированных и отстающих врачей
• Clog_RMSSDms.2_lag: HRV центрированных и отстающих пациентов

Вот моймодель:

library(dyadr)
library(dplyr)
library(nlme)
ctrl <- lmeControl(msMaxIter=10000,
                   MaxIter=100000,
                   msMaxEval=10000,
                   returnObject=TRUE,
                   niterEM=10000,
                  nlmStepMax=1000)

stab_infl_2int  <- lme(log_RMSSDms.2 ~ Doc_Pat.1 + 
Doc_Pat.1:Clog_RMSSDms.2_lag 
                       + Doc_Pat.1:Clog_RMSSDms.1_lag - 1, 
                       data = Final_data_moderators_1, 
                       random = ~ Doc_Pat.1 + Doc_Pat.1:Clog_RMSSDms.2_lag 
                       + Doc_Pat.1:Clog_RMSSDms.1_lag - 1|Dyad_ID,
                       correlation = corCompSymm(form = ~1|Dyad_ID/obs_ID),  
                       weights = varIdent(form = ~1|Doc_Pat.1),
                       na.action = na.omit,
                       control = ctrl)
smallsummary(stab_infl_2int)