PlotGeneNetworks
PlotGeneNetworks(Object = IPA.object3, labels=T, RRmin = 15)
Object a PathwayObject
labels Boolean show labels or not, or a vector with labels to be shown.
RRmin filter for different RR cutoff for what is an edge.
require(GeneSetCluster)
IPA.files <- c(system.file("extdata", "MM10.IPA.KO.uGvsMac.Canonical_pathways.xls", package = "GeneSetCluster"),
system.file("extdata", "MM10.IPA.WT.uGvsMac.Canonical_pathways.xls", package = "GeneSetCluster"),
system.file("extdata", "MM10.IPA.KO.uGvsMac.Functional_annotations.xls", package = "GeneSetCluster"),
system.file("extdata", "MM10.IPA.WT.uGvsMac.Functional_annotations.xls", package = "GeneSetCluster"))
canonical.files <- IPA.files[grep("Canonical", IPA.files)]
IPA.object1 <- LoadGeneSets(file_location = canonical.files, #where are the files
groupnames= c("KO", "WT"),#Names of the groups
P.cutoff = 1.3, #minumum cutoff if smaller than 1 it assumes normal pvalue, if larger than 1 it assumes a log10 palue
Mol.cutoff = 5,# amount of molecules interested in
Source = "IPA",#How was the data generated
type = "Canonical_Pathways",#What is the experiment e.g. canonical pathways, functional anotation
structure = "SYMBOL",#structure of the molecules e.g. genenames, ensembl_ID etc
seperator = ",")#How are the genes seperated
IPA.object2 <- CombineGeneSets(Object = IPA.object1)
IPA.object3 <- ClusterGeneSets(Object = IPA.object2,
clusters = 7,
method = "kmeans")
PlotGeneNetworks(Object = IPA.object3, labels=F, RRmin = 0)
PlotGeneNetworks(Object = IPA.object3, labels=T, RRmin = 0)
PlotGeneNetworks(Object = IPA.object3, labels=F, RRmin = 15)
PlotGeneNetworks(Object = IPA.object3, labels=T, RRmin = 15)
labels.idx <- c("Neuroinflammation Signaling Pathway", "Phagosome Maturation", "IL-4 Signaling")
PlotGeneNetworks(Object = IPA.object3, labels=labels.idx, RRmin = 15)
