¿Qué inspiró esta solicitud de lista de deseos?
Esta es una solicitud de Koji Uehara-sama de Nice Corporación.
Describa la característica o solución nueva o mejorada que le gustaría ver.
La prueba de Mantel para dos matrices de distancia se utiliza para comprobar la relación entre ellas. A continuación, se muestra un script de ejemplo (creado por mí).
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>¿Qué valor tendría para usted y/o para otros usuarios ofrecer esta nueva función?</STRONG> </P> <P>No lo sé.</P> <BR /> <P><STRONG>¿Es esta una solicitud de un cliente externo?</STRONG> </P> <P>Sí, Koji Uehara-sama de Nice Corporation solicitó esta función.</P> <BR /> <P><STRONG>¿Existe algún paquete de software competitivo que tenga esta característica o solución y que debamos conocer?</STRONG></P> <P>Creo que puedes hacer esto en R.</P> <BR /> <P><STRONG>¿Esto afectará una renovación, expansión o adopción por parte de un cliente? ¿Pone en riesgo una cuenta si no abordamos este problema?</STRONG></P> <P>No.</P> <BR />
Descubrí que la prueba de Mantel fue solicitada hace unos cuatro años.
https://jmpdiscovery.atlassian.net/browse/JMP-9834
El usuario puede utilizar JMP Pro (al menos ahora), y ahora tenemos la plataforma Distance Matrix en JMP Pro, por lo que parece el camino más corto para realizar la solicitud del usuario de que Distance Matrix soportará la prueba de Mantel.
