如何在交互因素中创建嵌套效应?你好,
我将使用拟合模型运行 RCBD。 我知道如何在单个因素内创建嵌套效果,但如何创建 ne相互作用因素内的特定效应?
前任。块(年 * 位置)
谢谢你!
回复:如何在交互因素中创建嵌套效果?你说, ”我还希望能够报告基因型中的个体是完全不同的。 所以我想使用 MS 而不是 REML 方法,因为 REML 不测试随机效应。 “ 随机效应被建模为方差。 REML 报告包括置信区间和 Wald 检验统计量。 您还可以检查个体的最佳线性无偏预测 (BLUP)。
你说, ”因此,当我完全按照您展示的方式进行操作时,个体被赋予唯一的 ID,基因型和 T 作为全因子效应,个体作为随机非相互作用块效应。 而且我得到了很多奇异点和缺失的 DF,就好像 JMP 在使用 MS 时试图将个体交叉到基因型和 T 一样。 (在 REML 中一切都很好)。 “ 这是 REML 的优势之一。 许多设计不适合旧的 MS 方法。
回复:如何在交互因素中创建嵌套效果?亲爱的马克,
非常感谢,这么快,在周末,你们太棒了。
是的,这几乎是我的设计,虽然每个人不是被测量 3 次,而是更像是 300 次(速度每秒记录 5 分钟,并且随着个人游泳、休息、再次游泳等而变化很大)。 当然这 300 次测量并不是完全独立的观察,但我觉得它们仍然可以用来估计个体之间的差异。 (两个相邻的测量值根本不是独立的,但经过几个游泳和休息周期后它们可能是独立的)。
是的,当然,我首先对每个人的速度进行平均,然后使用个人作为复制品。 但我也希望能够报告基因型中的个体是完全不同的。 所以我想使用 MS 而不是 REML 方法,因为 REML 不测试随机效应。
因此,当我完全按照您展示的方式进行操作时,个体被赋予唯一的 ID,基因型和 T 作为全因子效应,个体作为随机非相互作用块效应。 而且我得到了很多奇异点和缺失的 DF,就好像 JMP 在使用 MS 时试图将个体交叉到基因型和 T 一样。 (在 REML 中一切都很好)。
也许我无法估计同一模型中个体效应的重要性,应该在单独的分析中进行,比如分别针对每个 T 并嵌套在基因型中?
再次感谢您的建议!
回复:如何在交互因素中创建嵌套效果?这是一个简单的例子来说明我如何解释你的设计。

您包括 4 x 2 = 8 次治疗。 每个治疗对多个受试者进行重复(例如,5 次)。 每个受试者被观察多次(例如,3 次)。 这基本上是你所做的吗?
因此,您可以使用观察值的平均值作为每个人的响应。 这样,残差是单独的,您不需要模型中的项。 (您还可以创建观察值的方差并对其建模,以查看处理是否也影响变异性。)

或者,您可以使用个人观察结果。 您必须先堆叠单独的观察列,然后为个体的身份创建一个新列。 ID 列将作为随机效应添加。

新表如下所示:

所以模型更改为:

这样,个体贡献了随机效应,但它不是嵌套效应。 它就像一个随机块。
回复:如何在交互因素中创建嵌套效果?亲爱的马克,
感谢您 3 年前对 RCBD 设计的讨论。 我关注了这个讨论,但我仍然不确定如何指定我的模型,这实际上要简单得多。 它还包括嵌套在其他两个效果组合中的效果。 具体来说,我有 4 个基因型在 2 个温度下饲养,每个温度和表型组合重复测量 4-5 个个体(对每个个体的游泳速度进行多次观察)。 所以个体是嵌套在 genotypeXtemperature 中的随机效应。 当我按照你在上面的讨论中建议的那样指定具有唯一 ID 的个人并运行它时:
拟合模型(是( : v ) ,效果( :基因型, :吨, :吨* :基因型, :G _T_Ind &随机的) )
我迷路了 DF,没有任何评估。
必须有一种方法来指定个体在组合中是唯一的:吨* : 基因型,对吧?
提前致谢!
扬波尔斯基
回复:如何在交互因素中创建嵌套效果?你好马克,
感谢您为我比较检测能力。 我不想对待嵌套块(年*位置)效果固定,效果。 我想将其视为随机效果,但不知道如何通过嵌套来运行测试堵塞在的相互作用年和地点。但现在由于检测能力并不受简单命名的影响1~12块在 2 个地点和 2 年中,我将以这种方式进行分析。 再次感谢你!
回复:如何在交互因素中创建嵌套效果?我不知道您的响应的标准偏差或您试图检测的最小效果,但这种比较是公平的,因为我使用 1 表示标准偏差(预期 RMSE)和 1 表示参数值(带有因子水平编码,该值相当于标准偏差两倍的效应)在两种分析中。 下面是在模型中使用随机块效应的功效分析:

