使用 R 包 RSiena 创建 siena 数据时更改协变量的错误

Posted 2023-02-24

【中文标题】使用 R 包 RSiena 创建 siena 数据时更改协变量的错误

【英文标题】：Error of changing covariates when creating siena data by using R package RSiena

【发布时间】：2021-10-27 01:37:34

【问题描述】：

我的数据是在 2 个时间点从 200 人收集的，其中包括

A1：时间 1 的 200×200 网络矩阵（因变量） A2：时间 2 的 200×200 网络矩阵（因变量） C1：时间 1 的 200×200 网络矩阵 C2：时间 2 的 200×200 网络矩阵 S：200×1 向量 Y：200×1 向量 P：200×1 向量 N：200×1 向量

我使用的代码如下。

> A <- sienaDependent(array(c(A1, A2), dim=c(200, 200, 2)))
> C <- varDyadCovar(array(c(C1, C2), dim=c(200, 200, 2)))
> S <- coCovar(S)
> Y <- coCovar(Y)
> P <- sienaDependent(P, type = "behavior" )
> N <- sienaDependent(N, type = "behavior" )

当我尝试创建带有类似错误消息的 siena 数据对象时显示错误

> MyData <- sienaDataCreate(A, S, Y, C, P, N)

    Error in sienaDataCreate(A, S, Y, C, ：
    Changing covariates are not possible with only two waves

显示错误是因为我在时间 1 和时间 2 只有 2 个矩阵吗？ 我是否需要来自超过 3 个时间点的更多观察结果来创建变化的二元变量？ 谢谢。

【参考方案1】：

如您所料，您需要对同一个网络进行 3 次以上的观察才能获得有意义的影响效果结果。 RSiena 手册说明了各个协变量： “改变 协变量作为自变量，只有在有 3 个或更多观察值时才有意义 化时刻，因为对于 2 个观察时刻，常数和 改变协变量没有意义。”（RSiena 手册第 27 页） 一般来说，我建议您使用 RSiena 手册，尤其是当您刚接触 RSiena 时：http://www.stats.ox.ac.uk/~snijders/siena/RSiena_Manual.pdf

【讨论】：

我阅读了手册，然后解决了这个问题。还是谢谢！