M statistic commands: Interpoint distance distribution analysis
We implement the commands mstat and mtest to perform inference based on the M statistic, a statistic that can be used to compare the interpoint distance distribution across groups of observations. The analyses are based on the study of the interpoint distances between n points in a k-dimensional setting to produce a one-dimensional real-valued test statistic. The locations are distributed in a region of the plane. When we consider all (n 2) interpoint distances, the dependencies among them are difficult to express 2 analytically, but their distribution is informative, and the M statistic can be built to summarize one aspect of this information. The two commands can be used on a wide class of datasets to test the null hypothesis that two groups have the same (spatial) distribution. mstat and mtest return the exact M test statistic. Moreover, mtest executes a Monte Carlo–type permutation test, which returns the empirical p-value together with its confidence interval. This is the command to use in most situations, because the convergence of M to its asymptotic chi-squared distribution is slow. Both commands can be used to obtain graphical output of the empirical density function of the interpoint distance distributions in the two groups and the two- dimensional map of the n observations in the plane. The descriptions of the commands are accompanied by examples of applications with real and simulated data. We run the test on the Alt and Vach grave site dataset (Manjourides and Pagano, forthcoming, Statistics in Medicine) and reject the null hypothesis, in contradiction to other published analyses. We also show how to adapt the techniques to discrete datasets with more than one unit in each location. Finally, we report an extensive application on breast cancer data in Massachusetts; in the application, we show the compatibility of the M commands with Pisati's spmap package. Copyright 2011 by StataCorp LP.
Volume (Year): 11 (2011)
Issue (Month): 2 (June)
|Note:||to access software from within Stata, net describe http://www.stata-journal.com/software/sj11-2/st0228/|
|Contact details of provider:|| Web page: http://www.stata-journal.com/|
|Order Information:||Web: http://www.stata-journal.com/subscription.html|
References listed on IDEAS
Please report citation or reference errors to , or , if you are the registered author of the cited work, log in to your RePEc Author Service profile, click on "citations" and make appropriate adjustments.:
- Maurizio Pisati, 2004. "Simple thematic mapping," Stata Journal, StataCorp LP, vol. 4(4), pages 361-378, December.
- White, Laura Forsberg & Bonetti, Marco & Pagano, Marcello, 2009. "The choice of the number of bins for the M statistic," Computational Statistics & Data Analysis, Elsevier, vol. 53(10), pages 3640-3649, August.
- Scott Merryman, 2005. "USMAPS2: Stata module to provide US county map coordinates for tmap," Statistical Software Components S448404, Boston College Department of Economics, revised 24 Jan 2005.
When requesting a correction, please mention this item's handle: RePEc:tsj:stataj:v:11:y:2011:i:2:p:271-289. See general information about how to correct material in RePEc.
For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: (Christopher F. Baum)or (Lisa Gilmore)
If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.
If references are entirely missing, you can add them using this form.
If the full references list an item that is present in RePEc, but the system did not link to it, you can help with this form.
If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your profile, as there may be some citations waiting for confirmation.
Please note that corrections may take a couple of weeks to filter through the various RePEc services.