K-sample Anderson-Darling test

------=_Part_2606_8528153.1196260346456
Content-Type: text/plain; charset=ISO-8859-1
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

hi,

I have coded the K-sample Anderson-Darling test [*].

It tests if the populations of several (more than two) random samples
are drawned from the same distribution. This test is an extension of
the Kolmogorov-Smirnov test and it is known to be more sensitive
to discrepancies in the tails of the distribution.

The implementation follows exactly [*] and it only provides the "continuous"
distribution case. Note that the compensation given in [*] for the discrete case
(with ties), is not provided here.

adt.m execute three basic checks of ad.m

Ref: [*] K-sample Anderson-Darling tests, Schotz and Stephens
          Journal of American Statistical Association, p 918, 1987.

cheers,
eric.

------=_Part_2606_8528153.1196260346456
Content-Type: application/octet-stream; name=ad.m
Content-Transfer-Encoding: base64
X-Attachment-Id: f_f9jyhbfd0
Content-Disposition: attachment; filename=ad.m

IyMKIyMgQ29weXJpZ2h0IChDKSAyMDA3IEVyaWMgQ2hhc3NhbmRlLU1vdHRpbgojIwojIyBUaGlz
IHByb2dyYW0gbWF5IGJlIGZyZWVseSBjb3BpZWQgdW5kZXIgdGhlIHRlcm1zIG9mIHRoZSAKIyMg
R05VIEdlbmVyYWwgUHVibGljIExpY2Vuc2UgKEdQTCksIHZlcnNpb24gMiwgb3IgYXQgeW91ciAK
IyMgb3B0aW9uIGFueSBsYXRlciB2ZXJzaW9uLgojIwoKZnVuY3Rpb24gW3QsYWRzXSA9IGFkKFgs
YWxwaGEpCiMjIEstc2FtcGxlIEFuZGVyc29uLURhcmxpbmcgdGVzdAojIwojIyBBRCB0ZXN0cyBp
ZiB0aGUgcG9wdWxhdGlvbnMgb2Ygc2V2ZXJhbCAobW9yZSB0aGFuIHR3bykgcmFuZG9tIHNhbXBs
ZXMKIyMgYXJlIGRyYXduZWQgZnJvbSB0aGUgc2FtZSBkaXN0cmlidXRpb24uIFRoaXMgdGVzdCBp
cyBhbiBleHRlbnNpb24gb2YKIyMgdGhlIEtvbG1vZ29yb3YtU21pcm5vdiB0ZXN0IGFuZCBpdCBp
cyBrbm93biB0byBiZSBtdWNoIG1vcmUgc2Vuc2l0aXZlCiMjIHRvIGRpc2NyZXBhbmNpZXMgaW4g
dGhlIHRhaWxzIG9mIHRoZSBkaXN0cmlidXRpb24uCiMjCiMjIFVzYWdlOiBbdCxhZHNdID0gYWQo
WCxhbHBoYSkKIyMgCiMjIFg6IGNlbGwgYXJyYXkgb2Ygcm93IHZlY3RvcnMsIG9uZSBmb3IgZWFj
aCBzYW1wbGVzLgojIyBhbHBoYTogc2lnbmlmaWNhbmNlIGxldmVsCiMjCiMjIHQgaXMgdGhlIHJl
c3VsdCBvZiB0aGUgdGVzdCAoMDogbnVsbCBoeXBvdGhlc2lzLCB0aGUgcG9wdWxhdGlvbgojIyBh
cmUgZHJhd24gZnJvbSBpZGVudGljYWwgZGlzdHJpYnV0aW9uKS4KIyMgYWRzOiBzdGFuZGFyZGl6
ZWQgQW5kZXJzb24tRGFybGluZyBzdGF0aXN0aWMKIyMgCiMjCiMjIEV4YW1wbGU6IFg9Y2VsbCgy
LDEpOyBYezF9PXJhbmRuKDEsMzIpOyBYezJ9PXJhbmRuKDEsMzIpOyB0PWFkKFgsMC4wMSkKIyMK
IyMgUmVmOiBbKl0gSy1zYW1wbGUgQW5kZXJzb24tRGFybGluZyB0ZXN0cywgU2Nob3R6IGFuZCBT
dGVwaGVucwojIyAgICAgICAgICBKb3VybmFsIG9mIEFtZXJpY2FuIFN0YXRpc3RpY2FsIEFzc29j
aWF0aW9uLCBwIDkxOCwgMTk4Ny4KCmlmIChuYXJnaW4gPCAyKQogIGVycm9yKCJVc2FnZTogW3Qs
YWRzXSA9IGFkKFgsYWxwaGEpIik7CmVuZGlmOwoKaz1sZW5ndGgoWCk7CgppZiAoazwyKQogIGVy
cm9yKCJUaGUgbnVtYmVyIG9mIHNhbXBsZXMgc2hvdWxkIGJlIGF0IGxlYXN0IDIiKTsKZW5kaWY7
CgojIHBvb2xlZCBkYXRhIHNldAoKeD1bWHs6fV07Ck49bGVuZ3RoKHgpOwp4ID0gcmVzaGFwZSAo
eCwgTiwgMSk7CltzLGlkeF0gPSBzb3J0KHgpOwoKaj0oMTpOLTEpOwpkZW49ai4qZmxpcGxyKGop
OwoKcmVzPTAuMDsgbj0wOyBIPTAuMDsKCiMgY29tcHV0ZSBBRCBzdGF0aXN0aWMKCmZvciBpPTE6
aywKICB4aSA9IFtYe2l9XTsKICBuaSA9IGxlbmd0aCh4aSk7CiAgeGkgPSByZXNoYXBlICh4aSwg
bmksIDEpOwogIGNvdW50PXplcm9zKDEsTik7CiAgY291bnQgKGZpbmQgKChpZHggPiBuKSAmIChp
ZHggPD0gbituaSkpKSA9IDE7CiAgbT1jdW1zdW0oY291bnQoMTpOLTEpKTsKICBudW09IChOKm0g
LSBqKm5pKS5eMjsKICByZXMrPSBzdW0obnVtLi9kZW4pL25pOwogIG4rPW5pOwogIEgrPTEvbmk7
CmVuZGZvcjsKCnJlcy89TjsKCiMjIHJlc2NhbGUgd2l0aCB0aGVvcmV0aWNhbCBtZWFuIGFuZCB2
YXJpYW5jZSB1bmRlciBIXzAKCm09ay0xOwoKaD1zdW0oMS4vaik7Cmc9MC4wOwpmb3IgaT0xOk4t
MgogIGcrPXN1bSgxLi8oaSsxOk4tMSkpLyhOLWkpOwplbmRmb3I7CmE9ICg0KmctNikqKGstMSkg
KyAoMTAtNipnKSpIOwpiPSAoMipnLTQpKmteMiArIDgqaCprICsgKDIqZy0xNCpoLTQpKkggLSA4
KmgrNCpnLTY7CmM9ICg2KmgrMipnLTIpKmteMiArICg0KmgtNCpnKzYpKmsgKyAoMipoLTYpKkgg
KyA0Kmg7CmQ9KDIqaCs2KSprXjItNCpoKms7CnY9IChhKk5eMytiKk5eMitjKk4rZCkvKChOLTEp
KihOLTIpKihOLTMpKTsKCmFkcz0gKHJlcyAtIG0pL3NxcnQodik7CgojIyBjb21wdXRlIHRoZSBx
dWFudGlsZSwgaW50ZXJwb2xhdGUgZnJvbSB0aGUgdGFibGUgb2YgWypdCgpteV9hbHBoYSA9IFsw
LjI1ICAgIDAuMTAgICAgMC4wNSAgICAwLjAyNSAgIDAuMDFdOwpteV9iICAgICA9IFswLjY3NSAg
IDEuMjgxICAgMS42NDUgICAxLjk2MCAgIDIuMzI2OwogICAgICAgICAgIC0wLjI0NSAgIDAuMjUw
ICAgMC42NzggICAxLjE0OSAgIDEuODIyOwogICAgICAgICAgIC0wLjEwNSAgLTAuMzA1ICAtMC4z
NjIgIC0wLjM5MSAgLTAuMzk2XTsKCm15X3EgPSBteV9iKDEsOikrbXlfYigyLDopL3NxcnQobSkr
bXlfYigzLDopL207CgppZiAoYWxwaGE+PW1pbihteV9hbHBoYSkpJihhbHBoYTw9bWF4KG15X2Fs
cGhhKSkKICBxPWludGVycDEobXlfYWxwaGEsbXlfcSxhbHBoYSk7CmVsc2UKICBxPWludGVycDEo
bXlfYWxwaGEsbXlfcSxhbHBoYSwnZXh0cmFwJyk7CmVuZGlmOwoKJSBpcyB0aGUgc3RhdGlzdGlj
IHdpdGhpbiB0aGUgYWxwaGEgcXVhbnRpbGU/Cgp0PSBhZHMgPiBxOw==
------=_Part_2606_8528153.1196260346456
Content-Type: application/octet-stream; name=adt.m
Content-Transfer-Encoding: base64
X-Attachment-Id: f_f9jyhq5p1
Content-Disposition: attachment; filename=adt.m

