このウィッシュリストリクエストのきっかけは何ですか?
ナイス株式会社の上原浩治様からの依頼です。
希望する新しい機能や改善された機能、またはソリューションについて説明してください。
2つの距離行列のマンテル検定は、2つの距離行列間の関係性を確認する場合に用いられます。以下はサンプルスクリプト(私が作成したもの)です。
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>この新しい機能が提供されると、あなたや他のユーザーにとってどのような価値がもたらされるでしょうか?</STRONG> </P> <P>分かりません。</P> <BR /> <P><STRONG>これは外部顧客からのリクエストですか?</STRONG> </P> <P>はい、ナイス株式会社の上原浩治様からのご要望でした。</P> <BR /> <P><STRONG>私たちが知っておくべきこの機能またはソリューションを備えた競合ソフトウェア パッケージはありますか?</STRONG></P> <P>R でもこれができると思います。</P> <BR /> <P><STRONG>これは、顧客側での更新、拡張、または導入に影響しますか? また、これに対処しないとアカウントが危険にさらされますか?</STRONG></P> <P>いいえ。</P> <BR />
マンテルのテストが約 4 年前に要求されていたことがわかりました。
https://jmpdiscovery.atlassian.net/browse/JMP-9834
ユーザーは JMP Pro を使用できます (少なくとも現時点では)。また、JMP Pro には現在距離行列プラットフォームがあるため、距離行列で Mantel 検定をサポートするというユーザーの要求を実現するための最短ルートと思われます。