是什么促使我提出这个愿望清单请求的?这只是我的个人建议,并非用户建议。
请描述您希望看到的新的或改进的功能或解决方案。
我建议有一个功能来计算 Z*a 和 X*beta + Z+a 的 BLUP。
在当前的“响应筛选”REML中,如果使用R端REML,平台将无法输出BLUP或条件预测。理论上,我们无法在“y = X*beta + Z*a + e”中获得“a”的唯一BLUP。但我们可以计算Z*a和X*beta + Z*a的唯一BLUP。如果G = Z*Z,则Z*a的唯一BLUP为G*Inv(G + ve/vu*Identity(NRow(G))) * (y - X*b),其中b是beta的参数估计值。
要查看当前行为,请在 JMP19 上运行以下脚本。没有计算 BLUP 或条件预测的命令。
Names Default To Here(1);
Random Reset(555555);
n = 1000;
p = 20;
z = J(n, p, Random Integer(-1,1));
/* In JMP19.0, Z is always centered eveywhere. */
For(i = 1, i <= NCol(z), i++,
z[0, i] = z[0, i] - Mean(z[0, i])
);
a = J(p,1,Random Normal());
y = 3 + z*a + J(n,1,Random Normal());
/////////////////////////////////////////////////////////////////////////
dt = As Table(y||z);
Column(dt, 1)<
<SET name=""></SET>
<BR /> <P><STRONG>如果提供此新功能,对您和/或其他用户有何价值?</STRONG>在基因组选择中,我们需要BLUP来选择“更优”的个体。在当前的JMP中,我们无法根据“标记相关性”平台的结果计算BLUP。</P> <P><STRONG>这是外部客户请求吗?</STRONG>不是。</P> <P><STRONG>有没有哪个竞争性的软件包具备这个功能或解决方案,值得我们关注?</STRONG>我不知道。</P> <P><STRONG>这会影响客户的续约、扩展或采用吗?如果我们不解决这个问题,是否会危及账户安全?</STRONG>不会。</P> <BR />
我的条目中有错别字。
“X*beta + Z+a”不正确。正确的公式是“X*beta + Z*a”。