Создание условий взаимодействия вручную

Question

Я пытаюсь оценить панель с фиксированными эффектами с индивидуальными тенденциями времени, используя plm, и я столкнулся с той же проблемой, что и other people. Я более чем готов использовать обходное решение, описанное в связанном вопросе CrossValidated, но не могу понять, как создать необходимые столбцы фреймов данных.

То есть, у меня есть кадр данных формы

data.frame(date=rep(1:5,times=3),id=rep(1:3,each=5)) 

и хотел бы добавить к этому кадру данных столбца для каждого id, который называется date_idX, имеет такое же значение, как date для всех наблюдений где id==X и ноль в противном случае.

Конечно, любые элегантные решения моей проблемы также будут оценены.

Answer 1

> dfrm <- data.frame(date=rep(1:5,times=3),id=rep(1:3,each=5)) 
> 
> X <-3; dfrm$time_idX <- dfrm$date*(dfrm$id==X) 
> dfrm 
    date id time_idX 
1  1 1  0 
2  2 1  0 
3  3 1  0 
4  4 1  0 
5  5 1  0 
6  1 2  0 
7  2 2  0 
8  3 2  0 
9  4 2  0 
10 5 2  0 
11 1 3  1 
12 2 3  2 
13 3 3  3 
14 4 3  4 
15 5 3  5 

Я подозреваю, что вы действительно хотели сделать это в формуле регрессии. Для этого необходима функция I(). Это псевдо-код:

regfun(form = yield ~ I(date*(id==X)), data=dfrm) 

Я не гарантирующий это будет правильное решение проблемы использования PLM, но это метод, который должен работать с обычной регрессии. Вы должны отредактировать свой вопрос, чтобы включить правильный тестовый пример.