Y2xlYXIKCmFscGhhPTAuMDUKCnByaW50ZigidGVzdCAxIC0tLS0tLS0tLS0tLS0tLS0tLS0tLS0t
LS0tLS0tLS0tLS0tLS0tXG4iKTsKCngxPVsxIDIgM10KeDI9WzQgNSA2IDddClggPWNlbGwoMiwx
KTsKWHsxfT14MTsgWHsyfT14MjsKClt0LGFkc109YWQoWCxhbHBoYSkKCnByaW50Zigic3RhbmRh
cmRpemVkIHN0YXQgc2hvdWxkIGJlIDIuOTU2OCwgdGVzdCA9IDFcbiIpOwoKcHJpbnRmKCJ0ZXN0
IDIgLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS1cbiIpOwoKeDE9WzAuMDU3
NzU0IC0yLjM1MjAyOSAtMC44MjEzMDcgLTAuMjUzNTQ2XQp4Mj1bLTIuMDI1NDQgLTIuNDkyMzAg
MC42OTc1OCAtMC41MDcxMF0KWCA9Y2VsbCgyLDEpOwpYezF9PXgxOyBYezJ9PXgyOwoKW3QsYWRz
XT1hZChYLGFscGhhKQoKcHJpbnRmKCJzdGFuZGFyZGl6ZWQgc3RhdCBzaG91bGQgYmUgLTAuOTUy
NzMsIHRlc3QgPSAwXG4iKTsKCnByaW50ZigidGVzdCAzIC0tLS0tLS0tLS0tLS0tLS0tLS0tLS0t
LS0tLS0tLS0tLS0tLS0tXG4iKTsKCngxPXJhbmRuKDEsMzIpOwp4Mj1yYW5kbigxLDMyKTsKWCA9
Y2VsbCgyLDEpOwpYezF9PXgxOyBYezJ9PXgyOwoKW3QsYWRzXT1hZChYLGFscGhhKQoKcHJpbnRm
KCJ0ZXN0IHNob3VsZCBiZSAwXG4iKTsKCnByaW50ZigidGVzdCA0IC0tLS0tLS0tLS0tLS0tLS0t
LS0tLS0tLS0tLS0tLS0tLS0tLS0tXG4iKTsKCiMjIGZyb20gSy1zYW1wbGUgQW5kZXJzb24tRGFy
bGluZyB0ZXN0cywgU2Nob3R6IGFuZCBTdGVwaGVucwojIyAgICAgIEpvdXJuYWwgb2YgQW1lcmlj
YW4gU3RhdGlzdGljYWwgQXNzb2NpYXRpb24sIHAgOTIyLCAxOTg3LgoKYWxwaGE9MC4wNQoKeDE9
WzM4LjcgNDEuNSA0My44IDQ0LjUgNDUuNSA0Ni4wIDQ3LjcgNTguMF0KeDI9WzM5LjIgMzkuMyAz
OS43IDQxLjQgNDEuOCA0Mi45IDQzLjMgNDUuOF0KeDM9WzM0LjAgMzUuMCAzOS4wIDQwLjAgNDMu
MCA0My4wIDQ0LjAgNDUuMF0gCng0PVszNC4wIDM0LjggMzQuOCAzNS40IDM3LjIgMzcuOCA0MS4y
IDQyLjhdClg9Y2VsbCg0LDEpOwpYezF9PXgxOwpYezJ9PXgyOwpYezN9PXgzOwpYezR9PXg0OwoK
W3QsYWRzXT1hZChYLGFscGhhKQoKcHJpbnRmKCJzdGFuZGFyZGl6ZWQgc3RhdCBzaG91bGQgYmUg
NC40NDksIHRlc3QgPSAxXG4iKTsK
------=_Part_2606_8528153.1196260346456
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

_______________________________________________
Octave-sources mailing list
[email protected]
https://www.cae.wisc.edu/mailman/listinfo/octave-sources

------=_Part_2606_8528153.1196260346456--