# run_mrnet: Wrapper for MRNET method In dnapath: Differential Network Analysis using Gene Pathways

## Description

Conducts co-expression analysis using MRNET \insertCitemeyer07dnapath. Uses the implementation from the `minet` package \insertCitemeyer08dnapath. Can be used for the `network_inference` argument in `dnapath`.

## Usage

 `1` ```run_mrnet(x, estimator = "spearman", ...) ```

## Arguments

 `x` A n by p matrix of gene expression data (n samples and p genes). `estimator` Argument is passed into `build.mim`. `...` Additional arguments are ignored.

## Value

A p by p matrix of association scores.

## References

\insertRef

meyer07dnapath

\insertRef

meyer08dnapath

`run_aracne`, `run_bc3net`, `run_c3net`, `run_clr`, `run_corr`, `run_dwlasso`, `run_genie3`, `run_glasso`, `run_pcor`, and `run_silencer`
 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35``` ```data(meso) data(p53_pathways) # To create a short example, we subset on two pathways from the p53 pathway list, # and will only run 3 permutations for significance testing. pathway_list <- p53_pathways[c(8, 13)] n_perm <- 3 # Use this method to perform differential network analysis. # The parameters in run_mrnet() can be adjusted using the ... argument. # For example, the 'estimator' parameter can be specified as shown here. results <- dnapath(x = meso\$gene_expression, pathway_list = pathway_list, groups = meso\$groups, n_perm = n_perm, network_inference = run_mrnet, estimator = "spearman") summary(results) # The group-specific association matrices can be extracted using get_networks(). nw_list <- get_networks(results[[1]]) # Get networks for pathway 1. # nw_list has length 2 and contains the inferred networks for the two groups. # The gene names are the Entrezgene IDs from the original expression dataset. # Renaming the genes in the dnapath results to rename those in the networks. # NOTE: The temporary directory, tempdir(), is used in this example. In practice, # this argument can be removed or changed to an existing directory results <- rename_genes(results, to = "symbol", species = "human", dir_save = tempdir()) nw_list <- get_networks(results[[1]]) # The genes (columns) will have new names. # (Optional) Plot the network using SeqNet package (based on igraph plotting). # First rename entrezgene IDs into gene symbols. SeqNet::plot_network(nw_list[[1]]) ```