一键模拟功能可用于根据响应公式生成统计信息的分布,该响应公式必须包含随机函数(如“随机正态”)。这是一个文章有关如何在实验设计中使用一键式仿真进行功率分析的知识。
对于固定效果模型,公式可以很简单:Y = 3+ 3 * X + Random Normal()。如果在Y模型中添加随机分组因子,我首先创建一个分组列(group,下图)以指定各组的行所有权,然后创建另一列(G),其公式为随机生成每个分组的效果,稍后将估计为BLUP,公式变为:Y = 3 + 3 * X + G + Random Normal();其中G表示组效果-有关更多详细信息,请参见随附的JMP表文件。

但是,如果如上所述指定Y公式,则BLUP的仿真(通过下图所示的仿真函数)似乎不会在每次仿真时都对G进行评估,因为我认为,随机产生的群效应不会直接出现在Y公式中;随机组效应是由G列公式而不是由Y列公式生成的。

这样,如果我尝试使用Y进行模拟,则每次模拟时唯一要做的重新评估就是Y公式中随机误差项的重新评估,也就是说,我得到的是固定效果模型模拟。然后,如何为Y指定直接生成组效果的公式?我猜需要一些循环,但是我不知道该怎么做。
我认为您是对的,问题是没有为每个模拟重新评估G列... 可以通过将整个公式包含在一栏中来纠正这种情况,也许是这样的:
3 +:X * 3 + 如果 ( 行 ()== 1个 | :group!= 落后 (:group),r = 随机正态 (),r)+ 随机正态 ()
在此,为每个组重新定义变量“ r”。 假定数据表按“组”排序。 公式中的第一个随机法线对应于组效应,第二个是残差。
希望这可以帮助。
亚当