iDOS-package: Integrated Discovery of Oncogenic Signatures

Description Details Author(s) Examples

Description

A method to identify correlated changes on mRNA and DNA level. For details, see PMID: 27358048

Details

Package: iDOS
Type: Package
Version: 0.0.1
Date: 2014-11-14
License: GPL-2

Author(s)

Syed Haider, Francesca M Buffa
Maintainer: Syed Haider <Syed.Haider@oncology.ox.ac.uk>
MRC Weatherall Institute of Molecular Medicine
Department of Oncology
University of Oxford.

Examples

 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
36
37
38
39
40
41
42
43
# load test data
x <- get.test.data(data.types = c("mRNA.T", "mRNA.N", "CNA", "ann"));

# list of features to be assessed for differential expression
feature.ids <- rownames(x$mRNA.T$BLCA);

# get differentially expressed features
DE.results <- find.DE.features(
  exp.data.T = x$mRNA.T, 
  exp.data.N = x$mRNA.N, 
  feature.ids = feature.ids,
  test.name = "t.test"
  );

# get top features
top.features <- get.top.features(
  DE.features = cbind("FC" = DE.results[, 1], "P" = DE.results[, 2]),
  cna.data.fractions = x$CNA.fractions$BLCA, 
  mRNA.FC.up = 0.25, 
  mRNA.FC.down = 0.25, 
  mRNA.p = 0.05, 
  mRNA.top.n = NULL, 
  cna.fractions.gain = 0.2, 
  cna.fractions.loss = 0.2
  );

# temporary output directory
tmp.output.dir <- tempdir();

# estimate mRNA and CNA correlation using the pre-selected top features
correlated.features <- estimate.expression.cna.correlation(
  exp.data = x$mRNA.T$BLCA, 
  cna.data.log2 = x$CNA.log2$BLCA, 
  corr.threshold = 0.3, 
  corr.direction = "two.sided", 
  subtypes.metadata = list(
    "subtype.samples.list" = list("All" = colnames(x$mRNA.T$BLCA))
    ), 
  feature.ids = top.features, 
  cancer.type = "BLCA", 
  data.dir = paste(tmp.output.dir, "/data/BLCA/", sep = ""),
  graphs.dir = paste(tmp.output.dir, "/graphs/BLCA/", sep = "")
  );

iDOS documentation built on May 2, 2019, 2:28 p.m.