如何测试RSiena中的时间异质性？

Question

如果我有多个面板，如何测试模型参数的时间差异？

例如，我使用RSiena的库存数据集：

d <- sienaDataCreate(sienaDependent(array(c(s501, s502, s503), dim=c(50, 50, 3))))
e <- getEffects(d)
e <- includeEffects(e, transTrip)
ans <- siena07(sienaAlgorithmCreate(nsub=4, n3=1000), data=d, effects=e, batch=TRUE)

我想看看transTrip效果是否与第1期到第2期不同。

Answer 1

所需要的只是将ans传递给sienaTimeTest函数：

tt <- sienaTimeTest(ans)
summary(tt)

这为参数中是否存在任何时间异质性提供了一个很好的测试。