是什麼促使我提出這個願望清單請求的?這只是我的個人建議,並非用戶建議。
請描述您希望看到的新的或改進的功能或解決方案。
我建議有一個功能來計算 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”。