mutationsTCGA: Gather Mutations for TCGA Datasets

Description Usage Arguments Issues Note Author(s) See Also Examples

View source: R/mutationsTCGA.R

Description

Function gathers mutations over multiple TCGA datasets and extracts mutations and further informations about them for desired genes. See mutations.

Usage

1
2
mutationsTCGA(..., extract.cols = c("Hugo_Symbol", "Variant_Classification",
  "bcr_patient_barcode"), extract.names = TRUE, unique = TRUE)

Arguments

...

A data.frame or data.frames from TCGA study containing mutations information (RTCGA.mutations).

extract.cols

A character specifing the names of columns to be extracted with bcr_patient_barcode. If NULL all columns are returned.

extract.names

Logical, whether to extract names of passed data.frames in ....

unique

Should the outputed data be unique. By default it's TRUE.

Issues

If you have any problems, issues or think that something is missing or is not clear please post an issue on https://github.com/RTCGA/RTCGA/issues.

Note

Input data.frames should contain column bcr_patient_barcode if extract.cols is specified.

Author(s)

Marcin Kosinski, m.p.kosinski@gmail.com

See Also

RTCGA website http://rtcga.github.io/RTCGA/Visualizations.html.

Other RTCGA: RTCGA-package, boxplotTCGA, checkTCGA, convertTCGA, datasetsTCGA, downloadTCGA, expressionsTCGA, heatmapTCGA, infoTCGA, installTCGA, kmTCGA, pcaTCGA, readTCGA, survivalTCGA, theme_RTCGA

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
library(RTCGA.mutations)
library(dplyr)
mutationsTCGA(BRCA.mutations, OV.mutations) %>%
	filter(Hugo_Symbol == 'TP53') %>%
	filter(substr(bcr_patient_barcode, 14, 15) == "01") %>% # cancer tissue
	mutate(bcr_patient_barcode = substr(bcr_patient_barcode, 1, 12)) -> BRCA_OV.mutations

library(RTCGA.clinical)
survivalTCGA(BRCA.clinical, OV.clinical, extract.cols = "admin.disease_code") %>%
	rename(disease = admin.disease_code)-> BRCA_OV.clinical

BRCA_OV.clinical %>%
	left_join(BRCA_OV.mutations,
	by = "bcr_patient_barcode") %>%
	mutate(TP53 = ifelse(!is.na(Variant_Classification), "Mut",
 "WILDorNOINFO")) -> BRCA_OV.clinical_mutations

BRCA_OV.clinical_mutations %>%
	select(times, patient.vital_status, disease, TP53) -> BRCA_OV.2plot
kmTCGA(BRCA_OV.2plot, explanatory.names = c("TP53", "disease"),
			 break.time.by = 400, xlim = c(0,2000))

RTCGA documentation built on Nov. 8, 2020, 5:11 p.m.