问题标签 [panel-data]

我使用带有 plm 的线性面板模型创建了两个回归模型，以及使用带有 pglm 包的泊松的广义面板模型。

我现在想通过在一组散点图中绘制拟合值来以图形方式比较这两种拟合。最好沿着这些路线使用ggplot2：

我考虑过简单地使用 ggplot2's stat_smooth()，但是（也许不足为奇）它似乎无法识别我的数据的面板格式。手动提取预测值predict也似乎不适用于 pglm 模型。

如何在此图中叠加两个面板模型的预测值？

我查看了plm（面板模型的 R 包）如何实现 Breusch-Pagan 测试的随机效应，plmtest()并想知道它是否可以处理不平衡的面板。

对于不平衡的面板，我们需要另一个版本的 Breusch-Pagan 随机效应检验，正如 Baltagi/Li (1990) 给出的那样：具有不完整面板的误差分量模型的拉格朗日乘数检验，计量经济学评论，9:1, 103- 107，DOI：10.1080/07474939008800180。由于这篇论文有点难读，你也可以看看 STATA 是如何做到的：http: //www.stata.com/manuals13/xtxtregpostestimation.pdf

编辑 允许不平衡面板的修改测试现在位于 CRAN 的包中（从版本 1.6-4 开始）。

我有以下面板数据，其中一些人的观察比其他人多。

我想计算个人之间观察的相关性。R 应该忽略缺失的观察结果，例如

Mike 和 John 之间的相关性：只应考虑时间 1 到 4

John 和 Andrea 之间的相关性：仅时间 1 到 5 等。

实现这些结果的最佳方法是什么？

我一直在阅读有关reshape2R 和其他包中的包的文档，以将数据从宽格式重组为长格式，反之亦然。但是，我被困在这个问题上，因为我需要为时间创建一个新变量并按时间对我的激素测量值进行分组。我以前在 SPSS 中这样做过，尽管我正在切换到R与许多其他人一样的明显原因。我知道必须有一种简单的方法可以做到这一点R，但我无法弄清楚这一点。

该项目的数据来自一项纵向临床研究，其中在 5 个不同的时间点为每位患者测量了 20 种不同的激素（示例性数据如下 - a1 是第 1 次就诊时的激素“a”，a2 是激素“a” ' 在访问 2 等）。该研究共有 20 名患者，每个患者在电子表格 (id) 中都有唯一的标识符。激素数据（激素“a”、激素“b”等）在我的电子表格中以宽泛的形式排列如下：

当我之前在 SPSS 中执行此操作时，软件会提示我输入 id 变量以及变量名称以将纵向测量值折叠到其中。我会创建一个名为“visit”的新变量，它是从 1 到 5，对应于我对每种激素的 5 次测量。当我在 SPSS 中执行此操作时，它会以长格式创建一个新输出，如下所示：

我尝试过使用 reshape，并且该功能似乎可以工作，但是当我查看实际数据时，数字在宽格式和长格式之间混淆了。我必须做一些非常基本的错误，但我很难弄清楚。

我在为多维面板数据创建 xts 对象时遇到问题。数据如下所示：

在这篇文章中定位自己，我的 R 代码如下所示：

A 的摘要看起来像这样，对我来说这看起来不错：

但不知何故，周期性被搞砸了：

这在 axTicksByTime 函数中也可见：

当我在没有因子列的情况下输入数据时，一切正常，并且周期性正确显示为季度。

如何在 R 仍然识别数据格式的情况下正确输入我的多维数据？

我还没有找到这个问题的明确答案，所以希望有人能把我引向正确的方向！

我有一个嵌套数据框（面板数据），在多个个体中有多个观察结果。我想通过那些至少有 20 行数据的个人 (id) 对我的数据框进行子集化。

我尝试了以下方法：

但是，我仍然找到少于 20 行数据的个人。

任何人都可以提供解决方案吗？

提前致谢

我有不平衡的面板数据，其中包含一个二进制变量，指示事件是否发生。我想控制时间依赖性，所以我想创建一个变量来指示自上次事件以来经过的年数。数据按 dyad-year 组织。

这是一个可重现的示例，其中包含我想要实现的向量。谢谢！

˚

我有看起来像这样的面板数据：

随着更多的观察遵循相同的模式和更多与这个问题无关的变量。

我想建立一个 ID 处理样本，其中虚拟变量在 2010 年“切换”（当 year<2010 时为 0，当 year>=2010 时为 1）。在上面的示例数据中，1234 和 3456 将在样本中，而 2345 则不在。

我对 SAS 相当陌生，我想我对 CLASS 和 BY 语句还不够熟悉，无法弄清楚如何做到这一点。

到目前为止，我已经这样做了：

但这对数据的面板方面没有任何作用。我不知道如何做最后一步，只选择每年 trtmt_grp 为 1 的 ID。

感谢所有帮助！谢谢！

我现在一直在为我的时间序列横截面数据集苦苦挣扎，特别是在试图找到一种方法来定义每个国家和年份的一列的最大值时。我尝试了不同版本的 for 和 if/else 循环，但并没有真正成功。你能帮我找到任何线索吗？

对于我的数据结构，这将是一个可重复的小示例：

最后，我想创建一个新变量来说明给定年份和国家/地区百分比最高的主题：

最好我还会生成另一个变量，指定具有最高 perc 值的主题的确切百分比。

任何帮助将不胜感激。我发现的所有关于循环的教程都没有解决时间序列横截面问题......谢谢！

我有如下数据，我用函数创建了变量“B”：

而且我想要我的数据以这种方式并且不知道如何获得它。想到了一个“if-function”，但从未在 R 中使用过。

有人可以帮我吗？