TCGAWorkflow 1.3.0
Tiago Chedraoui Silva
| Snapshot Date: 2018-05-18 07:35:07 -0400 (Fri, 18 May 2018) |
| URL: https://git.bioconductor.org/packages/TCGAWorkflow |
| Branch: master |
| Last Commit: d7225d8 |
| Last Changed Date: 2018-04-30 10:10:33 -0400 (Mon, 30 Apr 2018) |
| malbec1 | Linux (Ubuntu 16.04.1 LTS) / x86_64 | OK | OK |  |
##############################################################################
##############################################################################
###
### Running command:
###
### /home/biocbuild/bbs-3.8-bioc/R/bin/R -q -e 'rmarkdown::render(".buildwebvig/TCGAWorkflow/TCGAWorkflow.Rmd", output_format="BiocStyle:::html_fragment")' && /home/biocbuild/bbs-3.8-bioc/R/bin/R -q -e 'knitr::purl(".buildwebvig/TCGAWorkflow/TCGAWorkflow.Rmd", ".buildwebvig/TCGAWorkflow/TCGAWorkflow.R")'
###
##############################################################################
##############################################################################
> rmarkdown::render(".buildwebvig/TCGAWorkflow/TCGAWorkflow.Rmd", output_format="BiocStyle:::html_fragment")
processing file: TCGAWorkflow.Rmd
|
| | 0%
|
| | 1%
inline R code fragments
|
|. | 1%
label: setup (with options)
List of 4
$ include: logi FALSE
$ echo : symbol F
$ warning: symbol F
$ message: symbol F
|
|. | 2%
ordinary text without R code
|
|.. | 3%
label: unnamed-chunk-1 (with options)
List of 4
$ echo : logi FALSE
$ results: chr "hide"
$ warning: logi FALSE
$ message: logi FALSE
inline R code fragments
|
|... | 4%
label: unnamed-chunk-2 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
ordinary text without R code
|
|.... | 5%
label: unnamed-chunk-3 (with options)
List of 4
$ eval : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
$ include: logi TRUE
|
|.... | 6%
ordinary text without R code
label: unnamed-chunk-4 (with options)
List of 4
$ eval : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
$ include: logi FALSE
|
|..... | 7%
ordinary text without R code
|
|..... | 8%
label: unnamed-chunk-5 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|...... | 8%
ordinary text without R code
|
|...... | 9%
label: unnamed-chunk-6 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|....... | 10%
ordinary text without R code
label: unnamed-chunk-7 (with options)
List of 2
$ eval : logi TRUE
$ include: logi TRUE
Loading required package: GenomicRanges
Loading required package: stats4
Loading required package: S4Vectors
Attaching package: 'S4Vectors'
The following object is masked from 'package:base':
expand.grid
Loading required package: IRanges
Loading required package: GenomeInfoDb
Loading required package: DelayedArray
Loading required package: matrixStats
Attaching package: 'matrixStats'
The following objects are masked from 'package:Biobase':
anyMissing, rowMedians
Loading required package: BiocParallel
Attaching package: 'DelayedArray'
The following objects are masked from 'package:matrixStats':
colMaxs, colMins, colRanges, rowMaxs, rowMins, rowRanges
The following objects are masked from 'package:base':
aperm, apply
|
|........ | 11%
ordinary text without R code
label: unnamed-chunk-8 (with options)
List of 2
$ eval : logi TRUE
$ include: logi TRUE
|
|........ | 12%
label: unnamed-chunk-9 (with options)
List of 3
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|......... | 13%
ordinary text without R code
label: unnamed-chunk-10 (with options)
List of 4
$ results: chr "hide"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
--------------------------------------
o GDCquery: Searching in GDC database
--------------------------------------
Genome of reference: hg38
--------------------------------------------
oo Accessing GDC. This might take a while...
--------------------------------------------
ooo Project: TCGA-GBM
--------------------
oo Filtering results
--------------------
ooo By barcode
----------------
oo Checking data
----------------
ooo Check if there are duplicated cases
ooo Check if there results for the query
-------------------
o Preparing output
-------------------
Downloading data for project TCGA-GBM
GDCdownload will download 2 files. A total of 60.227 KB
Downloading as: Fri_May_18_08_54_59_2018.tar.gz
To get the following information please change the clinical.info argument
=> new_tumor_events: new_tumor_event
=> drugs: drug
=> follow_ups: follow_up
=> radiations: radiation
Parsing follow up version: follow_up_v1.0
|
|.......... | 14%
label: unnamed-chunk-11 (with options)
List of 3
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|.......... | 15%
label: unnamed-chunk-12 (with options)
List of 4
$ results: chr "hide"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|........... | 15%
label: unnamed-chunk-13 (with options)
List of 3
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|........... | 16%
label: unnamed-chunk-14 (with options)
List of 4
$ results: chr "hide"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|............ | 17%
label: unnamed-chunk-15 (with options)
List of 3
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
label: unnamed-chunk-16 (with options)
List of 4
$ results: chr "hide"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|............ | 18%
label: unnamed-chunk-17 (with options)
List of 3
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|............. | 18%
ordinary text without R code
|
|............. | 19%
label: unnamed-chunk-18 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|.............. | 20%
ordinary text without R code
label: unnamed-chunk-19 (with options)
List of 3
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|............... | 21%
ordinary text without R code
|
|............... | 22%
label: unnamed-chunk-20 (with options)
List of 2
$ eval : logi TRUE
$ include: logi TRUE
gbm subtype information from:doi:10.1016/j.cell.2015.12.028
|
|................ | 22%
ordinary text without R code
|
|................ | 23%
label: unnamed-chunk-21 (with options)
List of 3
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|................ | 24%
ordinary text without R code
|
|................. | 24%
label: unnamed-chunk-22 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|................. | 25%
ordinary text without R code
|
|.................. | 25%
label: unnamed-chunk-23 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|.................. | 26%
ordinary text without R code
|
|................... | 27%
label: unnamed-chunk-24 (with options)
List of 4
$ eval : logi FALSE
$ message: logi FALSE
$ warning: logi FALSE
$ include: logi TRUE
ordinary text without R code
|
|.................... | 28%
label: unnamed-chunk-25 (with options)
List of 4
$ eval : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
$ include: logi TRUE
|
|.................... | 29%
ordinary text without R code
|
|..................... | 29%
label: unnamed-chunk-26 (with options)
List of 5
$ results: chr "hide"
$ eval : logi FALSE
$ message: logi FALSE
$ warning: logi FALSE
$ include: logi TRUE
|
|..................... | 30%
ordinary text without R code
|
|..................... | 31%
label: unnamed-chunk-27 (with options)
List of 5
$ results: chr "hide"
$ eval : logi FALSE
$ message: logi FALSE
$ warning: logi FALSE
$ include: logi TRUE
|
|...................... | 31%
ordinary text without R code
|
|...................... | 32%
label: unnamed-chunk-28 (with options)
List of 4
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
$ include: logi TRUE
Starting GBM
trying URL 'ftp://ftp.broadinstitute.org/pub/GISTIC2.0/hg19_support/CNV.hg19.bypos.111213.txt'
Content type 'unknown' length 358469 bytes (350 KB)
==================================================
Parsed with column specification:
cols(
ID = col_character(),
Chromosome = col_integer(),
Start = col_integer(),
End = col_integer(),
Flankingstart = col_integer(),
Flankingend = col_integer()
)
Loading Marker Informations
........................
Done
Loading Copy Number Data
................................................
Done
Performing Data Preprocessing
Done
Computing Discontinuity Matrix
.
Done
Computing Probability Distribution
....
Done
Assessing the Significance of Observed Data
..
Done
Writing GAIA_GBM_flt.txt.igv.gistic File for Integrative Genomics Viewer (IGV) Tool
..
Done
Running Homogeneous peel-off Algorithm With Significance Threshold of 0.25 and Homogenous Threshold of 0.12
..
Done
Writing Output File 'GAIA_GBM_flt.txt' Containing the Significant Regions
File 'GAIA_GBM_flt.txt' Saved
|
|....................... | 32%
ordinary text without R code
|
|....................... | 33%
label: unnamed-chunk-29 (with options)
List of 4
$ eval : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
$ include: logi FALSE
|
|........................ | 34%
ordinary text without R code
label: unnamed-chunk-30 (with options)
List of 4
$ eval : logi FALSE
$ message: logi FALSE
$ warning: logi FALSE
$ include: logi TRUE
|
|......................... | 35%
ordinary text without R code
|
|......................... | 36%
label: unnamed-chunk-31 (with options)
List of 4
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
$ include: logi TRUE
label: unnamed-chunk-32 (with options)
List of 4
$ results: chr "asis"
$ echo : logi FALSE
$ message: logi FALSE
$ warning: logi FALSE
|
|.......................... | 37%
ordinary text without R code
|
|.......................... | 38%
label: unnamed-chunk-33 (with options)
List of 5
$ results: chr "hide"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
$ eval : logi FALSE
|
|........................... | 38%
ordinary text without R code
|
|........................... | 39%
label: unnamed-chunk-34 (with options)
List of 4
$ results: chr "hide"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
silent variants: 38433
Summarizing..
Gene Summary..
NOTE: Possible FLAGS among top ten genes:
Checking clinical data..
Annotation missing for below samples in MAF
Done !
|
|............................ | 39%
ordinary text without R code
|
|............................ | 40%
label: unnamed-chunk-35 (with options)
List of 4
$ echo : logi TRUE
$ message : logi FALSE
$ warning : logi FALSE
$ fig.width: num 10
|
|............................. | 41%
ordinary text without R code
label: unnamed-chunk-36 (with options)
List of 5
$ echo : logi TRUE
$ message : logi FALSE
$ warning : logi FALSE
$ fig.height: num 10
$ fig.width : num 15
gbm subtype information from:doi:10.1016/j.cell.2015.12.028
lgg subtype information from:doi:10.1016/j.cell.2015.12.028
|
|............................. | 42%
ordinary text without R code
|
|.............................. | 43%
label: unnamed-chunk-37 (with options)
List of 3
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
Looking for clinical data in annoatation slot of MAF..
Number of mutated samples for given genes:
Median survival..
ordinary text without R code
|
|............................... | 44%
label: unnamed-chunk-38 (with options)
List of 5
$ results: chr "hide"
$ eval : logi FALSE
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|............................... | 45%
ordinary text without R code
|
|................................ | 45%
label: unnamed-chunk-39 (with options)
List of 4
$ results: chr "hide"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
========================================
circlize version 0.4.3
CRAN page: https://cran.r-project.org/package=circlize
Github page: https://github.com/jokergoo/circlize
Documentation: http://jokergoo.github.io/circlize_book/book/
If you use it in published research, please cite:
Gu, Z. circlize implements and enhances circular visualization
in R. Bioinformatics 2014.
========================================
Note: 2 points are out of plotting region in sector 'chr17', track '4'.
|
|................................ | 46%
ordinary text without R code
|
|................................. | 46%
label: unnamed-chunk-40 (with options)
List of 4
$ results: chr "asis"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
Note: 2 points are out of plotting region in sector 'chr17', track '4'.
Note: 2 points are out of plotting region in sector 'chr17', track '5'.
Note: 2 points are out of plotting region in sector 'chr17', track '5'.
Note: 1 point is out of plotting region in sector 'chr17', track '5'.
Note: 1 point is out of plotting region in sector 'chr17', track '5'.
|
|................................. | 47%
ordinary text without R code
|
|................................. | 48%
label: unnamed-chunk-41 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|.................................. | 48%
ordinary text without R code
|
|.................................. | 49%
label: unnamed-chunk-42 (with options)
List of 4
$ results: chr "asis"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
I Need about 11 seconds for this Complete Normalization Upper Quantile [Processing 80k elements /s]
Step 1 of 4: newSeqExpressionSet ...
Step 2 of 4: withinLaneNormalization ...
Step 3 of 4: betweenLaneNormalization ...
Step 4 of 4: .quantileNormalization ...
Batch correction skipped since no factors provided
----------------------- DEA -------------------------------
there are Cond1 type LGG in 20 samples
there are Cond2 type GBM in 20 samples
there are 14477 features as miRNA or genes
I Need about 19 seconds for this DEA. [Processing 30k elements /s]
----------------------- END DEA -------------------------------
|
|................................... | 50%
ordinary text without R code
label: unnamed-chunk-43 (with options)
List of 4
$ results: chr "asis"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|.................................... | 51%
ordinary text without R code
|
|.................................... | 52%
label: unnamed-chunk-44 (with options)
List of 6
$ results : chr "hide"
$ echo : logi TRUE
$ message : logi FALSE
$ warning : logi FALSE
$ fig.height: num 10
$ fig.width : num 10
|
|..................................... | 52%
label: unnamed-chunk-45 (with options)
List of 6
$ results : chr "hide"
$ echo : logi TRUE
$ message : logi FALSE
$ warning : logi FALSE
$ fig.height: num 10
$ fig.width : num 10
|
|..................................... | 53%
ordinary text without R code
|
|..................................... | 54%
label: unnamed-chunk-46 (with options)
List of 6
$ results : chr "hide"
$ echo : logi TRUE
$ message : logi FALSE
$ warning : logi FALSE
$ fig.height: num 10
$ fig.width : num 10
|
|...................................... | 54%
ordinary text without R code
|
|...................................... | 55%
label: unnamed-chunk-47 (with options)
List of 6
$ results : chr "hide"
$ echo : logi TRUE
$ message : logi FALSE
$ warning : logi FALSE
$ fig.height: num 12
$ fig.width : num 15
|
|....................................... | 55%
ordinary text without R code
|
|....................................... | 56%
label: unnamed-chunk-48 (with options)
List of 4
$ results: chr "asis"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
clusterProfiler v3.9.1 For help: https://guangchuangyu.github.io/software/clusterProfiler
If you use clusterProfiler in published research, please cite:
Guangchuang Yu, Li-Gen Wang, Yanyan Han, Qing-Yu He. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS: A Journal of Integrative Biology. 2012, 16(5):284-287.
Attaching package: 'clusterProfiler'
The following object is masked from 'package:DelayedArray':
simplify
Loading required package: org.Hs.eg.db
Loading required package: AnnotationDbi
'select()' returned 1:many mapping between keys and columns
Info: Downloading xml files for hsa05214, 1/1 pathways..
Info: Downloading png files for hsa05214, 1/1 pathways..
'select()' returned 1:1 mapping between keys and columns
Info: Working in directory /home/biocbuild/bbs-3.8-workflows/meat/.buildwebvig/TCGAWorkflow
Info: Writing image file hsa05214.pathview.png
Info: some node width is different from others, and hence adjusted!
|
|........................................ | 57%
ordinary text without R code
label: unnamed-chunk-49 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|......................................... | 58%
ordinary text without R code
|
|......................................... | 59%
label: unnamed-chunk-50 (with options)
List of 2
$ eval : logi TRUE
$ include: logi TRUE
Loading required package: igraph
Attaching package: 'igraph'
The following object is masked from 'package:clusterProfiler':
simplify
The following object is masked from 'package:circlize':
degree
The following objects are masked from 'package:DelayedArray':
path, simplify
The following object is masked from 'package:GenomicRanges':
union
The following object is masked from 'package:IRanges':
union
The following object is masked from 'package:S4Vectors':
union
The following objects are masked from 'package:BiocGenerics':
normalize, path, union
The following objects are masked from 'package:stats':
decompose, spectrum
The following object is masked from 'package:base':
union
ordinary text without R code
|
|.......................................... | 60%
label: unnamed-chunk-51 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|.......................................... | 61%
ordinary text without R code
|
|........................................... | 61%
label: unnamed-chunk-52 (with options)
List of 5
$ echo : logi TRUE
$ message : logi FALSE
$ warning : logi FALSE
$ fig.height: num 8
$ fig.width : num 8
==================== DATA Summary ====================
==================== END DATA Summary ====================
==================== T test results ====================
==================== END T test results ====================
|
|........................................... | 62%
ordinary text without R code
|
|............................................ | 62%
label: unnamed-chunk-53 (with options)
List of 4
$ results: chr "hide"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
Group1:TCGA-GBM
Group2:TCGA-LGG
Calculating the diference between the mean methylation of the groups...
Calculating the p-values of each probe...
|
|............................................ | 63%
ordinary text without R code
|
|............................................. | 64%
label: unnamed-chunk-54 (with options)
List of 4
$ results: chr "hide"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
Loading required package: grid
========================================
ComplexHeatmap version 1.19.0
Bioconductor page: http://bioconductor.org/packages/ComplexHeatmap/
Github page: https://github.com/jokergoo/ComplexHeatmap
Documentation: http://bioconductor.org/packages/ComplexHeatmap/
If you use it in published research, please cite:
Gu, Z. Complex heatmaps reveal patterns and correlations in multidimensional
genomic data. Bioinformatics 2016.
========================================
ordinary text without R code
|
|............................................. | 65%
label: unnamed-chunk-55 (with options)
List of 4
$ results: chr "asis"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
Loading required package: Biostrings
Loading required package: XVector
Attaching package: 'Biostrings'
The following object is masked from 'package:DelayedArray':
type
The following object is masked from 'package:base':
strsplit
Loading required package: BSgenome
Loading required package: rtracklayer
Attaching package: 'rtracklayer'
The following object is masked from 'package:igraph':
blocks
Loading required package: seqLogo
Loading required package: grImport
Loading required package: XML
Loading required package: MotIV
Attaching package: 'MotIV'
The following object is masked from 'package:rGADEM':
readPWMfile
The following object is masked from 'package:seqLogo':
makePWM
The following object is masked from 'package:igraph':
similarity
The following object is masked from 'package:stats':
filter
Loading required package: ade4
Attaching package: 'ade4'
The following object is masked from 'package:BSgenome':
score
The following object is masked from 'package:rtracklayer':
score
The following object is masked from 'package:Biostrings':
score
The following object is masked from 'package:GenomicRanges':
score
The following object is masked from 'package:IRanges':
score
The following object is masked from 'package:BiocGenerics':
score
|
|.............................................. | 66%
ordinary text without R code
label: unnamed-chunk-56 (with options)
List of 4
$ results: chr "asis"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|............................................... | 67%
ordinary text without R code
|
|............................................... | 68%
label: unnamed-chunk-57 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|................................................ | 68%
ordinary text without R code
|
|................................................ | 69%
label: unnamed-chunk-58 (with options)
List of 4
$ results: chr "asis"
$ echo : logi FALSE
$ message: logi FALSE
$ warning: logi FALSE
Found more than one class "textConnection" in cache; using the first, from namespace 'BiocGenerics'
Also defined by 'RJSONIO'
Found more than one class "textConnection" in cache; using the first, from namespace 'BiocGenerics'
Also defined by 'RJSONIO'
Found more than one class "textConnection" in cache; using the first, from namespace 'BiocGenerics'
Also defined by 'RJSONIO'
Downloading genome information (try:0) Using: Homo sapiens genes (GRCh37.p13)
o Fetching auxiliary information
oo Fetching probes genomic information
http://zwdzwd.io/InfiniumAnnotation/current/hm450/hm450.hg19.manifest.rds
trying URL 'http://zwdzwd.io/InfiniumAnnotation/current/hm450/hm450.hg19.manifest.rds'
Content type 'unknown' length 18546360 bytes (17.7 MB)
==================================================
downloaded 17.7 MB
oo Fetching TSS information
Downloading transcripts information. Using: Homo sapiens genes (GRCh37.p13)
o Mapping probes to nearest TSS
o Mapping results information
o Plotting figure
oo Adding names to genes
o Saving results
oo Saving significant results as: starburst_results.csv
oo It contains pair with changes both in the expression level of the nearest gene and in the DNA methylation level
|
|................................................. | 69%
ordinary text without R code
|
|................................................. | 70%
label: table4 (with options)
List of 4
$ echo : logi FALSE
$ message : logi FALSE
$ warnings: logi FALSE
$ results : chr "asis"
|
|................................................. | 71%
ordinary text without R code
|
|.................................................. | 71%
label: table5 (with options)
List of 4
$ echo : logi FALSE
$ message : logi FALSE
$ warnings: logi FALSE
$ results : chr "asis"
|
|.................................................. | 72%
ordinary text without R code
|
|................................................... | 73%
label: unnamed-chunk-59 (with options)
List of 5
$ results: chr "hide"
$ eval : logi TRUE
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
ChIPseeker v1.17.0 For help: https://guangchuangyu.github.io/ChIPseeker
If you use ChIPseeker in published research, please cite:
Guangchuang Yu, Li-Gen Wang, Qing-Yu He. ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics 2015, 31(14):2382-2383
ordinary text without R code
|
|.................................................... | 74%
label: unnamed-chunk-60 (with options)
List of 5
$ results: chr "hide"
$ eval : logi FALSE
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|.................................................... | 75%
ordinary text without R code
|
|..................................................... | 75%
label: unnamed-chunk-61 (with options)
List of 4
$ results: chr "hide"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|..................................................... | 76%
ordinary text without R code
|
|...................................................... | 76%
label: unnamed-chunk-62 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|...................................................... | 77%
ordinary text without R code
|
|...................................................... | 78%
label: unnamed-chunk-63 (with options)
List of 5
$ results : chr "asis"
$ echo : logi FALSE
$ fig.width: num 7
$ message : logi FALSE
$ warning : logi FALSE
|
|....................................................... | 78%
ordinary text without R code
|
|....................................................... | 79%
label: unnamed-chunk-64 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|........................................................ | 80%
ordinary text without R code
label: unnamed-chunk-65 (with options)
List of 4
$ results: chr "asis"
$ echo : logi FALSE
$ message: logi FALSE
$ warning: logi FALSE
Scale for 'x' is already present. Adding another scale for 'x', which will
replace the existing scale.
Scale for 'colour' is already present. Adding another scale for 'colour',
which will replace the existing scale.
|
|......................................................... | 81%
ordinary text without R code
|
|......................................................... | 82%
label: unnamed-chunk-66 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|.......................................................... | 82%
ordinary text without R code
|
|.......................................................... | 83%
label: unnamed-chunk-67 (with options)
List of 2
$ message: logi FALSE
$ warning: logi FALSE
Loading required package: ELMER.data
Downloading transcripts information. Using: Homo sapiens genes (GRCh37.p13)
Returning distal probes: 163576
Downloading genome information (try:0) Using: Homo sapiens genes (GRCh37.p13)
Loading from disk
Checking samples have both DNA methylation and Gene expression and they are in the same order...
Starting to add information to samples
=> Add clinical information to samples
Add FFPE information. More information at:
=> https://cancergenome.nih.gov/cancersselected/biospeccriteria
=> http://gdac.broadinstitute.org/runs/sampleReports/latest/FPPP_FFPE_Cases.html
=> Adding subtype information to samples
gbm subtype information from:doi:10.1016/j.cell.2015.12.028
=> Adding subtype information to samples
lgg subtype information from:doi:10.1016/j.cell.2015.12.028
Creating MultiAssayExperiment
MAE saved as: mae.rda
ordinary text without R code
|
|........................................................... | 84%
label: unnamed-chunk-68 (with options)
List of 2
$ warning: logi FALSE
$ results: chr "hide"
|
|........................................................... | 85%
ordinary text without R code
|
|............................................................ | 85%
label: unnamed-chunk-69 (with options)
List of 2
$ warning: logi FALSE
$ results: chr "hide"
Get diff methylated probes
ELMER will search for probes hypomethylated in group TCGA-GBM (n:10) compared to TCGA-LGG (n:10)
Number of probes: 2936
Saving results
Number of relevant probes found: 289
|
|............................................................ | 86%
label: unnamed-chunk-70 (with options)
List of 1
$ warning: logi FALSE
|
|............................................................. | 87%
ordinary text without R code
label: unnamed-chunk-71 (with options)
List of 2
$ warning: logi FALSE
$ results: chr "hide"
Get nearby genes
Searching for the 4 near genes
Get anti correlated probes-genes
Using pre-defined groups. U (unmethylated): TCGA-GBM, M (methylated): TCGA-LGG
-------------------
* Filtering probes
-------------------
For more information see function preAssociationProbeFiltering
Making sure we have at least 5% of beta values lesser than 0.3 and 5% of beta values greater 0.3.
Removing 2287 probes out of 2936
Calculating Pp (probe - gene) for all nearby genes
File created: elmer/hypo/getPair.hypo.all.pairs.statistic.csv
Downloading transcripts information. Using: Homo sapiens genes (GRCh37.p13)
Update the distance to gene to distance to the nearest TSS of the gene
|
|.............................................................. | 88%
label: unnamed-chunk-72 (with options)
List of 1
$ warning: logi FALSE
|
|.............................................................. | 89%
ordinary text without R code
label: unnamed-chunk-73 (with options)
List of 2
$ warning: logi FALSE
$ results: chr "hide"
Loading object: Probes.motif.hg19.450K
Retrieving TFClass family classification from ELMER.data.
Retrieving TFClass subfamily classification from ELMER.data.
------------------------------------
** Filtering motifs based on quality
------------------------------------
Number of enriched motifs with quality:
-----------
=> A: 15
=> B: 5
=> C: 6
=> D: 5
=> S: 0
-----------
Considering only motifs with quality from A up to DS: 31 motifs are enriched.
|
|............................................................... | 90%
label: unnamed-chunk-74 (with options)
List of 1
$ warning: logi FALSE
ordinary text without R code
|
|................................................................ | 91%
label: unnamed-chunk-75 (with options)
List of 2
$ warning: logi FALSE
$ results: chr "hide"
Using pre-defined groups. U (unmethylated): TCGA-GBM, M (methylated): TCGA-LGG
-------------------------------------------------------------------------------------------------------------------
** Downloading TF list from Lambert, Samuel A., et al. The human transcription factors. Cell 172.4 (2018): 650-665.
-------------------------------------------------------------------------------------------------------------------
Accessing TF families from TFClass database to indentify known potential TF
Retrieving TFClass family classification from ELMER.data.
Retrieving TFClass subfamily classification from ELMER.data.
Calculating the average methylation at all motif-adjacent probes
Performing Mann-Whitney U test
Finding potential TF and known potential TF
-----------------
** Creating plots
-----------------
TF rank plot highlighting TF in the same family (folder: elmer/hypo/TFrankPlot)
Retrieving TFClass family classification from ELMER.data.
Retrieving TFClass subfamily classification from ELMER.data.
Saving plot as: elmer/hypo/TFrankPlot/ETS1_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/KLF6_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/ELF5_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/ZN524_HUMAN.H11MO.0.D.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/FLI1_HUMAN.H11MO.1.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/EGR2_HUMAN.H11MO.1.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/FLI1_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/E2F4_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/ELF2_HUMAN.H11MO.0.C.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/E2F4_HUMAN.H11MO.1.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/PLAG1_HUMAN.H11MO.0.D.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/ZFP82_HUMAN.H11MO.0.C.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/GRHL1_HUMAN.H11MO.0.D.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/SP2_HUMAN.H11MO.1.B.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/PLAL1_HUMAN.H11MO.0.D.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/ERG_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/ZN341_HUMAN.H11MO.0.C.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/SP4_HUMAN.H11MO.1.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/ELK4_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/ZN467_HUMAN.H11MO.0.C.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/KLF14_HUMAN.H11MO.0.D.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/E2F3_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/ZN331_HUMAN.H11MO.0.C.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/FEV_HUMAN.H11MO.0.B.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/KLF15_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/KLF12_HUMAN.H11MO.0.C.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/MBD2_HUMAN.H11MO.0.B.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/EHF_HUMAN.H11MO.0.B.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/E2F1_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/ELF1_HUMAN.H11MO.0.A.TFrankPlot.pdf
Saving plot as: elmer/hypo/TFrankPlot/E2F7_HUMAN.H11MO.0.B.TFrankPlot.pdf
|
|................................................................ | 92%
ordinary text without R code
|
|................................................................. | 92%
label: unnamed-chunk-76 (with options)
List of 1
$ warning: logi FALSE
|
|................................................................. | 93%
ordinary text without R code
|
|.................................................................. | 94%
label: unnamed-chunk-77 (with options)
List of 6
$ echo : logi TRUE
$ message : logi FALSE
$ warnings : logi FALSE
$ results : chr "asis"
$ fig.height: num 6
$ fig.width : num 10
ordinary text without R code
|
|.................................................................. | 95%
label: unnamed-chunk-78 (with options)
List of 6
$ echo : logi TRUE
$ message : logi FALSE
$ warnings : logi FALSE
$ results : chr "asis"
$ fig.height: num 4
$ fig.width : num 15
|
|................................................................... | 96%
ordinary text without R code
label: unnamed-chunk-79 (with options)
List of 2
$ eval : logi FALSE
$ include: logi TRUE
|
|.................................................................... | 97%
label: unnamed-chunk-80 (with options)
List of 7
$ results : chr "asis"
$ echo : logi FALSE
$ message : logi FALSE
$ warning : logi FALSE
$ fig.align : chr "center"
$ fig.height: num 8
$ fig.width : num 8
Retrieving TFClass family classification from ELMER.data.
Retrieving TFClass subfamily classification from ELMER.data.
ordinary text without R code
|
|..................................................................... | 98%
label: unnamed-chunk-81 (with options)
List of 4
$ results: chr "asis"
$ echo : logi TRUE
$ message: logi FALSE
$ warning: logi FALSE
|
|..................................................................... | 99%
ordinary text without R code
|
|......................................................................| 99%
label: unnamed-chunk-82
|
|......................................................................| 100%
ordinary text without R code
output file: TCGAWorkflow.knit.md
/usr/bin/pandoc +RTS -K512m -RTS TCGAWorkflow.utf8.md --to html4 --from markdown+autolink_bare_uris+ascii_identifiers+tex_math_single_backslash+smart --output TCGAWorkflow.html --email-obfuscation none --self-contained --wrap preserve --standalone --section-divs --table-of-contents --toc-depth 3 --template /home/biocbuild/bbs-3.8-bioc/R/library/BiocStyle/resources/fragment.html --no-highlight --number-sections --include-in-header /tmp/RtmpVMEDM8/rmarkdown-str3530dc74ca0.html --mathjax --variable 'mathjax-url:https://mathjax.rstudio.com/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML' --filter /usr/bin/pandoc-citeproc
pandoc-citeproc: reference sham2014statistical not found
Output created: TCGAWorkflow.html
There were 50 or more warnings (use warnings() to see the first 50)
>
>
> knitr::purl(".buildwebvig/TCGAWorkflow/TCGAWorkflow.Rmd", ".buildwebvig/TCGAWorkflow/TCGAWorkflow.R")
processing file: .buildwebvig/TCGAWorkflow/TCGAWorkflow.Rmd
|
| | 0%
|
| | 1%
|
|. | 1%
|
|. | 2%
|
|.. | 3%
|
|.. | 4%
|
|... | 4%
|
|... | 5%
|
|.... | 6%
|
|..... | 7%
|
|..... | 8%
|
|...... | 9%
|
|...... | 10%
|
|....... | 10%
|
|....... | 11%
|
|........ | 12%
|
|........ | 13%
|
|......... | 13%
|
|......... | 14%
|
|.......... | 15%
|
|.......... | 16%
|
|........... | 17%
|
|............ | 18%
|
|............ | 19%
|
|............. | 20%
|
|.............. | 21%
|
|.............. | 22%
|
|............... | 23%
|
|............... | 24%
|
|................ | 24%
|
|................ | 25%
|
|................. | 25%
|
|................. | 26%
|
|................. | 27%
|
|.................. | 27%
|
|.................. | 28%
|
|................... | 29%
|
|................... | 30%
|
|.................... | 31%
|
|..................... | 32%
|
|...................... | 33%
|
|...................... | 34%
|
|....................... | 35%
|
|....................... | 36%
|
|........................ | 36%
|
|........................ | 37%
|
|........................ | 38%
|
|......................... | 38%
|
|......................... | 39%
|
|.......................... | 39%
|
|.......................... | 40%
|
|.......................... | 41%
|
|........................... | 41%
|
|........................... | 42%
|
|............................ | 43%
|
|............................. | 44%
|
|............................. | 45%
|
|.............................. | 46%
|
|............................... | 47%
|
|............................... | 48%
|
|................................ | 49%
|
|................................ | 50%
|
|................................. | 50%
|
|................................. | 51%
|
|.................................. | 52%
|
|.................................. | 53%
|
|................................... | 54%
|
|.................................... | 55%
|
|.................................... | 56%
|
|..................................... | 57%
|
|...................................... | 58%
|
|...................................... | 59%
|
|....................................... | 59%
|
|....................................... | 60%
|
|....................................... | 61%
|
|........................................ | 61%
|
|........................................ | 62%
|
|......................................... | 62%
|
|......................................... | 63%
|
|......................................... | 64%
|
|.......................................... | 64%
|
|.......................................... | 65%
|
|........................................... | 66%
|
|........................................... | 67%
|
|............................................ | 68%
|
|............................................. | 69%
|
|.............................................. | 70%
|
|.............................................. | 71%
|
|............................................... | 72%
|
|............................................... | 73%
|
|................................................ | 73%
|
|................................................ | 74%
|
|................................................ | 75%
|
|................................................. | 75%
|
|................................................. | 76%
|
|.................................................. | 76%
|
|.................................................. | 77%
|
|................................................... | 78%
|
|................................................... | 79%
|
|.................................................... | 80%
|
|..................................................... | 81%
|
|..................................................... | 82%
|
|...................................................... | 83%
|
|....................................................... | 84%
|
|....................................................... | 85%
|
|........................................................ | 86%
|
|........................................................ | 87%
|
|......................................................... | 87%
|
|......................................................... | 88%
|
|.......................................................... | 89%
|
|.......................................................... | 90%
|
|........................................................... | 90%
|
|........................................................... | 91%
|
|............................................................ | 92%
|
|............................................................ | 93%
|
|............................................................. | 94%
|
|.............................................................. | 95%
|
|.............................................................. | 96%
|
|............................................................... | 96%
|
|............................................................... | 97%
|
|................................................................ | 98%
|
|................................................................ | 99%
|
|.................................................................| 99%
|
|.................................................................| 100%
output file: .buildwebvig/TCGAWorkflow/TCGAWorkflow.R
[1] ".buildwebvig/TCGAWorkflow/TCGAWorkflow.R"
>
>