У меня есть набор данных панели из 216 человек, которые приняли участие в игре из 10 раундов, поэтому у меня 2160 наблюдений. Участников попросили внести вклад: C_it - это вклад всех i в t.
Все участники были частью группы, всего 18 групп ... группы независимы друг от друга, но C_it участников внутри одной и той же группы коррелируются.
Я хочу сделать регрессию gls в моих данных, включив в нее специфичные для субъекта c случайные эффекты и объединяя стандартные ошибки на уровне группы соответствия.
Я использовал пакет plm и функцию pggls
pdata = pdata.frame(perms, index=c("Idsubject","Period"))
attach(pdata)
test1 = pggls( contribution_it ~pfirst_dummy+tfirst_dummy+t , data = pdata, model = "random")
Группы определены в переменной в моем наборе данных, но я не знаю, как включить зависимость групп в мою регрессию.