data pulled: 2019-08-26
written: 2019-08-26
last ran: 2019-08-26
website: http://rpubs.com/navona/RISH_matching_MRP


General description

This code is used to identify matched controls, as required to test an implementation of the RISH harmonization method described in Karayumak et al., 2019. We use data from the SPINS study, and match 5 target sites (CMP, MRC, MRP, ZHH, ZHP) to the reference site (CMH).

In all cases, we aim to identify the gold-standard of n=20 matched ‘healthy control’ participants (or at least, the recommended n=16-18), on the basis of sex, age,IQ (here, established by WTAR), and handedness.


Data count

reference
targets
CMH CMP MRC MRP ZHH ZHP
female 16 6 9 10 5 17
male 23 3 13 11 7 17


Before matching


CMH MRP p test
n 39 21
age (mean (SD)) 27.08 (8.34) 34.29 (11.79) 0.008
sex = male (%) 23 (59.0) 11 (52.4) 0.827
wtar_std (mean (SD)) 115.33 (9.39) 109.19 (11.28) 0.028
males
CMH MRP p test
n 23 11
age (mean (SD)) 26.78 (8.23) 31.73 (11.95) 0.167
wtar_std (mean (SD)) 115.35 (9.20) 107.27 (13.20) 0.046
females
CMH MRP p test
n 16 10
age (mean (SD)) 27.50 (8.75) 37.10 (11.54) 0.024
wtar_std (mean (SD)) 115.31 (9.95) 111.30 (8.92) 0.309


After matching


CMH MRP p test
n 20 20
age (mean (SD)) 29.85 (9.24) 34.75 (11.89) 0.154
sex = male (%) 10 (50.0) 10 (50.0) 1.000
wtar_std (mean (SD)) 112.35 (9.49) 110.30 (10.33) 0.517
distance (mean (SD)) 0.37 (0.18) 0.46 (0.22) 0.179


List of matched pairs


males

SPN01_CMH_0044 SPN01_MRP_0088
SPN01_CMH_0008 SPN01_MRP_0162
SPN01_CMH_0040 SPN01_MRP_0100
SPN01_CMH_0152 SPN01_MRP_0152
SPN01_CMH_0113 SPN01_MRP_0144
SPN01_CMH_0030 SPN01_MRP_0095
SPN01_CMH_0073 SPN01_MRP_0136
SPN01_CMH_0123 SPN01_MRP_0125
SPN01_CMH_0065 SPN01_MRP_0157
SPN01_CMH_0031 SPN01_MRP_0102

females

SPN01_CMH_0066 SPN01_MRP_0122
SPN01_CMH_0007 SPN01_MRP_0139
SPN01_CMH_0023 SPN01_MRP_0082
SPN01_CMH_0134 SPN01_MRP_0101
SPN01_CMH_0005 SPN01_MRP_0081
SPN01_CMH_0026 SPN01_MRP_0123
SPN01_CMH_0016 SPN01_MRP_0076
SPN01_CMH_0194 SPN01_MRP_0156
SPN01_CMH_0020 SPN01_MRP_0128
SPN01_CMH_0125 SPN01_MRP_0127

Special notes

Because we could not match on the basis of handedness (established by EHI ≥ .5) for another site, we did not attempt to match on the basis of handedness here (nor will we for any other target site).