MUCH_merge {mergingTools}

MUCH Merge

Description

MUCH_merge merges HEMs coming from separate subexperiments into one single dataset. A Subexperiment is a measurement of a group of HEMs the size of the allowed PMCs by the machine. For MUCH you need to measure each HEM against the rest of them. For instance, let's say that one has 2 PMCs and 3 HEMs to measure. The input data coming from the experiments should look like this:

1	2	1	3	2	3
30	15	20	54	15	24
35	16	25	32	10	29
32	14	30	45	9	32

where the numbers on top are the codes for the HEMs on the T2080. In this case we have three subexperiments: subexp1 = (1, 2), subexp2 = (1, 3, 6), subexp3 = (2, 3); The data will be processed into:

H1	H2	H1	H3	H2	H3
30	15	20	54	15	24
35	16	25	32	10	29
32	14	30	45	9	32

The processing transforms the code of the HEMs to the reference name on the T2080 manual. Now with this data, MUCH_merge computes the correlation matrix of all HEMs and with it it constructs a multivariate Gaussian distribution (MVG). Then MUCH_merge uses the order statistics of the MVG to arrange the experimental data. Therefore the final input will look like this:

H1	H2	H3
30	15	45
35	16	54
32	14	32

Usage

MUCH_merge(splitted_data = NULL, n_runs = NULL, n_sims = NULL, dep_lvl = NULL)

Arguments

splitted_data

Sample data as list of dataframes.

n_runs

Number of rows for the output merged data

n_sims

Number of simulations for the multivariate Gaussian distribution to find the optimal merge

dep_lvl

Dependency level allowed. A number between 0 and 1 indicating the maximum correlation allowed between HEMs.

Details

Take into consideration that the input data has readings on the same HEM for multiple subexperiments. Therefore one must select the number of runs n_runs that they want the final output to have.

Value

Returns a merged dataframe as stated in the description

Examples

n_pmcs <- 6
data_much <- mergingTools::process_raw_experiments(data = data_much_raw_vignette,
                                                   n_pmcs = n_pmcs)
merged_data <- mergingTools::MUCH_merge(splitted_data = data_much,
                                        n_runs = 1000,
                                        n_sims = 10,
                                        dep_lvl = 0.85)