R 与以下 Stata 命令的等价物是什么？答案

【问题标题】：What is R's equivalent of the following Stata commands?R 与以下 Stata 命令的等价物是什么？
【发布时间】：2020-09-22 14:19:57
【问题描述】：

我最近遇到了以下 Stata 命令，我想使用 lm 和 lmer 将其转换为 R。但是我根本没有完全掌握Stata的语法，所以不知道如何将下面的Stata命令转换成R中的lm/lmer。

从互联网上做了一些研究，(1) 计算内部唯一效应。 xtreg 用于随机截取模型。所以我猜 (1) 和 (2) 都会转化为 lmer 函数。是不是和下面一样？

i(x2) 是否被指定为上述随机效应？

但是，查看guide here，模型公式如下，即通过减去组均值和加总来转换观察值。这个公式对我来说就像标准回归。还有其他资源说 xtreg 用于随机拦截模型。我现在真的很困惑。它应该是随机截距模型还是标准回归模型？我是否必须通过减去组均值并添加支持大均值来转换y 和x？

【问题讨论】：

不知道Stata不排除咨询stata.com/manuals/xtxtreg.pdf（通常的Google指向Stata 13手册中的条目，而该链接指向编写Stata时的最新版本16.1.)(
我读到了，但它并没有帮助我将 Stata 命令转换为 R。这就是我在这里问这个问题的原因。如果我根本不做任何挖掘，那我为什么会想出上面的lmer 代码？我在找人告诉我我的lmer 代码是否等同于 Stata 命令。 @尼克考克斯
我并没有指责你没有做任何研究。如果您理解并解释了 xtreg 命令在做什么，那么这可以帮助其他人理解什么应该是正确的 R 翻译，并且额外的解释不会造成任何伤害。或者，文档链接可以帮助您寻求帮助但可能不了解 Stata 的其他人。相反，如果您对我说，作为一个非常熟悉 Stata 的人，我可能不知道在 Stata 中与 R 中 xyzab 的等价物是什么，但是如果您解释了 xyzab 在统计上的作用或提供了文档链接，我会更有可能知道。
（我没有投反对票，也没有投票结束，但我猜这些行为是其他人有意为之的某种信号。如果 cmets 关于如何改进问题的问题是正确的，很难知道谁是对的可能会被误读。）
谢谢。我详细阐述了我的问题。希望这对我的困惑所在的社区有意义。

【解决方案1】：

您应该看看plm 包 - 它旨在估计所有不同类型的面板线性模型。

也就是说，你可以像上面那样做

## fixed effect model
lm(y ~ x1 + as.factor(x2), data=some_data)

## random effect model
lmer(y ~ x1 + (1|x2), data=some_data)

【讨论】：

【解决方案2】：

根据我的研究，我发现 xtreg 可用于拟合许多不同的混合或非混合模型，包括随机截距模型。

xtreg y x1, i(x2) fe 应该适合一个固定效应模型，其中特定于主题的效应被平均化，并且由此产生的边际模型对您的协变量具有总体平均效应。这与使用作为协变量的主题或集群标识符拟合常规回归模型相同（here 帖子提供了更多详细信息）。

【讨论】：