MAP.Matches(data, varname, t.cutoff = "98.00%", multiple = TRUE, perm = c("both", "columns", "labels")[1], nperm = 1000, test = c("t", "t.equalvar")[1], sig.col, sig.cutoff = 0.05)MetaArrayTRUE only paterrns with multiple '1' are used"labels" only class labels are permuted for statistical analysis (empirical significance), if "columns" only genes in each dataset are selected randomly, if "both" both class labels and genes are permuted and two p-values returned"t" then unequal variance t-test is used, if "t.equalvar" equal variance t-test is used"p.col.strong", "p.col.weak", "p.lab.strong", "p.lab.weak" , "col" refers to column permutations, "lab" to class labels, "weak" to soft match and "strong" to strong matchMAP.Matches.res containing
tests, 1 means the test statistics was higer than thresholdbin.matrix: number of selected genes in each dataset, genes with at least one 1 in pattern, probability of observing strong or soft match in the data n.strong and soft n.soft matches and number of genes involved n.sigdata(ColonData)
MAP.Matches(ColonData, "MSI", nperm = 100, sig.col="p.lab.strong")
Run the code above in your browser using DataLab