问题标签 [stata]

我是Stata和宏的新手。

我试图循环几个变量以从mlogit命令生成估计值，然后将它们保存在数据集中。那部分运作良好。

我遇到的问题是一个分类变量，我需要将其拆分为虚拟变量：

我以为通过设置SES为全局变量，循环会将其视为一组变量，但我错了。代码循环遍历每个变量，$SES因此我最终将每个虚拟变量trauma_main分别回归，这不是我想要的。

有没有办法“告诉”Stata 将虚拟变量视为一个块？此外，我知道我可以这样做i.SES并且使用它确实可以正常工作，但是使用的参考组不是我想要的。我已经用谷歌搜索了如何为类似的东西设置参考组i.var，但我没有找到任何有用的东西，可能是因为我使用了错误的搜索词。

提前感谢您的任何建议。

劣质煤

我目前正在设计一个调查系统（其中一个调查有很多问题，一个问题有很多答案，一个响应属于_用户、调查、问题和答案）。

我将在用户模型中拥有大量人口统计数据，并期望对各种问题等做出成千上万的回应。

例如，最终我们将要分析响应。80% 的男性喜欢香蕉，20% 的女性拥有福特之类的。

我正在研究 R、SAS 和 SPSS 等统计语言，我想知道我的数据是否需要以任何特定方式进行结构化才能被这些程序使用？还是他们都接受csv文件？

您在统计数据和构建数据模型方面有什么建议吗？

最后，SAS、SPSS 和 Stata 的成本是多少？

有谁知道用于解析/读取 .net 中的 Stata 文件的开源库？它是一种表格数据格式。

如果没有人看到可以重用的东西，我计划编写一个基于 IDataReader 的阅读器。

在不使用 Lumley 调查包的情况下，我非常感谢在 R 中指定概率权重的任何帮助。我正在使用 Imai 等人的中介包在 R 中进行中介分析，该软件包目前不支持 svyglm。

我目前正在运行的代码是：

但是，我不确定这是否正确地加权数据。原因是这段代码产生的标准错误与我在 Stata 中得到的不同。我正在运行的Stata代码是：

我想知道 R 中的权重选项是否可能不适用于逆概率权重，但我无法从文档、本论坛或其他地方确定这一点。如果我遗漏了什么，我真的很抱歉 - 我对 R 以及这个论坛都是新手。

预先感谢您的帮助。

我想在 Stata 中编写一个 TeX 文件，其中包含一个四舍五入到最接近 0.001 的 p 值。由于 sutex、outtab 等没有内置的后估计结果功能，我尝试了很多东西，比如

或者

有任何想法吗？

所以这是我的麻烦：我想将一个长格式的数据文件改成宽格式。但是，我没有唯一的“j”变量；长格式文件中的每条记录都有几个关键变量。

例如，我想这样做：

并将其变成这样：

但是，数据缺少 reshape 命令所需的后缀变量。Stata有没有办法自动生成这个后缀？

在d:\folder，我有很多Stata文件，如data_aa_1.dta、data_aa_2.dta、data_aa_3.dta、data_bb_1.dta、data_bb_2.dta、data_bb_3.dta、data_cc_1.dta....我想转换这些文件并获得与 R 中的 dta 文件一样多的数据帧。所以，我相信我必须遍历 c("aa"、"bb"、"cc") 和 c(1:3)。我尝试了类似以下的方法：

但是，这看起来是错误的——绝对是。

任何帮助将不胜感激。

谢谢！

假设你有宏

我使用引号是为了正确分组单词（MERCK SHARP DOHME 是一家公司，而不是三家不同的公司）。我正在尝试编写一个程序来遍历一个变量并在它具有 LabNames 字符串之一作为子字符串时替换它。

让我们从可以正常工作的代码部分开始。

在我的情况下，此代码根据需要进行 - 列出了 3 家不同的公司（不是 5 家）。但是，以下代码无法正确运行。3M 就坏了。

如果我们将此代码应用于

我们得到

您知道该怎么做才能使代码正确处理 3M 案例吗？

多年来，我对环境和气象变量（温度和湿度）进行了每小时测量的时间序列。从这些每小时值中，我想计算 24 小时运行平均值来创建曝光参数。为此，要求至少有 17 个每小时测量值可用，且连续缺失值不超过 6 小时。如果 24 小时内连续缺失超过 6 个小时值，则该特定日期的数据将设置为缺失。如何在 Stata 或 SAS 中实现这一点？

提前致谢

有谁知道支持固定效应、工具变量回归的 R 包，如xtivregstata（FE IV 回归）。是的，我可以只包含虚拟变量，但是当组数增加时，这将变得不可能。

谢谢！