DawnNormalize: DawnNormalize

Description Usage Arguments Value Examples

View source: R/DawnNormalize.R

Description

DawnNormalize takes a tumor and normal expression matrix and returns a standardized absolute differential expression matrix

Usage

1
DawnNormalize(tumorMat, normalMat)

Arguments

tumorMat,

a matrix representing the tumor expression

normalMat,

a matrix representing the normal expression

Value

the differential expression matrix, it is absolute standarized matrix

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
###using a small subset of the TCGA dataset and a small KEGG gene interaction network,
###We will show how to get DawnRank Results

library(DawnRank)

#load the mutation data
data(brcaExampleMutation)

#load the tumor expression data
data(brcaExampleTumorExpression)

#load the normal expression data
data(brcaExampleNormalExpression)

#load the pathway data
data(brcaExamplePathway)

#load the gold standard
data(goldStandard)

#normalize the tumor and normal data to get the differential expression
normalizedDawn<-DawnNormalize(tumorMat=brcaExampleTumorExpression,normalMat=brcaExampleNormalExpression)

#get the DawnRank Score
dawnRankScore<-DawnRank(adjMatrix=brcaExamplePathway,mutationMatrix=brcaExampleMutation,expressionMatrix=normalizedDawn, mu=3,goldStandard=goldStandard)

#look at the DawnRank scores for a few patients
dawnRankFrame<-dawnRankScore[[3]]
head(dawnRankFrame)

#get the aggregate DawnRank scores
aggregateDawnRankScore<-condorcetRanking(scoreMatrix=dawnRankScore[[2]],mutationMatrix=brcaExampleMutation)

#look at top 10 ranked genes
top10<-aggregateDawnRankScore[[2]][1:10]
top10

#get the individual cutoff for patient TCGA-A2-A04P
dawnRankFrame$isCGC<-dawnRankFrame$isGoldStandard
library(maxstat)
patspeccutoff(patient="TCGA-A2-A04P",ms=dawnRankFrame,default=95)->significance

MartinFXP/DawnRank documentation built on Sept. 10, 2021, 11:25 p.m.