仅针对某些变量在混合效应模型中缺少输出你好,
我试图模拟金属离子对蚯蚓不同代谢物的影响——即这是一个环境污染问题。 我有许多不同金属离子的浓度值,以及许多不同代谢物的浓度值,都基于每只蠕虫(个体)。
由于蠕虫是从许多不同的地点收集的,并且在地点内有明显的聚类 - 无论是代谢物相似性还是金属离子浓度(这一点并不奇怪,因为一些这些地点是专门采样的,因为它们被污染了)。
我试图将代谢物的影响建模为混合效应模型,使用“拟合模型”菜单,将“金属”作为固定效应,将“部位”和“部位 x 金属”作为随机效应(“金属”和“代谢物”) " 是连续变量,"site" 是分类变量)。
现在,这是我的问题。 对于某些代谢物,JMP 输出给出了我正在寻找的结果,金属离子对代谢物浓度影响的显着性检验(“固定效应检验”),但对于其他一些代谢物,它不计算值,只是给出“。”关于发生了什么的任何建议? 没有缺失值,所以样本和变量的数量是相同的——只是它给了我一些因变量的答案,而不是其他的。
谢谢你...
分析回复:仅针对某些变量在混合效应模型中缺少输出您可以将数据表作为附件或屏幕截图共享吗?
能否也分享一下JMP平台(如图)? 可能是您以混淆某些效果的方式指定了模型。
这些站点是否充当阻塞因素,具有随机效应? 这种解释是您所说的“站点内清晰的聚类”吗?
回复:仅针对某些变量在混合效应模型中缺少输出您可以将数据表作为附件或屏幕截图共享吗?
我将上传一个缩减版本 - 它与原始文件相同,但我刚刚删除了几乎所有的变量,只留下足以代表两次比较的内容。
(注意两个样本,第 21 行和第 22 行,被排除在分析之外。)
能否也分享一下JMP平台(如图)? 可能是您以混淆某些效果的方式指定了模型。
是的,我会试一试,谢谢 - 希望这就是你的意思,所选相关选项的屏幕截图?

注意,对我来说,显示的两个 Y 变量具有不同的结果,因为 JMP 成功地计算了一个的 P 值,而不是另一个。
这些站点是否充当阻塞因素,具有随机效应? 这种解释是您所说的“站点内清晰的聚类”吗?
这不是专门设计的田间试验——这些是来自土著居民的蠕虫,我们已经去那里挖掘了它们。 (为了更全面的参考:我们确实将站点分类为“控制”和“预测”,我也很想尝试对这种分类进行建模,但对于当前的问题,我正在尝试对潜在污染物(金属离子)基于每个蠕虫,而不是基于每个站点。)无论您是根据金属离子还是蠕虫代谢物浓度进行分类,这些站点都非常不同 - 例如,通过多变量方法(如 PCA)清晰地可视化。 我绝对不认为自己是使用混合效应模型的专家——我在这里的方法是基于与同事的讨论,如果我发现我正在尝试做一些不合理的事情,我不会感到震惊。 注意他们的评论之一是分配随机效应仅对分类变量真正有意义,但 JMP 确实允许将随机效应分配给连续变量:我应该假设这不是一个好主意吗?
总而言之,我可以得到一个有效的答案,使用 PCA 和 PLS 建模(基于站点的交叉验证保留集,使用与 JMP 不同的包 - 也许我应该考虑某种单变量重采样方法分析??),但我真的很想能够得到单变量情况的答案。原因有以下三个:
- 我希望能够在多变量分析的同时进行单变量分析。
- 即使我主要介绍多变量分析,能够提取关键变量作为潜在的生物标志物并更详细地分析这些也非常有用。
- 我希望能够了解发生了什么!
非常感谢 ...