Was hat mich zu dieser Wunschliste inspiriert? Dies ist nur mein persönlicher Vorschlag, nicht der eines Benutzers.
Bitte beschreiben Sie die neue oder verbesserte Funktion oder Lösung, die Sie gerne sehen würden.
Ich möchte eine Funktion zur Berechnung von BLUP für Z*a und X*beta + Z+a vorschlagen.
Wenn im aktuellen REML des Response Screenings R-seitiges REML verwendet wird, kann die Plattform weder BLUP noch bedingte Vorhersagen ausgeben. Theoretisch können wir kein eindeutiges BLUP für „a“ in „y = X*beta + Z*a + e“ erhalten. Aber wir können eindeutiges BLUP für Z*a und X*beta + Z*a berechnen. Wenn G = Z*Z`, dann ist eindeutiges BLUP für Z*a G*Inv(G + ve/vu*Identity(NRow(G))) * (y - X*b), wobei b die Parameterschätzung für Beta ist.
Um das aktuelle Verhalten zu sehen, führen Sie bitte das folgende Skript auf JMP19 aus. Es gibt keinen Befehl zum Berechnen von BLUP oder bedingter Vorhersage.
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>
<P> </P>
<P><STRONG>Welchen Wert hätte diese neue Funktion für Sie und/oder andere Benutzer?</STRONG> Bei der genomischen Selektion benötigen wir BLUP, um „bessere“ Individuen auszuwählen. Im aktuellen JMP können wir BLUPs nicht aus den Ergebnissen der Marker Relatedness-Plattform berechnen.</P>
<P><STRONG>Handelt es sich hierbei um eine externe Kundenanfrage?</STRONG> Nein.</P>
<P><STRONG>Gibt es ein konkurrierendes Softwarepaket mit dieser Funktion oder Lösung, das wir kennen sollten?</STRONG> Ich weiß es nicht.</P>
<P><STRONG>Wird dies Auswirkungen auf eine Erneuerung, Erweiterung oder Einführung bei einem Kunden haben und wird ein Konto dadurch gefährdet, wenn wir uns nicht darum kümmern?</STRONG> Nein</P>
<P> </P>In meinem Eintrag ist ein Tippfehler.
"X*beta + Z+a" ist falsch. Die richtige Formel lautet "X*beta + Z*a".
