问题标签 [stata]

我有“链”格式的数据，其中有受试者获得治疗“锁”和从每个“锁”招募的受试者或“链接”。因此，我的数据既宽又长——我如何编写 Stata .DTA 程序来重塑运行模型？我的数据是这样开始的

idlock idlink1 idlink2 ...

1 10 11 ...

2 20 21 ...

21 30 31 ...

并且链接可以稍后成为锁，但它仍然是原始锁链的一部分。因此，21 是链中以 1 开头的链接。每个新锁最多有 5 个链接（idlink1-idlink5）

我正在尝试做一些逻辑回归，并且听说logisticStata 中的命令可用于在 OR 尺度上呈现系数，logit在 ln(OR) 尺度上也可以用于同样的事情。R有什么等价的吗？我用过glm，输出只有ln(OR)，所以我需要自己计算，我担心我可能会误算一些东西。

这些是关于 R 统计编程的一些新手问题，我无法在网上找到答案。 我的数据框在下面的代码中标记为“eitc”。

1）一旦我加载了一个数据框，我想看看汇总统计数据。我已经使用了以下功能：

当满足某些条件时，如何在我的数据框上找到汇总统计信息。比如我想看变量“children”大于等于1时所有变量的汇总统计，等价的Stata代码为：

2）同样，当满足某些条件时，我如何找到特定的参数？例如，当“post93”变量等于 0 且“anykids”变量等于 1 时，我想找到变量“work”的平均值。等效的 Stata 代码是：

3）理想情况下，当我运行上面的汇总统计数据时，我想找出计算中包含多少观察值/符合标准。

4）当我读入我的数据框时，很高兴看到数据集中包含了多少观察值（也许有多少行有缺失值或“NA”）。

5）另外，我一直在使用以下代码创建虚拟变量。这是正确的方法还是有更有效的方法？

我正在尝试使用tabulateStata 中的命令来创建频率的时间序列。当我尝试合并tabulate每个日期之后的输出时，就会出现问题。tabulate当所讨论的变量的值不存在观察值时，将不包括 0 作为条目。例如，如果我想计算一个班级中 10 岁、11 岁和 12 岁的学生，如果只代表其中一个组，Stata 可能会输出 (8)，因此我们不知道 8 组是哪一组学生属于：可能是 (0,8,0) 或 (0,0,8)。

如果时间序列很短，这不是问题，因为“结果”窗口会显示哪些类别已表示或未表示。我的数据有更长的时间序列。有谁知道强制Stata在这些表格中包含零的解决方案/方法？我的代码的相关部分如下：

Stata 做得好的一件事是它构造新变量的方式（见下面的例子）。如何在 R 中做到这一点？

Stata的R是否有“％in％”的等效功能？

我想在 R 中进行瓦哈卡分解。我相信，它在例如劳动经济学中用于区分解释的方差与无法解释的方差。我一直无法在 R 中找到合适的解决方案，而且我不太愿意自己创建一个（我可能会搞砸）。

无论如何，这里简要说明该过程：

http://en.wikipedia.org/wiki/Ronald_Oaxaca

Stata 有幸为此提供了一个相当不错的软件包，但 Stata 对我来说并不容易获得。

www.stata.com/meeting/5german/SINNING_stata_presentation.pdf

请注意：我还在 R-help 上发布了一条消息，但没有得到回复。我希望也可以在此列表上发布。

在此先感谢，拉斯穆斯

编辑：我做了以下功能，这似乎产生了错误的答案（呃）。我尝试按照上面的 Stata 链接进行操作，但没有如我所愿：)

我正在监视 Stata/MP（Stata/SE 的多核版本）的 CPU 和内存使用情况，但我不是 Stata 程序员（更像是 Perl 人）。

任何人都可以发布一些代码，利用公共数据集在 Stata/MP 上产生足够的负载，以便在几分钟左右使用（甚至最大化）四个 CPU 内核？

如果您可以向我提供 .do 文件和 .dta 文件（或我可能需要的任何文件），我想我可以从那里获取。提前致谢！

这个问题纯粹出于求知欲。

相对经常浏览 Python 部分后，我看到了许多与此类似的问题，其中有人要求以编程方式定义全局变量。他们中的一些人意识到 的陷阱exec，而其他人则没有。

但是，我最近一直在 Stata 中进行编程，以下是常见的：

在 Stata 的说法中，N创建了一个具有名称的本地宏，并N计算为 100。在循环的每次迭代中foreach，从x1to的值x100被分配给本地宏x。然后，循环内的行将 的平方分配i给 的扩展x，这是一个与 结尾相同的局部宏i。也就是这个循环x4展开到4^2后，又x88展开到88^2。

在 Python 中，做类似事情的方法是：

然后squares[7]等于 7^2。

这是一个非常简单的例子。stata 宏还有很多其他用途。您可以将它们用作传递要评估的函数的一种方式，例如：

周围的双刻度线r将该宏扩展了两次，首先是r1/r2然后是regress/ corr，结果是运行线性回归，y作为因变量和x作为自变量，然后显示 和 之间的相关y性x。更复杂的东西是可能的。

我的问题基本上是，stata 是否属于更大的语言类别，其中变量赋值/评估采用这种“宏赋值/扩展”形式？任何解释为什么会这样设计一种语言，和/或其他语言中类似结构的例子的加分点。

我正在编写一个 Statado文件，如果用户不提供某些参数，我想提供默认值。为此，我想检查一个宏是否未定义。

我想出了一个hacky方法来做到这一点：

但我想知道是否有一种惯用的方法来做到这一点。