问题标签 [multilevel-analysis]

我有一个两级多级数据。下面是一个例子。我有两个参与者报告了他们两天的水果消费情况。

我想将数据框从 2 级数据重组为 3 级数据。我想 a) 为水果类型创建一个变量；b) 在一列中以天为单位的巢果消费量。所以数据框如下所示：

我正在 rstanarm 中创建一个贝叶斯多级模型，以通过人口控制来预测每年政党身份的满意度（使用年份作为第二级）。（我是 rstanarm 的新手，所以我的语法可能会关闭）。公式为：

响应变量“nhssat”采用李克特 5 点量表。因为我删除截距是为了预测实际的满意度（即 3.5），所以我的先验是系数将接近 3，并且绝对不会高于 5 且低于 1。

鉴于这是一个多级模型，我无法理解有关指定先验的小插曲。考虑到响应变量的李克特性质，我需要使用什么代码来指定更准确的先验？

我有随时间推移的年度数据（纵向数据），并对许多受试者进行重复测量。我认为我需要多级建模/回归来处理随着时间的推移对同一个人的肯定相关的测量集群。目前每年的数据都在单独的表格中。

我想知道 scikit-learn 中是否有一种内置方法，比如 LinearRegression()，它能够进行多级回归，其中 Level 1 是多年来的所有数据，Level 2 是集群上的受试者（每个受试者随时间测量的集群）。如果是这样，最好将纵向数据按长度排列（每个受试者随时间的测量都在一行中）或堆叠（每年的每个测量都是它自己的行）。

有没有办法做到这一点？

我不明白这两个模型有什么区别（我使用的是 lmer4 包）

当我比较这两个模型时，我得到了不同的自由度（模型 1 的自由度小于模型 2）。但我不确定 (0+..) 符号对随机效应有何影响，以及我如何在理论水平上决定使用哪种模型（除了比较 BIC 分数）。

我需要使用 PyMC3 拟合多级线性模型，我非常喜欢glm api，因为它提供了简洁性。我想问一下是否以及如何做到这一点。我发现的这篇博文提到：

glm() 还不能很好地处理分层模型

所以我有点怀疑这实际上是否可以做到，但它是几年前写的，所以这可能已经改变了。举个例子，下面是我想用 glm api 重写的模型

我认为它看起来像这样：

但是，在尝试堆叠系数时会在此处产生错误

我正在尝试mult.icc从数据集中计算多个 ICC。数据有组 ID 列和 9 个变量。所有变量都是数字。我已经运行了以下内容：

来自https://www.rdocumentation.org/packages/multilevel/versions/2.6/topics/mult.icc

并收到错误消息：

我附上了样本数据。

任何帮助将不胜感激。

我正在尝试mult.icc从数据集中计算多个 ICC。数据有组 ID 列和 9 个变量。所有变量都是数字。我已经运行了以下内容：

来自https://www.rdocumentation.org/packages/multilevel/versions/2.6/topics/mult.icc

并收到一条错误消息：

model.frame.default 中的错误（公式 = ~DV + grpid，数据 = ，：变量“DV”的类型（列表）无效

我附上了样本数据。

网上的资料很少。任何帮助将不胜感激。

在 R 会话中，我将有以下公式：

我将如何在数学公式中写这个？顺便问一下，它是一个随机斜率随机截距模型吗？

我目前正在尝试用 R 中的分层数据估计一个有序逻辑回归模型（我的因变量有三个类别，数据来自美国 50 个州；我只需要随机截距，没有随机斜率，但我也有一个交叉水平互动）。在 Stata 中，我会使用“meologit”命令，但在 R 中找不到等效命令。是否有命令可以让我计算这样的模型？我已经阅读了 ordinal 包中的 clmm 和 clmm2 以及 mixor 的内容，但据我了解，这些内容没有多级选项。或者，我尝试过 lmer 但我认为我不能指定我想在那里计算一个序数逻辑回归模型。

有人对我有什么建议吗？

非常感谢，祝您度过愉快的一周，R.Van

我正在尝试一种新颖的分层线性模型，但数据结构让我想知道这在 R 中是否可能。我之前对模型的尝试被错误地指定（哎呀），现在我不确定如何处理这工作。我在 HLM 的课程涵盖了多级模型和交叉分类模型，但不包括 3 级双交叉分类模型。

1级：

• 对二分法评分项目的反应。分类因变量，所以我想我将使用glmer(). （约 150 万个项目回复）

2级：

• 响应嵌套在项目中 - 一个项目将有许多响应（来自不同的人），但单个响应不会链接到多个项目。
• 响应也嵌套在测试实例中 - 一个测试实例将有许多响应 (50)，但一个响应不能链接到多个测试实例。
• 项目不与测试实例嵌套，测试实例也不嵌套在项目中。一个项目将出现在多个测试实例中（每次有人使用表格 A），并且一个测试实例将与多个项目相关（测试表单上的每个项目）。

3级：

• 项目嵌套在测试表单中 - 一个表单上可以有多个项目，但（在这种情况下）项目不能出现在多个表单上。
• 测试实例嵌套在人体内——一个人可以参与多个测试实例，但一个测试实例不能由多人执行。
• 测试实例也嵌套在位置内 - 一个位置可以有多个测试实例，但一个测试实例不能出现在多个位置
• 测试表单不嵌套在人员中，人员也不嵌套在测试表单中。一个人可以填写多个测试表，一个测试表可以由多个人填写。
• 人员不嵌套在位置中，位置也不嵌套在人员中 - 一个人可以在多个位置参加考试，而几个人可以在一个位置参加考试。
• 测试表单不嵌套在位置中，并且位置不嵌套在测试表单中 - 一个测试表单可以在多个位置使用；一个位置可用于管理许多测试表格。

我假设某些位置变量可能会对特定项目的表现产生影响，但我认为这将受到诸如参加测试的人的能力之类的影响。我有兴趣探索的位置、学生和项目级别的解释变量，例如噪音级别、GPA 和主题。

如果您有任何问题或建议，请告诉我。