xtdpdml：如何绘制交互图答案

【问题标题】：xtdpdml: how to draw interaction plotsxtdpdml：如何绘制交互图
【发布时间】：2021-04-14 04:35:16
【问题描述】：

我有一个 IV、一个 DV、一个主持人和两个控制变量。

面板数据与 5 年非常平衡。

有没有办法绘制 xtdpdml 结果的交互图？

我尝试了常用的方法，（iv：自变量，取值从 1 到 6；mo：调节器，取值从 1 到 6）

egen iv = rowmean (q1 q2 q3 q4 q5 q6)
egen dv = rowmean (qd1 qd2 qd3 qd4 qd5 qd6)
egen mo = rpwmean (qm1 qm2 qm3 qm4 qm5 qm6)
egen im = iv * dv

xtdpdml dv, inv(gender birthyear) pre(iv mo im) fiml
 
summarize mo

global moa = round(r(mean) + r(sd),0.1)

global mo = round(r(mean),0.1)

global mob = round(r(mean) - r(sd),0.1)

margins, at(iv=(1(1)6) mo=($moa $mo $mob))

并且错误消息显示为“iv ambiguous abbreviation”。

【问题讨论】：

编辑问题以澄清是一个非常好的主意。随意编辑问题和更改名称并没有太大帮助，并且使我的答案更难以理解。如果这是您的真实代码，则尝试使用rpwmean() 将失败，这是rowmean() 的拼写错误。这很容易，但除此之外我无能为力。 xtdpdml 是社区贡献的，但我从未使用过它。 margins是官方代码，但是我不是很熟练。
主要问题似乎与舍入无关，但为什么margins 不理解iv。您是否还有其他名称以iv 开头的变量？理想情况下，您需要向我们展示一个可重现的最小示例，这在这里很难，但如果没有它，这个问题可能不会得到很好的答案。
1) 对。我发现“您不应该尝试使用突出显示版本（例如预测、边距）执行大多数后估计命令，因为必要的信息可能不会存储在文件中；请改用完整版本。”来自 Williams, R.、Allison, P. D. 和 Moral-Benito, E. (2018)。使用最大似然和结构方程建模的线性动态面板数据估计。 Stata 杂志，18（2），293-326。但是，我不知道这个“完整版”是什么意思。 2) 不。我没有其他以 iv 开头的变量。
我没有阅读长篇论文并尝试使用该命令。也许直接通过电子邮件联系第一作者。
如果有任何更新，我会添加评论。谢谢。

标签： panel stata

【解决方案1】：

编辑 这个答案是对使用完全不同的变量名称的问题的先前版本的回应。它只解决了一个附带问题，即如何呈现四舍五入到所需小数位数的数字结果。

你的代码没有提到pnv，所以很难评论。

我明白你想做什么，但这是一种尴尬的方式，而且容易犯错。问题是Stata以二进制工作，只能精确地保存一些十进制小数0.01（0.01）0.99，特别是0.25 0.50 0.75。那是 99 分中的 3 分；其余的必须保持为近似值。其他 Stata 代码意味着您不会总是看到近似值，但不能保证。我会试试这个，这使得四舍五入变成了真正的格式问题：

summarize sns
local snsa : di %2.1f r(mean) + r(sd)
local sns: di %2.1f r(mean) 
local snsb: di %2.1f r(mean) - r(sd)

margins, at(priv=(0(1)6) sns=(`snsa' `sns' `snsb'))

我不能排除真正的问题出在代码的其他部分。

【讨论】：

我编辑了包含 xtdpdml 代码的问题。谢谢你的帮助！我错误地认为四舍五入不会造成任何问题。但我仍然遇到同样的错误。