Was hat Sie zu diesem Wunschzettel inspiriert?
Dies ist eine Anfrage von Koji Uehara-sama von der Nice Gesellschaft.
Bitte beschreiben Sie die neue oder verbesserte Funktion bzw. Lösung, die Sie sich wünschen.
Der Mantel-Test für zwei Distanzmatrizen wird verwendet, um die Beziehung zwischen zwei Distanzmatrizen zu überprüfen. Im Folgenden finden Sie ein Beispielskript (von mir erstellt).
Clear Symbols();
Names Default To Here(1);
///////////////////////////////////////////////////////
// Sample Data
Random Reset(11111111);
n = 8;
x1 = (1::n)`+J(n,1,Round(Random Normal(0,3),1));
x2 = (1::n)`+J(n,1,Round(Random Normal(0,3),1));
// The followings are two input ditance matrices.
d1 = Sqrt(Distance(x1,x1));
d2 = Sqrt(Distance(x2,x2));
///////////////////////////////////////////////////////
cacl permutation corr = Function({d1, d2, NSim=100000, seed=1111111}, {Default Local},
p = NCol(d1);
y1 = J(p*(p-1)/2,1,.);
y2 = J(p*(p-1)/2,1,.);
k = 0;
For(i = 1, i <= p - 1, i++,
For(j = i+1, j<=p, j++,
k++;
y1[k] = d1[i,j];
y2[k] = d2[i,j];
)
);
sr = Correlation(y1||y2)[1,2];
two tail = 0;
greater = 0;
smaller = 0;
eps = 1e-8;
r0 = J(NSim,1,.);
For(ii = 1, ii <= NSim, ii++,
Wait(0);
If(Mod(ii,5000)==0, Write("Simulation:"||Char(ii)||"\!N"));
Wait(0);
ind = Random Shuffle(1::NCol(d1));
d22 = d2[ind,ind];
y2 = J(p*(p-1)/2,1,.);
k = 0;
For(i = 1, i <= p - 1, i++,
For(j = i+1, j<=p, j++,
k++;
y2[k] = d22[i,j];
)
);
pr = Correlation(y1||y2)[1,2];
r0[ii] = pr;
two tail += (Abs(pr) >= Abs(sr)-eps);
greater += pr >= sr-eps;
smaller += pr <= sr+eps;
);
Eval List({two tail/NSim, greater/NSim, smaller/NSim, r0})
);
{two tailed p, one tailed greater p, one tailed smaller p, r0}
= cacl permutation corr(d1, d2, NSim=100000, seed=1111111);
Show(two tailed p, one tailed greater p, one tailed smaller p);
dt = As Table(r0,<
<COLUMN names=""></COLUMN>
<BR /> <P><STRONG>Welchen Nutzen hätte diese neue Funktion für Sie und/oder andere Nutzer?</STRONG> </P> <P>Ich weiß es nicht.</P> <BR /> <P><STRONG>Handelt es sich um eine externe Kundenanfrage?</STRONG> </P> <P>Ja, Koji Uehara-sama von der Nice Corporation hat diese Funktion angefordert.</P> <BR /> <P><STRONG>Gibt es eine Konkurrenzsoftware, die diese Funktion oder Lösung bietet und die wir kennen sollten?</STRONG></P> <P>Ich glaube, das geht mit R.</P> <BR /> <P><STRONG>Wird sich dies auf eine Vertragsverlängerung, -erweiterung oder -einführung bei einem Kunden auswirken, und birgt dies ein Risiko für ein Kundenkonto, wenn wir dieses Problem nicht angehen?</STRONG></P> <P>NEIN.</P> <BR />
Ich fand heraus, dass Mantels Test vor etwa vier Jahren angefordert wurde.
https://jmpdiscovery.atlassian.net/browse/JMP-9834
Der Benutzer kann (zumindest jetzt) JMP Pro verwenden, und wir haben jetzt die Distance Matrix Plattform in JMP Pro, daher scheint dies der kürzeste Weg zu sein, um den Wunsch des Benutzers zu erfüllen, dass die Distance Matrix den Mantel-Test unterstützt.
