find_node_perturbation {ceRNAnetsim} | R Documentation

## Calculates average expression changes of all (or specified) nodes except trigger and finds the perturbed node count for all (or specified) nodes in system.

### Description

Calculates average expression changes of all (or specified) nodes except trigger and finds the perturbed node count for all (or specified) nodes in system.

### Usage

    find_node_perturbation(input_graph, how = 2, cycle = 1, limit = 0, fast = 0)

### Arguments

`input_graph` | The graph object that was processed with priming_graph function.
`how` | The change of count (expression) of the given node in terms of fold change.
`cycle` | The iteration of simulation.
`limit` | The minimum fold change which can be taken into account for perturbation calculation on all nodes in terms of percentage.
`fast` | specifies percentage of affected target in target expression. For example, if fast = 1, the nodes that are affected from miRNA repression activity more than one percent of their expression is determined as subgraph.

### Details

find_node_perturbation calculates mean expression changes of elements after the change in the network in terms of percentage. It also calculates the number of nodes that have expression changes after the change occur in the network. The outputs of the function are the perturbation efficiency and perturbed count of nodes for each nodes.

### Value

It gives a tibble form dataset that includes node names, perturbation efficiency and perturbed count of nodes.

### Examples

    data('minsamp')
    data('midsamp')

     minsamp%>%
      priming_graph(competing_count = Competing_expression, miRNA_count = miRNA_expression)%>%
      find_node_perturbation()%>%
      select(name, perturbation_efficiency, perturbed_count)

     minsamp%>%
      priming_graph(competing_count = Competing_expression, miRNA_count = miRNA_expression,
        aff_factor = c(energy,seed_type), deg_factor = region)%>%
      find_node_perturbation(how = 3, cycle = 4)%>%
      select(name, perturbation_efficiency, perturbed_count)

     midsamp%>%
      priming_graph(competing_count = Gene_expression, miRNA_count = miRNA_expression)%>%
      find_node_perturbation(how = 2, cycle= 3, limit=1, fast = 5)%>%
      select(name, perturbation_efficiency, perturbed_count)

[Package ceRNAnetsim version 1.0.0 Index]