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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/plot_spectrum.R
\name{plot_spectrum}
\alias{plot_spectrum}
\title{Plot point mutation spectrum}
\usage{
plot_spectrum(
type_occurrences,
CT = FALSE,
by = NA,
indv_points = FALSE,
error_bars = c("95\%_CI", "stdev", "SEM", "none"),
colors = NA,
legend = TRUE,
condensed = FALSE
)
}
\arguments{
\item{type_occurrences}{Type occurrences matrix}
\item{CT}{Distinction between C>T at CpG and C>T at other
sites, default = FALSE}
\item{by}{Optional grouping variable}
\item{indv_points}{Whether to plot the individual samples
as points, default = FALSE}
\item{error_bars}{The type of error bars to plot.
* '95%_CI' for 95% Confidence intervals (default);
* 'stdev' for standard deviations;
* 'SEM' for the standard error of the mean (NOT recommended);
* 'none' Do not plot any error bars;}
\item{colors}{Optional color vector with 7 values}
\item{legend}{Plot legend, default = TRUE}
\item{condensed}{More condensed plotting format. Default = F.}
}
\value{
Spectrum plot
}
\description{
Plot point mutation spectrum
}
\examples{
## See the 'read_vcfs_as_granges()' example for how we obtained the
## following data:
vcfs <- readRDS(system.file("states/read_vcfs_as_granges_output.rds",
package = "MutationalPatterns"
))
## Load a reference genome.
ref_genome <- "BSgenome.Hsapiens.UCSC.hg19"
library(ref_genome, character.only = TRUE)
## Get the type occurrences for all VCF objects.
type_occurrences <- mut_type_occurrences(vcfs, ref_genome)
## Plot the point mutation spectrum over all samples
plot_spectrum(type_occurrences)
## Or with distinction of C>T at CpG sites
plot_spectrum(type_occurrences, CT = TRUE)
## You can also include individual sample points.
plot_spectrum(type_occurrences, CT = TRUE, indv_points = TRUE)
## You can also change the type of error bars
plot_spectrum(type_occurrences, error_bars = "stdev")
## Or plot spectrum per tissue
tissue <- c(
"colon", "colon", "colon",
"intestine", "intestine", "intestine",
"liver", "liver", "liver"
)
plot_spectrum(type_occurrences, by = tissue, CT = TRUE)
## Or plot the spectrum per sample. Error bars are set to 'none', because they can't be plotted.
plot_spectrum(type_occurrences, by = names(vcfs), CT = TRUE, error_bars = "none")
## Plot it in a more condensed manner,
## which is is ideal for publications.
plot_spectrum(type_occurrences,
by = names(vcfs),
CT = TRUE,
error_bars = "none",
condensed = TRUE)
## You can also set custom colors.
my_colors <- c(
"pink", "orange", "blue", "lightblue",
"green", "red", "purple"
)
## And use them in a plot.
plot_spectrum(type_occurrences,
CT = TRUE,
legend = TRUE,
colors = my_colors
)
}
\seealso{
\code{\link{read_vcfs_as_granges}},
\code{\link{mut_type_occurrences}}
}
|
