1 Introduction
Though there is no fundamental reason to separate the likelihood ratio testing and GSVA from the differential expression analyses, they are both tasks which reach back to the primary expression data, unlike GSEA. Therefore I chose to separate these two tasks into a their own document.
2 LRT
The likelihood ratio testing performed on the TMRC3 data is intended to look for shared patterns of expression across some known, constant factor. This may be time (are there patterns across the visits of each person), or cell type (do some genes act consistently across type). We could also ask this question of our two clinics, or indeed across the cure/fail samples.
2.1 Patterns across clinic
tc_clinical_filt <- normalize_expt(tc_clinical, filter = TRUE)## Removing 5633 low-count genes (14290 remaining).tc_lrt_clinic <- deseq_lrt(tc_clinical_filt, transform = "vst", interaction = FALSE,
interactor_column = "visitnumber",
interest_column = "clinic")## Warning in deseq_lrt(tc_clinical_filt, transform = "vst", interaction = FALSE, :
## The clinic should probably be a factor, set it with the 'factors' arg.## converting counts to integer mode## Warning in DESeqDataSet(se, design = design, ignoreRank): some variables in
## design formula are characters, converting to factors## estimating size factors## estimating dispersions## gene-wise dispersion estimates## mean-dispersion relationship## final dispersion estimates## fitting model and testing## -- replacing outliers and refitting for 169 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)## estimating dispersions## fitting model and testing## A large number of genes was given-- please, make sure this is not an error. Normally, only DE genes will be useful for this function.## Working with 6272 genes.## Working with 6269 genes after filtering: minc > 3## Joining, by = "merge"
## Joining, by = "merge"
summary(lrt_visit[["favorite_genes"]])## Error in h(simpleError(msg, call)): error in evaluating the argument 'object' in selecting a method for function 'summary': object 'lrt_visit' not foundt_clinical_filt <- normalize_expt(t_clinical, filter=TRUE)## Removing 5774 low-count genes (14149 remaining).lrt_visit <- deseq_lrt(t_clinical_filt, transform = "vst", interaction = FALSE,
interactor_column = "visitnumber",
interest_column = "finaloutcome")## converting counts to integer mode## estimating size factors## estimating dispersions## gene-wise dispersion estimates## mean-dispersion relationship## final dispersion estimates## fitting model and testing## -- replacing outliers and refitting for 120 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)## estimating dispersions## fitting model and testing## A large number of genes was given-- please, make sure this is not an error. Normally, only DE genes will be useful for this function.## Working with 5439 genes.## Working with 5435 genes after filtering: minc > 3## Joining, by = "merge"
## Joining, by = "merge"
summary(lrt_visit[["favorite_genes"]])## genes cluster
## Length:5435 Min. : 1.00
## Class :character 1st Qu.: 2.00
## Mode :character Median : 3.00
## Mean : 3.46
## 3rd Qu.: 3.00
## Max. :11.00written <- write_xlsx(data = as.data.frame(lrt_visit[["deseq_table"]]),
excel = glue::glue("excel/lrt_clinical_visit-v{ver}.xlsx"))
lrt_monocyte_visit <- deseq_lrt(t_visitcf_monocyte, transform = "vst",
interaction = FALSE,
interactor_column = "visitnumber",
interest_column = "finaloutcome")## converting counts to integer mode
## estimating size factors
## estimating dispersions
## gene-wise dispersion estimates
## mean-dispersion relationship
## final dispersion estimates
## fitting model and testing
## -- replacing outliers and refitting for 68 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)
## estimating dispersions
## fitting model and testing
## Working with 12 genes.
## Working with 12 genes after filtering: minc > 3
## Joining, by = "merge"Joining, by = "merge"
lrt_monocyte_visit$cluster_data$plot
lrt_monocyte_visit_v2 <- deseq_lrt(t_visitcf_monocyte, transform = "vst",
interaction = FALSE, minc = 1,
interactor_column = "visitnumber",
interest_column = "finaloutcome")## converting counts to integer mode
## estimating size factors
## estimating dispersions
## gene-wise dispersion estimates
## mean-dispersion relationship
## final dispersion estimates
## fitting model and testing
## -- replacing outliers and refitting for 68 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)
## estimating dispersions
## fitting model and testing
## Working with 12 genes.
## Working with 12 genes after filtering: minc > 1
## Joining, by = "merge"Joining, by = "merge"
2.2 Shared patterns across cell types
lrt_celltype_clinical_test <- deseq_lrt(tc_clinical, transform = "vst",
interactor_column = "typeofcells",
interest_column = "finaloutcome")## converting counts to integer mode## estimating size factors## estimating dispersions## gene-wise dispersion estimates## mean-dispersion relationship## final dispersion estimates## fitting model and testing## -- replacing outliers and refitting for 15 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)## estimating dispersions## fitting model and testing## Working with 552 genes.## Working with 551 genes after filtering: minc > 3## Joining, by = "merge"
## Joining, by = "merge"
hs_annot <- fData(hs_expt)
deseq_lrt_df <- merge(hs_annot, as.data.frame(lrt_celltype_clinical_test[["deseq_table"]]), all.y=TRUE,
by="row.names")
rownames(deseq_lrt_df) <- deseq_lrt_df[["Row.names"]]
deseq_lrt_df[["Row.names"]] <- NULL
written <- write_xlsx(data=deseq_lrt_df,
excel=glue::glue("excel/lrt_clinical_celltype-v{ver}.xlsx"))3 GSVA
3.1 Load some signatures
If one chooses, simple_gsva() can either load signatures from the GSVAdata package, which is a little old, or load an arbitrary set. load_gmt_signatures() provides a quick way to extract them from a gmt file.
broad_c7 <- load_gmt_signatures(signatures="reference/msigdb/c7.all.v7.5.1.entrez.gmt",
signature_category="c7")
broad_c2 <- load_gmt_signatures(signatures="reference/msigdb/c2.all.v7.5.1.entrez.gmt",
signature_category="c2")
broad_h <- load_gmt_signatures(signatures="reference/msigdb/h.all.v7.5.1.entrez.gmt",
signature_category="h")3.1.1 Clinical samples
3.1.1.1 Clinical C2
tc_celltype_gsva_c2 <- simple_gsva(
tc_valid, signatures=broad_c2,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 6353 rownames shared with the gsva result out of 6353.tc_celltype_gsva_c2_sig <- get_sig_gsva_categories(
tc_celltype_gsva_c2,
excel="analyses/3_cali_and_tumaco/GSVA/tc_valid_gsva_c2.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: failure_vs_cure. Adjust = BH## Limma step 6/6: 1/2: Creating table: cure. Adjust = BH## Limma step 6/6: 2/2: Creating table: failure. Adjust = BH## The factor cure has 122 rows.## The factor failure has 62 rows.## plot labels was not set and there are more than 100 samples, disabling it.## Testing each factor against the others.## Scoring cure against everything else.## Scoring failure against everything else.## Deleting the file analyses/3_cali_and_tumaco/GSVA/tc_valid_gsva_c2.xlsx before writing the tables.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.tc_celltype_gsva_c2_sig$subset_plottc_celltype_gsva_c2_sig$score_plot3.1.1.2 Valid C7
tc_celltype_gsva_c7 <- simple_gsva(
tc_valid, signatures=broad_c7,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c7")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 5209 rownames shared with the gsva result out of 5209.tc_celltype_gsva_c7_sig <- get_sig_gsva_categories(
tc_celltype_gsva_c7,
excel="analyses/3_cali_and_tumaco/GSVA/tc_valid_gsva_c7.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: failure_vs_cure. Adjust = BH## Limma step 6/6: 1/2: Creating table: cure. Adjust = BH## Limma step 6/6: 2/2: Creating table: failure. Adjust = BH## The factor cure has 122 rows.## The factor failure has 62 rows.## plot labels was not set and there are more than 100 samples, disabling it.## Testing each factor against the others.## Scoring cure against everything else.## Scoring failure against everything else.## Deleting the file analyses/3_cali_and_tumaco/GSVA/tc_valid_gsva_c7.xlsx before writing the tables.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.3.1.1.3 Valid H
tc_celltype_gsva_h <- simple_gsva(
tc_valid,
signatures=broad_h,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="h")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 50 rownames shared with the gsva result out of 50.tc_celltype_gsva_h_sig <- get_sig_gsva_categories(
tc_celltype_gsva_h,
excel="analyses/3_cali_and_tumaco/GSVA/tc_valid_gsva_h.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: failure_vs_cure. Adjust = BH## Limma step 6/6: 1/2: Creating table: cure. Adjust = BH## Limma step 6/6: 2/2: Creating table: failure. Adjust = BH## The factor cure has 122 rows.## The factor failure has 62 rows.## plot labels was not set and there are more than 100 samples, disabling it.## Testing each factor against the others.## Scoring cure against everything else.## Scoring failure against everything else.## Deleting the file analyses/3_cali_and_tumaco/GSVA/tc_valid_gsva_h.xlsx before writing the tables.3.1.2 Tumaco samples
3.1.2.1 Clinical C2
t_clinical_gsva_c2 <- simple_gsva(
t_clinical,
signatures=broad_c2,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c2")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 6353 rownames shared with the gsva result out of 6353.t_clinical_gsva_c2_sig <- get_sig_gsva_categories(
t_clinical_gsva_c2,
excel="analyses/4_tumaco/GSVA/t_clinical_gsva_c2.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: failure_vs_cure. Adjust = BH## Limma step 6/6: 1/2: Creating table: cure. Adjust = BH## Limma step 6/6: 2/2: Creating table: failure. Adjust = BH## The factor cure has 67 rows.## The factor failure has 56 rows.## plot labels was not set and there are more than 100 samples, disabling it.## Testing each factor against the others.## Scoring cure against everything else.## Scoring failure against everything else.## Deleting the file analyses/4_tumaco/GSVA/t_clinical_gsva_c2.xlsx before writing the tables.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.t_clinical_gsva_c2_sig$subset_plot
t_clinical_gsva_c2_sig$score_plot
3.1.2.2 Clinical C7
t_clinical_gsva_c7 <- simple_gsva(
t_clinical,
signatures=broad_c7,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c7")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 5209 rownames shared with the gsva result out of 5209.t_clinical_gsva_c7_sig <- get_sig_gsva_categories(
t_clinical_gsva_c7,
excel="analyses/4_tumaco/GSVA/t_clinical_gsva_c7.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: failure_vs_cure. Adjust = BH## Limma step 6/6: 1/2: Creating table: cure. Adjust = BH## Limma step 6/6: 2/2: Creating table: failure. Adjust = BH## The factor cure has 67 rows.## The factor failure has 56 rows.## plot labels was not set and there are more than 100 samples, disabling it.## Testing each factor against the others.## Scoring cure against everything else.## Scoring failure against everything else.## Deleting the file analyses/4_tumaco/GSVA/t_clinical_gsva_c7.xlsx before writing the tables.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.3.1.2.3 Clinical H
t_clinical_gsva_h <- simple_gsva(
t_clinical,
signatures=broad_h,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="h")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 50 rownames shared with the gsva result out of 50.t_clinical_gsva_h_sig <- get_sig_gsva_categories(
t_clinical_gsva_h,
excel="analyses/4_tumaco/GSVA/t_clinical_gsva_h.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: failure_vs_cure. Adjust = BH## Limma step 6/6: 1/2: Creating table: cure. Adjust = BH## Limma step 6/6: 2/2: Creating table: failure. Adjust = BH## The factor cure has 67 rows.## The factor failure has 56 rows.## plot labels was not set and there are more than 100 samples, disabling it.## Testing each factor against the others.## Scoring cure against everything else.## Scoring failure against everything else.## Deleting the file analyses/4_tumaco/GSVA/t_clinical_gsva_h.xlsx before writing the tables.3.1.2.4 Tumaco Biopsies c2
t_biopsy_gsva_c2 <- simple_gsva(
t_biopsies,
signatures=broad_c2,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c2")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 6353 rownames shared with the gsva result out of 6353.t_biopsy_gsva_c2_sig <- get_sig_gsva_categories(
t_biopsy_gsva_c2,
excel="analyses/4_tumaco/GSVA/t_biopsy_gsva_c2.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 9 rows.## The factor Tumaco_failure has 5 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.## Deleting the file analyses/4_tumaco/GSVA/t_biopsy_gsva_c2.xlsx before writing the tables.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.3.1.2.5 Tumaco Biopsies c7
t_biopsy_gsva_c7 <- simple_gsva(
t_biopsies,
signatures=broad_c7,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c7")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 5209 rownames shared with the gsva result out of 5209.t_biopsy_gsva_c7_sig <- get_sig_gsva_categories(
t_biopsy_gsva_c7,
excel="analyses/4_tumaco/GSVA/t_biopsy_gsva_c7.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 9 rows.## The factor Tumaco_failure has 5 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.## Deleting the file analyses/4_tumaco/GSVA/t_biopsy_gsva_c7.xlsx before writing the tables.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.3.1.2.6 Tumaco biopsies H
t_biopsy_gsva_h <- simple_gsva(
t_biopsies,
signatures=broad_h,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="h")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 50 rownames shared with the gsva result out of 50.t_biopsy_gsva_h_sig <- get_sig_gsva_categories(
t_biopsy_gsva_h,
excel="analyses/4_tumaco/GSVA/t_biopsy_gsva_h.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 9 rows.## The factor Tumaco_failure has 5 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.## Deleting the file analyses/4_tumaco/GSVA/t_biopsy_gsva_h.xlsx before writing the tables.3.1.2.7 Tumaco Eosinophils C7
t_eosinophil_gsva_c7 <- simple_gsva(
t_eosinophils,
signatures=broad_c7,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c7")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 5209 rownames shared with the gsva result out of 5209.t_eosinophil_gsva_c7_sig <- get_sig_gsva_categories(
t_eosinophil_gsva_c7,
excel="analyses/4_tumaco/GSVA/t_eosinophil_gsva_c7.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 17 rows.## The factor Tumaco_failure has 9 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.## Deleting the file analyses/4_tumaco/GSVA/t_eosinophil_gsva_c7.xlsx before writing the tables.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.3.1.2.8 Tumaco Eosinophils C2
t_eosinophil_gsva_c2 <- simple_gsva(
t_eosinophils,
signatures=broad_c2,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c2")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 6353 rownames shared with the gsva result out of 6353.t_eosinophil_gsva_c2_sig <- get_sig_gsva_categories(
t_eosinophil_gsva_c2,
excel="analyses/4_tumaco/GSVA/t_eosinophil_gsva_c2.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 17 rows.## The factor Tumaco_failure has 9 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.## Deleting the file analyses/4_tumaco/GSVA/t_eosinophil_gsva_c2.xlsx before writing the tables.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.3.1.2.9 Tumaco Eosinophils H
t_eosinophil_gsva_h <- simple_gsva(
t_eosinophils,
signatures=broad_h,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="h")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 50 rownames shared with the gsva result out of 50.t_eosinophil_gsva_h_sig <- get_sig_gsva_categories(
t_eosinophil_gsva_h,
excel="analyses/4_tumaco/GSVA/t_eosinophil_gsva_h.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 17 rows.## The factor Tumaco_failure has 9 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.## Deleting the file analyses/4_tumaco/GSVA/t_eosinophil_gsva_h.xlsx before writing the tables.3.1.2.10 Tumaco Monocytes C7
t_monocyte_gsva_c7 <- simple_gsva(
t_monocytes,
signatures=broad_c7,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c7")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 5209 rownames shared with the gsva result out of 5209.t_monocyte_gsva_c7_sig <- get_sig_gsva_categories(
t_monocyte_gsva_c7,
excel="analyses/4_tumaco/GSVA/t_monocyte_gsva_c7.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 21 rows.## The factor Tumaco_failure has 21 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.3.1.2.11 Tumaco Monocytes C2
t_monocyte_gsva_c2 <- simple_gsva(
t_monocytes,
signatures=broad_c2,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c2")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 6353 rownames shared with the gsva result out of 6353.t_monocyte_gsva_c2_sig <- get_sig_gsva_categories(
t_monocyte_gsva_c2,
excel="analyses/4_tumaco/GSVA/t_monocyte_gsva_c2.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 21 rows.## The factor Tumaco_failure has 21 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.3.1.2.12 Tumaco Monocytes H
t_monocyte_gsva_h <- simple_gsva(
t_monocytes,
signatures=broad_h,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="h")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 50 rownames shared with the gsva result out of 50.t_monocyte_gsva_h_sig <- get_sig_gsva_categories(
t_monocyte_gsva_h,
excel="analyses/4_tumaco/GSVA/t_monocyte_gsva_h.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 21 rows.## The factor Tumaco_failure has 21 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.3.1.2.13 Tumaco Neutrophils c7
t_neutrophil_gsva_c7 <- simple_gsva(
t_neutrophils,
signatures=broad_c7,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c7")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 5209 rownames shared with the gsva result out of 5209.t_neutrophil_gsva_c7_sig <- get_sig_gsva_categories(
t_neutrophil_gsva_c7,
excel="analyses/4_tumaco/GSVA/t_neutrophil_gsva_c7.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 20 rows.## The factor Tumaco_failure has 21 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.## Deleting the file analyses/4_tumaco/GSVA/t_neutrophil_gsva_c7.xlsx before writing the tables.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 84296|PLEKHA1,PL ... 20H1,KMT5B,51111 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 0750246|PRDM5,PR ... 9|OAS1,OAS1,4938 is truncated.
## Number of characters exeed the limit of 32767.3.1.2.14 Tumaco Neutrophils c2
t_neutrophil_gsva_c2 <- simple_gsva(
t_neutrophils,
signatures=broad_c2,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="c2")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 6353 rownames shared with the gsva result out of 6353.t_neutrophil_gsva_c2_sig <- get_sig_gsva_categories(
t_neutrophil_gsva_c2,
excel="analyses/4_tumaco/GSVA/t_neutrophil_gsva_c2.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 20 rows.## The factor Tumaco_failure has 21 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ENSG00000273641, ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2|39449_at,,|355 ... 14_at,UBA2,10054 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : t,NUFIP1,26747|2 ... x_at,EEF1A1,1915 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : NUP42,11097|1654 ... 4_at,UBE2D3,7323 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 203395_s_at,HES1 ... 4_s_at,CTBS,1486 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 16224|Arid4b,ARI ... gcs1,HMGCS1,3157 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 343|ENSG00000182 ... 499,FAM20C,56975 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : ,51283|1442393_a ... 11_s_at,DBP,1628 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 7|ENSG0000019736 ... 000288631,162699 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |219698_s_at,MET ... _at,TMEM97,27346 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 643180|218637_at ... 0_s_at,NREP,9315 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 701_s_at,CHCHD7, ... 0_at,SH2B3,10019 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 2804,KRTAP10-7,3 ... 8516,PRKACA,5566 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 101927051|204365 ... _at,EPS8L2,64787 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : |rc_AA957218_at, ... 8|rc_H33093_at,, is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : 824,HDGFL1,15415 ... 7495,TMED3,23423 is truncated.
## Number of characters exeed the limit of 32767.## Warning in wb$writeData(df = x, colNames = TRUE, sheet = sheet, startRow = startRow, : L2,27039|8108287 ... 7,CCDC187,399693 is truncated.
## Number of characters exeed the limit of 32767.3.1.2.15 Tumaco Neutrophils H
t_neutrophil_gsva_h <- simple_gsva(
t_neutrophils,
signatures=broad_h,
msig_xml="reference/msigdb/msigdb_v7.5.1.xml",
signature_category="h")## Converting the rownames() of the expressionset to ENTREZID.## 526 ENSEMBL ID's didn't have a matching ENTEREZ ID. Dropping them now.## Before conversion, the expressionset has 19923 entries.## After conversion, the expressionset has 19549 entries.## Adding annotations from reference/msigdb/msigdb_v7.5.1.xml.## Not subsetting the msigdb metadata, the wanted_meta argument was 'all'.## The downloaded msigdb contained 50 rownames shared with the gsva result out of 50.t_neutrophil_gsva_h_sig <- get_sig_gsva_categories(
t_neutrophil_gsva_h,
excel="analyses/4_tumaco/GSVA/t_neutrophil_gsva_h.xlsx")## Starting limma pairwise comparison.## libsize was not specified, this parameter has profound effects on limma's result.## Using the libsize from expt$libsize.## Limma step 1/6: choosing model.## Choosing the non-intercept containing model.## Assuming this data is similar to a micro array and not performign voom.## Limma step 3/6: running lmFit with method: ls.## Limma step 4/6: making and fitting contrasts with no intercept. (~ 0 + factors)## Limma step 5/6: Running eBayes with robust = FALSE and trend = FALSE.## Limma step 6/6: Writing limma outputs.## Limma step 6/6: 1/1: Creating table: Tumacofailure_vs_Tumacocure. Adjust = BH## Limma step 6/6: 1/2: Creating table: Tumacocure. Adjust = BH## Limma step 6/6: 2/2: Creating table: Tumacofailure. Adjust = BH## The factor Tumaco_cure has 20 rows.## The factor Tumaco_failure has 21 rows.## Testing each factor against the others.## Scoring Tumaco_cure against everything else.## Scoring Tumaco_failure against everything else.if (!isTRUE(get0("skip_load"))) {
pander::pander(sessionInfo())
message(paste0("This is hpgltools commit: ", get_git_commit()))
message(paste0("Saving to ", savefile))
tmp <- sm(saveme(filename=savefile))
}tmp <- loadme(filename=savefile)LS0tCnRpdGxlOiAiTC4gcGFuYW1lbnNpcyAyMDIyMDk6IExpa2VsaWhvb2QgUmF0aW8gVGVzdHMgYW5kIEdTVkEiCmF1dGhvcjogImF0YiBhYmVsZXdAZ21haWwuY29tIgpkYXRlOiAiYHIgU3lzLkRhdGUoKWAiCm91dHB1dDoKICBodG1sX2RvY3VtZW50OgogICAgY29kZV9kb3dubG9hZDogdHJ1ZQogICAgY29kZV9mb2xkaW5nOiBzaG93CiAgICBmaWdfY2FwdGlvbjogdHJ1ZQogICAgZmlnX2hlaWdodDogNwogICAgZmlnX3dpZHRoOiA3CiAgICBoaWdobGlnaHQ6IGRlZmF1bHQKICAgIGtlZXBfbWQ6IGZhbHNlCiAgICBtb2RlOiBzZWxmY29udGFpbmVkCiAgICBudW1iZXJfc2VjdGlvbnM6IHRydWUKICAgIHNlbGZfY29udGFpbmVkOiB0cnVlCiAgICB0aGVtZTogcmVhZGFibGUKICAgIHRvYzogdHJ1ZQogICAgdG9jX2Zsb2F0OgogICAgICBjb2xsYXBzZWQ6IGZhbHNlCiAgICAgIHNtb290aF9zY3JvbGw6IGZhbHNlCi0tLQoKPHN0eWxlPgpib2R5IC5tYWluLWNvbnRhaW5lciB7Cm1heC13aWR0aDogMTYwMHB4Owp9Cjwvc3R5bGU+CgpgYGB7ciBvcHRpb25zLCBpbmNsdWRlPUZBTFNFfQpsaWJyYXJ5KGhwZ2x0b29scykKbGlicmFyeShkcGx5cikKbGlicmFyeShmb3JjYXRzKQp0dCA8LSBzbShkZXZ0b29sczo6bG9hZF9hbGwoIn4vaHBnbHRvb2xzIikpCmtuaXRyOjpvcHRzX2tuaXQkc2V0KHByb2dyZXNzPVRSVUUsCiAgICAgICAgICAgICAgICAgICAgIHZlcmJvc2U9VFJVRSwKICAgICAgICAgICAgICAgICAgICAgd2lkdGg9OTAsCiAgICAgICAgICAgICAgICAgICAgIGVjaG89VFJVRSkKa25pdHI6Om9wdHNfY2h1bmskc2V0KGVycm9yPVRSVUUsCiAgICAgICAgICAgICAgICAgICAgICBmaWcud2lkdGg9OCwKICAgICAgICAgICAgICAgICAgICAgIGZpZy5oZWlnaHQ9OCwKICAgICAgICAgICAgICAgICAgICAgIGRwaT05NikKb2xkX29wdGlvbnMgPC0gb3B0aW9ucyhkaWdpdHM9NCwKICAgICAgICAgICAgICAgICAgICAgICBzdHJpbmdzQXNGYWN0b3JzPUZBTFNFLAogICAgICAgICAgICAgICAgICAgICAgIGtuaXRyLmR1cGxpY2F0ZS5sYWJlbD0iYWxsb3ciKQpnZ3Bsb3QyOjp0aGVtZV9zZXQoZ2dwbG90Mjo6dGhlbWVfYncoYmFzZV9zaXplPTEyKSkKdmVyIDwtICIyMDIyMDkiCnByZXZpb3VzX2ZpbGUgPC0gIiIKcnVuZGF0ZSA8LSBmb3JtYXQoU3lzLkRhdGUoKSwgZm9ybWF0PSIlWSVtJWQiKQoKIyN0bXAgPC0gdHJ5KHNtKGxvYWRtZShmaWxlbmFtZT1nc3ViKHBhdHRlcm49IlxcLlJtZCIsIHJlcGxhY2U9IlxcLnJkYVxcLnh6IiwgeD1wcmV2aW91c19maWxlKSkpKQpybWRfZmlsZSA8LSBnbHVlOjpnbHVlKCJ0bXJjM19scnRfZ3N2YV97dmVyfS5SbWQiKQpzYXZlZmlsZSA8LSBnc3ViKHBhdHRlcm49IlxcLlJtZCIsIHJlcGxhY2U9IlxcLnJkYVxcLnh6IiwgeD1ybWRfZmlsZSkKbG9hZGVkIDwtIGxvYWQoZmlsZT1nbHVlOjpnbHVlKCJyZGEvdG1yYzNfZGF0YV9zdHJ1Y3R1cmVzLXZ7dmVyfS5yZGEiKSkKYGBgCgojIEludHJvZHVjdGlvbgoKVGhvdWdoIHRoZXJlIGlzIG5vIGZ1bmRhbWVudGFsIHJlYXNvbiB0byBzZXBhcmF0ZSB0aGUgbGlrZWxpaG9vZCByYXRpbwp0ZXN0aW5nIGFuZCBHU1ZBIGZyb20gdGhlIGRpZmZlcmVudGlhbCBleHByZXNzaW9uIGFuYWx5c2VzLCB0aGV5IGFyZQpib3RoIHRhc2tzIHdoaWNoIHJlYWNoIGJhY2sgdG8gdGhlIHByaW1hcnkgZXhwcmVzc2lvbiBkYXRhLCB1bmxpa2UKR1NFQS4gIFRoZXJlZm9yZSBJIGNob3NlIHRvIHNlcGFyYXRlIHRoZXNlIHR3byB0YXNrcyBpbnRvIGEgdGhlaXIgb3duCmRvY3VtZW50LgoKIyBMUlQKClRoZSBsaWtlbGlob29kIHJhdGlvIHRlc3RpbmcgcGVyZm9ybWVkIG9uIHRoZSBUTVJDMyBkYXRhIGlzIGludGVuZGVkCnRvIGxvb2sgZm9yIHNoYXJlZCBwYXR0ZXJucyBvZiBleHByZXNzaW9uIGFjcm9zcyBzb21lIGtub3duLCBjb25zdGFudApmYWN0b3IuICBUaGlzIG1heSBiZSB0aW1lIChhcmUgdGhlcmUgcGF0dGVybnMgYWNyb3NzIHRoZSB2aXNpdHMgb2YKZWFjaCBwZXJzb24pLCBvciBjZWxsIHR5cGUgKGRvIHNvbWUgZ2VuZXMgYWN0IGNvbnNpc3RlbnRseSBhY3Jvc3MKdHlwZSkuICBXZSBjb3VsZCBhbHNvIGFzayB0aGlzIHF1ZXN0aW9uIG9mIG91ciB0d28gY2xpbmljcywgb3IgaW5kZWVkCmFjcm9zcyB0aGUgY3VyZS9mYWlsIHNhbXBsZXMuCgojIyBQYXR0ZXJucyBhY3Jvc3MgY2xpbmljCgpgYGB7ciBscnRfY2xpbmljc30KdGNfY2xpbmljYWxfZmlsdCA8LSBub3JtYWxpemVfZXhwdCh0Y19jbGluaWNhbCwgZmlsdGVyID0gVFJVRSkKCnRjX2xydF9jbGluaWMgPC0gZGVzZXFfbHJ0KHRjX2NsaW5pY2FsX2ZpbHQsIHRyYW5zZm9ybSA9ICJ2c3QiLCBpbnRlcmFjdGlvbiA9IEZBTFNFLAogICAgICAgICAgICAgICAgICAgICAgICAgICBpbnRlcmFjdG9yX2NvbHVtbiA9ICJ2aXNpdG51bWJlciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgIGludGVyZXN0X2NvbHVtbiA9ICJjbGluaWMiKQpzdW1tYXJ5KGxydF92aXNpdFtbImZhdm9yaXRlX2dlbmVzIl1dKQpgYGAKCmBgYHtyIGxydF92aXNpdH0KdF9jbGluaWNhbF9maWx0IDwtIG5vcm1hbGl6ZV9leHB0KHRfY2xpbmljYWwsIGZpbHRlcj1UUlVFKQoKbHJ0X3Zpc2l0IDwtIGRlc2VxX2xydCh0X2NsaW5pY2FsX2ZpbHQsIHRyYW5zZm9ybSA9ICJ2c3QiLCBpbnRlcmFjdGlvbiA9IEZBTFNFLAogICAgICAgICAgICAgICAgICAgICAgIGludGVyYWN0b3JfY29sdW1uID0gInZpc2l0bnVtYmVyIiwKICAgICAgICAgICAgICAgICAgICAgICBpbnRlcmVzdF9jb2x1bW4gPSAiZmluYWxvdXRjb21lIikKc3VtbWFyeShscnRfdmlzaXRbWyJmYXZvcml0ZV9nZW5lcyJdXSkKd3JpdHRlbiA8LSB3cml0ZV94bHN4KGRhdGEgPSBhcy5kYXRhLmZyYW1lKGxydF92aXNpdFtbImRlc2VxX3RhYmxlIl1dKSwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvbHJ0X2NsaW5pY2FsX3Zpc2l0LXZ7dmVyfS54bHN4IikpCgpscnRfbW9ub2N5dGVfdmlzaXQgPC0gZGVzZXFfbHJ0KHRfdmlzaXRjZl9tb25vY3l0ZSwgdHJhbnNmb3JtID0gInZzdCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgaW50ZXJhY3Rpb24gPSBGQUxTRSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBpbnRlcmFjdG9yX2NvbHVtbiA9ICJ2aXNpdG51bWJlciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgaW50ZXJlc3RfY29sdW1uID0gImZpbmFsb3V0Y29tZSIpCmxydF9tb25vY3l0ZV92aXNpdCRjbHVzdGVyX2RhdGEkcGxvdAoKbHJ0X21vbm9jeXRlX3Zpc2l0X3YyIDwtIGRlc2VxX2xydCh0X3Zpc2l0Y2ZfbW9ub2N5dGUsIHRyYW5zZm9ybSA9ICJ2c3QiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGludGVyYWN0aW9uID0gRkFMU0UsIG1pbmMgPSAxLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGludGVyYWN0b3JfY29sdW1uID0gInZpc2l0bnVtYmVyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBpbnRlcmVzdF9jb2x1bW4gPSAiZmluYWxvdXRjb21lIikKYGBgCgojIyBTaGFyZWQgcGF0dGVybnMgYWNyb3NzIGNlbGwgdHlwZXMKCmBgYHtyIGxydF9jZWxsdHlwZX0KbHJ0X2NlbGx0eXBlX2NsaW5pY2FsX3Rlc3QgPC0gZGVzZXFfbHJ0KHRjX2NsaW5pY2FsLCB0cmFuc2Zvcm0gPSAidnN0IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGludGVyYWN0b3JfY29sdW1uID0gInR5cGVvZmNlbGxzIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGludGVyZXN0X2NvbHVtbiA9ICJmaW5hbG91dGNvbWUiKQoKaHNfYW5ub3QgPC0gZkRhdGEoaHNfZXhwdCkKZGVzZXFfbHJ0X2RmIDwtIG1lcmdlKGhzX2Fubm90LCBhcy5kYXRhLmZyYW1lKGxydF9jZWxsdHlwZV9jbGluaWNhbF90ZXN0W1siZGVzZXFfdGFibGUiXV0pLCBhbGwueT1UUlVFLAogICAgICAgICAgICAgICAgICAgICAgYnk9InJvdy5uYW1lcyIpCnJvd25hbWVzKGRlc2VxX2xydF9kZikgPC0gZGVzZXFfbHJ0X2RmW1siUm93Lm5hbWVzIl1dCmRlc2VxX2xydF9kZltbIlJvdy5uYW1lcyJdXSA8LSBOVUxMCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhPWRlc2VxX2xydF9kZiwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsPWdsdWU6OmdsdWUoImV4Y2VsL2xydF9jbGluaWNhbF9jZWxsdHlwZS12e3Zlcn0ueGxzeCIpKQpgYGAKCiMgR1NWQQoKIyMgTG9hZCBzb21lIHNpZ25hdHVyZXMKCklmIG9uZSBjaG9vc2VzLCBzaW1wbGVfZ3N2YSgpIGNhbiBlaXRoZXIgbG9hZCBzaWduYXR1cmVzIGZyb20gdGhlCkdTVkFkYXRhIHBhY2thZ2UsIHdoaWNoIGlzIGEgbGl0dGxlIG9sZCwgb3IgbG9hZCBhbiBhcmJpdHJhcnkgc2V0Lgpsb2FkX2dtdF9zaWduYXR1cmVzKCkgcHJvdmlkZXMgYSBxdWljayB3YXkgdG8gZXh0cmFjdCB0aGVtIGZyb20gYSBnbXQKZmlsZS4KCgpgYGB7ciBsb2FkX3NpZ25hdHVyZXN9CmJyb2FkX2M3IDwtIGxvYWRfZ210X3NpZ25hdHVyZXMoc2lnbmF0dXJlcz0icmVmZXJlbmNlL21zaWdkYi9jNy5hbGwudjcuNS4xLmVudHJlei5nbXQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHNpZ25hdHVyZV9jYXRlZ29yeT0iYzciKQpicm9hZF9jMiA8LSBsb2FkX2dtdF9zaWduYXR1cmVzKHNpZ25hdHVyZXM9InJlZmVyZW5jZS9tc2lnZGIvYzIuYWxsLnY3LjUuMS5lbnRyZXouZ210IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzaWduYXR1cmVfY2F0ZWdvcnk9ImMyIikKYnJvYWRfaCA8LSBsb2FkX2dtdF9zaWduYXR1cmVzKHNpZ25hdHVyZXM9InJlZmVyZW5jZS9tc2lnZGIvaC5hbGwudjcuNS4xLmVudHJlei5nbXQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc2lnbmF0dXJlX2NhdGVnb3J5PSJoIikKYGBgCgojIyMgQ2xpbmljYWwgc2FtcGxlcwoKIyMjIyBDbGluaWNhbCBDMgoKYGBge3IgZ3N2YSwgZmlnLnNob3c9ImhpZGUifQp0Y19jZWxsdHlwZV9nc3ZhX2MyIDwtIHNpbXBsZV9nc3ZhKAogICAgdGNfdmFsaWQsIHNpZ25hdHVyZXM9YnJvYWRfYzIsCiAgICBtc2lnX3htbD0icmVmZXJlbmNlL21zaWdkYi9tc2lnZGJfdjcuNS4xLnhtbCIpCnRjX2NlbGx0eXBlX2dzdmFfYzJfc2lnIDwtIGdldF9zaWdfZ3N2YV9jYXRlZ29yaWVzKAogICAgdGNfY2VsbHR5cGVfZ3N2YV9jMiwKICAgIGV4Y2VsPSJhbmFseXNlcy8zX2NhbGlfYW5kX3R1bWFjby9HU1ZBL3RjX3ZhbGlkX2dzdmFfYzIueGxzeCIpCnRjX2NlbGx0eXBlX2dzdmFfYzJfc2lnJHN1YnNldF9wbG90CnRjX2NlbGx0eXBlX2dzdmFfYzJfc2lnJHNjb3JlX3Bsb3QKYGBgCgojIyMjIFZhbGlkICBDNwoKYGBge3IgZ3N2YV92YWxpZF9jN30KdGNfY2VsbHR5cGVfZ3N2YV9jNyA8LSBzaW1wbGVfZ3N2YSgKICAgIHRjX3ZhbGlkLCBzaWduYXR1cmVzPWJyb2FkX2M3LAogICAgbXNpZ194bWw9InJlZmVyZW5jZS9tc2lnZGIvbXNpZ2RiX3Y3LjUuMS54bWwiLAogICAgc2lnbmF0dXJlX2NhdGVnb3J5PSJjNyIpCnRjX2NlbGx0eXBlX2dzdmFfYzdfc2lnIDwtIGdldF9zaWdfZ3N2YV9jYXRlZ29yaWVzKAogICAgdGNfY2VsbHR5cGVfZ3N2YV9jNywKICAgIGV4Y2VsPSJhbmFseXNlcy8zX2NhbGlfYW5kX3R1bWFjby9HU1ZBL3RjX3ZhbGlkX2dzdmFfYzcueGxzeCIpCmBgYAoKIyMjIyBWYWxpZCBICgpgYGB7ciBnc3ZhX3ZhbGlkX2h9CnRjX2NlbGx0eXBlX2dzdmFfaCA8LSBzaW1wbGVfZ3N2YSgKICAgIHRjX3ZhbGlkLAogICAgc2lnbmF0dXJlcz1icm9hZF9oLAogICAgbXNpZ194bWw9InJlZmVyZW5jZS9tc2lnZGIvbXNpZ2RiX3Y3LjUuMS54bWwiLAogICAgc2lnbmF0dXJlX2NhdGVnb3J5PSJoIikKdGNfY2VsbHR5cGVfZ3N2YV9oX3NpZyA8LSBnZXRfc2lnX2dzdmFfY2F0ZWdvcmllcygKICAgIHRjX2NlbGx0eXBlX2dzdmFfaCwKICAgIGV4Y2VsPSJhbmFseXNlcy8zX2NhbGlfYW5kX3R1bWFjby9HU1ZBL3RjX3ZhbGlkX2dzdmFfaC54bHN4IikKYGBgCgojIyMgVHVtYWNvIHNhbXBsZXMKCiMjIyMgQ2xpbmljYWwgQzIKCmBgYHtyIGdzdmFfdHVtYWNvX2NsaW5pY2FsX2MyfQp0X2NsaW5pY2FsX2dzdmFfYzIgPC0gc2ltcGxlX2dzdmEoCiAgICB0X2NsaW5pY2FsLAogICAgc2lnbmF0dXJlcz1icm9hZF9jMiwKICAgIG1zaWdfeG1sPSJyZWZlcmVuY2UvbXNpZ2RiL21zaWdkYl92Ny41LjEueG1sIiwKICAgIHNpZ25hdHVyZV9jYXRlZ29yeT0iYzIiKQp0X2NsaW5pY2FsX2dzdmFfYzJfc2lnIDwtIGdldF9zaWdfZ3N2YV9jYXRlZ29yaWVzKAogICAgdF9jbGluaWNhbF9nc3ZhX2MyLAogICAgZXhjZWw9ImFuYWx5c2VzLzRfdHVtYWNvL0dTVkEvdF9jbGluaWNhbF9nc3ZhX2MyLnhsc3giKQp0X2NsaW5pY2FsX2dzdmFfYzJfc2lnJHN1YnNldF9wbG90CnRfY2xpbmljYWxfZ3N2YV9jMl9zaWckc2NvcmVfcGxvdApgYGAKCiMjIyMgQ2xpbmljYWwgQzcKCmBgYHtyIHRfZ3N2YV9jbGluaWNhbF9jN30KdF9jbGluaWNhbF9nc3ZhX2M3IDwtIHNpbXBsZV9nc3ZhKAogICAgdF9jbGluaWNhbCwKICAgIHNpZ25hdHVyZXM9YnJvYWRfYzcsCiAgICBtc2lnX3htbD0icmVmZXJlbmNlL21zaWdkYi9tc2lnZGJfdjcuNS4xLnhtbCIsCiAgICBzaWduYXR1cmVfY2F0ZWdvcnk9ImM3IikKdF9jbGluaWNhbF9nc3ZhX2M3X3NpZyA8LSBnZXRfc2lnX2dzdmFfY2F0ZWdvcmllcygKICAgIHRfY2xpbmljYWxfZ3N2YV9jNywKICAgIGV4Y2VsPSJhbmFseXNlcy80X3R1bWFjby9HU1ZBL3RfY2xpbmljYWxfZ3N2YV9jNy54bHN4IikKYGBgCgojIyMjIENsaW5pY2FsIEgKCmBgYHtyIHRfY2xpbmljYWxfZ3N2YV9ofQp0X2NsaW5pY2FsX2dzdmFfaCA8LSBzaW1wbGVfZ3N2YSgKICAgIHRfY2xpbmljYWwsCiAgICBzaWduYXR1cmVzPWJyb2FkX2gsCiAgICBtc2lnX3htbD0icmVmZXJlbmNlL21zaWdkYi9tc2lnZGJfdjcuNS4xLnhtbCIsCiAgICBzaWduYXR1cmVfY2F0ZWdvcnk9ImgiKQp0X2NsaW5pY2FsX2dzdmFfaF9zaWcgPC0gZ2V0X3NpZ19nc3ZhX2NhdGVnb3JpZXMoCiAgICB0X2NsaW5pY2FsX2dzdmFfaCwKICAgIGV4Y2VsPSJhbmFseXNlcy80X3R1bWFjby9HU1ZBL3RfY2xpbmljYWxfZ3N2YV9oLnhsc3giKQpgYGAKCiMjIyMgVHVtYWNvIEJpb3BzaWVzIGMyCgpgYGB7ciBnc3ZhX3R1bWFjb19iaW9wc3lfYzJ9CnRfYmlvcHN5X2dzdmFfYzIgPC0gc2ltcGxlX2dzdmEoCiAgICB0X2Jpb3BzaWVzLAogICAgc2lnbmF0dXJlcz1icm9hZF9jMiwKICAgIG1zaWdfeG1sPSJyZWZlcmVuY2UvbXNpZ2RiL21zaWdkYl92Ny41LjEueG1sIiwKICAgIHNpZ25hdHVyZV9jYXRlZ29yeT0iYzIiKQp0X2Jpb3BzeV9nc3ZhX2MyX3NpZyA8LSBnZXRfc2lnX2dzdmFfY2F0ZWdvcmllcygKICAgIHRfYmlvcHN5X2dzdmFfYzIsCiAgICBleGNlbD0iYW5hbHlzZXMvNF90dW1hY28vR1NWQS90X2Jpb3BzeV9nc3ZhX2MyLnhsc3giKQpgYGAKCiMjIyMgVHVtYWNvIEJpb3BzaWVzIGM3CgpgYGB7ciBnc3ZhX3R1bWFjb19iaW9wc3lfYzd9CnRfYmlvcHN5X2dzdmFfYzcgPC0gc2ltcGxlX2dzdmEoCiAgICB0X2Jpb3BzaWVzLAogICAgc2lnbmF0dXJlcz1icm9hZF9jNywKICAgIG1zaWdfeG1sPSJyZWZlcmVuY2UvbXNpZ2RiL21zaWdkYl92Ny41LjEueG1sIiwKICAgIHNpZ25hdHVyZV9jYXRlZ29yeT0iYzciKQp0X2Jpb3BzeV9nc3ZhX2M3X3NpZyA8LSBnZXRfc2lnX2dzdmFfY2F0ZWdvcmllcygKICAgIHRfYmlvcHN5X2dzdmFfYzcsCiAgICBleGNlbD0iYW5hbHlzZXMvNF90dW1hY28vR1NWQS90X2Jpb3BzeV9nc3ZhX2M3Lnhsc3giKQpgYGAKCiMjIyMgVHVtYWNvIGJpb3BzaWVzIEgKCmBgYHtyIGdzdmFfdHVtYWNvX2Jpb3BzaWVzX2h9CnRfYmlvcHN5X2dzdmFfaCA8LSBzaW1wbGVfZ3N2YSgKICAgIHRfYmlvcHNpZXMsCiAgICBzaWduYXR1cmVzPWJyb2FkX2gsCiAgICBtc2lnX3htbD0icmVmZXJlbmNlL21zaWdkYi9tc2lnZGJfdjcuNS4xLnhtbCIsCiAgICBzaWduYXR1cmVfY2F0ZWdvcnk9ImgiKQp0X2Jpb3BzeV9nc3ZhX2hfc2lnIDwtIGdldF9zaWdfZ3N2YV9jYXRlZ29yaWVzKAogICAgdF9iaW9wc3lfZ3N2YV9oLAogICAgZXhjZWw9ImFuYWx5c2VzLzRfdHVtYWNvL0dTVkEvdF9iaW9wc3lfZ3N2YV9oLnhsc3giKQpgYGAKCiMjIyMgVHVtYWNvIEVvc2lub3BoaWxzIEM3CgpgYGB7ciBnc3ZhX3R1bWFjb19lb3Npbm9waGlsX2M3fQp0X2Vvc2lub3BoaWxfZ3N2YV9jNyA8LSBzaW1wbGVfZ3N2YSgKICAgIHRfZW9zaW5vcGhpbHMsCiAgICBzaWduYXR1cmVzPWJyb2FkX2M3LAogICAgbXNpZ194bWw9InJlZmVyZW5jZS9tc2lnZGIvbXNpZ2RiX3Y3LjUuMS54bWwiLAogICAgc2lnbmF0dXJlX2NhdGVnb3J5PSJjNyIpCnRfZW9zaW5vcGhpbF9nc3ZhX2M3X3NpZyA8LSBnZXRfc2lnX2dzdmFfY2F0ZWdvcmllcygKICAgIHRfZW9zaW5vcGhpbF9nc3ZhX2M3LAogICAgZXhjZWw9ImFuYWx5c2VzLzRfdHVtYWNvL0dTVkEvdF9lb3Npbm9waGlsX2dzdmFfYzcueGxzeCIpCmBgYAoKIyMjIyBUdW1hY28gRW9zaW5vcGhpbHMgQzIKCmBgYHtyIGdzdmFfdHVtYWNvX2Vvc2lub3BoaWxfYzJ9CnRfZW9zaW5vcGhpbF9nc3ZhX2MyIDwtIHNpbXBsZV9nc3ZhKAogICAgdF9lb3Npbm9waGlscywKICAgIHNpZ25hdHVyZXM9YnJvYWRfYzIsCiAgICBtc2lnX3htbD0icmVmZXJlbmNlL21zaWdkYi9tc2lnZGJfdjcuNS4xLnhtbCIsCiAgICBzaWduYXR1cmVfY2F0ZWdvcnk9ImMyIikKdF9lb3Npbm9waGlsX2dzdmFfYzJfc2lnIDwtIGdldF9zaWdfZ3N2YV9jYXRlZ29yaWVzKAogICAgdF9lb3Npbm9waGlsX2dzdmFfYzIsCiAgICBleGNlbD0iYW5hbHlzZXMvNF90dW1hY28vR1NWQS90X2Vvc2lub3BoaWxfZ3N2YV9jMi54bHN4IikKYGBgCgojIyMjIFR1bWFjbyBFb3Npbm9waGlscyBICgpgYGB7ciBnc3ZhX3R1bWFjb19lb3Npbm9waGlsX2h9CnRfZW9zaW5vcGhpbF9nc3ZhX2ggPC0gc2ltcGxlX2dzdmEoCiAgICB0X2Vvc2lub3BoaWxzLAogICAgc2lnbmF0dXJlcz1icm9hZF9oLAogICAgbXNpZ194bWw9InJlZmVyZW5jZS9tc2lnZGIvbXNpZ2RiX3Y3LjUuMS54bWwiLAogICAgc2lnbmF0dXJlX2NhdGVnb3J5PSJoIikKdF9lb3Npbm9waGlsX2dzdmFfaF9zaWcgPC0gZ2V0X3NpZ19nc3ZhX2NhdGVnb3JpZXMoCiAgICB0X2Vvc2lub3BoaWxfZ3N2YV9oLAogICAgZXhjZWw9ImFuYWx5c2VzLzRfdHVtYWNvL0dTVkEvdF9lb3Npbm9waGlsX2dzdmFfaC54bHN4IikKYGBgCgojIyMjIFR1bWFjbyBNb25vY3l0ZXMgQzcKCmBgYHtyIGdzdmFfdHVtYWNvX21vbm9jeXRlX2M3fQp0X21vbm9jeXRlX2dzdmFfYzcgPC0gc2ltcGxlX2dzdmEoCiAgICB0X21vbm9jeXRlcywKICAgIHNpZ25hdHVyZXM9YnJvYWRfYzcsCiAgICBtc2lnX3htbD0icmVmZXJlbmNlL21zaWdkYi9tc2lnZGJfdjcuNS4xLnhtbCIsCiAgICBzaWduYXR1cmVfY2F0ZWdvcnk9ImM3IikKdF9tb25vY3l0ZV9nc3ZhX2M3X3NpZyA8LSBnZXRfc2lnX2dzdmFfY2F0ZWdvcmllcygKICAgIHRfbW9ub2N5dGVfZ3N2YV9jNywKICAgIGV4Y2VsPSJhbmFseXNlcy80X3R1bWFjby9HU1ZBL3RfbW9ub2N5dGVfZ3N2YV9jNy54bHN4IikKYGBgCgojIyMjIFR1bWFjbyBNb25vY3l0ZXMgQzIKCmBgYHtyIGdzdmFfdHVtYWNvX21vbm9jeXRlX2MyfQp0X21vbm9jeXRlX2dzdmFfYzIgPC0gc2ltcGxlX2dzdmEoCiAgICB0X21vbm9jeXRlcywKICAgIHNpZ25hdHVyZXM9YnJvYWRfYzIsCiAgICBtc2lnX3htbD0icmVmZXJlbmNlL21zaWdkYi9tc2lnZGJfdjcuNS4xLnhtbCIsCiAgICBzaWduYXR1cmVfY2F0ZWdvcnk9ImMyIikKdF9tb25vY3l0ZV9nc3ZhX2MyX3NpZyA8LSBnZXRfc2lnX2dzdmFfY2F0ZWdvcmllcygKICAgIHRfbW9ub2N5dGVfZ3N2YV9jMiwKICAgIGV4Y2VsPSJhbmFseXNlcy80X3R1bWFjby9HU1ZBL3RfbW9ub2N5dGVfZ3N2YV9jMi54bHN4IikKYGBgCgojIyMjIFR1bWFjbyBNb25vY3l0ZXMgSAoKYGBge3IgZ3N2YV9tb25vY3l0ZV90dW1hY29faH0KdF9tb25vY3l0ZV9nc3ZhX2ggPC0gc2ltcGxlX2dzdmEoCiAgICB0X21vbm9jeXRlcywKICAgIHNpZ25hdHVyZXM9YnJvYWRfaCwKICAgIG1zaWdfeG1sPSJyZWZlcmVuY2UvbXNpZ2RiL21zaWdkYl92Ny41LjEueG1sIiwKICAgIHNpZ25hdHVyZV9jYXRlZ29yeT0iaCIpCnRfbW9ub2N5dGVfZ3N2YV9oX3NpZyA8LSBnZXRfc2lnX2dzdmFfY2F0ZWdvcmllcygKICAgIHRfbW9ub2N5dGVfZ3N2YV9oLAogICAgZXhjZWw9ImFuYWx5c2VzLzRfdHVtYWNvL0dTVkEvdF9tb25vY3l0ZV9nc3ZhX2gueGxzeCIpCmBgYAoKIyMjIyBUdW1hY28gTmV1dHJvcGhpbHMgYzcKCmBgYHtyIGdzdmFfbmV1dHJvcGhpbHNfdHVtYWNvX2M3fQp0X25ldXRyb3BoaWxfZ3N2YV9jNyA8LSBzaW1wbGVfZ3N2YSgKICAgIHRfbmV1dHJvcGhpbHMsCiAgICBzaWduYXR1cmVzPWJyb2FkX2M3LAogICAgbXNpZ194bWw9InJlZmVyZW5jZS9tc2lnZGIvbXNpZ2RiX3Y3LjUuMS54bWwiLAogICAgc2lnbmF0dXJlX2NhdGVnb3J5PSJjNyIpCnRfbmV1dHJvcGhpbF9nc3ZhX2M3X3NpZyA8LSBnZXRfc2lnX2dzdmFfY2F0ZWdvcmllcygKICAgIHRfbmV1dHJvcGhpbF9nc3ZhX2M3LAogICAgZXhjZWw9ImFuYWx5c2VzLzRfdHVtYWNvL0dTVkEvdF9uZXV0cm9waGlsX2dzdmFfYzcueGxzeCIpCmBgYAoKIyMjIyBUdW1hY28gTmV1dHJvcGhpbHMgYzIKCmBgYHtyIGdzdmFfdHVtYWNvX25ldXRyb3BoaWxfYzJ9CnRfbmV1dHJvcGhpbF9nc3ZhX2MyIDwtIHNpbXBsZV9nc3ZhKAogICAgdF9uZXV0cm9waGlscywKICAgIHNpZ25hdHVyZXM9YnJvYWRfYzIsCiAgICBtc2lnX3htbD0icmVmZXJlbmNlL21zaWdkYi9tc2lnZGJfdjcuNS4xLnhtbCIsCiAgICBzaWduYXR1cmVfY2F0ZWdvcnk9ImMyIikKdF9uZXV0cm9waGlsX2dzdmFfYzJfc2lnIDwtIGdldF9zaWdfZ3N2YV9jYXRlZ29yaWVzKAogICAgdF9uZXV0cm9waGlsX2dzdmFfYzIsCiAgICBleGNlbD0iYW5hbHlzZXMvNF90dW1hY28vR1NWQS90X25ldXRyb3BoaWxfZ3N2YV9jMi54bHN4IikKYGBgCgojIyMjIFR1bWFjbyBOZXV0cm9waGlscyBICgpgYGB7ciB0dW1hY29fZ3N2YV9uZXV0cm9waGlsX2h9CnRfbmV1dHJvcGhpbF9nc3ZhX2ggPC0gc2ltcGxlX2dzdmEoCiAgICB0X25ldXRyb3BoaWxzLAogICAgc2lnbmF0dXJlcz1icm9hZF9oLAogICAgbXNpZ194bWw9InJlZmVyZW5jZS9tc2lnZGIvbXNpZ2RiX3Y3LjUuMS54bWwiLAogICAgc2lnbmF0dXJlX2NhdGVnb3J5PSJoIikKdF9uZXV0cm9waGlsX2dzdmFfaF9zaWcgPC0gZ2V0X3NpZ19nc3ZhX2NhdGVnb3JpZXMoCiAgICB0X25ldXRyb3BoaWxfZ3N2YV9oLAogICAgZXhjZWw9ImFuYWx5c2VzLzRfdHVtYWNvL0dTVkEvdF9uZXV0cm9waGlsX2dzdmFfaC54bHN4IikKYGBgCgpgYGB7ciBzYXZlbWUsIGV2YWw9RkFMU0V9CmlmICghaXNUUlVFKGdldDAoInNraXBfbG9hZCIpKSkgewogIHBhbmRlcjo6cGFuZGVyKHNlc3Npb25JbmZvKCkpCiAgbWVzc2FnZShwYXN0ZTAoIlRoaXMgaXMgaHBnbHRvb2xzIGNvbW1pdDogIiwgZ2V0X2dpdF9jb21taXQoKSkpCiAgbWVzc2FnZShwYXN0ZTAoIlNhdmluZyB0byAiLCBzYXZlZmlsZSkpCiAgdG1wIDwtIHNtKHNhdmVtZShmaWxlbmFtZT1zYXZlZmlsZSkpCn0KYGBgCgpgYGB7ciBsb2FkbWVfYWZ0ZXIsIGV2YWw9RkFMU0V9CnRtcCA8LSBsb2FkbWUoZmlsZW5hbWU9c2F2ZWZpbGUpCmBgYAo=