R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(RTCGA)
Welcome to the RTCGA (version: 1.2.2).
> png(filename="/home/ddbj/snapshot/RGM3/R_BC/result/RTCGA/kmTCGA.Rd_%03d_medium.png", width=480, height=480)
> ### Name: kmTCGA
> ### Title: Plot Kaplan-Meier Estimates of Survival Curves for Survival Data
> ### Aliases: kmTCGA
>
> ### ** Examples
>
>
> ## Extracting Survival Data
> library(RTCGA.clinical)
> survivalTCGA(BRCA.clinical, OV.clinical, extract.cols = "admin.disease_code") -> BRCAOV.survInfo
>
> # first munge data, then extract survival info
> library(dplyr)
Attaching package: 'dplyr'
The following objects are masked from 'package:stats':
filter, lag
The following objects are masked from 'package:base':
intersect, setdiff, setequal, union
> BRCA.clinical %>%
+ filter(patient.drugs.drug.therapy_types.therapy_type %in%
+ c("chemotherapy", "hormone therapy")) %>%
+ rename(therapy = patient.drugs.drug.therapy_types.therapy_type) %>%
+ survivalTCGA(extract.cols = c("therapy")) -> BRCA.survInfo.chemo
>
> # first extract survival info, then munge data
> survivalTCGA(BRCA.clinical,
+ extract.cols = c("patient.drugs.drug.therapy_types.therapy_type")) %>%
+ filter(patient.drugs.drug.therapy_types.therapy_type %in%
+ c("chemotherapy", "hormone therapy")) %>%
+ rename(therapy = patient.drugs.drug.therapy_types.therapy_type) -> BRCA.survInfo.chemo
>
> ## Kaplan-Meier Survival Curves
> kmTCGA(BRCAOV.survInfo, explanatory.names = "admin.disease_code", pval = TRUE)
>
> kmTCGA(BRCAOV.survInfo, explanatory.names = "admin.disease_code", main = "",
+ xlim = c(0,4000))
>
> kmTCGA(BRCA.survInfo.chemo, explanatory.names = "therapy", xlim = c(0, 3000), conf.int = FALSE)
>
>
>
>
>
>
> dev.off()
null device
1
>