Description Usage Arguments Value References See Also Examples

This function tests whether covariates in a treatment group and a matched control group are balanced in observational studies through the minimum spanning tree constructed on the subjects.

1 |

`distM` |
The distance matrix for the pooled observations (pooled over the treated subjects and the matched controls). If there are n treated subjects and n matched controls, then this distance matrix is a 2n by 2n matrix with the [i,j] element the distance between observation i and observation j. What distance to use is decided by users. Some simple choices are the Euclidean distance, L1 distance, and mahalanobis distance. |

`treated.index` |
The subject indices of the treated subjects. The subjects are ordered in the same way as for calculating the distance matrix, distM. |

`perm` |
The number of permutations performed to calculate the p-value of the test. The default value is 0, which means the permutation is not performed and only approximate p-value based on asymptotic theory is provided. Doing permutation could be time consuming, so be cautious if you want to set this value to be larger than 10,000. |

`k` |
Set as positive integer values, indicates k-MST is used. |

`discrete.correction` |
When this is set as TRUE (recommended), a continuation correction is done for computing the asymptotic p-value to account for the discrete nature of the statistic. |

`test.stat.Z` |
The standardized test statistic (ZR in the reference paper. |

`pval.appr` |
The approximated p-value based on asymptotic theory. |

`pval.perm` |
The permutation p-value when argument 'perm' is positive. |

Chen, H. and Small, D. (2019) New multivariate tests for assessing covariate balance in matched observational studies.

1 2 3 4 5 6 7 8 | ```
## A snippet of the smoking example in the reference paper.
## smoking.rda contains a 300 by 300 distance matrix, smokingDist.
## The indices of the treated subjects are 1:150.
data(smoking)
CrossMST(smokingDist, 1:150)
## Uncomment the following line to get permutation p-value with 1,000 permutations.
# CrossMST(smokingDist, 1:150, perm=1000)
``` |

```
$test.stat.R
[1] 88
$test.stat.Z
[1] 2.379564
$pval.appr
[1] 0.01712261
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.