(模型中具有随机块效应的 RCBD)
这是使用您的模型的功率分析:

(模型中具有固定块效应的 RCBD)
所以你可以看到,对于感兴趣的效果,功效本质上是相同的。 您的选择应该归结为您如何看待块效应。 随机效应似乎适合我。 为什么要将 Block 建模为固定效应? 我并不是说这个问题是一个挑战。 这只是需要考虑的事情。
回复:如何在交互因素中创建嵌套效果?使用这两种方法检测品种差异的能力是否存在差异?
谢谢你!!
回复:如何在交互因素中创建嵌套效果?通过给每个块一个唯一的级别 (1-12),您可以间接指示嵌套。 如果您在四个环境中使用了 instread 1-3,那么您将不得不显式地嵌套这些块。
回复:如何在交互因素中创建嵌套效果?标记,
非常感谢您花费时间和精力模拟数据。
我试图治疗嵌套块 (Y*L)作为随机效应,不固定。
我想知道这两种方法在区分因素差异和减少现场差异影响方面的差异是什么? 我如何证明 myslef 没有在统计上嵌套块?
再次感谢你 !
回复:如何在交互因素中创建嵌套效果?我做了一个定制设计,应该等同于你描述的 RCBD。 (我用 Monte Carlo 模拟组成了响应值。)它包括具有随机块效应的模型(表脚本)。
我附上它供您检查。
回复:如何在交互因素中创建嵌套效果?谢谢你的澄清,这真的很有帮助。
所以可以肯定的是,我明白了,你在位置 (2) 和年份 (2) 的每个组合中都有所有 10 个变体,一个环境。 每个品种一式三份进行测试,并在每个环境中安排为 3 个复合块。 所以你有 10 x 2 x 2 x 3 = 120 个观察值,排列成 12 个块?
在我看来,模型应该将块视为没有嵌套的随机效果。 不指示 block = 1、2 或 3,而是在数据表中将其指示为 1 到 12。 选择模型中的块效应,点击属性旁边的红色小三角,然后选择随机效应。
我不明白为什么您会将块视为固定效果或需要根据您的描述进行嵌套。
回复:如何在交互因素中创建嵌套效果?你好马克,
这三个因素是:10 个品种、2 年、2 个地点(每个环境中有 3 个街区,即每年和地点)
这是我针对一个环境的 RCBD 现场设计:我正在尝试在(年*位置)中下一个块效果

回复:如何在交互因素中创建嵌套效果?可能的模型不止一种,很容易误解您的设计和模型。 请提供以下资料以作澄清?
- 你的 RCBD 是什么? 即,什么是因素、因素水平、处理和区组。 特别是,嵌套在您的设计中是如何发生的?
- 你想要的线性模型是什么?
对于线性模型,为清楚起见,请使用此命名法:
Y = I + X1 + X2[X1] + X1*X2[X1] + B
我指出有两个因素,X1(例如,站点)和 X2(例如,线路)。 它们已作为两个主效应和一个交互效应输入。 X2 中的更改嵌套在 X1 级别中。 B 是区块(例如,Day)。
(顺便说一句,我不保证我的例子是可估计的,而只是为了展示命名法是如何工作的!)
回复:如何在交互因素中创建嵌套效果?嗨,马克,
我知道如何创建交互项,但如何将“块”嵌套在 ( Year*Location ) 中?
谢谢。
回复:如何在交互因素中创建嵌套效果?您只需先选择交互中涉及的列,然后单击交叉。