вычислить абсолютное расстояние/взвешенное скользящее среднее. Stata или R

Question

Я хочу, чтобы захватить в переменной (называемой «абсолютное сопротивление») абсолютное значение расстояния от числовой переменной (называемой «v1») между данным объектом (идентифицированным с var «id», всего 50 предметов) в течение определенного года (с указанием «длительности», «всего» 13 лет), которое подвержено заболеванию (испытать «событие», зафиксированное зависимым/исходным двоичным var, называемым «больным», кодируется 1, когда субъект испытывает событие в данном году, 0 в противном случае) и значение «v1» всех предметов, которые заболели (испытуемые, которые ранее переживали событие).

Ценность испытуемых, которые ранее заболевали, - это значение предмета «v1» во время переживания события (в год, когда испытуемые заболели). Значение «v1» во время падения больного представляет собой взвешенную скользящую среднюю, в которой последний пациент, заболевший больным, имеет тот же вес, что и все другие больные предметы. Таким образом, за каждый год в наборе данных, я хочу, чтобы вычислить расстояние («absolutedistance») для каждого здорового человека, используя следующую формулу:

absolutedistance formula

«больной» может быть испытано более чем один раз, но для простоты сосредоточимся на наборе данных, в котором нас интересует время, пока первый «больной» (общ. после первого события не будет «подвергнут цензуре» и не будет проанализирован, данные будут проанализированы с выдержкой/длительность).

Теперь все усложняется, если учесть, что некоторые испытуемые испытали событие (заболели) в том же году (давайте обозначим их как «связи»). «Связи» вносят вклад в «среднее значение v1 всех других испытуемых, которые заболели до», со «средним значением v1 этих связей». Например, три пациента заболели в 2000 году. Эти «связи» будут способствовать «среднему показателю v1 всех других пациентов, которые заболевали до», начиная с 2001 года, со средним значением v1 этих трех субъектов в 2000 году.

я оценю комментарии о том, как вычислить это в Stata или R.

для неэффективного, подверженного ошибкам ручного обходного для расчета этого первенствует см мой пост: excel workaround

Answer 1

вы можете wan't к загляните в rollapply из пакета зоопарка. Например:

z <- 1:20 
rollapply(z, 5, mean)