У меня есть следующие данные, для которых я пытаюсь вменять недостающие значения для двух конечных переменных y и z. План представляет собой трехуровневый иерархический план, который включает повторяющиеся измерения (уровень 1) отдельных лиц (уровень 2) в кластерах (уровень 3) через дискретные интервалы времени в течение исследования.
id time y z cluster intervention
1 1 0.5 0.23 1 0
1 2 . 0.11 1 1
1 3 . . 1 1
2 1 0.15 . 2 0
2 2 0.05 0.05 2 0
2 3 . . 2 1
3 1 0.90 0.90 1 0
3 2 0.23 0.81 1 1
3 3 0.22 0.22 1 1
* 1003 data, рекомендуется сначала преобразовать данные из длинного в широкий формат (
https://www.stata.com/support/faqs/s...and-mi-impute/) в соответствии со стратегией 3: использовать многомерную нормальную модель для одновременного вменения всех кластеров. Мне удалось это успешно сделать, вот так
reshape wide y z intervention, i(id) j(time) string
id y3 z3 intervention3 y2 z2 intervention2 y1 z1 intervention1 cluster
1 . . 1 . .11 1 .5 .23 0 1
2 . . 1 .05 .05 0 .15 . 0 2
3 .45 . 1 .23 .81 1 .9 .9 0 1
Мои вопросы следующие:
Ссылка показывает пример в стратегии 3 двухуровневая модель, а у меня трехуровневая. Было бы в принципе правильным изменить форму данных дважды, а затем вменять? Если да, то как мне изменить форму во второй раз?
После того, как я изменил форму данных, очевидно, что каждая переменная в наборе данных (кроме идентификатора и кластера) была умножена на 3, чтобы отразить три момента времени (например, y1, y2, y3). Код для моего вменения должен выглядеть так.
mi set wide
mi register imputed y1 y2 y3 z1 z2 z3
mi impute mvn y1 y2 y3 z1 z2 z3 intervention1 intervention2 intervention3 cluster id, add(100) noisily --Fixed effects for cluster and id added here because they are added as random effects parameters in the analytic model
mi reshape long y z intervention, i(id) j(time) string
mi estimate: mixed y intervention time ||cluster: ||id:
mi estimate: mixed z intervention time ||cluster: ||id:
Теоретически модель вменения и модель аналити c должны содержать одни и те же переменные (включая зависимую переменную). Однако мне нужно включить время в качестве переменной с фиксированным эффектом в мою аналитическую модель c с оценкой mi, приведенной выше, которая, очевидно, была «удалена» в процессе преобразования. Как мне это согласовать?
Также обратите внимание, что в аналитической модели c я ввел одну переменную для y, вмешательства и времени, но в модели вменения (mi impute) у меня есть каждая из этих переменных. переменные повторяются 3 раза (например, y1 y2, y3). Разрешено ли это?
Можно ли дважды запустить команду оценки mi, поскольку у меня две конечные переменные?
Stata версии 15.1.
Спасибо!