1 Annotation
We take the annotation data from ensembl’s biomart instance. The genome which was used to map the data was hg38 revision 100. My default when using biomart is to load the data from 1 year before the current date.
hs_annot <- sm(load_biomart_annotations(year = "2020"))
hs_annot <- hs_annot[["annotation"]]
hs_annot[["transcript"]] <- paste0(rownames(hs_annot), ".", hs_annot[["version"]])
rownames(hs_annot) <- make.names(hs_annot[["ensembl_gene_id"]], unique = TRUE)
tx_gene_map <- hs_annot[, c("transcript", "ensembl_gene_id")]
summary(hs_annot)## ensembl_transcript_id ensembl_gene_id version transcript_version
## Length:227921 Length:227921 Min. : 1.0 Min. : 1.00
## Class :character Class :character 1st Qu.: 6.0 1st Qu.: 1.00
## Mode :character Mode :character Median :12.0 Median : 1.00
## Mean :10.7 Mean : 3.08
## 3rd Qu.:16.0 3rd Qu.: 5.00
## Max. :29.0 Max. :17.00
##
## hgnc_symbol description gene_biotype cds_length
## Length:227921 Length:227921 Length:227921 Min. : 3
## Class :character Class :character Class :character 1st Qu.: 357
## Mode :character Mode :character Mode :character Median : 694
## Mean : 1139
## 3rd Qu.: 1446
## Max. :107976
## NA's :127343
## chromosome_name strand start_position end_position
## Length:227921 Length:227921 Min. :5.77e+02 Min. :6.47e+02
## Class :character Class :character 1st Qu.:3.11e+07 1st Qu.:3.12e+07
## Mode :character Mode :character Median :6.04e+07 Median :6.06e+07
## Mean :7.41e+07 Mean :7.42e+07
## 3rd Qu.:1.09e+08 3rd Qu.:1.09e+08
## Max. :2.49e+08 Max. :2.49e+08
##
## transcript
## Length:227921
## Class :character
## Mode :character
##
##
##
## hs_go <- sm(load_biomart_go()[["go"]])
hs_length <- hs_annot[, c("ensembl_gene_id", "cds_length")]
colnames(hs_length) <- c("ID", "length")2 Introduction
This document is intended to provide an overview of TMRC3 samples which have been sequenced. It includes some plots and analyses showing the relationships among the samples as well as some differential analyses when possible.
3 Sample Estimation
3.1 Generate expressionsets
The sample sheet is copied from our shared online sheet and updated with each release of sequencing data.
samplesheet <- "sample_sheets/tmrc3_samples_20210601.xlsx"3.1.1 Hisat2 expressionsets
The first thing to note is the large range in coverage. There are multiple samples with coverage which is too low to use. These will be removed shortly.
In the following block I immediately exclude any non-coding reads as well.
## Create the expressionset and immediately pass it to a filter
## removing the non protein coding genes.
sanitize_columns <- c("visitnumber", "clinicaloutcome", "donor",
"typeofcells", "clinicalpresentation",
"condition", "batch")
hs_expt <- create_expt(samplesheet,
file_column = "hg38100hisatfile",
savefile = glue::glue("rda/hs_expt_all-v{ver}.rda"),
gene_info = hs_annot) %>%
exclude_genes_expt(column = "gene_biotype", method = "keep",
pattern = "protein_coding", meta_column = "ncrna_lost") %>%
sanitize_expt_metadata(columns = sanitize_columns) %>%
set_expt_factors(columns = sanitize_columns, class = "factor")## Reading the sample metadata.## Dropped 71 rows from the sample metadata because they were blank.## The sample definitions comprises: 173 rows(samples) and 75 columns(metadata fields).## Warning in create_expt(samplesheet, file_column = "hg38100hisatfile", savefile =
## glue::glue("rda/hs_expt_all-v{ver}.rda"), : Some samples were removed when cross
## referencing the samples against the count data.## Matched 21452 annotations and counts.## Bringing together the count matrix and gene information.## Some annotations were lost in merging, setting them to 'undefined'.## The final expressionset has 21481 rows and 138 columns.## Before removal, there were 21481 genes, now there are 19941.## There are 13 samples which kept less than 90 percent counts.## TMRC30015 TMRC30017 TMRC30019 TMRC30044 TMRC30045 TMRC30097 TMRC30075 TMRC30087
## 79.24 85.72 89.75 80.34 73.33 89.90 86.97 83.63
## TMRC30101 TMRC30104 TMRC30114 TMRC30131 TMRC30073
## 88.41 80.29 87.62 86.82 89.26levels(pData(hs_expt[["expressionset"]])[["visitnumber"]]) <- list(
'0' = "notapplicable", '1' = 1, '2' = 2, '3' = 3)Split this data into CDS and lncRNA. Oh crap in order to do that I need to recount the data. Running now (20210518)
## lnc_expt <- create_expt(samplesheet,
## file_column = "hg38100lncfile",
## gene_info = hs_annot)3.1.1.1 Initial metrics
Once the data was loaded, there are a couple of metrics which may be plotted immediately.
nonzero <- plot_nonzero(hs_expt)
nonzero$plot## Warning: ggrepel: 111 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
ncrna_lost_df <- as.data.frame(pData(hs_expt)[["ncrna_lost"]])
rownames(ncrna_lost_df) <- rownames(pData(hs_expt))
colnames(ncrna_lost_df) <- "ncrna_lost"
tmpdf <- merge(nonzero$table, ncrna_lost_df, by = "row.names")
rownames(tmpdf) <- tmpdf[["Row.names"]]
tmpdf[["Row.names"]] <- NULL
ggplot(tmpdf, aes(x=ncrna_lost, y=nonzero_genes)) +
ggplot2::geom_point() +
ggplot2::ggtitle("Nonzero genes with respect to percent counts
lost when ncRNA was removed.")
Najib doesn’t want this plot, but I am using it to check new samples, so will hide it from general use.
libsize <- plot_libsize(hs_expt)
libsize$plot
3.2 Minimum coverage sample filtering
I arbitrarily chose 11,000 non-zero genes as a minimum. We may want this to be higher.
hs_valid <- subset_expt(hs_expt, nonzero = 11000)## The samples (and read coverage) removed when filtering 11000 non-zero genes are:## TMRC30010 TMRC30050 TMRC30052
## 52471 808149 3087347## subset_expt(): There were 138, now there are 135 samples.valid_write <- sm(write_expt(hs_valid, excel = glue("excel/hs_valid-v{ver}.xlsx")))4 Project Aims
The project seeks to determine the relationship of the innate immune response and inflammatory signaling to the clinical outcome of antileishmanial drug treatment. We will test the hypothesis that the profile of innate immune cell activation and their dynamics through the course of treatment differ between CL patients with prospectively determined therapeutic cure or failure.
This will be achieved through the characterization of the in vivo dynamics of blood-derived monocyte, neutrophil and eosinophil transcriptome before, during and at the end of treatment in CL patients. Cell-type specific transcriptomes, composite signatures and time-response expression profiles will be contrasted among patients with therapeutic cure or failure.
4.1 Preparation
To address these, I added to the end of the sample sheet columns named ‘condition’, ‘batch’, ‘donor’, and ‘time’. These are filled in with shorthand values according to the above.
4.2 Global view
Before addressing the questions explicitly by subsetting the data, I want to get a look at the samples as they are.
new_names <- pData(hs_valid)[["samplename"]]
hs_valid <- hs_valid %>%
set_expt_batches(fact = "cellssource") %>%
set_expt_conditions(fact = "typeofcells") %>%
set_expt_samplenames(newnames = new_names)
all_norm <- sm(normalize_expt(hs_valid, transform = "log2", norm = "quant",
convert = "cpm", filter = TRUE))
all_pca <- plot_pca(all_norm, plot_labels = FALSE,
plot_title = "PCA - Cell type", size_column = "visitnumber")
pp(file = glue("images/tmrc3_pca_nolabels-v{ver}.png"), image = all_pca$plot)
write.csv(all_pca$table, file = "coords/hs_donor_pca_coords.csv")
plot_corheat(all_norm, plot_title = "Heirarchical clustering:
cell types")$plot
4.3 Examine samples relevant to clinical outcome
Now let us consider only the samples for which we have a clinical outcome. These fall primarily into either ‘cured’ or ‘failed’, but some people have not yet returned to the clinic after the first or second visit. These are deemed ‘lost’.
hs_clinical <- hs_valid %>%
set_expt_conditions(fact = "clinicaloutcome") %>%
set_expt_batches(fact = "typeofcells") %>%
subset_expt(subset = "typeofcells!='pbmcs'&typeofcells!='macrophages'")## subset_expt(): There were 135, now there are 117 samples.chosen_colors <- c("#D95F02", "#7570B3", "#1B9E77", "#FF0000", "#FF0000")
names(chosen_colors) <- c("cure", "failure", "lost", "null", "notapplicable")
hs_clinical <- set_expt_colors(hs_clinical, colors = chosen_colors)
newnames <- make.names(pData(hs_clinical)[["samplename"]], unique = TRUE)
hs_clinical <- set_expt_samplenames(hs_clinical, newnames = newnames)
hs_clinical_norm <- sm(normalize_expt(hs_clinical, filter = TRUE, transform = "log2",
convert = "cpm", norm = "quant"))
clinical_pca <- plot_pca(hs_clinical_norm, plot_labels = FALSE,
size_column = "visitnumber", cis = NULL,
plot_title = "PCA - clinical samples")
pp(file = glue("images/all_clinical_nobatch_pca-v{ver}.png"), image = clinical_pca$plot,
height = 8, width = 20)
4.3.1 Repeat without the biopsy samples
hs_clinical_nobiop <- hs_clinical %>%
subset_expt(subset = "typeofcells!='biopsy'")## subset_expt(): There were 117, now there are 75 samples.hs_clinical_nobiop_norm <- sm(normalize_expt(hs_clinical_nobiop, filter = TRUE, transform = "log2",
convert = "cpm", norm = "quant"))
clinical_nobiop_pca <- plot_pca(hs_clinical_nobiop_norm, plot_labels = FALSE, cis = NULL,
plot_title = "PCA - clinical samples without biopsies")
pp(file = glue("images/all_clinical_nobiop_nobatch_pca-v{ver}.png"),
image = clinical_nobiop_pca$plot)
4.3.2 Attempt to correct for the surrogate variables
At this time we have two primary data structures of interest: hs_clinical and hs_clinical_nobiop
hs_clinical_nb <- normalize_expt(hs_clinical, filter = TRUE, batch = "svaseq",
transform = "log2", convert = "cpm")## Removing 5312 low-count genes (14629 remaining).## batch_counts: Before batch/surrogate estimation, 108221 entries are x==0: 6%.## batch_counts: Before batch/surrogate estimation, 314017 entries are 0<x<1: 18%.## Setting 22060 low elements to zero.## transform_counts: Found 22060 values equal to 0, adding 1 to the matrix.clinical_batch_pca <- plot_pca(hs_clinical_nb, plot_labels = FALSE, cis = NULL,
size_column = "visitnumber", plot_title = "PCA - clinical samples")
clinical_batch_pca$plot
hs_clinical_nobiop_nb <- sm(normalize_expt(hs_clinical_nobiop, filter = TRUE, batch = "svaseq",
transform = "log2", convert = "cpm"))
clinical_nobiop_batch_pca <- plot_pca(hs_clinical_nobiop_nb,
plot_title = "PCA - clinical samples without biopsies",
plot_labels = FALSE)
pp(file = "images/clinical_batch.png", image = clinical_nobiop_batch_pca$plot)
test <- plot_pca(hs_clinical_nobiop_nb, size_column = "visitnumber",
plot_title = "PCA - clinical samples without biopsies",
plot_labels = FALSE)
test$plot
clinical_nobiop_batch_tsne <- plot_tsne(hs_clinical_nobiop_nb,
plot_title = "tSNE - clinical samples without biopsies",
plot_labels = FALSE)
clinical_nobiop_batch_tsne$plot
4.3.2.1 Look at remaining variance with variancePartition
test <- simple_varpart(hs_clinical_nobiop)##
## Total:100 stest$partition_plot
4.4 Perform DE of the clinical samples cure vs. fail
individual_celltypes <- subset_expt(hs_clinical_nobiop, subset="condition!='lost'")## subset_expt(): There were 75, now there are 60 samples.hs_clinic_de <- sm(all_pairwise(individual_celltypes, model_batch = "svaseq", filter = TRUE))
hs_clinic_table <- sm(combine_de_tables(
hs_clinic_de,
excel = glue::glue("excel/individual_celltypes_table-v{ver}.xlsx")))
hs_clinic_sig <- sm(extract_significant_genes(
hs_clinic_table,
excel = glue::glue("excel/individual_celltypes_sig-v{ver}.xlsx")))
hs_clinic_sig[["summary_df"]]## limma_V1 limma_V2 edger_V1 edger_V2 deseq_V1 deseq_V2 ebseq_V1 ebseq_V2
## 1 169 207 260 289 227 323 63 194
## basic_V1 basic_V2
## 1 34 9hs_clinic_de[["comparison"]][["heat"]]## NULL4.4.1 Perform LRT with the clinical samples
I am not sure if we have enough samples across the three visit to completely work as well as we would like, but there is only 1 way to find out! Now that I think about it, one thing which might be awesome is to use cell type as an interacting factor…
4.4.1.1 With biopsy samples
I figure this might be a place where the biopsy samples might prove useful.
clinical_nolost <- subset_expt(hs_clinical, subset="condition!='lost'")## subset_expt(): There were 117, now there are 100 samples.lrt_visit_clinical_test <- deseq_lrt(clinical_nolost, transform = "vst",
interactor_column = "visitnumber",
interest_column = "clinicaloutcome")## 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 231 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)## estimating dispersions## fitting model and testing## Working with 2 genes.## Warning: `distinct_()` was deprecated in dplyr 0.7.0.
## Please use `distinct()` instead.
## See vignette('programming') for more help## Working with 0 genes after filtering: minc > 3## Error: object 'cutoff' not foundlrt_visit_clinical_test[["favorite_genes"]]## Error in eval(expr, envir, enclos): object 'lrt_visit_clinical_test' not foundlrt_celltype_clinical_test <- deseq_lrt(clinical_nolost, transform = "vst",
interactor_column = "typeofcells",
interest_column = "clinicaloutcome")## 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 57 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)## estimating dispersions## fitting model and testing## Working with 882 genes.## Working with 880 genes after filtering: minc > 3## Joining, by = "merge"
## Joining, by = "merge"
lrt_celltype_clinical_test[["favorite_genes"]]## genes cluster
## ENSG00000001497 ENSG00000001497 1
## ENSG00000003137 ENSG00000003137 2
## ENSG00000004455 ENSG00000004455 3
## ENSG00000004846 ENSG00000004846 4
## ENSG00000005020 ENSG00000005020 5
## ENSG00000005175 ENSG00000005175 6
## ENSG00000005486 ENSG00000005486 7
## ENSG00000005801 ENSG00000005801 1
## ENSG00000007047 ENSG00000007047 5
## ENSG00000007392 ENSG00000007392 8
## ENSG00000008056 ENSG00000008056 9
## ENSG00000008086 ENSG00000008086 10
## ENSG00000010270 ENSG00000010270 11
## ENSG00000010278 ENSG00000010278 8
## ENSG00000010704 ENSG00000010704 11
## ENSG00000010803 ENSG00000010803 12
## ENSG00000011028 ENSG00000011028 2
## ENSG00000011426 ENSG00000011426 4
## ENSG00000013306 ENSG00000013306 13
## ENSG00000013583 ENSG00000013583 14
## ENSG00000017427 ENSG00000017427 15
## ENSG00000019144 ENSG00000019144 2
## ENSG00000020256 ENSG00000020256 8
## ENSG00000020577 ENSG00000020577 2
## ENSG00000020633 ENSG00000020633 16
## ENSG00000021355 ENSG00000021355 9
## ENSG00000028203 ENSG00000028203 1
## ENSG00000030066 ENSG00000030066 17
## ENSG00000032444 ENSG00000032444 16
## ENSG00000035141 ENSG00000035141 3
## ENSG00000035720 ENSG00000035720 17
## ENSG00000038945 ENSG00000038945 2
## ENSG00000039560 ENSG00000039560 2
## ENSG00000042088 ENSG00000042088 1
## ENSG00000044459 ENSG00000044459 13
## ENSG00000048828 ENSG00000048828 5
## ENSG00000050820 ENSG00000050820 18
## ENSG00000055130 ENSG00000055130 13
## ENSG00000057704 ENSG00000057704 19
## ENSG00000060138 ENSG00000060138 14
## ENSG00000060982 ENSG00000060982 2
## ENSG00000061918 ENSG00000061918 3
## ENSG00000062038 ENSG00000062038 4
## ENSG00000062598 ENSG00000062598 5
## ENSG00000064763 ENSG00000064763 10
## ENSG00000064933 ENSG00000064933 3
## ENSG00000065150 ENSG00000065150 1
## ENSG00000065809 ENSG00000065809 7
## ENSG00000065989 ENSG00000065989 12
## ENSG00000066427 ENSG00000066427 5
## ENSG00000066455 ENSG00000066455 15
## ENSG00000066651 ENSG00000066651 2
## ENSG00000067533 ENSG00000067533 1
## ENSG00000068305 ENSG00000068305 16
## ENSG00000069020 ENSG00000069020 15
## ENSG00000069345 ENSG00000069345 7
## ENSG00000069998 ENSG00000069998 13
## ENSG00000070269 ENSG00000070269 10
## ENSG00000072657 ENSG00000072657 20
## ENSG00000073150 ENSG00000073150 5
## ENSG00000073350 ENSG00000073350 8
## ENSG00000073417 ENSG00000073417 12
## ENSG00000073737 ENSG00000073737 10
## ENSG00000073754 ENSG00000073754 2
## ENSG00000074211 ENSG00000074211 4
## ENSG00000074657 ENSG00000074657 2
## ENSG00000074935 ENSG00000074935 3
## ENSG00000075643 ENSG00000075643 1
## ENSG00000076242 ENSG00000076242 3
## ENSG00000076716 ENSG00000076716 2
## ENSG00000077616 ENSG00000077616 17
## ENSG00000077782 ENSG00000077782 2
## ENSG00000077935 ENSG00000077935 13
## ENSG00000078098 ENSG00000078098 15
## ENSG00000078589 ENSG00000078589 22
## ENSG00000078596 ENSG00000078596 3
## ENSG00000079215 ENSG00000079215 1
## ENSG00000079482 ENSG00000079482 1
## ENSG00000079616 ENSG00000079616 8
## ENSG00000081386 ENSG00000081386 3
## ENSG00000081913 ENSG00000081913 9
## ENSG00000082516 ENSG00000082516 3
## ENSG00000083457 ENSG00000083457 12
## ENSG00000083828 ENSG00000083828 19
## ENSG00000084207 ENSG00000084207 13
## ENSG00000085871 ENSG00000085871 7
## ENSG00000085999 ENSG00000085999 1
## ENSG00000086730 ENSG00000086730 9
## ENSG00000087116 ENSG00000087116 1
## ENSG00000087269 ENSG00000087269 1
## ENSG00000088543 ENSG00000088543 3
## ENSG00000088726 ENSG00000088726 3
## ENSG00000088826 ENSG00000088826 1
## ENSG00000088827 ENSG00000088827 14
## ENSG00000088992 ENSG00000088992 16
## ENSG00000089012 ENSG00000089012 18
## ENSG00000089127 ENSG00000089127 14
## ENSG00000089818 ENSG00000089818 5
## ENSG00000090013 ENSG00000090013 14
## ENSG00000090674 ENSG00000090674 7
## ENSG00000091409 ENSG00000091409 18
## ENSG00000091972 ENSG00000091972 15
## ENSG00000092445 ENSG00000092445 1
## ENSG00000094841 ENSG00000094841 1
## ENSG00000095002 ENSG00000095002 18
## ENSG00000095585 ENSG00000095585 2
## ENSG00000099260 ENSG00000099260 18
## ENSG00000099783 ENSG00000099783 16
## ENSG00000100036 ENSG00000100036 13
## ENSG00000100065 ENSG00000100065 3
## ENSG00000100124 ENSG00000100124 15
## ENSG00000100138 ENSG00000100138 1
## ENSG00000100281 ENSG00000100281 16
## ENSG00000100288 ENSG00000100288 1
## ENSG00000100292 ENSG00000100292 14
## ENSG00000100335 ENSG00000100335 23
## ENSG00000100354 ENSG00000100354 10
## ENSG00000100376 ENSG00000100376 13
## ENSG00000100413 ENSG00000100413 23
## ENSG00000100558 ENSG00000100558 1
## ENSG00000100596 ENSG00000100596 16
## ENSG00000100767 ENSG00000100767 6
## ENSG00000100842 ENSG00000100842 15
## ENSG00000101000 ENSG00000101000 4
## ENSG00000101096 ENSG00000101096 17
## ENSG00000101188 ENSG00000101188 11
## ENSG00000101189 ENSG00000101189 23
## ENSG00000101255 ENSG00000101255 9
## ENSG00000101384 ENSG00000101384 3
## ENSG00000101665 ENSG00000101665 23
## ENSG00000101844 ENSG00000101844 18
## ENSG00000102007 ENSG00000102007 5
## ENSG00000102048 ENSG00000102048 1
## ENSG00000102221 ENSG00000102221 1
## ENSG00000102290 ENSG00000102290 13
## ENSG00000102384 ENSG00000102384 18
## ENSG00000102409 ENSG00000102409 3
## ENSG00000102572 ENSG00000102572 5
## ENSG00000102738 ENSG00000102738 18
## ENSG00000102967 ENSG00000102967 11
## ENSG00000102978 ENSG00000102978 9
## ENSG00000103047 ENSG00000103047 13
## ENSG00000103056 ENSG00000103056 17
## ENSG00000103257 ENSG00000103257 4
## ENSG00000103274 ENSG00000103274 13
## ENSG00000103319 ENSG00000103319 8
## ENSG00000103740 ENSG00000103740 18
## ENSG00000104133 ENSG00000104133 10
## ENSG00000104213 ENSG00000104213 2
## ENSG00000104218 ENSG00000104218 1
## ENSG00000104321 ENSG00000104321 15
## ENSG00000104689 ENSG00000104689 13
## ENSG00000104738 ENSG00000104738 3
## ENSG00000104936 ENSG00000104936 2
## ENSG00000104969 ENSG00000104969 16
## ENSG00000104980 ENSG00000104980 13
## ENSG00000104998 ENSG00000104998 13
## ENSG00000105146 ENSG00000105146 4
## ENSG00000105220 ENSG00000105220 12
## ENSG00000105245 ENSG00000105245 20
## ENSG00000105278 ENSG00000105278 6
## ENSG00000105281 ENSG00000105281 13
## ENSG00000105366 ENSG00000105366 17
## ENSG00000105499 ENSG00000105499 2
## ENSG00000105576 ENSG00000105576 3
## ENSG00000105612 ENSG00000105612 14
## ENSG00000105649 ENSG00000105649 10
## ENSG00000105656 ENSG00000105656 5
## ENSG00000105668 ENSG00000105668 15
## ENSG00000105928 ENSG00000105928 2
## ENSG00000105967 ENSG00000105967 14
## ENSG00000105997 ENSG00000105997 6
## ENSG00000106100 ENSG00000106100 22
## ENSG00000106366 ENSG00000106366 1
## ENSG00000106617 ENSG00000106617 16
## ENSG00000106635 ENSG00000106635 16
## ENSG00000106638 ENSG00000106638 13
## ENSG00000106804 ENSG00000106804 14
## ENSG00000107130 ENSG00000107130 4
## ENSG00000107185 ENSG00000107185 5
## ENSG00000107281 ENSG00000107281 1
## ENSG00000107438 ENSG00000107438 8
## ENSG00000107485 ENSG00000107485 1
## ENSG00000107719 ENSG00000107719 2
## ENSG00000107731 ENSG00000107731 3
## ENSG00000107816 ENSG00000107816 1
## ENSG00000107819 ENSG00000107819 23
## ENSG00000108219 ENSG00000108219 16
## ENSG00000108389 ENSG00000108389 10
## ENSG00000108821 ENSG00000108821 15
## ENSG00000108963 ENSG00000108963 13
## ENSG00000109113 ENSG00000109113 13
## ENSG00000109182 ENSG00000109182 4
## ENSG00000109685 ENSG00000109685 8
## ENSG00000109917 ENSG00000109917 23
## ENSG00000109919 ENSG00000109919 1
## ENSG00000109971 ENSG00000109971 2
## ENSG00000110092 ENSG00000110092 3
## ENSG00000110171 ENSG00000110171 23
## ENSG00000110200 ENSG00000110200 12
## ENSG00000110455 ENSG00000110455 12
## ENSG00000110660 ENSG00000110660 1
## ENSG00000110665 ENSG00000110665 24
## ENSG00000110851 ENSG00000110851 17
## ENSG00000111144 ENSG00000111144 11
## ENSG00000111145 ENSG00000111145 9
## ENSG00000111331 ENSG00000111331 11
## ENSG00000111335 ENSG00000111335 14
## ENSG00000111424 ENSG00000111424 9
## ENSG00000111640 ENSG00000111640 14
## ENSG00000111799 ENSG00000111799 15
## ENSG00000111801 ENSG00000111801 2
## ENSG00000111850 ENSG00000111850 6
## ENSG00000112033 ENSG00000112033 5
## ENSG00000112130 ENSG00000112130 1
## ENSG00000112208 ENSG00000112208 18
## ENSG00000112299 ENSG00000112299 11
## ENSG00000112715 ENSG00000112715 16
## ENSG00000113273 ENSG00000113273 24
## ENSG00000113360 ENSG00000113360 3
## ENSG00000113649 ENSG00000113649 13
## ENSG00000113716 ENSG00000113716 12
## ENSG00000114439 ENSG00000114439 10
## ENSG00000114735 ENSG00000114735 8
## ENSG00000114737 ENSG00000114737 16
## ENSG00000114861 ENSG00000114861 6
## ENSG00000115155 ENSG00000115155 11
## ENSG00000115414 ENSG00000115414 2
## ENSG00000115423 ENSG00000115423 1
## ENSG00000115598 ENSG00000115598 2
## ENSG00000115648 ENSG00000115648 15
## ENSG00000115718 ENSG00000115718 9
## ENSG00000115738 ENSG00000115738 23
## ENSG00000115761 ENSG00000115761 1
## ENSG00000115762 ENSG00000115762 9
## ENSG00000115904 ENSG00000115904 18
## ENSG00000115919 ENSG00000115919 14
## ENSG00000115963 ENSG00000115963 15
## ENSG00000116127 ENSG00000116127 8
## ENSG00000116132 ENSG00000116132 15
## ENSG00000116141 ENSG00000116141 3
## ENSG00000116157 ENSG00000116157 3
## ENSG00000116478 ENSG00000116478 24
## ENSG00000116514 ENSG00000116514 5
## ENSG00000116678 ENSG00000116678 1
## ENSG00000116688 ENSG00000116688 5
## ENSG00000116874 ENSG00000116874 1
## ENSG00000116898 ENSG00000116898 2
## ENSG00000116984 ENSG00000116984 8
## ENSG00000117226 ENSG00000117226 2
## ENSG00000117228 ENSG00000117228 20
## ENSG00000117245 ENSG00000117245 15
## ENSG00000117298 ENSG00000117298 5
## ENSG00000117448 ENSG00000117448 1
## ENSG00000117481 ENSG00000117481 1
## ENSG00000117501 ENSG00000117501 15
## ENSG00000117597 ENSG00000117597 1
## ENSG00000117724 ENSG00000117724 3
## ENSG00000118004 ENSG00000118004 4
## ENSG00000118680 ENSG00000118680 6
## ENSG00000118785 ENSG00000118785 15
## ENSG00000118849 ENSG00000118849 18
## ENSG00000118946 ENSG00000118946 18
## ENSG00000119121 ENSG00000119121 19
## ENSG00000119698 ENSG00000119698 6
## ENSG00000119865 ENSG00000119865 15
## ENSG00000119915 ENSG00000119915 2
## ENSG00000120008 ENSG00000120008 13
## ENSG00000120053 ENSG00000120053 3
## ENSG00000120162 ENSG00000120162 2
## ENSG00000120254 ENSG00000120254 1
## ENSG00000120278 ENSG00000120278 2
## ENSG00000120675 ENSG00000120675 2
## ENSG00000120696 ENSG00000120696 6
## ENSG00000120868 ENSG00000120868 10
## ENSG00000121057 ENSG00000121057 1
## ENSG00000121067 ENSG00000121067 10
## ENSG00000121152 ENSG00000121152 3
## ENSG00000121210 ENSG00000121210 10
## ENSG00000121410 ENSG00000121410 13
## ENSG00000121807 ENSG00000121807 14
## ENSG00000121966 ENSG00000121966 24
## ENSG00000122483 ENSG00000122483 22
## ENSG00000122729 ENSG00000122729 6
## ENSG00000123219 ENSG00000123219 3
## ENSG00000123384 ENSG00000123384 14
## ENSG00000123607 ENSG00000123607 14
## ENSG00000123838 ENSG00000123838 11
## ENSG00000124006 ENSG00000124006 3
## ENSG00000124145 ENSG00000124145 1
## ENSG00000124216 ENSG00000124216 14
## ENSG00000124357 ENSG00000124357 5
## ENSG00000124615 ENSG00000124615 4
## ENSG00000124766 ENSG00000124766 1
## ENSG00000124780 ENSG00000124780 18
## ENSG00000124785 ENSG00000124785 2
## ENSG00000125247 ENSG00000125247 13
## ENSG00000125434 ENSG00000125434 12
## ENSG00000125650 ENSG00000125650 2
## ENSG00000125733 ENSG00000125733 2
## ENSG00000125735 ENSG00000125735 19
## ENSG00000125753 ENSG00000125753 5
## ENSG00000125779 ENSG00000125779 10
## ENSG00000125821 ENSG00000125821 1
## ENSG00000125863 ENSG00000125863 2
## ENSG00000125952 ENSG00000125952 24
## ENSG00000125968 ENSG00000125968 13
## ENSG00000126457 ENSG00000126457 1
## ENSG00000126870 ENSG00000126870 3
## ENSG00000128245 ENSG00000128245 13
## ENSG00000128274 ENSG00000128274 2
## ENSG00000128731 ENSG00000128731 8
## ENSG00000129071 ENSG00000129071 19
## ENSG00000129484 ENSG00000129484 1
## ENSG00000129493 ENSG00000129493 22
## ENSG00000129675 ENSG00000129675 24
## ENSG00000130167 ENSG00000130167 5
## ENSG00000130203 ENSG00000130203 1
## ENSG00000130208 ENSG00000130208 2
## ENSG00000130300 ENSG00000130300 2
## ENSG00000130479 ENSG00000130479 16
## ENSG00000130529 ENSG00000130529 1
## ENSG00000130559 ENSG00000130559 25
## ENSG00000130638 ENSG00000130638 1
## ENSG00000130649 ENSG00000130649 1
## ENSG00000130725 ENSG00000130725 5
## ENSG00000130768 ENSG00000130768 24
## ENSG00000131019 ENSG00000131019 1
## ENSG00000131042 ENSG00000131042 5
## ENSG00000131238 ENSG00000131238 14
## ENSG00000131409 ENSG00000131409 2
## ENSG00000131508 ENSG00000131508 9
## ENSG00000131828 ENSG00000131828 3
## ENSG00000132256 ENSG00000132256 20
## ENSG00000132305 ENSG00000132305 1
## ENSG00000132530 ENSG00000132530 11
## ENSG00000132746 ENSG00000132746 3
## ENSG00000132792 ENSG00000132792 12
## ENSG00000132965 ENSG00000132965 19
## ENSG00000133030 ENSG00000133030 8
## ENSG00000133103 ENSG00000133103 3
## ENSG00000133116 ENSG00000133116 2
## ENSG00000133997 ENSG00000133997 19
## ENSG00000134057 ENSG00000134057 18
## ENSG00000134107 ENSG00000134107 7
## ENSG00000134207 ENSG00000134207 13
## ENSG00000134245 ENSG00000134245 2
## ENSG00000134318 ENSG00000134318 22
## ENSG00000134333 ENSG00000134333 1
## ENSG00000134369 ENSG00000134369 1
## ENSG00000134453 ENSG00000134453 12
## ENSG00000134460 ENSG00000134460 17
## ENSG00000134463 ENSG00000134463 7
## ENSG00000134686 ENSG00000134686 5
## ENSG00000134697 ENSG00000134697 1
## ENSG00000134759 ENSG00000134759 3
## ENSG00000134809 ENSG00000134809 2
## ENSG00000134824 ENSG00000134824 3
## ENSG00000135047 ENSG00000135047 2
## ENSG00000135077 ENSG00000135077 14
## ENSG00000135083 ENSG00000135083 19
## ENSG00000135094 ENSG00000135094 13
## ENSG00000135116 ENSG00000135116 22
## ENSG00000135205 ENSG00000135205 10
## ENSG00000135245 ENSG00000135245 18
## ENSG00000135362 ENSG00000135362 24
## ENSG00000135363 ENSG00000135363 11
## ENSG00000135372 ENSG00000135372 1
## ENSG00000135503 ENSG00000135503 12
## ENSG00000135723 ENSG00000135723 5
## ENSG00000135823 ENSG00000135823 9
## ENSG00000135916 ENSG00000135916 8
## ENSG00000135924 ENSG00000135924 17
## ENSG00000135929 ENSG00000135929 11
## ENSG00000136011 ENSG00000136011 15
## ENSG00000136045 ENSG00000136045 3
## ENSG00000136068 ENSG00000136068 8
## ENSG00000136235 ENSG00000136235 2
## ENSG00000136319 ENSG00000136319 3
## ENSG00000136631 ENSG00000136631 1
## ENSG00000136717 ENSG00000136717 14
## ENSG00000136732 ENSG00000136732 2
## ENSG00000136807 ENSG00000136807 16
## ENSG00000136830 ENSG00000136830 13
## ENSG00000136840 ENSG00000136840 1
## ENSG00000136874 ENSG00000136874 20
## ENSG00000136930 ENSG00000136930 16
## ENSG00000136932 ENSG00000136932 7
## ENSG00000136938 ENSG00000136938 23
## ENSG00000137094 ENSG00000137094 2
## ENSG00000137166 ENSG00000137166 23
## ENSG00000137265 ENSG00000137265 2
## ENSG00000137312 ENSG00000137312 5
## ENSG00000137547 ENSG00000137547 1
## ENSG00000137571 ENSG00000137571 4
## ENSG00000137628 ENSG00000137628 20
## ENSG00000137672 ENSG00000137672 17
## ENSG00000137673 ENSG00000137673 15
## ENSG00000137807 ENSG00000137807 18
## ENSG00000137812 ENSG00000137812 3
## ENSG00000137959 ENSG00000137959 11
## ENSG00000137965 ENSG00000137965 11
## ENSG00000138061 ENSG00000138061 14
## ENSG00000138180 ENSG00000138180 3
## ENSG00000138246 ENSG00000138246 20
## ENSG00000138442 ENSG00000138442 3
## ENSG00000138646 ENSG00000138646 20
## ENSG00000138709 ENSG00000138709 1
## ENSG00000138760 ENSG00000138760 2
## ENSG00000138764 ENSG00000138764 10
## ENSG00000139112 ENSG00000139112 5
## ENSG00000139160 ENSG00000139160 6
## ENSG00000139211 ENSG00000139211 2
## ENSG00000139514 ENSG00000139514 1
## ENSG00000139597 ENSG00000139597 14
## ENSG00000139684 ENSG00000139684 1
## ENSG00000139722 ENSG00000139722 5
## ENSG00000140297 ENSG00000140297 15
## ENSG00000140525 ENSG00000140525 18
## ENSG00000140564 ENSG00000140564 7
## ENSG00000140577 ENSG00000140577 23
## ENSG00000140650 ENSG00000140650 2
## ENSG00000140859 ENSG00000140859 4
## ENSG00000141298 ENSG00000141298 5
## ENSG00000141574 ENSG00000141574 11
## ENSG00000141576 ENSG00000141576 1
## ENSG00000141837 ENSG00000141837 6
## ENSG00000142528 ENSG00000142528 20
## ENSG00000142621 ENSG00000142621 2
## ENSG00000142627 ENSG00000142627 1
## ENSG00000142910 ENSG00000142910 1
## ENSG00000143079 ENSG00000143079 2
## ENSG00000143127 ENSG00000143127 3
## ENSG00000143178 ENSG00000143178 10
## ENSG00000143416 ENSG00000143416 15
## ENSG00000143420 ENSG00000143420 19
## ENSG00000143458 ENSG00000143458 8
## ENSG00000143493 ENSG00000143493 3
## ENSG00000143498 ENSG00000143498 18
## ENSG00000143669 ENSG00000143669 10
## ENSG00000143801 ENSG00000143801 3
## ENSG00000143845 ENSG00000143845 1
## ENSG00000143847 ENSG00000143847 25
## ENSG00000143851 ENSG00000143851 7
## ENSG00000144331 ENSG00000144331 15
## ENSG00000144580 ENSG00000144580 1
## ENSG00000144815 ENSG00000144815 25
## ENSG00000144908 ENSG00000144908 18
## ENSG00000145040 ENSG00000145040 2
## ENSG00000145191 ENSG00000145191 1
## ENSG00000145244 ENSG00000145244 20
## ENSG00000145246 ENSG00000145246 2
## ENSG00000145247 ENSG00000145247 3
## ENSG00000145348 ENSG00000145348 3
## ENSG00000145362 ENSG00000145362 2
## ENSG00000145375 ENSG00000145375 3
## ENSG00000145416 ENSG00000145416 23
## ENSG00000145685 ENSG00000145685 14
## ENSG00000145721 ENSG00000145721 6
## ENSG00000145817 ENSG00000145817 2
## ENSG00000145945 ENSG00000145945 1
## ENSG00000146021 ENSG00000146021 17
## ENSG00000146205 ENSG00000146205 22
## ENSG00000146243 ENSG00000146243 15
## ENSG00000146281 ENSG00000146281 6
## ENSG00000146416 ENSG00000146416 6
## ENSG00000146918 ENSG00000146918 1
## ENSG00000147138 ENSG00000147138 3
## ENSG00000147174 ENSG00000147174 6
## ENSG00000147257 ENSG00000147257 1
## ENSG00000147408 ENSG00000147408 6
## ENSG00000147454 ENSG00000147454 5
## ENSG00000147614 ENSG00000147614 15
## ENSG00000147647 ENSG00000147647 1
## ENSG00000148219 ENSG00000148219 1
## ENSG00000148225 ENSG00000148225 18
## ENSG00000148248 ENSG00000148248 23
## ENSG00000148334 ENSG00000148334 13
## ENSG00000148335 ENSG00000148335 12
## ENSG00000148606 ENSG00000148606 1
## ENSG00000148690 ENSG00000148690 3
## ENSG00000148737 ENSG00000148737 11
## ENSG00000148814 ENSG00000148814 1
## ENSG00000149418 ENSG00000149418 23
## ENSG00000149599 ENSG00000149599 2
## ENSG00000149633 ENSG00000149633 18
## ENSG00000149679 ENSG00000149679 16
## ENSG00000149972 ENSG00000149972 6
## ENSG00000150048 ENSG00000150048 20
## ENSG00000150051 ENSG00000150051 2
## ENSG00000150347 ENSG00000150347 2
## ENSG00000150556 ENSG00000150556 4
## ENSG00000150756 ENSG00000150756 18
## ENSG00000150990 ENSG00000150990 13
## ENSG00000151150 ENSG00000151150 18
## ENSG00000151320 ENSG00000151320 2
## ENSG00000151353 ENSG00000151353 1
## ENSG00000151490 ENSG00000151490 1
## ENSG00000151503 ENSG00000151503 3
## ENSG00000151576 ENSG00000151576 11
## ENSG00000151689 ENSG00000151689 8
## ENSG00000151692 ENSG00000151692 24
## ENSG00000151693 ENSG00000151693 1
## ENSG00000151778 ENSG00000151778 3
## ENSG00000151789 ENSG00000151789 2
## ENSG00000151790 ENSG00000151790 4
## ENSG00000152056 ENSG00000152056 2
## ENSG00000152061 ENSG00000152061 10
## ENSG00000152137 ENSG00000152137 1
## ENSG00000152804 ENSG00000152804 16
## ENSG00000152952 ENSG00000152952 3
## ENSG00000153823 ENSG00000153823 13
## ENSG00000153976 ENSG00000153976 4
## ENSG00000153982 ENSG00000153982 18
## ENSG00000154240 ENSG00000154240 4
## ENSG00000154265 ENSG00000154265 17
## ENSG00000154277 ENSG00000154277 18
## ENSG00000154305 ENSG00000154305 10
## ENSG00000154447 ENSG00000154447 14
## ENSG00000154451 ENSG00000154451 20
## ENSG00000154589 ENSG00000154589 11
## ENSG00000154760 ENSG00000154760 17
## ENSG00000155016 ENSG00000155016 1
## ENSG00000155158 ENSG00000155158 2
## ENSG00000155189 ENSG00000155189 3
## ENSG00000155252 ENSG00000155252 14
## ENSG00000155275 ENSG00000155275 25
## ENSG00000155363 ENSG00000155363 11
## ENSG00000155380 ENSG00000155380 1
## ENSG00000155906 ENSG00000155906 8
## ENSG00000156030 ENSG00000156030 5
## ENSG00000156042 ENSG00000156042 10
## ENSG00000156049 ENSG00000156049 17
## ENSG00000156239 ENSG00000156239 3
## ENSG00000156398 ENSG00000156398 2
## ENSG00000156711 ENSG00000156711 19
## ENSG00000156802 ENSG00000156802 8
## ENSG00000156804 ENSG00000156804 1
## ENSG00000156970 ENSG00000156970 3
## ENSG00000157036 ENSG00000157036 1
## ENSG00000157227 ENSG00000157227 2
## ENSG00000157456 ENSG00000157456 1
## ENSG00000157654 ENSG00000157654 3
## ENSG00000157657 ENSG00000157657 8
## ENSG00000157933 ENSG00000157933 7
## ENSG00000157985 ENSG00000157985 2
## ENSG00000158062 ENSG00000158062 14
## ENSG00000158373 ENSG00000158373 6
## ENSG00000158402 ENSG00000158402 4
## ENSG00000158406 ENSG00000158406 6
## ENSG00000158710 ENSG00000158710 5
## ENSG00000158715 ENSG00000158715 12
## ENSG00000158856 ENSG00000158856 8
## ENSG00000159216 ENSG00000159216 12
## ENSG00000159261 ENSG00000159261 4
## ENSG00000159307 ENSG00000159307 3
## ENSG00000160094 ENSG00000160094 12
## ENSG00000160113 ENSG00000160113 3
## ENSG00000160789 ENSG00000160789 1
## ENSG00000160791 ENSG00000160791 2
## ENSG00000160856 ENSG00000160856 2
## ENSG00000160932 ENSG00000160932 11
## ENSG00000161640 ENSG00000161640 16
## ENSG00000161835 ENSG00000161835 12
## ENSG00000161960 ENSG00000161960 3
## ENSG00000162104 ENSG00000162104 13
## ENSG00000162129 ENSG00000162129 13
## ENSG00000162139 ENSG00000162139 13
## ENSG00000162367 ENSG00000162367 8
## ENSG00000162377 ENSG00000162377 1
## ENSG00000162390 ENSG00000162390 8
## ENSG00000162627 ENSG00000162627 2
## ENSG00000162669 ENSG00000162669 9
## ENSG00000162722 ENSG00000162722 12
## ENSG00000162757 ENSG00000162757 15
## ENSG00000162877 ENSG00000162877 18
## ENSG00000162928 ENSG00000162928 3
## ENSG00000163116 ENSG00000163116 6
## ENSG00000163191 ENSG00000163191 9
## ENSG00000163328 ENSG00000163328 20
## ENSG00000163399 ENSG00000163399 13
## ENSG00000163406 ENSG00000163406 8
## ENSG00000163449 ENSG00000163449 2
## ENSG00000163466 ENSG00000163466 5
## ENSG00000163513 ENSG00000163513 19
## ENSG00000163521 ENSG00000163521 14
## ENSG00000163534 ENSG00000163534 2
## ENSG00000163624 ENSG00000163624 2
## ENSG00000163666 ENSG00000163666 2
## ENSG00000163702 ENSG00000163702 13
## ENSG00000163879 ENSG00000163879 15
## ENSG00000163995 ENSG00000163995 2
## ENSG00000164124 ENSG00000164124 13
## ENSG00000164125 ENSG00000164125 14
## ENSG00000164136 ENSG00000164136 14
## ENSG00000164187 ENSG00000164187 18
## ENSG00000164309 ENSG00000164309 15
## ENSG00000164379 ENSG00000164379 2
## ENSG00000164440 ENSG00000164440 2
## ENSG00000164442 ENSG00000164442 22
## ENSG00000164466 ENSG00000164466 8
## ENSG00000164741 ENSG00000164741 2
## ENSG00000164818 ENSG00000164818 1
## ENSG00000164985 ENSG00000164985 10
## ENSG00000165028 ENSG00000165028 8
## ENSG00000165259 ENSG00000165259 13
## ENSG00000165280 ENSG00000165280 9
## ENSG00000165527 ENSG00000165527 16
## ENSG00000165685 ENSG00000165685 13
## ENSG00000165733 ENSG00000165733 1
## ENSG00000165804 ENSG00000165804 1
## ENSG00000165819 ENSG00000165819 24
## ENSG00000165891 ENSG00000165891 17
## ENSG00000165943 ENSG00000165943 7
## ENSG00000165949 ENSG00000165949 2
## ENSG00000165966 ENSG00000165966 4
## ENSG00000166123 ENSG00000166123 1
## ENSG00000166164 ENSG00000166164 14
## ENSG00000166199 ENSG00000166199 1
## ENSG00000166257 ENSG00000166257 18
## ENSG00000166432 ENSG00000166432 24
## ENSG00000166448 ENSG00000166448 2
## ENSG00000166484 ENSG00000166484 16
## ENSG00000166503 ENSG00000166503 3
## ENSG00000166508 ENSG00000166508 19
## ENSG00000166510 ENSG00000166510 15
## ENSG00000166526 ENSG00000166526 17
## ENSG00000166557 ENSG00000166557 8
## ENSG00000166592 ENSG00000166592 14
## ENSG00000166801 ENSG00000166801 20
## ENSG00000166928 ENSG00000166928 14
## ENSG00000166949 ENSG00000166949 12
## ENSG00000167193 ENSG00000167193 9
## ENSG00000167291 ENSG00000167291 2
## ENSG00000167315 ENSG00000167315 1
## ENSG00000167543 ENSG00000167543 8
## ENSG00000167562 ENSG00000167562 10
## ENSG00000167566 ENSG00000167566 9
## ENSG00000167703 ENSG00000167703 5
## ENSG00000167772 ENSG00000167772 1
## ENSG00000167925 ENSG00000167925 16
## ENSG00000167994 ENSG00000167994 2
## ENSG00000167995 ENSG00000167995 5
## ENSG00000168016 ENSG00000168016 10
## ENSG00000168209 ENSG00000168209 19
## ENSG00000168256 ENSG00000168256 14
## ENSG00000168273 ENSG00000168273 13
## ENSG00000168421 ENSG00000168421 17
## ENSG00000168439 ENSG00000168439 2
## ENSG00000168569 ENSG00000168569 1
## ENSG00000168679 ENSG00000168679 17
## ENSG00000168795 ENSG00000168795 3
## ENSG00000168944 ENSG00000168944 17
## ENSG00000168994 ENSG00000168994 2
## ENSG00000169047 ENSG00000169047 1
## ENSG00000169239 ENSG00000169239 23
## ENSG00000169432 ENSG00000169432 20
## ENSG00000169860 ENSG00000169860 1
## ENSG00000169908 ENSG00000169908 15
## ENSG00000169946 ENSG00000169946 15
## ENSG00000170027 ENSG00000170027 2
## ENSG00000170037 ENSG00000170037 23
## ENSG00000170234 ENSG00000170234 10
## ENSG00000170312 ENSG00000170312 18
## ENSG00000170370 ENSG00000170370 18
## ENSG00000170473 ENSG00000170473 1
## ENSG00000170677 ENSG00000170677 13
## ENSG00000170989 ENSG00000170989 8
## ENSG00000171100 ENSG00000171100 24
## ENSG00000171262 ENSG00000171262 3
## ENSG00000171365 ENSG00000171365 14
## ENSG00000171530 ENSG00000171530 1
## ENSG00000171604 ENSG00000171604 14
## ENSG00000171617 ENSG00000171617 12
## ENSG00000171729 ENSG00000171729 14
## ENSG00000171812 ENSG00000171812 14
## ENSG00000171877 ENSG00000171877 4
## ENSG00000172123 ENSG00000172123 22
## ENSG00000172159 ENSG00000172159 11
## ENSG00000172331 ENSG00000172331 6
## ENSG00000172403 ENSG00000172403 15
## ENSG00000172594 ENSG00000172594 1
## ENSG00000172716 ENSG00000172716 2
## ENSG00000173124 ENSG00000173124 18
## ENSG00000173156 ENSG00000173156 2
## ENSG00000173166 ENSG00000173166 14
## ENSG00000173391 ENSG00000173391 13
## ENSG00000173457 ENSG00000173457 1
## ENSG00000173511 ENSG00000173511 3
## ENSG00000173530 ENSG00000173530 14
## ENSG00000173818 ENSG00000173818 23
## ENSG00000173917 ENSG00000173917 1
## ENSG00000174177 ENSG00000174177 12
## ENSG00000174197 ENSG00000174197 8
## ENSG00000174371 ENSG00000174371 3
## ENSG00000174705 ENSG00000174705 2
## ENSG00000174989 ENSG00000174989 8
## ENSG00000175130 ENSG00000175130 7
## ENSG00000175445 ENSG00000175445 3
## ENSG00000175544 ENSG00000175544 3
## ENSG00000175691 ENSG00000175691 3
## ENSG00000176014 ENSG00000176014 13
## ENSG00000176046 ENSG00000176046 15
## ENSG00000176105 ENSG00000176105 3
## ENSG00000176125 ENSG00000176125 18
## ENSG00000176531 ENSG00000176531 1
## ENSG00000176903 ENSG00000176903 3
## ENSG00000177119 ENSG00000177119 13
## ENSG00000177294 ENSG00000177294 20
## ENSG00000177311 ENSG00000177311 18
## ENSG00000177469 ENSG00000177469 3
## ENSG00000177479 ENSG00000177479 23
## ENSG00000177989 ENSG00000177989 11
## ENSG00000178209 ENSG00000178209 23
## ENSG00000178338 ENSG00000178338 1
## ENSG00000178573 ENSG00000178573 1
## ENSG00000178700 ENSG00000178700 2
## ENSG00000178741 ENSG00000178741 13
## ENSG00000178896 ENSG00000178896 9
## ENSG00000179044 ENSG00000179044 20
## ENSG00000179262 ENSG00000179262 23
## ENSG00000179409 ENSG00000179409 1
## ENSG00000179476 ENSG00000179476 15
## ENSG00000179528 ENSG00000179528 1
## ENSG00000180376 ENSG00000180376 8
## ENSG00000180543 ENSG00000180543 1
## ENSG00000181045 ENSG00000181045 14
## ENSG00000181350 ENSG00000181350 1
## ENSG00000181392 ENSG00000181392 22
## ENSG00000181523 ENSG00000181523 14
## ENSG00000181873 ENSG00000181873 1
## ENSG00000182158 ENSG00000182158 8
## ENSG00000182173 ENSG00000182173 8
## ENSG00000182240 ENSG00000182240 17
## ENSG00000182263 ENSG00000182263 4
## ENSG00000182378 ENSG00000182378 8
## ENSG00000182504 ENSG00000182504 15
## ENSG00000182557 ENSG00000182557 12
## ENSG00000182704 ENSG00000182704 2
## ENSG00000182885 ENSG00000182885 19
## ENSG00000182901 ENSG00000182901 6
## ENSG00000182952 ENSG00000182952 20
## ENSG00000182986 ENSG00000182986 3
## ENSG00000183044 ENSG00000183044 10
## ENSG00000183087 ENSG00000183087 1
## ENSG00000183091 ENSG00000183091 6
## ENSG00000183117 ENSG00000183117 14
## ENSG00000183283 ENSG00000183283 5
## ENSG00000183431 ENSG00000183431 7
## ENSG00000183578 ENSG00000183578 2
## ENSG00000183617 ENSG00000183617 13
## ENSG00000183763 ENSG00000183763 3
## ENSG00000183785 ENSG00000183785 22
## ENSG00000183801 ENSG00000183801 18
## ENSG00000183853 ENSG00000183853 2
## ENSG00000183955 ENSG00000183955 7
## ENSG00000184226 ENSG00000184226 15
## ENSG00000184384 ENSG00000184384 11
## ENSG00000184402 ENSG00000184402 10
## ENSG00000184661 ENSG00000184661 18
## ENSG00000185033 ENSG00000185033 11
## ENSG00000185043 ENSG00000185043 7
## ENSG00000185561 ENSG00000185561 4
## ENSG00000185669 ENSG00000185669 5
## ENSG00000185686 ENSG00000185686 4
## ENSG00000185736 ENSG00000185736 19
## ENSG00000185825 ENSG00000185825 9
## ENSG00000185842 ENSG00000185842 6
## ENSG00000185920 ENSG00000185920 18
## ENSG00000186197 ENSG00000186197 2
## ENSG00000186645 ENSG00000186645 10
## ENSG00000186818 ENSG00000186818 14
## ENSG00000186951 ENSG00000186951 12
## ENSG00000187097 ENSG00000187097 8
## ENSG00000187105 ENSG00000187105 17
## ENSG00000187164 ENSG00000187164 13
## ENSG00000187231 ENSG00000187231 11
## ENSG00000187554 ENSG00000187554 11
## ENSG00000187566 ENSG00000187566 1
## ENSG00000188037 ENSG00000188037 9
## ENSG00000188343 ENSG00000188343 15
## ENSG00000188452 ENSG00000188452 22
## ENSG00000188636 ENSG00000188636 13
## ENSG00000188921 ENSG00000188921 10
## ENSG00000188938 ENSG00000188938 9
## ENSG00000189013 ENSG00000189013 2
## ENSG00000189067 ENSG00000189067 6
## ENSG00000189077 ENSG00000189077 5
## ENSG00000189159 ENSG00000189159 19
## ENSG00000189195 ENSG00000189195 20
## ENSG00000189337 ENSG00000189337 6
## ENSG00000196072 ENSG00000196072 5
## ENSG00000196123 ENSG00000196123 14
## ENSG00000196126 ENSG00000196126 13
## ENSG00000196199 ENSG00000196199 10
## ENSG00000196230 ENSG00000196230 1
## ENSG00000196247 ENSG00000196247 6
## ENSG00000196305 ENSG00000196305 1
## ENSG00000196369 ENSG00000196369 11
## ENSG00000196378 ENSG00000196378 12
## ENSG00000196405 ENSG00000196405 8
## ENSG00000196526 ENSG00000196526 3
## ENSG00000196652 ENSG00000196652 6
## ENSG00000196872 ENSG00000196872 6
## ENSG00000196923 ENSG00000196923 5
## ENSG00000197008 ENSG00000197008 3
## ENSG00000197093 ENSG00000197093 6
## ENSG00000197283 ENSG00000197283 8
## ENSG00000197302 ENSG00000197302 18
## ENSG00000197451 ENSG00000197451 13
## ENSG00000197461 ENSG00000197461 3
## ENSG00000197603 ENSG00000197603 3
## ENSG00000197782 ENSG00000197782 18
## ENSG00000197959 ENSG00000197959 18
## ENSG00000198178 ENSG00000198178 20
## ENSG00000198286 ENSG00000198286 1
## ENSG00000198336 ENSG00000198336 18
## ENSG00000198435 ENSG00000198435 3
## ENSG00000198502 ENSG00000198502 13
## ENSG00000198690 ENSG00000198690 10
## ENSG00000198700 ENSG00000198700 13
## ENSG00000198795 ENSG00000198795 15
## ENSG00000198805 ENSG00000198805 13
## ENSG00000198829 ENSG00000198829 18
## ENSG00000198959 ENSG00000198959 2
## ENSG00000203546 ENSG00000203546 1
## ENSG00000204103 ENSG00000204103 14
## ENSG00000204388 ENSG00000204388 2
## ENSG00000204389 ENSG00000204389 6
## ENSG00000204520 ENSG00000204520 6
## ENSG00000204568 ENSG00000204568 1
## ENSG00000204869 ENSG00000204869 15
## ENSG00000204909 ENSG00000204909 6
## ENSG00000204991 ENSG00000204991 18
## ENSG00000205038 ENSG00000205038 3
## ENSG00000205060 ENSG00000205060 1
## ENSG00000205090 ENSG00000205090 10
## ENSG00000205362 ENSG00000205362 15
## ENSG00000205639 ENSG00000205639 12
## ENSG00000205730 ENSG00000205730 14
## ENSG00000211448 ENSG00000211448 15
## ENSG00000213588 ENSG00000213588 3
## ENSG00000213719 ENSG00000213719 11
## ENSG00000213859 ENSG00000213859 3
## ENSG00000214226 ENSG00000214226 2
## ENSG00000214706 ENSG00000214706 1
## ENSG00000214872 ENSG00000214872 5
## ENSG00000215788 ENSG00000215788 1
## ENSG00000219481 ENSG00000219481 8
## ENSG00000235750 ENSG00000235750 20
## ENSG00000239920 ENSG00000239920 10
## ENSG00000240445 ENSG00000240445 4
## ENSG00000244165 ENSG00000244165 13
## ENSG00000244405 ENSG00000244405 2
## ENSG00000248405 ENSG00000248405 20
## ENSG00000248993 ENSG00000248993 13
## ENSG00000249242 ENSG00000249242 1
## ENSG00000254979 ENSG00000254979 23
## ENSG00000256043 ENSG00000256043 1
## ENSG00000256235 ENSG00000256235 12
## ENSG00000257335 ENSG00000257335 5
## ENSG00000259330 ENSG00000259330 23
## ENSG00000261652 ENSG00000261652 18
## ENSG00000261832 ENSG00000261832 24
## ENSG00000263001 ENSG00000263001 25
## ENSG00000263528 ENSG00000263528 13
## ENSG00000263715 ENSG00000263715 3
## ENSG00000264364 ENSG00000264364 19
## ENSG00000268182 ENSG00000268182 8
## ENSG00000268861 ENSG00000268861 24
## ENSG00000273802 ENSG00000273802 6
## ENSG00000274810 ENSG00000274810 18
## ENSG00000275895 ENSG00000275895 7
## ENSG00000276085 ENSG00000276085 6
## ENSG00000277117 ENSG00000277117 13
## ENSG00000282988 ENSG00000282988 10
## ENSG00000283149 ENSG00000283149 13
## ENSG00000283977 ENSG00000283977 16
## ENSG00000285253 ENSG00000285253 1
## ENSG00000285708 ENSG00000285708 104.4.2 Look at only the differential genes
A good suggestion from Theresa was to examine only the most variant genes from failure vs. cure and see how they change the clustering/etc results. This is my attempt to address this query.
hs_clinic_topn <- sm(extract_significant_genes(hs_clinic_table, n = 100))
table <- "failure_vs_cure"
wanted <- rbind(hs_clinic_topn[["deseq"]][["ups"]][[table]],
hs_clinic_topn[["deseq"]][["downs"]][[table]])
small_expt <- exclude_genes_expt(hs_clinical_nobiop, ids = rownames(wanted), method = "keep")## Before removal, there were 19941 genes, now there are 200.## There are 75 samples which kept less than 90 percent counts.## X1017n1 X1017m1 X1034n1 X1034n2 X1034m2 X1034m2. X2052e1 X2052m2
## 0.4192 0.3754 2.4868 2.9008 1.2743 1.2381 0.5331 0.6164
## X2052n2 X2052m3 X2052n3 X2065m1 X2065n1 X2066m1 X2066n1 X2065m2
## 1.4080 0.7016 1.1583 1.0719 2.8747 0.3940 0.7222 0.5819
## X2065n2 X2065e2 X2066m2 X2066n2 X2066e2 X2068m1 X2068n1 X2068e1
## 0.4933 0.8420 0.4952 0.6676 0.8413 0.5006 0.5287 1.0548
## X2072m1 X2072n1 X2072e1 X2071m1 X2071n1 X2073m1 X2073n1 X2073e1
## 0.4799 0.4190 0.6492 0.7265 1.6127 0.6570 1.7804 0.7564
## X2068m2 X2068n2 X2068e2 X2072m2 X2072n2 X2072e2 X2073m2 X2073n2
## 0.3427 0.4429 0.6802 0.3992 0.5430 0.5317 1.0397 2.8228
## X2073e2 X2066m3 X2066n3 X2065e3 X2068m3 X2068n3 X2068e3 X2072m3
## 0.7690 0.4304 0.5995 0.4657 0.4092 0.4483 0.6113 0.7292
## X2072n3 X2072e3 X2073m3 X2073n3 X2073e3 X2162m1 X2162n1 X2162e1
## 1.9625 0.7389 0.5093 0.6651 0.4507 0.4594 0.4140 0.7208
## X2167m1 X2168n1 X2168e1 X2168m2 X2168n2 X2168e2 X2167m2 X2167n3
## 0.6211 2.3153 0.9847 0.9395 2.3681 0.8613 0.3830 1.0254
## X2167e3 X2168m3 X2168n3 X2168e3 X2172n1 X2172e1 X1168m1 X1168n1
## 1.1411 1.3360 3.7049 1.0558 0.3235 0.3720 0.7320 1.4522
## X1168m2 X1168e2 X1168n3
## 0.7185 0.5010 1.0138small_norm <- sm(normalize_expt(small_expt, transform = "log2", convert = "cpm",
norm = "quant", filter = TRUE))
plot_pca(small_norm)$plot## Warning: ggrepel: 15 unlabeled data points (too many overlaps). Consider
## increasing max.overlaps
small_nb <- normalize_expt(small_expt, transform = "log2", convert = "cpm",
batch = "svaseq", norm = "quant", filter = TRUE)## Warning in normalize_expt(small_expt, transform = "log2", convert = "cpm", :
## Quantile normalization and sva do not always play well together.## Removing 0 low-count genes (200 remaining).## batch_counts: Before batch/surrogate estimation, 2 entries are x==0: 0%.## batch_counts: Before batch/surrogate estimation, 526 entries are 0<x<1: 4%.## Setting 58 low elements to zero.## transform_counts: Found 58 values equal to 0, adding 1 to the matrix.plot_pca(small_nb)$plot
## DESeq2 MA plot of failure / cure
hs_clinic_table[["plots"]][["failure_vs_cure"]][["deseq_ma_plots"]]$plot
## DESeq2 Volcano plot of failure / cure
hs_clinic_table[["plots"]][["failure_vs_cure"]][["deseq_vol_plots"]]$plot
4.4.3 g:Profiler results using the significant up and down genes
ups <- hs_clinic_sig[["deseq"]][["ups"]][[1]]
downs <- hs_clinic_sig[["deseq"]][["downs"]][[1]]
hs_clinic_gprofiler_ups <- simple_gprofiler(ups)## Performing gProfiler GO search of 227 genes against hsapiens.## GO search found 71 hits.## Performing gProfiler KEGG search of 227 genes against hsapiens.## KEGG search found 5 hits.## Performing gProfiler REAC search of 227 genes against hsapiens.## REAC search found 9 hits.## Performing gProfiler MI search of 227 genes against hsapiens.## MI search found 0 hits.## Performing gProfiler TF search of 227 genes against hsapiens.## TF search found 43 hits.## Performing gProfiler CORUM search of 227 genes against hsapiens.## CORUM search found 0 hits.## Performing gProfiler HP search of 227 genes against hsapiens.## HP search found 0 hits.hs_clinic_gprofiler_ups[["pvalue_plots"]][["bpp_plot_over"]]
hs_clinic_gprofiler_ups[["pvalue_plots"]][["mfp_plot_over"]]
hs_clinic_gprofiler_ups[["pvalue_plots"]][["reactome_plot_over"]]
##hs_try2 <- simple_gprofiler2(ups)
hs_clinic_gprofiler_downs <- simple_gprofiler(downs)## Performing gProfiler GO search of 323 genes against hsapiens.## GO search found 106 hits.## Performing gProfiler KEGG search of 323 genes against hsapiens.## KEGG search found 3 hits.## Performing gProfiler REAC search of 323 genes against hsapiens.## REAC search found 5 hits.## Performing gProfiler MI search of 323 genes against hsapiens.## MI search found 0 hits.## Performing gProfiler TF search of 323 genes against hsapiens.## TF search found 40 hits.## Performing gProfiler CORUM search of 323 genes against hsapiens.## CORUM search found 0 hits.## Performing gProfiler HP search of 323 genes against hsapiens.## HP search found 3 hits.hs_clinic_gprofiler_downs[["pvalue_plots"]][["bpp_plot_over"]]
hs_clinic_gprofiler_downs[["pvalue_plots"]][["mfp_plot_over"]]
hs_clinic_gprofiler_downs[["pvalue_plots"]][["reactome_plot_over"]]
4.5 Perform GSVA on the clinical samples
hs_celltype_gsva_c2 <- sm(simple_gsva(individual_celltypes))
hs_celltype_gsva_c2_sig <- sm(get_sig_gsva_categories(
hs_celltype_gsva_c2,
excel = "excel/individual_celltypes_gsva_c2.xlsx"))
broad_c7 <- GSEABase::getGmt("reference/msigdb/c7.all.v7.2.entrez.gmt",
collectionType = GSEABase::BroadCollection(category = "c7"),
geneIdType = GSEABase::EntrezIdentifier())
hs_celltype_gsva_c7 <- sm(simple_gsva(individual_celltypes, signatures = broad_c7,
msig_xml = "reference/msigdb_v7.2.xml", cores = 10))
hs_celltype_gsva_c7_sig <- sm(get_sig_gsva_categories(
hs_celltype_gsva_c7,
excel = "excel/individual_celltypes_gsva_c7.xlsx"))4.5.1 Print some plots of the GSVA outputs
## The raw heatmap of the C2 values
hs_celltype_gsva_c2_sig[["raw_plot"]]
## The 'significance' scores of the C2 values
hs_celltype_gsva_c2_sig[["score_plot"]]
## The subset of scores for categories deemed significantly different.
hs_celltype_gsva_c2_sig[["subset_plot"]]
## The raw heatmap of the C7 values
hs_celltype_gsva_c7_sig[["raw_plot"]]
## The 'significance' scores of the C7 values
hs_celltype_gsva_c7_sig[["score_plot"]]
## The subset of scores for categories deemed significantly different.
hs_celltype_gsva_c7_sig[["subset_plot"]]
5 Individual Cell types
The following blocks split the samples into a few groups by sample type and look at the distributions between them.
5.1 Implementation details
Get top/bottom n genes for each cell type, using clinical outcome as the factor of interest. For the moment, use sva for the DE analysis. Provide cpms for the top/bottom n genes.
Start with top/bottom 200. Perform default logFC and p-value as well.
5.2 Monocytes
5.2.1 Evaluate Monocyte samples
mono <- subset_expt(hs_valid, subset = "typeofcells=='monocytes'") %>%
set_expt_conditions(fact = "clinicaloutcome") %>%
set_expt_batches(fact = "donor") %>%
set_expt_colors(colors = chosen_colors)## subset_expt(): There were 135, now there are 27 samples.## FIXME set_expt_colors should speak up if there are mismatches here!!!
save_result <- save(mono, file = "rda/monocyte_expt.rda")
mono_norm <- normalize_expt(mono, convert = "cpm", filter = TRUE,
transform = "log2", norm = "quant")## Removing 8915 low-count genes (11026 remaining).## transform_counts: Found 9 values equal to 0, adding 1 to the matrix.plt <- plot_pca(mono_norm, plot_labels = FALSE)$plot
pp(file = glue("images/mono_pca_normalized-v{ver}.pdf"), image = plt)
mono_nb <- normalize_expt(mono, convert = "cpm", filter = TRUE,
transform = "log2", batch = "svaseq")## Removing 8915 low-count genes (11026 remaining).## batch_counts: Before batch/surrogate estimation, 1404 entries are x==0: 0%.## batch_counts: Before batch/surrogate estimation, 17611 entries are 0<x<1: 6%.## Setting 556 low elements to zero.## transform_counts: Found 556 values equal to 0, adding 1 to the matrix.plt <- plot_pca(mono_nb, plot_labels = FALSE)$plot
pp(file = glue("images/mono_pca_normalized_batch-v{ver}.pdf"), image = plt)
5.2.2 DE of Monocyte samples
5.2.2.1 Without sva
mono_de <- sm(all_pairwise(mono, model_batch = FALSE, filter = TRUE))mono_tables <- sm(combine_de_tables(
mono_de, keepers = keepers,
excel = glue::glue("excel/monocyte_clinical_all_tables-v{ver}.xlsx")))
written <- write_xlsx(data = mono_tables[["data"]][[1]],
excel = glue::glue("excel/monocyte_clinical_table-v{ver}.xlsx"))
mono_sig <- sm(extract_significant_genes(mono_tables, according_to = "deseq"))
written <- write_xlsx(data = mono_sig[["deseq"]][["ups"]][[1]],
excel = glue::glue("excel/monocyte_clinical_sigup-v{ver}.xlsx"))
written <- write_xlsx(data = mono_sig[["deseq"]][["downs"]][[1]],
excel = glue::glue("excel/monocyte_clinical_sigdown-v{ver}.xlsx"))
mono_pct_sig <- sm(extract_significant_genes(mono_tables, n = 200,
lfc = NULL, p = NULL, according_to = "deseq"))
written <- write_xlsx(data = mono_pct_sig[["deseq"]][["ups"]][[1]],
excel = glue::glue("excel/monocyte_clinical_sigup_pct-v{ver}.xlsx"))
written <- write_xlsx(data = mono_pct_sig[["deseq"]][["downs"]][[1]],
excel = glue::glue("excel/monocyte_clinical_sigdown_pct-v{ver}.xlsx"))
mono_sig$summary_df## data frame with 0 columns and 1 row## Print out a table of the cpm values for other explorations.
mono_cpm <- sm(normalize_expt(mono, convert = "cpm"))
written <- write_xlsx(data = exprs(mono_cpm),
excel = glue::glue("excel/monocyte_cpm_before_batch-v{ver}.xlsx"))
mono_bcpm <- sm(normalize_expt(mono, filter = TRUE, convert = "cpm", batch = "svaseq"))
written <- write_xlsx(data = exprs(mono_bcpm),
excel = glue::glue("excel/monocyte_cpm_after_batch-v{ver}.xlsx"))5.2.2.2 With sva
mono_de_sva <- sm(all_pairwise(mono, model_batch = "svaseq", filter = TRUE))mono_tables_sva <- sm(combine_de_tables(
mono_de_sva, keepers = keepers,
excel = glue::glue("excel/monocyte_clinical_all_tables_sva-v{ver}.xlsx")))
mono_sig_sva <- sm(extract_significant_genes(
mono_tables_sva,
excel = glue::glue("excel/monocyte_clinical_sig_tables_sva-v{ver}.xlsx"),
according_to = "deseq"))5.2.2.3 Monocyte DE plots
First print out the DE plots without and then with sva estimates.
## DESeq2 MA plot of failure / cure
mono_tables[["plots"]][["fail_vs_cure"]][["deseq_ma_plots"]]$plot
## DESeq2 Volcano plot of failure / cure
mono_tables[["plots"]][["fail_vs_cure"]][["deseq_vol_plots"]]$plot
## DESeq2 MA plot of failure / cure with svaseq
mono_tables_sva[["plots"]][["fail_vs_cure"]][["deseq_ma_plots"]]$plot
## DESeq2 Volcano plot of failure / cure with svaseq
mono_tables_sva[["plots"]][["fail_vs_cure"]][["deseq_vol_plots"]]$plot
5.2.2.4 Monocyte ontology search
ups <- mono_sig[["deseq"]][["ups"]][["fail_vs_cure"]]
downs <- mono_sig[["deseq"]][["downs"]][["fail_vs_cure"]]
mono_up_gp <- simple_gprofiler(sig_genes = ups)## Performing gProfiler GO search of 103 genes against hsapiens.## GO search found 22 hits.## Performing gProfiler KEGG search of 103 genes against hsapiens.## KEGG search found 4 hits.## Performing gProfiler REAC search of 103 genes against hsapiens.## REAC search found 4 hits.## Performing gProfiler MI search of 103 genes against hsapiens.## MI search found 0 hits.## Performing gProfiler TF search of 103 genes against hsapiens.## TF search found 13 hits.## Performing gProfiler CORUM search of 103 genes against hsapiens.## CORUM search found 1 hits.## Performing gProfiler HP search of 103 genes against hsapiens.## HP search found 19 hits.mono_up_gp[["pvalue_plots"]][["bpp_plot_over"]]
mono_up_gp[["pvalue_plots"]][["mfp_plot_over"]]
mono_up_gp[["pvalue_plots"]][["reactome_plot_over"]]
mono_down_gp <- simple_gprofiler(sig_genes = downs)## Performing gProfiler GO search of 265 genes against hsapiens.## GO search found 97 hits.## Performing gProfiler KEGG search of 265 genes against hsapiens.## KEGG search found 2 hits.## Performing gProfiler REAC search of 265 genes against hsapiens.## REAC search found 6 hits.## Performing gProfiler MI search of 265 genes against hsapiens.## MI search found 0 hits.## Performing gProfiler TF search of 265 genes against hsapiens.## TF search found 136 hits.## Performing gProfiler CORUM search of 265 genes against hsapiens.## CORUM search found 0 hits.## Performing gProfiler HP search of 265 genes against hsapiens.## HP search found 0 hits.mono_down_gp[["pvalue_plots"]][["bpp_plot_over"]]
mono_down_gp[["pvalue_plots"]][["mfp_plot_over"]]
mono_down_gp[["pvalue_plots"]][["reactome_plot_over"]]
5.2.2.5 Monocyte MSigDB query
broad_c7 <- GSEABase::getGmt("reference/msigdb/c7.all.v7.2.entrez.gmt",
collectionType = GSEABase::BroadCollection(category = "c7"),
geneIdType = GSEABase::EntrezIdentifier())
mono_up_goseq_msig <- goseq_msigdb(sig_genes = ups, signatures = broad_c7,
signature_category = "c7", length_db = hs_length)## Before conversion: 103, after conversion: 104.## Before conversion: 227921, after conversion: 35341.## Found 99 go_db genes and 104 length_db genes out of 104.mono_down_goseq_msig <- goseq_msigdb(sig_genes = downs, signatures = broad_c7,
signature_category = "c7", length_db = hs_length)## Before conversion: 265, after conversion: 264.## Before conversion: 227921, after conversion: 35341.## Found 254 go_db genes and 264 length_db genes out of 264.5.2.2.6 Plot of similar experiments
## Monocyte genes with increased expression in the failed samples
## share genes with the following experiments
mono_up_goseq_msig[["pvalue_plots"]][["mfp_plot_over"]]
## Monocyte genes with increased expression in the cured samples
## share genes with the following experiments
mono_down_goseq_msig[["pvalue_plots"]][["mfp_plot_over"]]
5.2.3 Evaluate Neutrophil samples
neut <- subset_expt(hs_valid, subset = "typeofcells=='neutrophils'") %>%
set_expt_conditions(fact = "clinicaloutcome") %>%
set_expt_batches(fact = "donor") %>%
set_expt_colors(colors = chosen_colors)## subset_expt(): There were 135, now there are 28 samples.save_result <- save(neut, file = "rda/neutrophil_expt.rda")
neut_norm <- sm(normalize_expt(neut, convert = "cpm", filter = TRUE, transform = "log2"))
plt <- plot_pca(neut_norm, plot_labels = FALSE)$plot
pp(file = glue("images/neut_pca_normalized-v{ver}.pdf"), image = plt)
neut_nb <- sm(normalize_expt(neut, convert = "cpm", filter = TRUE,
transform = "log2", batch = "svaseq"))
plt <- plot_pca(neut_nb, plot_labels = FALSE)$plot
pp(file = glue("images/neut_pca_normalized_svaseq-v{ver}.pdf"), image = plt)
5.2.4 DE of Netrophil samples
5.2.4.1 Without sva
neut_de <- sm(all_pairwise(neut, model_batch = FALSE, filter = TRUE))neut_tables <- sm(combine_de_tables(
neut_de, keepers = keepers,
excel = glue::glue("excel/neutrophil_clinical_all_tables-v{ver}.xlsx")))
written <- write_xlsx(data = neut_tables[["data"]][[1]],
excel = glue::glue("excel/neutrophil_clinical_table-v{ver}.xlsx"))
neut_sig <- sm(extract_significant_genes(neut_tables, according_to = "deseq"))
written <- write_xlsx(data = neut_sig[["deseq"]][["ups"]][[1]],
excel = glue::glue("excel/neutrophil_clinical_sigup-v{ver}.xlsx"))
written <- write_xlsx(data = neut_sig[["deseq"]][["downs"]][[1]],
excel = glue::glue("excel/neutrophil_clinical_sigdown-v{ver}.xlsx"))
neut_pct_sig <- sm(extract_significant_genes(neut_tables, n = 200, lfc = NULL,
p = NULL, according_to = "deseq"))
written <- write_xlsx(data = neut_sig[["deseq"]][["ups"]][[1]],
excel = glue::glue("excel/neutrophil_clinical_sigup_pct-v{ver}.xlsx"))
written <- write_xlsx(data = neut_sig[["deseq"]][["downs"]][[1]],
excel = glue::glue("excel/neutrophil_clinical_sigdown_pct-v{ver}.xlsx"))
neut_cpm <- sm(normalize_expt(neut, convert = "cpm"))
written <- write_xlsx(data = exprs(neut_cpm),
excel = glue::glue("excel/neutrophil_cpm_before_batch-v{ver}.xlsx"))
neut_bcpm <- sm(normalize_expt(neut, filter = TRUE, batch = "svaseq", convert = "cpm"))
written <- write_xlsx(data = exprs(neut_bcpm),
excel = glue::glue("excel/neutrophil_cpm_after_batch-v{ver}.xlsx"))5.2.4.2 With sva
neut_de_sva <- sm(all_pairwise(neut, model_batch = "svaseq", filter = TRUE))neut_tables_sva <- sm(combine_de_tables(
neut_de_sva, keepers = keepers,
excel = glue::glue("excel/neutrophil_clinical_all_tables_sva-v{ver}.xlsx")))
neut_sig_sva <- sm(extract_significant_genes(
neut_tables_sva,
excel = glue::glue("excel/neutrophil_clinical_sig_tables_sva-v{ver}.xlsx"),
according_to = "deseq"))5.2.4.3 Neutrophil DE plots
## DESeq2 MA plot of failure / cure
neut_tables[["plots"]][["fail_vs_cure"]][["deseq_ma_plots"]]$plot
## DESeq2 Volcano plot of failure / cure
neut_tables[["plots"]][["fail_vs_cure"]][["deseq_vol_plots"]]$plot
## DESeq2 MA plot of failure / cure with sva
neut_tables_sva[["plots"]][["fail_vs_cure"]][["deseq_ma_plots"]]$plot
## DESeq2 Volcano plot of failure / cure with sva
neut_tables_sva[["plots"]][["fail_vs_cure"]][["deseq_vol_plots"]]$plot
5.2.4.4 Neutrophil ontology search
ups <- neut_sig[["deseq"]][["ups"]][["fail_vs_cure"]]
downs <- neut_sig[["deseq"]][["downs"]][["fail_vs_cure"]]
neut_up_gp <- simple_gprofiler(sig_genes = ups)## Performing gProfiler GO search of 122 genes against hsapiens.## GO search found 0 hits.## Performing gProfiler KEGG search of 122 genes against hsapiens.## KEGG search found 1 hits.## Performing gProfiler REAC search of 122 genes against hsapiens.## REAC search found 1 hits.## Performing gProfiler MI search of 122 genes against hsapiens.## MI search found 0 hits.## Performing gProfiler TF search of 122 genes against hsapiens.## TF search found 6 hits.## Performing gProfiler CORUM search of 122 genes against hsapiens.## CORUM search found 0 hits.## Performing gProfiler HP search of 122 genes against hsapiens.## HP search found 1 hits.neut_up_gp[["pvalue_plots"]][["bpp_plot_over"]]## NULLneut_up_gp[["pvalue_plots"]][["mfp_plot_over"]]## NULLneut_up_gp[["pvalue_plots"]][["reactome_plot_over"]]
neut_down_gp <- simple_gprofiler(downs)## Performing gProfiler GO search of 137 genes against hsapiens.## GO search found 2 hits.## Performing gProfiler KEGG search of 137 genes against hsapiens.## KEGG search found 0 hits.## Performing gProfiler REAC search of 137 genes against hsapiens.## REAC search found 0 hits.## Performing gProfiler MI search of 137 genes against hsapiens.## MI search found 0 hits.## Performing gProfiler TF search of 137 genes against hsapiens.## TF search found 6 hits.## Performing gProfiler CORUM search of 137 genes against hsapiens.## CORUM search found 0 hits.## Performing gProfiler HP search of 137 genes against hsapiens.## HP search found 5 hits.neut_down_gp[["pvalue_plots"]][["bpp_plot_over"]]
neut_down_gp[["pvalue_plots"]][["mfp_plot_over"]]## NULLneut_down_gp[["pvalue_plots"]][["reactome_plot_over"]]## NULL5.2.4.5 Neutrophil GSVA query
neut_up_goseq_msig <- goseq_msigdb(sig_genes = ups, signatures = broad_c7,
signature_category = "c7", length_db = hs_length)## Before conversion: 122, after conversion: 121.## Before conversion: 227921, after conversion: 35341.## Found 114 go_db genes and 121 length_db genes out of 121.neut_down_goseq_msig <- goseq_msigdb(sig_genes = downs, signatures = broad_c7,
signature_category = "c7", length_db = hs_length)## Before conversion: 137, after conversion: 133.## Before conversion: 227921, after conversion: 35341.## Found 127 go_db genes and 133 length_db genes out of 133.5.2.4.6 Plot of similar experiments
## Neutrophil genes with increased expression in the failed samples
## share genes with the following experiments
neut_up_goseq_msig[["pvalue_plots"]][["mfp_plot_over"]]
## Neutrophil genes with increased expression in the cured samples
## share genes with the following experiments
neut_down_goseq_msig[["pvalue_plots"]][["mfp_plot_over"]]
5.3 Eosinophils
5.3.1 Evaluate Eosinophil samples
eo <- subset_expt(hs_valid, subset = "typeofcells=='eosinophils'") %>%
set_expt_conditions(fact = "clinicaloutcome") %>%
set_expt_batches(fact = "donor") %>%
set_expt_colors(colors = chosen_colors)## subset_expt(): There were 135, now there are 20 samples.save_result <- save(eo, file = "rda/eosinophil_expt.rda")
eo_norm <- sm(normalize_expt(eo, convert = "cpm", transform = "log2",
norm = "quant", filter = TRUE))
plt <- plot_pca(eo_norm, plot_labels = FALSE)$plot
pp(file = glue("images/eo_pca_normalized-v{ver}.pdf"), image = plt)
eo_nb <- sm(normalize_expt(eo, convert = "cpm", transform = "log2",
filter = TRUE, batch = "svaseq"))
plot_pca(eo_nb)$plot
5.3.2 DE of Eosinophil samples
5.3.2.1 Withouth sva
eo_de <- sm(all_pairwise(eo, model_batch = FALSE, filter = TRUE))eo_tables <- sm(combine_de_tables(
eo_de, keepers = keepers,
excel = glue::glue("excel/eosinophil_clinical_all_tables-v{ver}.xlsx")))
written <- write_xlsx(data = eo_tables[["data"]][[1]],
excel = glue::glue("excel/eosinophil_clinical_table-v{ver}.xlsx"))
eo_sig <- sm(extract_significant_genes(eo_tables, according_to = "deseq"))
written <- write_xlsx(data = eo_sig[["deseq"]][["ups"]][[1]],
excel = glue::glue("excel/eosinophil_clinical_sigup-v{ver}.xlsx"))
written <- write_xlsx(data = eo_sig[["deseq"]][["downs"]][[1]],
excel = glue::glue("excel/eosinophil_clinical_sigdown-v{ver}.xlsx"))
eo_pct_sig <- sm(extract_significant_genes(eo_tables, n = 200,
lfc = NULL, p = NULL, according_to = "deseq"))
written <- write_xlsx(data = eo_pct_sig[["deseq"]][["ups"]][[1]],
excel = glue::glue("excel/eosinophil_clinical_sigup_pct-v{ver}.xlsx"))
written <- write_xlsx(data = eo_pct_sig[["deseq"]][["downs"]][[1]],
excel = glue::glue("excel/eosinophil_clinical_sigdown_pct-v{ver}.xlsx"))
eo_cpm <- sm(normalize_expt(eo, convert = "cpm"))
written <- write_xlsx(data = exprs(eo_cpm),
excel = glue::glue("excel/eosinophil_cpm_before_batch-v{ver}.xlsx"))
eo_bcpm <- sm(normalize_expt(eo, filter = TRUE, batch = "svaseq", convert = "cpm"))
written <- write_xlsx(data = exprs(eo_bcpm),
excel = glue::glue("excel/eosinophil_cpm_after_batch-v{ver}.xlsx"))5.3.2.2 With sva
eo_de_sva <- sm(all_pairwise(eo, model_batch = "svaseq", filter = TRUE))eo_tables_sva <- sm(combine_de_tables(
eo_de_sva, keepers = keepers,
excel = glue::glue("excel/eosinophil_clinical_all_tables_sva-v{ver}.xlsx")))
eo_sig_sva <- sm(extract_significant_genes(
eo_tables_sva,
excel = glue::glue("excel/eosinophil_clinical_sig_tables_sva-v{ver}.xlsx"),
according_to = "deseq"))5.3.2.3 Eosinophil DE plots
## DESeq2 MA plot of failure / cure
eo_tables[["plots"]][["fail_vs_cure"]][["deseq_ma_plots"]]$plot
## DESeq2 Volcano plot of failure / cure
eo_tables[["plots"]][["fail_vs_cure"]][["deseq_vol_plots"]]$plot
## DESeq2 MA plot of failure / cure with sva
eo_tables_sva[["plots"]][["fail_vs_cure"]][["deseq_ma_plots"]]$plot
## DESeq2 Volcano plot of failure / cure with sva
eo_tables_sva[["plots"]][["fail_vs_cure"]][["deseq_vol_plots"]]$plot
5.3.2.4 Eosinophil ontology search
ups <- eo_sig[["deseq"]][["ups"]][["fail_vs_cure"]]
downs <- eo_sig[["deseq"]][["downs"]][["fail_vs_cure"]]
eo_up_gp <- simple_gprofiler(sig_genes = ups)## Performing gProfiler GO search of 103 genes against hsapiens.## GO search found 97 hits.## Performing gProfiler KEGG search of 103 genes against hsapiens.## KEGG search found 9 hits.## Performing gProfiler REAC search of 103 genes against hsapiens.## REAC search found 9 hits.## Performing gProfiler MI search of 103 genes against hsapiens.## MI search found 0 hits.## Performing gProfiler TF search of 103 genes against hsapiens.## TF search found 50 hits.## Performing gProfiler CORUM search of 103 genes against hsapiens.## CORUM search found 0 hits.## Performing gProfiler HP search of 103 genes against hsapiens.## HP search found 1 hits.eo_up_gp[["pvalue_plots"]][["bpp_plot_over"]]
eo_up_gp[["pvalue_plots"]][["mfp_plot_over"]]
eo_up_gp[["pvalue_plots"]][["reactome_plot_over"]]
eo_down_gp <- simple_gprofiler(downs)## Performing gProfiler GO search of 40 genes against hsapiens.## GO search found 0 hits.## Performing gProfiler KEGG search of 40 genes against hsapiens.## KEGG search found 0 hits.## Performing gProfiler REAC search of 40 genes against hsapiens.## REAC search found 0 hits.## Performing gProfiler MI search of 40 genes against hsapiens.## MI search found 0 hits.## Performing gProfiler TF search of 40 genes against hsapiens.## TF search found 0 hits.## Performing gProfiler CORUM search of 40 genes against hsapiens.## CORUM search found 0 hits.## Performing gProfiler HP search of 40 genes against hsapiens.## HP search found 0 hits.eo_down_gp[["pvalue_plots"]][["bpp_plot_over"]]## NULLeo_down_gp[["pvalue_plots"]][["mfp_plot_over"]]## NULLeo_down_gp[["pvalue_plots"]][["reactome_plot_over"]]## NULL5.3.2.5 Eosinophil MSigDB query
eo_up_goseq_msig <- goseq_msigdb(sig_genes = ups, signatures = broad_c7,
signature_category = "c7", length_db = hs_length)## Before conversion: 103, after conversion: 103.## Before conversion: 227921, after conversion: 35341.## Found 101 go_db genes and 103 length_db genes out of 103.eo_down_goseq_msig <- goseq_msigdb(sig_genes = downs, signatures = broad_c7,
signature_category = "c7", length_db = hs_length)## Before conversion: 40, after conversion: 39.## Before conversion: 227921, after conversion: 35341.## Found 34 go_db genes and 39 length_db genes out of 39.5.3.2.6 Plot of similar experiments
## Eosinophil genes with increased expression in the failed samples
## share genes with the following experiments
eo_up_goseq_msig[["pvalue_plots"]][["mfp_plot_over"]]
## Eosinophil genes with increased expression in the cured samples
## share genes with the following experiments
eo_down_goseq_msig[["pvalue_plots"]][["mfp_plot_over"]]
5.4 Biopsies
5.4.1 Evaluate Biopsy samples
biop <- subset_expt(hs_valid, subset = "typeofcells=='biopsy'") %>%
set_expt_conditions(fact = "clinicaloutcome") %>%
set_expt_batches(fact = "donor") %>%
set_expt_colors(colors = chosen_colors)## subset_expt(): There were 135, now there are 42 samples.save_result <- save(biop, file = "rda/biopsy_expt.rda")
biop_norm <- normalize_expt(biop, filter = TRUE, convert = "cpm",
transform = "log2", norm = "quant")## Removing 5799 low-count genes (14142 remaining).## transform_counts: Found 5 values equal to 0, adding 1 to the matrix.plt <- plot_pca(biop_norm, plot_labels = FALSE)$plot
pp(file = glue("images/biop_pca_normalized-v{ver}.pdf"), image = plt)
biop_nb <- sm(normalize_expt(biop, convert = "cpm", filter = TRUE,
transform = "log2", batch = "svaseq"))
plt <- plot_pca(biop_nb, plot_labels = FALSE)$plot
pp(file = glue("images/biop_pca_normalized_svaseq-v{ver}.pdf"), image = plt)
5.4.2 DE of Biopsy samples
5.4.2.1 Without sva
biop_de <- sm(all_pairwise(biop, model_batch = FALSE, filter = TRUE))biop_tables <- combine_de_tables(biop_de, keepers = keepers,
excel = glue::glue("excel/biopsy_clinical_all_tables-v{ver}.xlsx"))## Deleting the file excel/biopsy_clinical_all_tables-v202106.xlsx before writing the tables.written <- write_xlsx(data = biop_tables[["data"]][[1]],
excel = glue::glue("excel/biopsy_clinical_table-v{ver}.xlsx"))
biop_sig <- extract_significant_genes(biop_tables, according_to = "deseq")
##written <- write_xlsx(data = biop_sig[["deseq"]][["ups"]][[1]],
## excel = glue::glue("excel/biopsy_clinical_sigup-v{ver}.xlsx"))
written <- write_xlsx(data = biop_sig[["deseq"]][["downs"]][[1]],
excel = glue::glue("excel/biopsy_clinical_sigdown-v{ver}.xlsx"))
biop_pct_sig <- extract_significant_genes(biop_tables, n = 200, lfc = NULL, p = NULL, according_to = "deseq")## Getting the top and bottom 200 genes.written <- write_xlsx(data = biop_pct_sig[["deseq"]][["ups"]][[1]],
excel = glue::glue("excel/biopsy_clinical_sigup_pct-v{ver}.xlsx"))
written <- write_xlsx(data = biop_pct_sig[["deseq"]][["downs"]][[1]],
excel = glue::glue("excel/biopsy_clinical_sigdown_pct-v{ver}.xlsx"))
biop_cpm <- sm(normalize_expt(biop, convert = "cpm"))
written <- write_xlsx(data = exprs(biop_cpm),
excel = glue::glue("excel/biopsy_cpm_before_batch-v{ver}.xlsx"))
biop_bcpm <- sm(normalize_expt(biop, filter = TRUE, batch = "svaseq", convert = "cpm"))
written <- write_xlsx(data = exprs(biop_bcpm),
excel = glue::glue("excel/biopsy_cpm_after_batch-v{ver}.xlsx"))5.4.2.2 with sva
biop_de_sva <- sm(all_pairwise(biop, model_batch = "svaseq", filter = TRUE))biop_tables_sva <- sm(combine_de_tables(
biop_de_sva, keepers = keepers,
excel = glue::glue("excel/biopsy_clinical_all_tables_sva-v{ver}.xlsx")))
biop_sig_sva <- sm(extract_significant_genes(
biop_tables_sva,
excel = glue::glue("excel/biopsy_clinical_sig_tables_sva-v{ver}.xlsx"),
according_to = "deseq"))5.4.2.3 Biopsy DE plots
## DESeq2 MA plot of failure / cure
biop_tables[["plots"]][["fail_vs_cure"]][["deseq_ma_plots"]]$plot
## DESeq2 Volcano plot of failure / cure
biop_tables[["plots"]][["fail_vs_cure"]][["deseq_vol_plots"]]$plot
## DESeq2 MA plot of failure / cure
biop_tables_sva[["plots"]][["fail_vs_cure"]][["deseq_ma_plots"]]$plot
## DESeq2 Volcano plot of failure / cure
biop_tables_sva[["plots"]][["fail_vs_cure"]][["deseq_vol_plots"]]$plot
7 Monocytes by visit
- Can you please share with us a PCA (SVA and non-SVA) of the monocytes of the TMRC3 project, but labeling them based on the visit (V1, V2, V3)?
- Can you please share DE lists of V1 vs V2, V1 vs V3, V1 vs. V2+V3 and V2 vs V3?
visit_colors <- chosen_colors <- c("#D95F02", "#7570B3", "#1B9E77")
names(visit_colors) <- c(1, 2, 3)
mono_visit <- subset_expt(hs_valid, subset = "typeofcells=='monocytes'") %>%
set_expt_conditions(fact = "visitnumber") %>%
set_expt_batches(fact = "clinicaloutcome") %>%
set_expt_colors(colors = chosen_colors)## subset_expt(): There were 135, now there are 27 samples.mono_visit_norm <- normalize_expt(mono_visit, filter = TRUE, norm = "quant", convert = "cpm",
transform = "log2")## Removing 8915 low-count genes (11026 remaining).## transform_counts: Found 9 values equal to 0, adding 1 to the matrix.mono_visit_pca <- plot_pca(mono_visit_norm)
pp(file = "images/monocyte_by_visit.png", image = mono_visit_pca$plot)
mono_visit_nb <- normalize_expt(mono_visit, filter = TRUE, convert = "cpm",
batch = "svaseq", transform = "log2")## Removing 8915 low-count genes (11026 remaining).## batch_counts: Before batch/surrogate estimation, 1404 entries are x==0: 0%.## batch_counts: Before batch/surrogate estimation, 17611 entries are 0<x<1: 6%.## Setting 441 low elements to zero.## transform_counts: Found 441 values equal to 0, adding 1 to the matrix.mono_visit_nb_pca <- plot_pca(mono_visit_nb)
pp(file = "images/monocyte_by_visit_nb.png", image = mono_visit_nb_pca$plot)
table(pData(mono_visit_norm)$batch)##
## cure failure lost
## 9 12 6keepers <- list(
"second_vs_first" = c("c2", "c1"),
"third_vs_second" = c("c3", "c2"),
"third_vs_first" = c("c3", "c1"))
mono_visit_de <- all_pairwise(mono_visit, model_batch = "svaseq", filter = TRUE)## batch_counts: Before batch/surrogate estimation, 1404 entries are x==0: 0%.## Plotting a PCA before surrogate/batch inclusion.## Using svaseq to visualize before/after batch inclusion.## Performing a test normalization with: raw## Removing 0 low-count genes (11026 remaining).## batch_counts: Before batch/surrogate estimation, 1404 entries are x==0: 0%.## batch_counts: Before batch/surrogate estimation, 17611 entries are 0<x<1: 6%.## Setting 441 low elements to zero.## transform_counts: Found 441 values equal to 0, adding 1 to the matrix.## Finished running DE analyses, collecting outputs.## Comparing analyses.mono_visit_tables <- combine_de_tables(
mono_visit_de,
keepers = keepers,
excel = glue::glue("excel/mono_visit_tables-v{ver}.xlsx"))new_factor <- as.character(pData(mono_visit)[["visitnumber"]])
not_one_idx <- new_factor != 1
new_factor[not_one_idx] <- "not_1"
mono_one_vs <- set_expt_conditions(mono_visit, new_factor)
mono_one_vs_de <- all_pairwise(mono_one_vs, model_batch = "svaseq", filter = TRUE)## batch_counts: Before batch/surrogate estimation, 1404 entries are x==0: 0%.## Plotting a PCA before surrogate/batch inclusion.## Using svaseq to visualize before/after batch inclusion.## Performing a test normalization with: raw## Removing 0 low-count genes (11026 remaining).## batch_counts: Before batch/surrogate estimation, 1404 entries are x==0: 0%.## batch_counts: Before batch/surrogate estimation, 17611 entries are 0<x<1: 6%.## Setting 436 low elements to zero.## transform_counts: Found 436 values equal to 0, adding 1 to the matrix.## Finished running DE analyses, collecting outputs.## Comparing analyses.mono_one_vs_tables <- combine_de_tables(
mono_one_vs_de,
excel = glue::glue("excel/mono_one_vs_tables-v{ver}.xlsx"))8 Test TSP
In writing the following, I quickly realized that tspair was not joking when it said it is intended for small numbers of genes. For a full expressionset of human data it is struggling. I like the idea, it may prove worth while to spend some time optimizing the package so that it is more usable.
expt <- hs_clinical_nobiop
simple_tsp <- function(expt, column = "condition") {
facts <- levels(as.factor(pData(expt)[[column]]))
retlist <- list()
if (length(facts) < 2) {
stop("This requires factors with at least 2 levels.")
} else if (length(facts) == 2) {
retlist <- simple_tsp_pair(expt, column = column)
} else {
for (first in 1:(length(facts) - 1)) {
for (second in 2:(length(facts))) {
if (first < second) {
name <- glue::glue("{facts[first]}_vs_{facts[second]}")
message("Starting ", name, ".")
substring <- glue::glue("{column}=='{facts[first]}'|{column}=='{facts[second]}'")
subby <- subset_expt(expt, subset=as.character(substring))
retlist[[name]] <- simple_tsp_pair(subby, column = column)
}
}
}
}
}
simple_tsp_pair <- function(subby, column = "condition", repetitions = 50) {
tsp_input <- subby[["expressionset"]]
tsp_output <- tspcalc(tsp_input, column)
tsp_scores <- tspsig(tsp_input, column, B = repetitions)
}
tsp1 <- tspcalc(tsp_input, "condition")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))
}## If you wish to reproduce this exact build of hpgltools, invoke the following:## > git clone http://github.com/abelew/hpgltools.git## > git reset 72947fcc6afe09da22d71967059edd84e3063341## This is hpgltools commit: Tue Jun 1 15:57:56 2021 -0400: 72947fcc6afe09da22d71967059edd84e3063341## Saving to tmrc3_02sample_estimation_v202106.rda.xztmp <- loadme(filename = savefile)LS0tCnRpdGxlOiAiVE1SQzMgQ29tcHJlaGVuc2l2ZSBEYXRhIEFuYWx5c2lzOiAyMDIxMDYiCmF1dGhvcjogImF0YiBhYmVsZXdAZ21haWwuY29tIgpkYXRlOiAiYHIgU3lzLkRhdGUoKWAiCm91dHB1dDoKIGh0bWxfZG9jdW1lbnQ6CiAgY29kZV9kb3dubG9hZDogdHJ1ZQogIGNvZGVfZm9sZGluZzogc2hvdwogIGZpZ19jYXB0aW9uOiB0cnVlCiAgZmlnX2hlaWdodDogNwogIGZpZ193aWR0aDogNwogIGhpZ2hsaWdodDogZGVmYXVsdAogIGtlZXBfbWQ6IGZhbHNlCiAgbW9kZTogc2VsZmNvbnRhaW5lZAogIG51bWJlcl9zZWN0aW9uczogdHJ1ZQogIHNlbGZfY29udGFpbmVkOiB0cnVlCiAgdGhlbWU6IHJlYWRhYmxlCiAgdG9jOiB0cnVlCiAgdG9jX2Zsb2F0OgogICBjb2xsYXBzZWQ6IGZhbHNlCiAgIHNtb290aF9zY3JvbGw6IGZhbHNlCi0tLQoKPHN0eWxlPgogIGJvZHkgLm1haW4tY29udGFpbmVyIHsKICAgIG1heC13aWR0aDogMTYwMHB4OwogIH0KPC9zdHlsZT4KCmBgYHtyIG9wdGlvbnMsIGluY2x1ZGUgPSBGQUxTRX0KbGlicmFyeShocGdsdG9vbHMpCnR0IDwtIHNtKGRldnRvb2xzOjpsb2FkX2FsbCgifi9ocGdsdG9vbHMiKSkKa25pdHI6Om9wdHNfa25pdCRzZXQocHJvZ3Jlc3MgPSBUUlVFLAogICAgICAgICAgICAgICAgICAgICB2ZXJib3NlID0gVFJVRSwKICAgICAgICAgICAgICAgICAgICAgd2lkdGggPSAxMjAsCiAgICAgICAgICAgICAgICAgICAgIGVjaG8gPSBUUlVFKQprbml0cjo6b3B0c19jaHVuayRzZXQoZXJyb3IgPSBUUlVFLAogICAgICAgICAgICAgICAgICAgICAgZmlnLndpZHRoID0gMTIsCiAgICAgICAgICAgICAgICAgICAgICBmaWcuaGVpZ2h0ID0gMTIsCiAgICAgICAgICAgICAgICAgICAgICBkcGkgPSA5NikKb2xkX29wdGlvbnMgPC0gb3B0aW9ucyhkaWdpdHMgPSA0LAogICAgICAgICAgICAgICAgICAgICAgIHN0cmluZ3NBc0ZhY3RvcnMgPSBGQUxTRSwKICAgICAgICAgICAgICAgICAgICAgICBrbml0ci5kdXBsaWNhdGUubGFiZWwgPSAiYWxsb3ciKQpnZ3Bsb3QyOjp0aGVtZV9zZXQoZ2dwbG90Mjo6dGhlbWVfYncoYmFzZV9zaXplID0gMTIpKQp2ZXIgPC0gIjIwMjEwNiIKcnVuZGF0ZSA8LSBmb3JtYXQoU3lzLkRhdGUoKSwgZm9ybWF0ID0gIiVZJW0lZCIpCgpybWRfZmlsZSA8LSBnbHVlOjpnbHVlKCJ0bXJjM18wMnNhbXBsZV9lc3RpbWF0aW9uX3Z7dmVyfS5SbWQiKQpzYXZlZmlsZSA8LSBnc3ViKHBhdHRlcm4gPSAiXFwuUm1kIiwgcmVwbGFjZSA9ICJcXC5yZGFcXC54eiIsIHggPSBybWRfZmlsZSkKYGBgCgojIEFubm90YXRpb24KCldlIHRha2UgdGhlIGFubm90YXRpb24gZGF0YSBmcm9tIGVuc2VtYmwncyBiaW9tYXJ0IGluc3RhbmNlLiAgVGhlIGdlbm9tZSB3aGljaAp3YXMgdXNlZCB0byBtYXAgdGhlIGRhdGEgd2FzIGhnMzggcmV2aXNpb24gMTAwLiAgTXkgZGVmYXVsdCB3aGVuIHVzaW5nIGJpb21hcnQgaXMKdG8gbG9hZCB0aGUgZGF0YSBmcm9tIDEgeWVhciBiZWZvcmUgdGhlIGN1cnJlbnQgZGF0ZS4KCmBgYHtyIGhzX2Fubm90fQpoc19hbm5vdCA8LSBzbShsb2FkX2Jpb21hcnRfYW5ub3RhdGlvbnMoeWVhciA9ICIyMDIwIikpCmhzX2Fubm90IDwtIGhzX2Fubm90W1siYW5ub3RhdGlvbiJdXQpoc19hbm5vdFtbInRyYW5zY3JpcHQiXV0gPC0gcGFzdGUwKHJvd25hbWVzKGhzX2Fubm90KSwgIi4iLCBoc19hbm5vdFtbInZlcnNpb24iXV0pCnJvd25hbWVzKGhzX2Fubm90KSA8LSBtYWtlLm5hbWVzKGhzX2Fubm90W1siZW5zZW1ibF9nZW5lX2lkIl1dLCB1bmlxdWUgPSBUUlVFKQp0eF9nZW5lX21hcCA8LSBoc19hbm5vdFssIGMoInRyYW5zY3JpcHQiLCAiZW5zZW1ibF9nZW5lX2lkIildCgpzdW1tYXJ5KGhzX2Fubm90KQpgYGAKCmBgYHtyIGhzX2dvfQpoc19nbyA8LSBzbShsb2FkX2Jpb21hcnRfZ28oKVtbImdvIl1dKQpoc19sZW5ndGggPC0gaHNfYW5ub3RbLCBjKCJlbnNlbWJsX2dlbmVfaWQiLCAiY2RzX2xlbmd0aCIpXQpjb2xuYW1lcyhoc19sZW5ndGgpIDwtIGMoIklEIiwgImxlbmd0aCIpCmBgYAoKIyBJbnRyb2R1Y3Rpb24KClRoaXMgZG9jdW1lbnQgaXMgaW50ZW5kZWQgdG8gcHJvdmlkZSBhbiBvdmVydmlldyBvZiBUTVJDMyBzYW1wbGVzIHdoaWNoIGhhdmUKYmVlbiBzZXF1ZW5jZWQuICBJdCBpbmNsdWRlcyBzb21lIHBsb3RzIGFuZCBhbmFseXNlcyBzaG93aW5nIHRoZSByZWxhdGlvbnNoaXBzCmFtb25nIHRoZSBzYW1wbGVzIGFzIHdlbGwgYXMgc29tZSBkaWZmZXJlbnRpYWwgYW5hbHlzZXMgd2hlbiBwb3NzaWJsZS4KCiMgU2FtcGxlIEVzdGltYXRpb24KCiMjIEdlbmVyYXRlIGV4cHJlc3Npb25zZXRzCgpUaGUgc2FtcGxlIHNoZWV0IGlzIGNvcGllZCBmcm9tIG91ciBzaGFyZWQgb25saW5lIHNoZWV0IGFuZCB1cGRhdGVkIHdpdGggZWFjaCByZWxlYXNlCm9mIHNlcXVlbmNpbmcgZGF0YS4KCmBgYHtyIHNhbXBsZXNoZWV0fQpzYW1wbGVzaGVldCA8LSAic2FtcGxlX3NoZWV0cy90bXJjM19zYW1wbGVzXzIwMjEwNjAxLnhsc3giCmBgYAoKIyMjIEhpc2F0MiBleHByZXNzaW9uc2V0cwoKVGhlIGZpcnN0IHRoaW5nIHRvIG5vdGUgaXMgdGhlIGxhcmdlIHJhbmdlIGluIGNvdmVyYWdlLiAgVGhlcmUgYXJlIG11bHRpcGxlCnNhbXBsZXMgd2l0aCBjb3ZlcmFnZSB3aGljaCBpcyB0b28gbG93IHRvIHVzZS4gIFRoZXNlIHdpbGwgYmUgcmVtb3ZlZCBzaG9ydGx5LgoKSW4gdGhlIGZvbGxvd2luZyBibG9jayBJIGltbWVkaWF0ZWx5IGV4Y2x1ZGUgYW55IG5vbi1jb2RpbmcgcmVhZHMgYXMgd2VsbC4KCmBgYHtyIGFsbF9uZXdfaGlzYXQyfQojIyBDcmVhdGUgdGhlIGV4cHJlc3Npb25zZXQgYW5kIGltbWVkaWF0ZWx5IHBhc3MgaXQgdG8gYSBmaWx0ZXIKIyMgcmVtb3ZpbmcgdGhlIG5vbiBwcm90ZWluIGNvZGluZyBnZW5lcy4Kc2FuaXRpemVfY29sdW1ucyA8LSBjKCJ2aXNpdG51bWJlciIsICJjbGluaWNhbG91dGNvbWUiLCAiZG9ub3IiLAogICAgICAgICAgICAgICAgICAgICAgInR5cGVvZmNlbGxzIiwgImNsaW5pY2FscHJlc2VudGF0aW9uIiwKICAgICAgICAgICAgICAgICAgICAgICJjb25kaXRpb24iLCAiYmF0Y2giKQpoc19leHB0IDwtIGNyZWF0ZV9leHB0KHNhbXBsZXNoZWV0LAogICAgICAgICAgICAgICAgICAgICAgIGZpbGVfY29sdW1uID0gImhnMzgxMDBoaXNhdGZpbGUiLAogICAgICAgICAgICAgICAgICAgICAgIHNhdmVmaWxlID0gZ2x1ZTo6Z2x1ZSgicmRhL2hzX2V4cHRfYWxsLXZ7dmVyfS5yZGEiKSwKICAgICAgICAgICAgICAgICAgICAgICBnZW5lX2luZm8gPSBoc19hbm5vdCkgJT4lCiAgZXhjbHVkZV9nZW5lc19leHB0KGNvbHVtbiA9ICJnZW5lX2Jpb3R5cGUiLCBtZXRob2QgPSAia2VlcCIsCiAgICAgICAgICAgICAgICAgICAgIHBhdHRlcm4gPSAicHJvdGVpbl9jb2RpbmciLCBtZXRhX2NvbHVtbiA9ICJuY3JuYV9sb3N0IikgJT4lCiAgc2FuaXRpemVfZXhwdF9tZXRhZGF0YShjb2x1bW5zID0gc2FuaXRpemVfY29sdW1ucykgJT4lCiAgc2V0X2V4cHRfZmFjdG9ycyhjb2x1bW5zID0gc2FuaXRpemVfY29sdW1ucywgY2xhc3MgPSAiZmFjdG9yIikKCmxldmVscyhwRGF0YShoc19leHB0W1siZXhwcmVzc2lvbnNldCJdXSlbWyJ2aXNpdG51bWJlciJdXSkgPC0gbGlzdCgKICAgICcwJyA9ICJub3RhcHBsaWNhYmxlIiwgJzEnID0gMSwgJzInID0gMiwgJzMnID0gMykKYGBgCgpTcGxpdCB0aGlzIGRhdGEgaW50byBDRFMgYW5kIGxuY1JOQS4gIE9oIGNyYXAgaW4gb3JkZXIgdG8gZG8gdGhhdCBJIG5lZWQgdG8gcmVjb3VudCB0aGUgZGF0YS4KUnVubmluZyBub3cgKDIwMjEwNTE4KQoKYGBge3IgbG5jX2Nkc30KIyMgbG5jX2V4cHQgPC0gY3JlYXRlX2V4cHQoc2FtcGxlc2hlZXQsCiMjICAgICAgICAgICAgICAgICAgICAgICAgIGZpbGVfY29sdW1uID0gImhnMzgxMDBsbmNmaWxlIiwKIyMgICAgICAgICAgICAgICAgICAgICAgICAgZ2VuZV9pbmZvID0gaHNfYW5ub3QpCmBgYAoKIyMjIyBJbml0aWFsIG1ldHJpY3MKCk9uY2UgdGhlIGRhdGEgd2FzIGxvYWRlZCwgdGhlcmUgYXJlIGEgY291cGxlIG9mIG1ldHJpY3Mgd2hpY2ggbWF5IGJlIHBsb3R0ZWQgaW1tZWRpYXRlbHkuCgpgYGB7ciBpbml0aWFsX21ldHJpY3N9Cm5vbnplcm8gPC0gcGxvdF9ub256ZXJvKGhzX2V4cHQpCm5vbnplcm8kcGxvdAoKbmNybmFfbG9zdF9kZiA8LSBhcy5kYXRhLmZyYW1lKHBEYXRhKGhzX2V4cHQpW1sibmNybmFfbG9zdCJdXSkKcm93bmFtZXMobmNybmFfbG9zdF9kZikgPC0gcm93bmFtZXMocERhdGEoaHNfZXhwdCkpCmNvbG5hbWVzKG5jcm5hX2xvc3RfZGYpIDwtICJuY3JuYV9sb3N0IgoKdG1wZGYgPC0gbWVyZ2Uobm9uemVybyR0YWJsZSwgbmNybmFfbG9zdF9kZiwgYnkgPSAicm93Lm5hbWVzIikKcm93bmFtZXModG1wZGYpIDwtIHRtcGRmW1siUm93Lm5hbWVzIl1dCnRtcGRmW1siUm93Lm5hbWVzIl1dIDwtIE5VTEwKCmdncGxvdCh0bXBkZiwgYWVzKHg9bmNybmFfbG9zdCwgeT1ub256ZXJvX2dlbmVzKSkgKwogIGdncGxvdDI6Omdlb21fcG9pbnQoKSArCiAgZ2dwbG90Mjo6Z2d0aXRsZSgiTm9uemVybyBnZW5lcyB3aXRoIHJlc3BlY3QgdG8gcGVyY2VudCBjb3VudHMgCmxvc3Qgd2hlbiBuY1JOQSB3YXMgcmVtb3ZlZC4iKQpgYGAKCk5hamliIGRvZXNuJ3Qgd2FudCB0aGlzIHBsb3QsIGJ1dCBJIGFtIHVzaW5nIGl0IHRvIGNoZWNrIG5ldyBzYW1wbGVzLApzbyB3aWxsIGhpZGUgaXQgZnJvbSBnZW5lcmFsIHVzZS4KCmBgYHtyIGxpYnNpemV9CmxpYnNpemUgPC0gcGxvdF9saWJzaXplKGhzX2V4cHQpCmxpYnNpemUkcGxvdApgYGAKCiMjIE1pbmltdW0gY292ZXJhZ2Ugc2FtcGxlIGZpbHRlcmluZwoKSSBhcmJpdHJhcmlseSBjaG9zZSAxMSwwMDAgbm9uLXplcm8gZ2VuZXMgYXMgYSBtaW5pbXVtLiAgV2UgbWF5CndhbnQgdGhpcyB0byBiZSBoaWdoZXIuCgpgYGB7ciBoaXNhdDJfd3JpdGUsIGZpZy5zaG93ID0gImhpZGUifQpoc192YWxpZCA8LSBzdWJzZXRfZXhwdChoc19leHB0LCBub256ZXJvID0gMTEwMDApCgp2YWxpZF93cml0ZSA8LSBzbSh3cml0ZV9leHB0KGhzX3ZhbGlkLCBleGNlbCA9IGdsdWUoImV4Y2VsL2hzX3ZhbGlkLXZ7dmVyfS54bHN4IikpKQpgYGAKCiMgUHJvamVjdCBBaW1zCgpUaGUgcHJvamVjdCBzZWVrcyB0byBkZXRlcm1pbmUgdGhlIHJlbGF0aW9uc2hpcCBvZiB0aGUgaW5uYXRlIGltbXVuZSByZXNwb25zZQphbmQgaW5mbGFtbWF0b3J5IHNpZ25hbGluZyB0byB0aGUgY2xpbmljYWwgb3V0Y29tZSBvZiBhbnRpbGVpc2htYW5pYWwgZHJ1Zwp0cmVhdG1lbnQuIFdlIHdpbGwgdGVzdCB0aGUgaHlwb3RoZXNpcyB0aGF0IHRoZSBwcm9maWxlIG9mIGlubmF0ZSBpbW11bmUgY2VsbAphY3RpdmF0aW9uIGFuZCB0aGVpciBkeW5hbWljcyB0aHJvdWdoIHRoZSBjb3Vyc2Ugb2YgdHJlYXRtZW50IGRpZmZlciBiZXR3ZWVuIENMCnBhdGllbnRzIHdpdGggcHJvc3BlY3RpdmVseSBkZXRlcm1pbmVkIHRoZXJhcGV1dGljIGN1cmUgb3IgZmFpbHVyZS4KClRoaXMgd2lsbCBiZSBhY2hpZXZlZCB0aHJvdWdoIHRoZSBjaGFyYWN0ZXJpemF0aW9uIG9mIHRoZSBpbiB2aXZvIGR5bmFtaWNzIG9mCmJsb29kLWRlcml2ZWQgbW9ub2N5dGUsIG5ldXRyb3BoaWwgYW5kIGVvc2lub3BoaWwgdHJhbnNjcmlwdG9tZSBiZWZvcmUsIGR1cmluZwphbmQgYXQgdGhlIGVuZCBvZiB0cmVhdG1lbnQgaW4gQ0wgcGF0aWVudHMuIENlbGwtdHlwZSBzcGVjaWZpYyB0cmFuc2NyaXB0b21lcywKY29tcG9zaXRlIHNpZ25hdHVyZXMgYW5kIHRpbWUtcmVzcG9uc2UgZXhwcmVzc2lvbiBwcm9maWxlcyB3aWxsIGJlIGNvbnRyYXN0ZWQKYW1vbmcgcGF0aWVudHMgd2l0aCB0aGVyYXBldXRpYyBjdXJlIG9yIGZhaWx1cmUuCgojIyBQcmVwYXJhdGlvbgoKVG8gYWRkcmVzcyB0aGVzZSwgSSBhZGRlZCB0byB0aGUgZW5kIG9mIHRoZSBzYW1wbGUgc2hlZXQgY29sdW1ucyBuYW1lZAonY29uZGl0aW9uJywgJ2JhdGNoJywgJ2Rvbm9yJywgYW5kICd0aW1lJy4gIFRoZXNlIGFyZSBmaWxsZWQgaW4gd2l0aCBzaG9ydGhhbmQKdmFsdWVzIGFjY29yZGluZyB0byB0aGUgYWJvdmUuCgojIyBHbG9iYWwgdmlldwoKQmVmb3JlIGFkZHJlc3NpbmcgdGhlIHF1ZXN0aW9ucyBleHBsaWNpdGx5IGJ5IHN1YnNldHRpbmcgdGhlIGRhdGEsIEkgd2FudCB0byBnZXQKYSBsb29rIGF0IHRoZSBzYW1wbGVzIGFzIHRoZXkgYXJlLgoKYGBge3IgcHJlX3F1ZXN0aW9uc30KbmV3X25hbWVzIDwtIHBEYXRhKGhzX3ZhbGlkKVtbInNhbXBsZW5hbWUiXV0KaHNfdmFsaWQgPC0gaHNfdmFsaWQgJT4lCiAgc2V0X2V4cHRfYmF0Y2hlcyhmYWN0ID0gImNlbGxzc291cmNlIikgJT4lCiAgc2V0X2V4cHRfY29uZGl0aW9ucyhmYWN0ID0gInR5cGVvZmNlbGxzIikgJT4lCiAgc2V0X2V4cHRfc2FtcGxlbmFtZXMobmV3bmFtZXMgPSBuZXdfbmFtZXMpCgphbGxfbm9ybSA8LSBzbShub3JtYWxpemVfZXhwdChoc192YWxpZCwgdHJhbnNmb3JtID0gImxvZzIiLCBub3JtID0gInF1YW50IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29udmVydCA9ICJjcG0iLCBmaWx0ZXIgPSBUUlVFKSkKCmFsbF9wY2EgPC0gcGxvdF9wY2EoYWxsX25vcm0sIHBsb3RfbGFiZWxzID0gRkFMU0UsCiAgICAgICAgICAgICAgICAgICAgcGxvdF90aXRsZSA9ICJQQ0EgLSBDZWxsIHR5cGUiLCBzaXplX2NvbHVtbiA9ICJ2aXNpdG51bWJlciIpCnBwKGZpbGUgPSBnbHVlKCJpbWFnZXMvdG1yYzNfcGNhX25vbGFiZWxzLXZ7dmVyfS5wbmciKSwgaW1hZ2UgPSBhbGxfcGNhJHBsb3QpCgp3cml0ZS5jc3YoYWxsX3BjYSR0YWJsZSwgZmlsZSA9ICJjb29yZHMvaHNfZG9ub3JfcGNhX2Nvb3Jkcy5jc3YiKQpwbG90X2NvcmhlYXQoYWxsX25vcm0sIHBsb3RfdGl0bGUgPSAiSGVpcmFyY2hpY2FsIGNsdXN0ZXJpbmc6CiAgICAgICAgIGNlbGwgdHlwZXMiKSRwbG90CmBgYAoKIyMgRXhhbWluZSBzYW1wbGVzIHJlbGV2YW50IHRvIGNsaW5pY2FsIG91dGNvbWUKCk5vdyBsZXQgdXMgY29uc2lkZXIgb25seSB0aGUgc2FtcGxlcyBmb3Igd2hpY2ggd2UgaGF2ZSBhIGNsaW5pY2FsIG91dGNvbWUuClRoZXNlIGZhbGwgcHJpbWFyaWx5IGludG8gZWl0aGVyICdjdXJlZCcgb3IgJ2ZhaWxlZCcsIGJ1dCBzb21lIHBlb3BsZSBoYXZlIG5vdAp5ZXQgcmV0dXJuZWQgdG8gdGhlIGNsaW5pYyBhZnRlciB0aGUgZmlyc3Qgb3Igc2Vjb25kIHZpc2l0LiAgVGhlc2UgYXJlIGRlZW1lZAonbG9zdCcuCgpgYGB7ciBhbGxfY2xpbmljYWx9CmhzX2NsaW5pY2FsIDwtIGhzX3ZhbGlkICU+JQogIHNldF9leHB0X2NvbmRpdGlvbnMoZmFjdCA9ICJjbGluaWNhbG91dGNvbWUiKSAlPiUKICBzZXRfZXhwdF9iYXRjaGVzKGZhY3QgPSAidHlwZW9mY2VsbHMiKSAlPiUKICBzdWJzZXRfZXhwdChzdWJzZXQgPSAidHlwZW9mY2VsbHMhPSdwYm1jcycmdHlwZW9mY2VsbHMhPSdtYWNyb3BoYWdlcyciKQoKY2hvc2VuX2NvbG9ycyA8LSBjKCIjRDk1RjAyIiwgIiM3NTcwQjMiLCAiIzFCOUU3NyIsICIjRkYwMDAwIiwgIiNGRjAwMDAiKQpuYW1lcyhjaG9zZW5fY29sb3JzKSA8LSBjKCJjdXJlIiwgImZhaWx1cmUiLCAibG9zdCIsICJudWxsIiwgIm5vdGFwcGxpY2FibGUiKQpoc19jbGluaWNhbCA8LSBzZXRfZXhwdF9jb2xvcnMoaHNfY2xpbmljYWwsIGNvbG9ycyA9IGNob3Nlbl9jb2xvcnMpCgpuZXduYW1lcyA8LSBtYWtlLm5hbWVzKHBEYXRhKGhzX2NsaW5pY2FsKVtbInNhbXBsZW5hbWUiXV0sIHVuaXF1ZSA9IFRSVUUpCmhzX2NsaW5pY2FsIDwtIHNldF9leHB0X3NhbXBsZW5hbWVzKGhzX2NsaW5pY2FsLCBuZXduYW1lcyA9IG5ld25hbWVzKQoKaHNfY2xpbmljYWxfbm9ybSA8LSBzbShub3JtYWxpemVfZXhwdChoc19jbGluaWNhbCwgZmlsdGVyID0gVFJVRSwgdHJhbnNmb3JtID0gImxvZzIiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbnZlcnQgPSAiY3BtIiwgbm9ybSA9ICJxdWFudCIpKQpjbGluaWNhbF9wY2EgPC0gcGxvdF9wY2EoaHNfY2xpbmljYWxfbm9ybSwgcGxvdF9sYWJlbHMgPSBGQUxTRSwKICAgICAgICAgICAgICAgICAgICAgICAgIHNpemVfY29sdW1uID0gInZpc2l0bnVtYmVyIiwgY2lzID0gTlVMTCwKICAgICAgICAgICAgICAgICAgICAgICAgIHBsb3RfdGl0bGUgPSAiUENBIC0gY2xpbmljYWwgc2FtcGxlcyIpCnBwKGZpbGUgPSBnbHVlKCJpbWFnZXMvYWxsX2NsaW5pY2FsX25vYmF0Y2hfcGNhLXZ7dmVyfS5wbmciKSwgaW1hZ2UgPSBjbGluaWNhbF9wY2EkcGxvdCwKICAgaGVpZ2h0ID0gOCwgd2lkdGggPSAyMCkKYGBgCgojIyMgUmVwZWF0IHdpdGhvdXQgdGhlIGJpb3BzeSBzYW1wbGVzCgpgYGB7ciBpYmlkX25vYmlvcHN5fQpoc19jbGluaWNhbF9ub2Jpb3AgPC0gaHNfY2xpbmljYWwgJT4lCiAgc3Vic2V0X2V4cHQoc3Vic2V0ID0gInR5cGVvZmNlbGxzIT0nYmlvcHN5JyIpCgpoc19jbGluaWNhbF9ub2Jpb3Bfbm9ybSA8LSBzbShub3JtYWxpemVfZXhwdChoc19jbGluaWNhbF9ub2Jpb3AsIGZpbHRlciA9IFRSVUUsIHRyYW5zZm9ybSA9ICJsb2cyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29udmVydCA9ICJjcG0iLCBub3JtID0gInF1YW50IikpCmNsaW5pY2FsX25vYmlvcF9wY2EgPC0gcGxvdF9wY2EoaHNfY2xpbmljYWxfbm9iaW9wX25vcm0sIHBsb3RfbGFiZWxzID0gRkFMU0UsIGNpcyA9IE5VTEwsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgcGxvdF90aXRsZSA9ICJQQ0EgLSBjbGluaWNhbCBzYW1wbGVzIHdpdGhvdXQgYmlvcHNpZXMiKQpwcChmaWxlID0gZ2x1ZSgiaW1hZ2VzL2FsbF9jbGluaWNhbF9ub2Jpb3Bfbm9iYXRjaF9wY2Etdnt2ZXJ9LnBuZyIpLAogICBpbWFnZSA9IGNsaW5pY2FsX25vYmlvcF9wY2EkcGxvdCkKYGBgCgojIyMgQXR0ZW1wdCB0byBjb3JyZWN0IGZvciB0aGUgc3Vycm9nYXRlIHZhcmlhYmxlcwoKQXQgdGhpcyB0aW1lIHdlIGhhdmUgdHdvIHByaW1hcnkgZGF0YSBzdHJ1Y3R1cmVzIG9mIGludGVyZXN0OiBoc19jbGluaWNhbCBhbmQgaHNfY2xpbmljYWxfbm9iaW9wCgpgYGB7ciBjbGluaWNhbF9zdmF9CmhzX2NsaW5pY2FsX25iIDwtIG5vcm1hbGl6ZV9leHB0KGhzX2NsaW5pY2FsLCBmaWx0ZXIgPSBUUlVFLCBiYXRjaCA9ICJzdmFzZXEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB0cmFuc2Zvcm0gPSAibG9nMiIsIGNvbnZlcnQgPSAiY3BtIikKY2xpbmljYWxfYmF0Y2hfcGNhIDwtIHBsb3RfcGNhKGhzX2NsaW5pY2FsX25iLCBwbG90X2xhYmVscyA9IEZBTFNFLCBjaXMgPSBOVUxMLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc2l6ZV9jb2x1bW4gPSAidmlzaXRudW1iZXIiLCBwbG90X3RpdGxlID0gIlBDQSAtIGNsaW5pY2FsIHNhbXBsZXMiKQpjbGluaWNhbF9iYXRjaF9wY2EkcGxvdAoKaHNfY2xpbmljYWxfbm9iaW9wX25iIDwtIHNtKG5vcm1hbGl6ZV9leHB0KGhzX2NsaW5pY2FsX25vYmlvcCwgZmlsdGVyID0gVFJVRSwgYmF0Y2ggPSAic3Zhc2VxIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHRyYW5zZm9ybSA9ICJsb2cyIiwgY29udmVydCA9ICJjcG0iKSkKY2xpbmljYWxfbm9iaW9wX2JhdGNoX3BjYSA8LSBwbG90X3BjYShoc19jbGluaWNhbF9ub2Jpb3BfbmIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgcGxvdF90aXRsZSA9ICJQQ0EgLSBjbGluaWNhbCBzYW1wbGVzIHdpdGhvdXQgYmlvcHNpZXMiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHBsb3RfbGFiZWxzID0gRkFMU0UpCnBwKGZpbGUgPSAiaW1hZ2VzL2NsaW5pY2FsX2JhdGNoLnBuZyIsIGltYWdlID0gY2xpbmljYWxfbm9iaW9wX2JhdGNoX3BjYSRwbG90KQp0ZXN0IDwtIHBsb3RfcGNhKGhzX2NsaW5pY2FsX25vYmlvcF9uYiwgc2l6ZV9jb2x1bW4gPSAidmlzaXRudW1iZXIiLAogICAgICAgICAgICAgICAgIHBsb3RfdGl0bGUgPSAiUENBIC0gY2xpbmljYWwgc2FtcGxlcyB3aXRob3V0IGJpb3BzaWVzIiwKICAgICAgICAgICAgICAgICBwbG90X2xhYmVscyA9IEZBTFNFKQp0ZXN0JHBsb3QKCmNsaW5pY2FsX25vYmlvcF9iYXRjaF90c25lIDwtIHBsb3RfdHNuZShoc19jbGluaWNhbF9ub2Jpb3BfbmIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBwbG90X3RpdGxlID0gInRTTkUgLSBjbGluaWNhbCBzYW1wbGVzIHdpdGhvdXQgYmlvcHNpZXMiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgcGxvdF9sYWJlbHMgPSBGQUxTRSkKY2xpbmljYWxfbm9iaW9wX2JhdGNoX3RzbmUkcGxvdApgYGAKCiMjIyMgTG9vayBhdCByZW1haW5pbmcgdmFyaWFuY2Ugd2l0aCB2YXJpYW5jZVBhcnRpdGlvbgoKYGBge3IgdmFyaWFuY2VfcGFydGl0aW9ufQp0ZXN0IDwtIHNpbXBsZV92YXJwYXJ0KGhzX2NsaW5pY2FsX25vYmlvcCkKdGVzdCRwYXJ0aXRpb25fcGxvdApgYGAKCiMjIFBlcmZvcm0gREUgb2YgdGhlIGNsaW5pY2FsIHNhbXBsZXMgY3VyZSB2cy4gZmFpbAoKYGBge3IgY2xpbmljYWxfZGUsIGZpZy5zaG93PSJoaWRlIn0KaW5kaXZpZHVhbF9jZWxsdHlwZXMgPC0gc3Vic2V0X2V4cHQoaHNfY2xpbmljYWxfbm9iaW9wLCBzdWJzZXQ9ImNvbmRpdGlvbiE9J2xvc3QnIikKaHNfY2xpbmljX2RlIDwtIHNtKGFsbF9wYWlyd2lzZShpbmRpdmlkdWFsX2NlbGx0eXBlcywgbW9kZWxfYmF0Y2ggPSAic3Zhc2VxIiwgZmlsdGVyID0gVFJVRSkpCgpoc19jbGluaWNfdGFibGUgPC0gc20oY29tYmluZV9kZV90YWJsZXMoCiAgICBoc19jbGluaWNfZGUsCiAgICBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL2luZGl2aWR1YWxfY2VsbHR5cGVzX3RhYmxlLXZ7dmVyfS54bHN4IikpKQoKaHNfY2xpbmljX3NpZyA8LSBzbShleHRyYWN0X3NpZ25pZmljYW50X2dlbmVzKAogICAgaHNfY2xpbmljX3RhYmxlLAogICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9pbmRpdmlkdWFsX2NlbGx0eXBlc19zaWctdnt2ZXJ9Lnhsc3giKSkpCgpoc19jbGluaWNfc2lnW1sic3VtbWFyeV9kZiJdXQpgYGAKCmBgYHtyIGRlX2hlYXRtYXB9CmhzX2NsaW5pY19kZVtbImNvbXBhcmlzb24iXV1bWyJoZWF0Il1dCmBgYAoKIyMjIFBlcmZvcm0gTFJUIHdpdGggdGhlIGNsaW5pY2FsIHNhbXBsZXMKCkkgYW0gbm90IHN1cmUgaWYgd2UgaGF2ZSBlbm91Z2ggc2FtcGxlcyBhY3Jvc3MgdGhlIHRocmVlIHZpc2l0IHRvCmNvbXBsZXRlbHkgd29yayBhcyB3ZWxsIGFzIHdlIHdvdWxkIGxpa2UsIGJ1dCB0aGVyZSBpcyBvbmx5IDEgd2F5IHRvCmZpbmQgb3V0ISAgTm93IHRoYXQgSSB0aGluayBhYm91dCBpdCwgb25lIHRoaW5nIHdoaWNoIG1pZ2h0IGJlIGF3ZXNvbWUKaXMgdG8gdXNlIGNlbGwgdHlwZSBhcyBhbiBpbnRlcmFjdGluZyBmYWN0b3IuLi4KCiMjIyMgV2l0aCBiaW9wc3kgc2FtcGxlcwoKSSBmaWd1cmUgdGhpcyBtaWdodCBiZSBhIHBsYWNlIHdoZXJlIHRoZSBiaW9wc3kgc2FtcGxlcyBtaWdodCBwcm92ZSB1c2VmdWwuCgpgYGB7ciBscnRfdGVzdH0KY2xpbmljYWxfbm9sb3N0IDwtIHN1YnNldF9leHB0KGhzX2NsaW5pY2FsLCBzdWJzZXQ9ImNvbmRpdGlvbiE9J2xvc3QnIikKbHJ0X3Zpc2l0X2NsaW5pY2FsX3Rlc3QgPC0gZGVzZXFfbHJ0KGNsaW5pY2FsX25vbG9zdCwgdHJhbnNmb3JtID0gInZzdCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBpbnRlcmFjdG9yX2NvbHVtbiA9ICJ2aXNpdG51bWJlciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBpbnRlcmVzdF9jb2x1bW4gPSAiY2xpbmljYWxvdXRjb21lIikKbHJ0X3Zpc2l0X2NsaW5pY2FsX3Rlc3RbWyJmYXZvcml0ZV9nZW5lcyJdXQoKbHJ0X2NlbGx0eXBlX2NsaW5pY2FsX3Rlc3QgPC0gZGVzZXFfbHJ0KGNsaW5pY2FsX25vbG9zdCwgdHJhbnNmb3JtID0gInZzdCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBpbnRlcmFjdG9yX2NvbHVtbiA9ICJ0eXBlb2ZjZWxscyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBpbnRlcmVzdF9jb2x1bW4gPSAiY2xpbmljYWxvdXRjb21lIikKbHJ0X2NlbGx0eXBlX2NsaW5pY2FsX3Rlc3RbWyJmYXZvcml0ZV9nZW5lcyJdXQpgYGAKCiMjIyBMb29rIGF0IG9ubHkgdGhlIGRpZmZlcmVudGlhbCBnZW5lcwoKQSBnb29kIHN1Z2dlc3Rpb24gZnJvbSBUaGVyZXNhIHdhcyB0byBleGFtaW5lIG9ubHkgdGhlIG1vc3QgdmFyaWFudApnZW5lcyBmcm9tIGZhaWx1cmUgdnMuIGN1cmUgYW5kIHNlZSBob3cgdGhleSBjaGFuZ2UgdGhlIGNsdXN0ZXJpbmcvZXRjCnJlc3VsdHMuICBUaGlzIGlzIG15IGF0dGVtcHQgdG8gYWRkcmVzcyB0aGlzIHF1ZXJ5LgoKYGBge3Igc21hbGxfZXhwbG9yZX0KaHNfY2xpbmljX3RvcG4gPC0gc20oZXh0cmFjdF9zaWduaWZpY2FudF9nZW5lcyhoc19jbGluaWNfdGFibGUsIG4gPSAxMDApKQp0YWJsZSA8LSAiZmFpbHVyZV92c19jdXJlIgp3YW50ZWQgPC0gcmJpbmQoaHNfY2xpbmljX3RvcG5bWyJkZXNlcSJdXVtbInVwcyJdXVtbdGFibGVdXSwKICAgICAgICAgICAgICAgIGhzX2NsaW5pY190b3BuW1siZGVzZXEiXV1bWyJkb3ducyJdXVtbdGFibGVdXSkKCnNtYWxsX2V4cHQgPC0gZXhjbHVkZV9nZW5lc19leHB0KGhzX2NsaW5pY2FsX25vYmlvcCwgaWRzID0gcm93bmFtZXMod2FudGVkKSwgbWV0aG9kID0gImtlZXAiKQpzbWFsbF9ub3JtIDwtIHNtKG5vcm1hbGl6ZV9leHB0KHNtYWxsX2V4cHQsIHRyYW5zZm9ybSA9ICJsb2cyIiwgY29udmVydCA9ICJjcG0iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG5vcm0gPSAicXVhbnQiLCBmaWx0ZXIgPSBUUlVFKSkKcGxvdF9wY2Eoc21hbGxfbm9ybSkkcGxvdAoKc21hbGxfbmIgPC0gbm9ybWFsaXplX2V4cHQoc21hbGxfZXhwdCwgdHJhbnNmb3JtID0gImxvZzIiLCBjb252ZXJ0ID0gImNwbSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgIGJhdGNoID0gInN2YXNlcSIsIG5vcm0gPSAicXVhbnQiLCBmaWx0ZXIgPSBUUlVFKQpwbG90X3BjYShzbWFsbF9uYikkcGxvdApgYGAKCmBgYHtyIGNsaW5pY2FsX3Bsb3R9CiMjIERFU2VxMiBNQSBwbG90IG9mIGZhaWx1cmUgLyBjdXJlCmhzX2NsaW5pY190YWJsZVtbInBsb3RzIl1dW1siZmFpbHVyZV92c19jdXJlIl1dW1siZGVzZXFfbWFfcGxvdHMiXV0kcGxvdAoKIyMgREVTZXEyIFZvbGNhbm8gcGxvdCBvZiBmYWlsdXJlIC8gY3VyZQpoc19jbGluaWNfdGFibGVbWyJwbG90cyJdXVtbImZhaWx1cmVfdnNfY3VyZSJdXVtbImRlc2VxX3ZvbF9wbG90cyJdXSRwbG90CmBgYAoKIyMjIGc6UHJvZmlsZXIgcmVzdWx0cyB1c2luZyB0aGUgc2lnbmlmaWNhbnQgdXAgYW5kIGRvd24gZ2VuZXMKCmBgYHtyIHBlcmZvcm1fZ3Byb2ZpbGVyfQp1cHMgPC0gaHNfY2xpbmljX3NpZ1tbImRlc2VxIl1dW1sidXBzIl1dW1sxXV0KZG93bnMgPC0gaHNfY2xpbmljX3NpZ1tbImRlc2VxIl1dW1siZG93bnMiXV1bWzFdXQoKaHNfY2xpbmljX2dwcm9maWxlcl91cHMgPC0gc2ltcGxlX2dwcm9maWxlcih1cHMpCmhzX2NsaW5pY19ncHJvZmlsZXJfdXBzW1sicHZhbHVlX3Bsb3RzIl1dW1siYnBwX3Bsb3Rfb3ZlciJdXQpoc19jbGluaWNfZ3Byb2ZpbGVyX3Vwc1tbInB2YWx1ZV9wbG90cyJdXVtbIm1mcF9wbG90X292ZXIiXV0KaHNfY2xpbmljX2dwcm9maWxlcl91cHNbWyJwdmFsdWVfcGxvdHMiXV1bWyJyZWFjdG9tZV9wbG90X292ZXIiXV0KCiMjaHNfdHJ5MiA8LSBzaW1wbGVfZ3Byb2ZpbGVyMih1cHMpCgpoc19jbGluaWNfZ3Byb2ZpbGVyX2Rvd25zIDwtIHNpbXBsZV9ncHJvZmlsZXIoZG93bnMpCmhzX2NsaW5pY19ncHJvZmlsZXJfZG93bnNbWyJwdmFsdWVfcGxvdHMiXV1bWyJicHBfcGxvdF9vdmVyIl1dCmhzX2NsaW5pY19ncHJvZmlsZXJfZG93bnNbWyJwdmFsdWVfcGxvdHMiXV1bWyJtZnBfcGxvdF9vdmVyIl1dCmhzX2NsaW5pY19ncHJvZmlsZXJfZG93bnNbWyJwdmFsdWVfcGxvdHMiXV1bWyJyZWFjdG9tZV9wbG90X292ZXIiXV0KYGBgCgojIyBQZXJmb3JtIEdTVkEgb24gdGhlIGNsaW5pY2FsIHNhbXBsZXMKCmBgYHtyIGdzdmEsIGZpZy5zaG93ID0gImhpZGUifQpoc19jZWxsdHlwZV9nc3ZhX2MyIDwtIHNtKHNpbXBsZV9nc3ZhKGluZGl2aWR1YWxfY2VsbHR5cGVzKSkKaHNfY2VsbHR5cGVfZ3N2YV9jMl9zaWcgPC0gc20oZ2V0X3NpZ19nc3ZhX2NhdGVnb3JpZXMoCiAgICBoc19jZWxsdHlwZV9nc3ZhX2MyLAogICAgZXhjZWwgPSAiZXhjZWwvaW5kaXZpZHVhbF9jZWxsdHlwZXNfZ3N2YV9jMi54bHN4IikpCgpicm9hZF9jNyA8LSBHU0VBQmFzZTo6Z2V0R210KCJyZWZlcmVuY2UvbXNpZ2RiL2M3LmFsbC52Ny4yLmVudHJlei5nbXQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbGxlY3Rpb25UeXBlID0gR1NFQUJhc2U6OkJyb2FkQ29sbGVjdGlvbihjYXRlZ29yeSA9ICJjNyIpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIGdlbmVJZFR5cGUgPSBHU0VBQmFzZTo6RW50cmV6SWRlbnRpZmllcigpKQpoc19jZWxsdHlwZV9nc3ZhX2M3IDwtIHNtKHNpbXBsZV9nc3ZhKGluZGl2aWR1YWxfY2VsbHR5cGVzLCBzaWduYXR1cmVzID0gYnJvYWRfYzcsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbXNpZ194bWwgPSAicmVmZXJlbmNlL21zaWdkYl92Ny4yLnhtbCIsIGNvcmVzID0gMTApKQpoc19jZWxsdHlwZV9nc3ZhX2M3X3NpZyA8LSBzbShnZXRfc2lnX2dzdmFfY2F0ZWdvcmllcygKICAgIGhzX2NlbGx0eXBlX2dzdmFfYzcsCiAgICBleGNlbCA9ICJleGNlbC9pbmRpdmlkdWFsX2NlbGx0eXBlc19nc3ZhX2M3Lnhsc3giKSkKYGBgCgojIyMgUHJpbnQgc29tZSBwbG90cyBvZiB0aGUgR1NWQSBvdXRwdXRzCgpgYGB7ciBnc3ZhX3Bsb3RzfQojIyBUaGUgcmF3IGhlYXRtYXAgb2YgdGhlIEMyIHZhbHVlcwpoc19jZWxsdHlwZV9nc3ZhX2MyX3NpZ1tbInJhd19wbG90Il1dCiMjIFRoZSAnc2lnbmlmaWNhbmNlJyBzY29yZXMgb2YgdGhlIEMyIHZhbHVlcwpoc19jZWxsdHlwZV9nc3ZhX2MyX3NpZ1tbInNjb3JlX3Bsb3QiXV0KIyMgVGhlIHN1YnNldCBvZiBzY29yZXMgZm9yIGNhdGVnb3JpZXMgZGVlbWVkIHNpZ25pZmljYW50bHkgZGlmZmVyZW50Lgpoc19jZWxsdHlwZV9nc3ZhX2MyX3NpZ1tbInN1YnNldF9wbG90Il1dCgojIyBUaGUgcmF3IGhlYXRtYXAgb2YgdGhlIEM3IHZhbHVlcwpoc19jZWxsdHlwZV9nc3ZhX2M3X3NpZ1tbInJhd19wbG90Il1dCiMjIFRoZSAnc2lnbmlmaWNhbmNlJyBzY29yZXMgb2YgdGhlIEM3IHZhbHVlcwpoc19jZWxsdHlwZV9nc3ZhX2M3X3NpZ1tbInNjb3JlX3Bsb3QiXV0KIyMgVGhlIHN1YnNldCBvZiBzY29yZXMgZm9yIGNhdGVnb3JpZXMgZGVlbWVkIHNpZ25pZmljYW50bHkgZGlmZmVyZW50Lgpoc19jZWxsdHlwZV9nc3ZhX2M3X3NpZ1tbInN1YnNldF9wbG90Il1dCmBgYAoKIyBJbmRpdmlkdWFsIENlbGwgdHlwZXMKClRoZSBmb2xsb3dpbmcgYmxvY2tzIHNwbGl0IHRoZSBzYW1wbGVzIGludG8gYSBmZXcgZ3JvdXBzIGJ5IHNhbXBsZSB0eXBlIGFuZCBsb29rCmF0IHRoZSBkaXN0cmlidXRpb25zIGJldHdlZW4gdGhlbS4KCiMjIEltcGxlbWVudGF0aW9uIGRldGFpbHMKCkdldCB0b3AvYm90dG9tIG4gZ2VuZXMgZm9yIGVhY2ggY2VsbCB0eXBlLCB1c2luZyBjbGluaWNhbCBvdXRjb21lIGFzIHRoZSBmYWN0b3Igb2YgaW50ZXJlc3QuCkZvciB0aGUgbW9tZW50LCB1c2Ugc3ZhIGZvciB0aGUgREUgYW5hbHlzaXMuClByb3ZpZGUgY3BtcyBmb3IgdGhlIHRvcC9ib3R0b20gbiBnZW5lcy4KClN0YXJ0IHdpdGggdG9wL2JvdHRvbSAyMDAuClBlcmZvcm0gZGVmYXVsdCBsb2dGQyBhbmQgcC12YWx1ZSBhcyB3ZWxsLgoKIyMjIFNoYXJlZCBjb250cmFzdHMKCkhlcmUgaXMgdGhlIGNvbnRyYXN0IHdlIHdpbGwgdXNlIHRocm91Z2hwdXQsIEkgYW0gbGVhdmluZyBvcGVuIHRoZSBvcHRpb24gdG8gYWRkIG1vcmUuCgpgYGB7ciBrZWVwZXJzfQprZWVwZXJzIDwtIGxpc3QoCiAgImZhaWxfdnNfY3VyZSIgPSBjKCJmYWlsdXJlIiwgImN1cmUiKSkKYGBgCgojIyBNb25vY3l0ZXMKCiMjIyBFdmFsdWF0ZSBNb25vY3l0ZSBzYW1wbGVzCgpgYGB7ciBtb25vY3l0ZXN9Cm1vbm8gPC0gc3Vic2V0X2V4cHQoaHNfdmFsaWQsIHN1YnNldCA9ICJ0eXBlb2ZjZWxscz09J21vbm9jeXRlcyciKSAlPiUKICBzZXRfZXhwdF9jb25kaXRpb25zKGZhY3QgPSAiY2xpbmljYWxvdXRjb21lIikgJT4lCiAgc2V0X2V4cHRfYmF0Y2hlcyhmYWN0ID0gImRvbm9yIikgJT4lCiAgc2V0X2V4cHRfY29sb3JzKGNvbG9ycyA9IGNob3Nlbl9jb2xvcnMpCiMjIEZJWE1FIHNldF9leHB0X2NvbG9ycyBzaG91bGQgc3BlYWsgdXAgaWYgdGhlcmUgYXJlIG1pc21hdGNoZXMgaGVyZSEhIQoKc2F2ZV9yZXN1bHQgPC0gc2F2ZShtb25vLCBmaWxlID0gInJkYS9tb25vY3l0ZV9leHB0LnJkYSIpCm1vbm9fbm9ybSA8LSBub3JtYWxpemVfZXhwdChtb25vLCBjb252ZXJ0ID0gImNwbSIsIGZpbHRlciA9IFRSVUUsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICB0cmFuc2Zvcm0gPSAibG9nMiIsIG5vcm0gPSAicXVhbnQiKQpwbHQgPC0gcGxvdF9wY2EobW9ub19ub3JtLCBwbG90X2xhYmVscyA9IEZBTFNFKSRwbG90CnBwKGZpbGUgPSBnbHVlKCJpbWFnZXMvbW9ub19wY2Ffbm9ybWFsaXplZC12e3Zlcn0ucGRmIiksIGltYWdlID0gcGx0KQoKbW9ub19uYiA8LSBub3JtYWxpemVfZXhwdChtb25vLCBjb252ZXJ0ID0gImNwbSIsIGZpbHRlciA9IFRSVUUsCiAgICAgICAgICAgICAgICAgICAgICAgICAgdHJhbnNmb3JtID0gImxvZzIiLCBiYXRjaCA9ICJzdmFzZXEiKQpwbHQgPC0gcGxvdF9wY2EobW9ub19uYiwgcGxvdF9sYWJlbHMgPSBGQUxTRSkkcGxvdApwcChmaWxlID0gZ2x1ZSgiaW1hZ2VzL21vbm9fcGNhX25vcm1hbGl6ZWRfYmF0Y2gtdnt2ZXJ9LnBkZiIpLCBpbWFnZSA9IHBsdCkKYGBgCgojIyMgREUgb2YgTW9ub2N5dGUgc2FtcGxlcwoKIyMjIyBXaXRob3V0IHN2YQoKYGBge3IgZGVfbW9ub2N5dGUsIGZpZy5zaG93ID0gImhpZGUifQptb25vX2RlIDwtIHNtKGFsbF9wYWlyd2lzZShtb25vLCBtb2RlbF9iYXRjaCA9IEZBTFNFLCBmaWx0ZXIgPSBUUlVFKSkKbW9ub190YWJsZXMgPC0gc20oY29tYmluZV9kZV90YWJsZXMoCiAgICBtb25vX2RlLCBrZWVwZXJzID0ga2VlcGVycywKICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvbW9ub2N5dGVfY2xpbmljYWxfYWxsX3RhYmxlcy12e3Zlcn0ueGxzeCIpKSkKd3JpdHRlbiA8LSB3cml0ZV94bHN4KGRhdGEgPSBtb25vX3RhYmxlc1tbImRhdGEiXV1bWzFdXSwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvbW9ub2N5dGVfY2xpbmljYWxfdGFibGUtdnt2ZXJ9Lnhsc3giKSkKbW9ub19zaWcgPC0gc20oZXh0cmFjdF9zaWduaWZpY2FudF9nZW5lcyhtb25vX3RhYmxlcywgYWNjb3JkaW5nX3RvID0gImRlc2VxIikpCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhID0gbW9ub19zaWdbWyJkZXNlcSJdXVtbInVwcyJdXVtbMV1dLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9tb25vY3l0ZV9jbGluaWNhbF9zaWd1cC12e3Zlcn0ueGxzeCIpKQp3cml0dGVuIDwtIHdyaXRlX3hsc3goZGF0YSA9IG1vbm9fc2lnW1siZGVzZXEiXV1bWyJkb3ducyJdXVtbMV1dLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9tb25vY3l0ZV9jbGluaWNhbF9zaWdkb3duLXZ7dmVyfS54bHN4IikpCgptb25vX3BjdF9zaWcgPC0gc20oZXh0cmFjdF9zaWduaWZpY2FudF9nZW5lcyhtb25vX3RhYmxlcywgbiA9IDIwMCwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbGZjID0gTlVMTCwgcCA9IE5VTEwsIGFjY29yZGluZ190byA9ICJkZXNlcSIpKQp3cml0dGVuIDwtIHdyaXRlX3hsc3goZGF0YSA9IG1vbm9fcGN0X3NpZ1tbImRlc2VxIl1dW1sidXBzIl1dW1sxXV0sCiAgICAgICAgICAgICAgICAgICAgICBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL21vbm9jeXRlX2NsaW5pY2FsX3NpZ3VwX3BjdC12e3Zlcn0ueGxzeCIpKQp3cml0dGVuIDwtIHdyaXRlX3hsc3goZGF0YSA9IG1vbm9fcGN0X3NpZ1tbImRlc2VxIl1dW1siZG93bnMiXV1bWzFdXSwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvbW9ub2N5dGVfY2xpbmljYWxfc2lnZG93bl9wY3Qtdnt2ZXJ9Lnhsc3giKSkKbW9ub19zaWckc3VtbWFyeV9kZgoKIyMgUHJpbnQgb3V0IGEgdGFibGUgb2YgdGhlIGNwbSB2YWx1ZXMgZm9yIG90aGVyIGV4cGxvcmF0aW9ucy4KbW9ub19jcG0gPC0gc20obm9ybWFsaXplX2V4cHQobW9ubywgY29udmVydCA9ICJjcG0iKSkKd3JpdHRlbiA8LSB3cml0ZV94bHN4KGRhdGEgPSBleHBycyhtb25vX2NwbSksCiAgICAgICAgICAgICAgICAgICAgICBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL21vbm9jeXRlX2NwbV9iZWZvcmVfYmF0Y2gtdnt2ZXJ9Lnhsc3giKSkKbW9ub19iY3BtIDwtIHNtKG5vcm1hbGl6ZV9leHB0KG1vbm8sIGZpbHRlciA9IFRSVUUsIGNvbnZlcnQgPSAiY3BtIiwgYmF0Y2ggPSAic3Zhc2VxIikpCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhID0gZXhwcnMobW9ub19iY3BtKSwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvbW9ub2N5dGVfY3BtX2FmdGVyX2JhdGNoLXZ7dmVyfS54bHN4IikpCmBgYAoKIyMjIyBXaXRoIHN2YQoKYGBge3IgZGVfbW9ub19zdmEsIGZpZy5zaG93ID0gImhpZGUifQptb25vX2RlX3N2YSA8LSBzbShhbGxfcGFpcndpc2UobW9ubywgbW9kZWxfYmF0Y2ggPSAic3Zhc2VxIiwgZmlsdGVyID0gVFJVRSkpCm1vbm9fdGFibGVzX3N2YSA8LSBzbShjb21iaW5lX2RlX3RhYmxlcygKICAgIG1vbm9fZGVfc3ZhLCBrZWVwZXJzID0ga2VlcGVycywKICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvbW9ub2N5dGVfY2xpbmljYWxfYWxsX3RhYmxlc19zdmEtdnt2ZXJ9Lnhsc3giKSkpCm1vbm9fc2lnX3N2YSA8LSBzbShleHRyYWN0X3NpZ25pZmljYW50X2dlbmVzKAogICAgbW9ub190YWJsZXNfc3ZhLAogICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9tb25vY3l0ZV9jbGluaWNhbF9zaWdfdGFibGVzX3N2YS12e3Zlcn0ueGxzeCIpLAogICAgYWNjb3JkaW5nX3RvID0gImRlc2VxIikpCmBgYAoKIyMjIyBNb25vY3l0ZSBERSBwbG90cwoKRmlyc3QgcHJpbnQgb3V0IHRoZSBERSBwbG90cyB3aXRob3V0IGFuZCB0aGVuIHdpdGggc3ZhIGVzdGltYXRlcy4KCmBgYHtyIG1vbm9fZGVfcGxvdHN9CiMjIERFU2VxMiBNQSBwbG90IG9mIGZhaWx1cmUgLyBjdXJlCm1vbm9fdGFibGVzW1sicGxvdHMiXV1bWyJmYWlsX3ZzX2N1cmUiXV1bWyJkZXNlcV9tYV9wbG90cyJdXSRwbG90CgojIyBERVNlcTIgVm9sY2FubyBwbG90IG9mIGZhaWx1cmUgLyBjdXJlCm1vbm9fdGFibGVzW1sicGxvdHMiXV1bWyJmYWlsX3ZzX2N1cmUiXV1bWyJkZXNlcV92b2xfcGxvdHMiXV0kcGxvdAoKIyMgREVTZXEyIE1BIHBsb3Qgb2YgZmFpbHVyZSAvIGN1cmUgd2l0aCBzdmFzZXEKbW9ub190YWJsZXNfc3ZhW1sicGxvdHMiXV1bWyJmYWlsX3ZzX2N1cmUiXV1bWyJkZXNlcV9tYV9wbG90cyJdXSRwbG90CgojIyBERVNlcTIgVm9sY2FubyBwbG90IG9mIGZhaWx1cmUgLyBjdXJlIHdpdGggc3Zhc2VxCm1vbm9fdGFibGVzX3N2YVtbInBsb3RzIl1dW1siZmFpbF92c19jdXJlIl1dW1siZGVzZXFfdm9sX3Bsb3RzIl1dJHBsb3QKYGBgCgojIyMjIE1vbm9jeXRlIG9udG9sb2d5IHNlYXJjaAoKYGBge3IgbW9ub19ncHJvZmlsZXJ9CnVwcyA8LSBtb25vX3NpZ1tbImRlc2VxIl1dW1sidXBzIl1dW1siZmFpbF92c19jdXJlIl1dCmRvd25zIDwtIG1vbm9fc2lnW1siZGVzZXEiXV1bWyJkb3ducyJdXVtbImZhaWxfdnNfY3VyZSJdXQoKbW9ub191cF9ncCA8LSBzaW1wbGVfZ3Byb2ZpbGVyKHNpZ19nZW5lcyA9IHVwcykKbW9ub191cF9ncFtbInB2YWx1ZV9wbG90cyJdXVtbImJwcF9wbG90X292ZXIiXV0KbW9ub191cF9ncFtbInB2YWx1ZV9wbG90cyJdXVtbIm1mcF9wbG90X292ZXIiXV0KbW9ub191cF9ncFtbInB2YWx1ZV9wbG90cyJdXVtbInJlYWN0b21lX3Bsb3Rfb3ZlciJdXQoKbW9ub19kb3duX2dwIDwtIHNpbXBsZV9ncHJvZmlsZXIoc2lnX2dlbmVzID0gZG93bnMpCm1vbm9fZG93bl9ncFtbInB2YWx1ZV9wbG90cyJdXVtbImJwcF9wbG90X292ZXIiXV0KbW9ub19kb3duX2dwW1sicHZhbHVlX3Bsb3RzIl1dW1sibWZwX3Bsb3Rfb3ZlciJdXQptb25vX2Rvd25fZ3BbWyJwdmFsdWVfcGxvdHMiXV1bWyJyZWFjdG9tZV9wbG90X292ZXIiXV0KYGBgCgojIyMjIE1vbm9jeXRlIE1TaWdEQiBxdWVyeQoKYGBge3IgbXNpZ19tb25vX2dvc2VxLCBmaWcuc2hvdyA9ICJoaWRlIn0KYnJvYWRfYzcgPC0gR1NFQUJhc2U6OmdldEdtdCgicmVmZXJlbmNlL21zaWdkYi9jNy5hbGwudjcuMi5lbnRyZXouZ210IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb2xsZWN0aW9uVHlwZSA9IEdTRUFCYXNlOjpCcm9hZENvbGxlY3Rpb24oY2F0ZWdvcnkgPSAiYzciKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICBnZW5lSWRUeXBlID0gR1NFQUJhc2U6OkVudHJleklkZW50aWZpZXIoKSkKCm1vbm9fdXBfZ29zZXFfbXNpZyA8LSBnb3NlcV9tc2lnZGIoc2lnX2dlbmVzID0gdXBzLCBzaWduYXR1cmVzID0gYnJvYWRfYzcsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc2lnbmF0dXJlX2NhdGVnb3J5ID0gImM3IiwgbGVuZ3RoX2RiID0gaHNfbGVuZ3RoKQoKbW9ub19kb3duX2dvc2VxX21zaWcgPC0gZ29zZXFfbXNpZ2RiKHNpZ19nZW5lcyA9IGRvd25zLCBzaWduYXR1cmVzID0gYnJvYWRfYzcsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzaWduYXR1cmVfY2F0ZWdvcnkgPSAiYzciLCBsZW5ndGhfZGIgPSBoc19sZW5ndGgpCmBgYAoKIyMjIyBQbG90IG9mIHNpbWlsYXIgZXhwZXJpbWVudHMKCmBgYHtyIG1zaWdfcGxvdHN9CiMjIE1vbm9jeXRlIGdlbmVzIHdpdGggaW5jcmVhc2VkIGV4cHJlc3Npb24gaW4gdGhlIGZhaWxlZCBzYW1wbGVzCiMjIHNoYXJlIGdlbmVzIHdpdGggdGhlIGZvbGxvd2luZyBleHBlcmltZW50cwptb25vX3VwX2dvc2VxX21zaWdbWyJwdmFsdWVfcGxvdHMiXV1bWyJtZnBfcGxvdF9vdmVyIl1dCgojIyBNb25vY3l0ZSBnZW5lcyB3aXRoIGluY3JlYXNlZCBleHByZXNzaW9uIGluIHRoZSBjdXJlZCBzYW1wbGVzCiMjIHNoYXJlIGdlbmVzIHdpdGggdGhlIGZvbGxvd2luZyBleHBlcmltZW50cwptb25vX2Rvd25fZ29zZXFfbXNpZ1tbInB2YWx1ZV9wbG90cyJdXVtbIm1mcF9wbG90X292ZXIiXV0KYGBgCgojIyMgRXZhbHVhdGUgTmV1dHJvcGhpbCBzYW1wbGVzCgpgYGB7ciBuZXV0cm9waGlsc30KbmV1dCA8LSBzdWJzZXRfZXhwdChoc192YWxpZCwgc3Vic2V0ID0gInR5cGVvZmNlbGxzPT0nbmV1dHJvcGhpbHMnIikgJT4lCiAgc2V0X2V4cHRfY29uZGl0aW9ucyhmYWN0ID0gImNsaW5pY2Fsb3V0Y29tZSIpICU+JQogIHNldF9leHB0X2JhdGNoZXMoZmFjdCA9ICJkb25vciIpICU+JQogIHNldF9leHB0X2NvbG9ycyhjb2xvcnMgPSBjaG9zZW5fY29sb3JzKQoKc2F2ZV9yZXN1bHQgPC0gc2F2ZShuZXV0LCBmaWxlID0gInJkYS9uZXV0cm9waGlsX2V4cHQucmRhIikKbmV1dF9ub3JtIDwtIHNtKG5vcm1hbGl6ZV9leHB0KG5ldXQsIGNvbnZlcnQgPSAiY3BtIiwgZmlsdGVyID0gVFJVRSwgdHJhbnNmb3JtID0gImxvZzIiKSkKcGx0IDwtIHBsb3RfcGNhKG5ldXRfbm9ybSwgcGxvdF9sYWJlbHMgPSBGQUxTRSkkcGxvdApwcChmaWxlID0gZ2x1ZSgiaW1hZ2VzL25ldXRfcGNhX25vcm1hbGl6ZWQtdnt2ZXJ9LnBkZiIpLCBpbWFnZSA9IHBsdCkKCm5ldXRfbmIgPC0gc20obm9ybWFsaXplX2V4cHQobmV1dCwgY29udmVydCA9ICJjcG0iLCBmaWx0ZXIgPSBUUlVFLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIHRyYW5zZm9ybSA9ICJsb2cyIiwgYmF0Y2ggPSAic3Zhc2VxIikpCnBsdCA8LSBwbG90X3BjYShuZXV0X25iLCBwbG90X2xhYmVscyA9IEZBTFNFKSRwbG90CnBwKGZpbGUgPSBnbHVlKCJpbWFnZXMvbmV1dF9wY2Ffbm9ybWFsaXplZF9zdmFzZXEtdnt2ZXJ9LnBkZiIpLCBpbWFnZSA9IHBsdCkKYGBgCgojIyMgREUgb2YgTmV0cm9waGlsIHNhbXBsZXMKCiMjIyMgV2l0aG91dCBzdmEKCmBgYHtyIG5ldXRyb3BoaWxfZGUsIGZpZy5zaG93ID0gImhpZGUifQpuZXV0X2RlIDwtIHNtKGFsbF9wYWlyd2lzZShuZXV0LCBtb2RlbF9iYXRjaCA9IEZBTFNFLCBmaWx0ZXIgPSBUUlVFKSkKbmV1dF90YWJsZXMgPC0gc20oY29tYmluZV9kZV90YWJsZXMoCiAgICBuZXV0X2RlLCBrZWVwZXJzID0ga2VlcGVycywKICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvbmV1dHJvcGhpbF9jbGluaWNhbF9hbGxfdGFibGVzLXZ7dmVyfS54bHN4IikpKQp3cml0dGVuIDwtIHdyaXRlX3hsc3goZGF0YSA9IG5ldXRfdGFibGVzW1siZGF0YSJdXVtbMV1dLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9uZXV0cm9waGlsX2NsaW5pY2FsX3RhYmxlLXZ7dmVyfS54bHN4IikpCm5ldXRfc2lnIDwtIHNtKGV4dHJhY3Rfc2lnbmlmaWNhbnRfZ2VuZXMobmV1dF90YWJsZXMsIGFjY29yZGluZ190byA9ICJkZXNlcSIpKQp3cml0dGVuIDwtIHdyaXRlX3hsc3goZGF0YSA9IG5ldXRfc2lnW1siZGVzZXEiXV1bWyJ1cHMiXV1bWzFdXSwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvbmV1dHJvcGhpbF9jbGluaWNhbF9zaWd1cC12e3Zlcn0ueGxzeCIpKQp3cml0dGVuIDwtIHdyaXRlX3hsc3goZGF0YSA9IG5ldXRfc2lnW1siZGVzZXEiXV1bWyJkb3ducyJdXVtbMV1dLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9uZXV0cm9waGlsX2NsaW5pY2FsX3NpZ2Rvd24tdnt2ZXJ9Lnhsc3giKSkKCm5ldXRfcGN0X3NpZyA8LSBzbShleHRyYWN0X3NpZ25pZmljYW50X2dlbmVzKG5ldXRfdGFibGVzLCBuID0gMjAwLCBsZmMgPSBOVUxMLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBwID0gTlVMTCwgYWNjb3JkaW5nX3RvID0gImRlc2VxIikpCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhID0gbmV1dF9zaWdbWyJkZXNlcSJdXVtbInVwcyJdXVtbMV1dLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9uZXV0cm9waGlsX2NsaW5pY2FsX3NpZ3VwX3BjdC12e3Zlcn0ueGxzeCIpKQp3cml0dGVuIDwtIHdyaXRlX3hsc3goZGF0YSA9IG5ldXRfc2lnW1siZGVzZXEiXV1bWyJkb3ducyJdXVtbMV1dLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9uZXV0cm9waGlsX2NsaW5pY2FsX3NpZ2Rvd25fcGN0LXZ7dmVyfS54bHN4IikpCm5ldXRfY3BtIDwtIHNtKG5vcm1hbGl6ZV9leHB0KG5ldXQsIGNvbnZlcnQgPSAiY3BtIikpCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhID0gZXhwcnMobmV1dF9jcG0pLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9uZXV0cm9waGlsX2NwbV9iZWZvcmVfYmF0Y2gtdnt2ZXJ9Lnhsc3giKSkKbmV1dF9iY3BtIDwtIHNtKG5vcm1hbGl6ZV9leHB0KG5ldXQsIGZpbHRlciA9IFRSVUUsIGJhdGNoID0gInN2YXNlcSIsIGNvbnZlcnQgPSAiY3BtIikpCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhID0gZXhwcnMobmV1dF9iY3BtKSwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvbmV1dHJvcGhpbF9jcG1fYWZ0ZXJfYmF0Y2gtdnt2ZXJ9Lnhsc3giKSkKYGBgCgojIyMjIFdpdGggc3ZhCgpgYGB7ciBkZV9uZXV0X3N2YSwgZmlnLnNob3cgPSAiaGlkZSJ9Cm5ldXRfZGVfc3ZhIDwtIHNtKGFsbF9wYWlyd2lzZShuZXV0LCBtb2RlbF9iYXRjaCA9ICJzdmFzZXEiLCBmaWx0ZXIgPSBUUlVFKSkKbmV1dF90YWJsZXNfc3ZhIDwtIHNtKGNvbWJpbmVfZGVfdGFibGVzKAogICAgbmV1dF9kZV9zdmEsIGtlZXBlcnMgPSBrZWVwZXJzLAogICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9uZXV0cm9waGlsX2NsaW5pY2FsX2FsbF90YWJsZXNfc3ZhLXZ7dmVyfS54bHN4IikpKQpuZXV0X3NpZ19zdmEgPC0gc20oZXh0cmFjdF9zaWduaWZpY2FudF9nZW5lcygKICAgIG5ldXRfdGFibGVzX3N2YSwKICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvbmV1dHJvcGhpbF9jbGluaWNhbF9zaWdfdGFibGVzX3N2YS12e3Zlcn0ueGxzeCIpLAogICAgYWNjb3JkaW5nX3RvID0gImRlc2VxIikpCmBgYAoKIyMjIyBOZXV0cm9waGlsIERFIHBsb3RzCgpgYGB7ciBuZXV0X2RlX3Bsb3RzfQojIyBERVNlcTIgTUEgcGxvdCBvZiBmYWlsdXJlIC8gY3VyZQpuZXV0X3RhYmxlc1tbInBsb3RzIl1dW1siZmFpbF92c19jdXJlIl1dW1siZGVzZXFfbWFfcGxvdHMiXV0kcGxvdAoKIyMgREVTZXEyIFZvbGNhbm8gcGxvdCBvZiBmYWlsdXJlIC8gY3VyZQpuZXV0X3RhYmxlc1tbInBsb3RzIl1dW1siZmFpbF92c19jdXJlIl1dW1siZGVzZXFfdm9sX3Bsb3RzIl1dJHBsb3QKCiMjIERFU2VxMiBNQSBwbG90IG9mIGZhaWx1cmUgLyBjdXJlIHdpdGggc3ZhCm5ldXRfdGFibGVzX3N2YVtbInBsb3RzIl1dW1siZmFpbF92c19jdXJlIl1dW1siZGVzZXFfbWFfcGxvdHMiXV0kcGxvdAoKIyMgREVTZXEyIFZvbGNhbm8gcGxvdCBvZiBmYWlsdXJlIC8gY3VyZSB3aXRoIHN2YQpuZXV0X3RhYmxlc19zdmFbWyJwbG90cyJdXVtbImZhaWxfdnNfY3VyZSJdXVtbImRlc2VxX3ZvbF9wbG90cyJdXSRwbG90CmBgYAoKIyMjIyBOZXV0cm9waGlsIG9udG9sb2d5IHNlYXJjaAoKYGBge3IgbmV1dF9ncH0KdXBzIDwtIG5ldXRfc2lnW1siZGVzZXEiXV1bWyJ1cHMiXV1bWyJmYWlsX3ZzX2N1cmUiXV0KZG93bnMgPC0gbmV1dF9zaWdbWyJkZXNlcSJdXVtbImRvd25zIl1dW1siZmFpbF92c19jdXJlIl1dCgpuZXV0X3VwX2dwIDwtIHNpbXBsZV9ncHJvZmlsZXIoc2lnX2dlbmVzID0gdXBzKQpuZXV0X3VwX2dwW1sicHZhbHVlX3Bsb3RzIl1dW1siYnBwX3Bsb3Rfb3ZlciJdXQpuZXV0X3VwX2dwW1sicHZhbHVlX3Bsb3RzIl1dW1sibWZwX3Bsb3Rfb3ZlciJdXQpuZXV0X3VwX2dwW1sicHZhbHVlX3Bsb3RzIl1dW1sicmVhY3RvbWVfcGxvdF9vdmVyIl1dCgpuZXV0X2Rvd25fZ3AgPC0gc2ltcGxlX2dwcm9maWxlcihkb3ducykKbmV1dF9kb3duX2dwW1sicHZhbHVlX3Bsb3RzIl1dW1siYnBwX3Bsb3Rfb3ZlciJdXQpuZXV0X2Rvd25fZ3BbWyJwdmFsdWVfcGxvdHMiXV1bWyJtZnBfcGxvdF9vdmVyIl1dCm5ldXRfZG93bl9ncFtbInB2YWx1ZV9wbG90cyJdXVtbInJlYWN0b21lX3Bsb3Rfb3ZlciJdXQpgYGAKCiMjIyMgTmV1dHJvcGhpbCBHU1ZBIHF1ZXJ5CgpgYGB7ciBtc2lnX25ldXRfZ29zZXEsIGZpZy5zaG93ID0gImhpZGUifQpuZXV0X3VwX2dvc2VxX21zaWcgPC0gZ29zZXFfbXNpZ2RiKHNpZ19nZW5lcyA9IHVwcywgc2lnbmF0dXJlcyA9IGJyb2FkX2M3LAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHNpZ25hdHVyZV9jYXRlZ29yeSA9ICJjNyIsIGxlbmd0aF9kYiA9IGhzX2xlbmd0aCkKCm5ldXRfZG93bl9nb3NlcV9tc2lnIDwtIGdvc2VxX21zaWdkYihzaWdfZ2VuZXMgPSBkb3ducywgc2lnbmF0dXJlcyA9IGJyb2FkX2M3LAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc2lnbmF0dXJlX2NhdGVnb3J5ID0gImM3IiwgbGVuZ3RoX2RiID0gaHNfbGVuZ3RoKQpgYGAKCiMjIyMgUGxvdCBvZiBzaW1pbGFyIGV4cGVyaW1lbnRzCgpgYGB7ciBtc2lnX3Bsb3RzX25ldXR9CiMjIE5ldXRyb3BoaWwgZ2VuZXMgd2l0aCBpbmNyZWFzZWQgZXhwcmVzc2lvbiBpbiB0aGUgZmFpbGVkIHNhbXBsZXMKIyMgc2hhcmUgZ2VuZXMgd2l0aCB0aGUgZm9sbG93aW5nIGV4cGVyaW1lbnRzCm5ldXRfdXBfZ29zZXFfbXNpZ1tbInB2YWx1ZV9wbG90cyJdXVtbIm1mcF9wbG90X292ZXIiXV0KCiMjIE5ldXRyb3BoaWwgZ2VuZXMgd2l0aCBpbmNyZWFzZWQgZXhwcmVzc2lvbiBpbiB0aGUgY3VyZWQgc2FtcGxlcwojIyBzaGFyZSBnZW5lcyB3aXRoIHRoZSBmb2xsb3dpbmcgZXhwZXJpbWVudHMKbmV1dF9kb3duX2dvc2VxX21zaWdbWyJwdmFsdWVfcGxvdHMiXV1bWyJtZnBfcGxvdF9vdmVyIl1dCmBgYAoKIyMgRW9zaW5vcGhpbHMKCiMjIyBFdmFsdWF0ZSBFb3Npbm9waGlsIHNhbXBsZXMKCmBgYHtyIGVvc2lub3BoaWxzfQplbyA8LSBzdWJzZXRfZXhwdChoc192YWxpZCwgc3Vic2V0ID0gInR5cGVvZmNlbGxzPT0nZW9zaW5vcGhpbHMnIikgJT4lCiAgc2V0X2V4cHRfY29uZGl0aW9ucyhmYWN0ID0gImNsaW5pY2Fsb3V0Y29tZSIpICU+JQogIHNldF9leHB0X2JhdGNoZXMoZmFjdCA9ICJkb25vciIpICU+JQogIHNldF9leHB0X2NvbG9ycyhjb2xvcnMgPSBjaG9zZW5fY29sb3JzKQoKc2F2ZV9yZXN1bHQgPC0gc2F2ZShlbywgZmlsZSA9ICJyZGEvZW9zaW5vcGhpbF9leHB0LnJkYSIpCmVvX25vcm0gPC0gc20obm9ybWFsaXplX2V4cHQoZW8sIGNvbnZlcnQgPSAiY3BtIiwgdHJhbnNmb3JtID0gImxvZzIiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIG5vcm0gPSAicXVhbnQiLCBmaWx0ZXIgPSBUUlVFKSkKcGx0IDwtIHBsb3RfcGNhKGVvX25vcm0sIHBsb3RfbGFiZWxzID0gRkFMU0UpJHBsb3QKcHAoZmlsZSA9IGdsdWUoImltYWdlcy9lb19wY2Ffbm9ybWFsaXplZC12e3Zlcn0ucGRmIiksIGltYWdlID0gcGx0KQoKZW9fbmIgPC0gc20obm9ybWFsaXplX2V4cHQoZW8sIGNvbnZlcnQgPSAiY3BtIiwgdHJhbnNmb3JtID0gImxvZzIiLAogICAgICAgICAgICAgICAgICAgICAgICAgICBmaWx0ZXIgPSBUUlVFLCBiYXRjaCA9ICJzdmFzZXEiKSkKcGxvdF9wY2EoZW9fbmIpJHBsb3QKYGBgCgojIyMgREUgb2YgRW9zaW5vcGhpbCBzYW1wbGVzCgojIyMjIFdpdGhvdXRoIHN2YQoKYGBge3IgZW9zaW5vcGhpbF9kZSwgZmlnLnNob3cgPSAiaGlkZSJ9CmVvX2RlIDwtIHNtKGFsbF9wYWlyd2lzZShlbywgbW9kZWxfYmF0Y2ggPSBGQUxTRSwgZmlsdGVyID0gVFJVRSkpCmVvX3RhYmxlcyA8LSBzbShjb21iaW5lX2RlX3RhYmxlcygKICAgIGVvX2RlLCBrZWVwZXJzID0ga2VlcGVycywKICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvZW9zaW5vcGhpbF9jbGluaWNhbF9hbGxfdGFibGVzLXZ7dmVyfS54bHN4IikpKQp3cml0dGVuIDwtIHdyaXRlX3hsc3goZGF0YSA9IGVvX3RhYmxlc1tbImRhdGEiXV1bWzFdXSwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvZW9zaW5vcGhpbF9jbGluaWNhbF90YWJsZS12e3Zlcn0ueGxzeCIpKQplb19zaWcgPC0gc20oZXh0cmFjdF9zaWduaWZpY2FudF9nZW5lcyhlb190YWJsZXMsIGFjY29yZGluZ190byA9ICJkZXNlcSIpKQp3cml0dGVuIDwtIHdyaXRlX3hsc3goZGF0YSA9IGVvX3NpZ1tbImRlc2VxIl1dW1sidXBzIl1dW1sxXV0sCiAgICAgICAgICAgICAgICAgICAgICBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL2Vvc2lub3BoaWxfY2xpbmljYWxfc2lndXAtdnt2ZXJ9Lnhsc3giKSkKd3JpdHRlbiA8LSB3cml0ZV94bHN4KGRhdGEgPSBlb19zaWdbWyJkZXNlcSJdXVtbImRvd25zIl1dW1sxXV0sCiAgICAgICAgICAgICAgICAgICAgICBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL2Vvc2lub3BoaWxfY2xpbmljYWxfc2lnZG93bi12e3Zlcn0ueGxzeCIpKQoKZW9fcGN0X3NpZyA8LSBzbShleHRyYWN0X3NpZ25pZmljYW50X2dlbmVzKGVvX3RhYmxlcywgbiA9IDIwMCwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGxmYyA9IE5VTEwsIHAgPSBOVUxMLCBhY2NvcmRpbmdfdG8gPSAiZGVzZXEiKSkKd3JpdHRlbiA8LSB3cml0ZV94bHN4KGRhdGEgPSBlb19wY3Rfc2lnW1siZGVzZXEiXV1bWyJ1cHMiXV1bWzFdXSwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvZW9zaW5vcGhpbF9jbGluaWNhbF9zaWd1cF9wY3Qtdnt2ZXJ9Lnhsc3giKSkKd3JpdHRlbiA8LSB3cml0ZV94bHN4KGRhdGEgPSBlb19wY3Rfc2lnW1siZGVzZXEiXV1bWyJkb3ducyJdXVtbMV1dLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9lb3Npbm9waGlsX2NsaW5pY2FsX3NpZ2Rvd25fcGN0LXZ7dmVyfS54bHN4IikpCgplb19jcG0gPC0gc20obm9ybWFsaXplX2V4cHQoZW8sIGNvbnZlcnQgPSAiY3BtIikpCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhID0gZXhwcnMoZW9fY3BtKSwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvZW9zaW5vcGhpbF9jcG1fYmVmb3JlX2JhdGNoLXZ7dmVyfS54bHN4IikpCmVvX2JjcG0gPC0gc20obm9ybWFsaXplX2V4cHQoZW8sIGZpbHRlciA9IFRSVUUsIGJhdGNoID0gInN2YXNlcSIsIGNvbnZlcnQgPSAiY3BtIikpCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhID0gZXhwcnMoZW9fYmNwbSksCiAgICAgICAgICAgICAgICAgICAgICBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL2Vvc2lub3BoaWxfY3BtX2FmdGVyX2JhdGNoLXZ7dmVyfS54bHN4IikpCmBgYAoKIyMjIyBXaXRoIHN2YQoKYGBge3IgZGVfZW9fc3ZhLCBmaWcuc2hvdyA9ICJoaWRlIn0KZW9fZGVfc3ZhIDwtIHNtKGFsbF9wYWlyd2lzZShlbywgbW9kZWxfYmF0Y2ggPSAic3Zhc2VxIiwgZmlsdGVyID0gVFJVRSkpCmVvX3RhYmxlc19zdmEgPC0gc20oY29tYmluZV9kZV90YWJsZXMoCiAgICBlb19kZV9zdmEsIGtlZXBlcnMgPSBrZWVwZXJzLAogICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9lb3Npbm9waGlsX2NsaW5pY2FsX2FsbF90YWJsZXNfc3ZhLXZ7dmVyfS54bHN4IikpKQplb19zaWdfc3ZhIDwtIHNtKGV4dHJhY3Rfc2lnbmlmaWNhbnRfZ2VuZXMoCiAgICBlb190YWJsZXNfc3ZhLAogICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9lb3Npbm9waGlsX2NsaW5pY2FsX3NpZ190YWJsZXNfc3ZhLXZ7dmVyfS54bHN4IiksCiAgICBhY2NvcmRpbmdfdG8gPSAiZGVzZXEiKSkKYGBgCgojIyMjIEVvc2lub3BoaWwgREUgcGxvdHMKCmBgYHtyIGVvX2RlX3Bsb3RzfQojIyBERVNlcTIgTUEgcGxvdCBvZiBmYWlsdXJlIC8gY3VyZQplb190YWJsZXNbWyJwbG90cyJdXVtbImZhaWxfdnNfY3VyZSJdXVtbImRlc2VxX21hX3Bsb3RzIl1dJHBsb3QKCiMjIERFU2VxMiBWb2xjYW5vIHBsb3Qgb2YgZmFpbHVyZSAvIGN1cmUKZW9fdGFibGVzW1sicGxvdHMiXV1bWyJmYWlsX3ZzX2N1cmUiXV1bWyJkZXNlcV92b2xfcGxvdHMiXV0kcGxvdAoKIyMgREVTZXEyIE1BIHBsb3Qgb2YgZmFpbHVyZSAvIGN1cmUgd2l0aCBzdmEKZW9fdGFibGVzX3N2YVtbInBsb3RzIl1dW1siZmFpbF92c19jdXJlIl1dW1siZGVzZXFfbWFfcGxvdHMiXV0kcGxvdAoKIyMgREVTZXEyIFZvbGNhbm8gcGxvdCBvZiBmYWlsdXJlIC8gY3VyZSB3aXRoIHN2YQplb190YWJsZXNfc3ZhW1sicGxvdHMiXV1bWyJmYWlsX3ZzX2N1cmUiXV1bWyJkZXNlcV92b2xfcGxvdHMiXV0kcGxvdApgYGAKCiMjIyMgRW9zaW5vcGhpbCBvbnRvbG9neSBzZWFyY2gKCmBgYHtyIGVvX2dwcm9maWxlcn0KdXBzIDwtIGVvX3NpZ1tbImRlc2VxIl1dW1sidXBzIl1dW1siZmFpbF92c19jdXJlIl1dCmRvd25zIDwtIGVvX3NpZ1tbImRlc2VxIl1dW1siZG93bnMiXV1bWyJmYWlsX3ZzX2N1cmUiXV0KCmVvX3VwX2dwIDwtIHNpbXBsZV9ncHJvZmlsZXIoc2lnX2dlbmVzID0gdXBzKQplb191cF9ncFtbInB2YWx1ZV9wbG90cyJdXVtbImJwcF9wbG90X292ZXIiXV0KZW9fdXBfZ3BbWyJwdmFsdWVfcGxvdHMiXV1bWyJtZnBfcGxvdF9vdmVyIl1dCmVvX3VwX2dwW1sicHZhbHVlX3Bsb3RzIl1dW1sicmVhY3RvbWVfcGxvdF9vdmVyIl1dCgplb19kb3duX2dwIDwtIHNpbXBsZV9ncHJvZmlsZXIoZG93bnMpCmVvX2Rvd25fZ3BbWyJwdmFsdWVfcGxvdHMiXV1bWyJicHBfcGxvdF9vdmVyIl1dCmVvX2Rvd25fZ3BbWyJwdmFsdWVfcGxvdHMiXV1bWyJtZnBfcGxvdF9vdmVyIl1dCmVvX2Rvd25fZ3BbWyJwdmFsdWVfcGxvdHMiXV1bWyJyZWFjdG9tZV9wbG90X292ZXIiXV0KYGBgCgojIyMjIEVvc2lub3BoaWwgTVNpZ0RCIHF1ZXJ5CgpgYGB7ciBtc2lnX2VvX2dvc2VxLCBmaWcuc2hvdyA9ICJoaWRlIn0KZW9fdXBfZ29zZXFfbXNpZyA8LSBnb3NlcV9tc2lnZGIoc2lnX2dlbmVzID0gdXBzLCBzaWduYXR1cmVzID0gYnJvYWRfYzcsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHNpZ25hdHVyZV9jYXRlZ29yeSA9ICJjNyIsIGxlbmd0aF9kYiA9IGhzX2xlbmd0aCkKCmVvX2Rvd25fZ29zZXFfbXNpZyA8LSBnb3NlcV9tc2lnZGIoc2lnX2dlbmVzID0gZG93bnMsIHNpZ25hdHVyZXMgPSBicm9hZF9jNywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzaWduYXR1cmVfY2F0ZWdvcnkgPSAiYzciLCBsZW5ndGhfZGIgPSBoc19sZW5ndGgpCmBgYAoKIyMjIyBQbG90IG9mIHNpbWlsYXIgZXhwZXJpbWVudHMKCmBgYHtyIG1zaWdfcGxvdHNfZW99CiMjIEVvc2lub3BoaWwgZ2VuZXMgd2l0aCBpbmNyZWFzZWQgZXhwcmVzc2lvbiBpbiB0aGUgZmFpbGVkIHNhbXBsZXMKIyMgc2hhcmUgZ2VuZXMgd2l0aCB0aGUgZm9sbG93aW5nIGV4cGVyaW1lbnRzCmVvX3VwX2dvc2VxX21zaWdbWyJwdmFsdWVfcGxvdHMiXV1bWyJtZnBfcGxvdF9vdmVyIl1dCgojIyBFb3Npbm9waGlsIGdlbmVzIHdpdGggaW5jcmVhc2VkIGV4cHJlc3Npb24gaW4gdGhlIGN1cmVkIHNhbXBsZXMKIyMgc2hhcmUgZ2VuZXMgd2l0aCB0aGUgZm9sbG93aW5nIGV4cGVyaW1lbnRzCmVvX2Rvd25fZ29zZXFfbXNpZ1tbInB2YWx1ZV9wbG90cyJdXVtbIm1mcF9wbG90X292ZXIiXV0KYGBgCgojIyBCaW9wc2llcwoKIyMjIEV2YWx1YXRlIEJpb3BzeSBzYW1wbGVzCgpgYGB7ciBiaW9wc2llc30KYmlvcCA8LSBzdWJzZXRfZXhwdChoc192YWxpZCwgc3Vic2V0ID0gInR5cGVvZmNlbGxzPT0nYmlvcHN5JyIpICU+JQogIHNldF9leHB0X2NvbmRpdGlvbnMoZmFjdCA9ICJjbGluaWNhbG91dGNvbWUiKSAlPiUKICBzZXRfZXhwdF9iYXRjaGVzKGZhY3QgPSAiZG9ub3IiKSAlPiUKICBzZXRfZXhwdF9jb2xvcnMoY29sb3JzID0gY2hvc2VuX2NvbG9ycykKCnNhdmVfcmVzdWx0IDwtIHNhdmUoYmlvcCwgZmlsZSA9ICJyZGEvYmlvcHN5X2V4cHQucmRhIikKYmlvcF9ub3JtIDwtIG5vcm1hbGl6ZV9leHB0KGJpb3AsIGZpbHRlciA9IFRSVUUsIGNvbnZlcnQgPSAiY3BtIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgIHRyYW5zZm9ybSA9ICJsb2cyIiwgbm9ybSA9ICJxdWFudCIpCnBsdCA8LSBwbG90X3BjYShiaW9wX25vcm0sIHBsb3RfbGFiZWxzID0gRkFMU0UpJHBsb3QKcHAoZmlsZSA9IGdsdWUoImltYWdlcy9iaW9wX3BjYV9ub3JtYWxpemVkLXZ7dmVyfS5wZGYiKSwgaW1hZ2UgPSBwbHQpCgpiaW9wX25iIDwtIHNtKG5vcm1hbGl6ZV9leHB0KGJpb3AsIGNvbnZlcnQgPSAiY3BtIiwgZmlsdGVyID0gVFJVRSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICB0cmFuc2Zvcm0gPSAibG9nMiIsIGJhdGNoID0gInN2YXNlcSIpKQpwbHQgPC0gcGxvdF9wY2EoYmlvcF9uYiwgcGxvdF9sYWJlbHMgPSBGQUxTRSkkcGxvdApwcChmaWxlID0gZ2x1ZSgiaW1hZ2VzL2Jpb3BfcGNhX25vcm1hbGl6ZWRfc3Zhc2VxLXZ7dmVyfS5wZGYiKSwgaW1hZ2UgPSBwbHQpCmBgYAoKIyMjIERFIG9mIEJpb3BzeSBzYW1wbGVzCgojIyMjIFdpdGhvdXQgc3ZhCgpgYGB7ciBkZV9iaW9wc3ksIGZpZy5zaG93ID0gImhpZGUifQpiaW9wX2RlIDwtIHNtKGFsbF9wYWlyd2lzZShiaW9wLCBtb2RlbF9iYXRjaCA9IEZBTFNFLCBmaWx0ZXIgPSBUUlVFKSkKYmlvcF90YWJsZXMgPC0gY29tYmluZV9kZV90YWJsZXMoYmlvcF9kZSwga2VlcGVycyA9IGtlZXBlcnMsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvYmlvcHN5X2NsaW5pY2FsX2FsbF90YWJsZXMtdnt2ZXJ9Lnhsc3giKSkKd3JpdHRlbiA8LSB3cml0ZV94bHN4KGRhdGEgPSBiaW9wX3RhYmxlc1tbImRhdGEiXV1bWzFdXSwKICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvYmlvcHN5X2NsaW5pY2FsX3RhYmxlLXZ7dmVyfS54bHN4IikpCmJpb3Bfc2lnIDwtIGV4dHJhY3Rfc2lnbmlmaWNhbnRfZ2VuZXMoYmlvcF90YWJsZXMsIGFjY29yZGluZ190byA9ICJkZXNlcSIpCiMjd3JpdHRlbiA8LSB3cml0ZV94bHN4KGRhdGEgPSBiaW9wX3NpZ1tbImRlc2VxIl1dW1sidXBzIl1dW1sxXV0sCiMjICAgICAgICAgICAgICAgICAgICAgIGV4Y2VsID0gZ2x1ZTo6Z2x1ZSgiZXhjZWwvYmlvcHN5X2NsaW5pY2FsX3NpZ3VwLXZ7dmVyfS54bHN4IikpCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhID0gYmlvcF9zaWdbWyJkZXNlcSJdXVtbImRvd25zIl1dW1sxXV0sCiAgICAgICAgICAgICAgICAgICAgICBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL2Jpb3BzeV9jbGluaWNhbF9zaWdkb3duLXZ7dmVyfS54bHN4IikpCmJpb3BfcGN0X3NpZyA8LSBleHRyYWN0X3NpZ25pZmljYW50X2dlbmVzKGJpb3BfdGFibGVzLCBuID0gMjAwLCBsZmMgPSBOVUxMLCBwID0gTlVMTCwgYWNjb3JkaW5nX3RvID0gImRlc2VxIikKd3JpdHRlbiA8LSB3cml0ZV94bHN4KGRhdGEgPSBiaW9wX3BjdF9zaWdbWyJkZXNlcSJdXVtbInVwcyJdXVtbMV1dLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9iaW9wc3lfY2xpbmljYWxfc2lndXBfcGN0LXZ7dmVyfS54bHN4IikpCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhID0gYmlvcF9wY3Rfc2lnW1siZGVzZXEiXV1bWyJkb3ducyJdXVtbMV1dLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9iaW9wc3lfY2xpbmljYWxfc2lnZG93bl9wY3Qtdnt2ZXJ9Lnhsc3giKSkKCmJpb3BfY3BtIDwtIHNtKG5vcm1hbGl6ZV9leHB0KGJpb3AsIGNvbnZlcnQgPSAiY3BtIikpCndyaXR0ZW4gPC0gd3JpdGVfeGxzeChkYXRhID0gZXhwcnMoYmlvcF9jcG0pLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9iaW9wc3lfY3BtX2JlZm9yZV9iYXRjaC12e3Zlcn0ueGxzeCIpKQpiaW9wX2JjcG0gPC0gc20obm9ybWFsaXplX2V4cHQoYmlvcCwgZmlsdGVyID0gVFJVRSwgYmF0Y2ggPSAic3Zhc2VxIiwgY29udmVydCA9ICJjcG0iKSkKd3JpdHRlbiA8LSB3cml0ZV94bHN4KGRhdGEgPSBleHBycyhiaW9wX2JjcG0pLAogICAgICAgICAgICAgICAgICAgICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9iaW9wc3lfY3BtX2FmdGVyX2JhdGNoLXZ7dmVyfS54bHN4IikpCmBgYAoKIyMjIyB3aXRoIHN2YQoKYGBge3IgZGVfYmlvcHN5X3N2YSwgZmlnLnNob3cgPSAiaGlkZSJ9CmJpb3BfZGVfc3ZhIDwtIHNtKGFsbF9wYWlyd2lzZShiaW9wLCBtb2RlbF9iYXRjaCA9ICJzdmFzZXEiLCBmaWx0ZXIgPSBUUlVFKSkKYmlvcF90YWJsZXNfc3ZhIDwtIHNtKGNvbWJpbmVfZGVfdGFibGVzKAogICAgYmlvcF9kZV9zdmEsIGtlZXBlcnMgPSBrZWVwZXJzLAogICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9iaW9wc3lfY2xpbmljYWxfYWxsX3RhYmxlc19zdmEtdnt2ZXJ9Lnhsc3giKSkpCmJpb3Bfc2lnX3N2YSA8LSBzbShleHRyYWN0X3NpZ25pZmljYW50X2dlbmVzKAogICAgYmlvcF90YWJsZXNfc3ZhLAogICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9iaW9wc3lfY2xpbmljYWxfc2lnX3RhYmxlc19zdmEtdnt2ZXJ9Lnhsc3giKSwKICAgIGFjY29yZGluZ190byA9ICJkZXNlcSIpKQpgYGAKCiMjIyMgQmlvcHN5IERFIHBsb3RzCgpgYGB7ciBiaW9wX2RlX3Bsb3RzfQojIyBERVNlcTIgTUEgcGxvdCBvZiBmYWlsdXJlIC8gY3VyZQpiaW9wX3RhYmxlc1tbInBsb3RzIl1dW1siZmFpbF92c19jdXJlIl1dW1siZGVzZXFfbWFfcGxvdHMiXV0kcGxvdAoKIyMgREVTZXEyIFZvbGNhbm8gcGxvdCBvZiBmYWlsdXJlIC8gY3VyZQpiaW9wX3RhYmxlc1tbInBsb3RzIl1dW1siZmFpbF92c19jdXJlIl1dW1siZGVzZXFfdm9sX3Bsb3RzIl1dJHBsb3QKCiMjIERFU2VxMiBNQSBwbG90IG9mIGZhaWx1cmUgLyBjdXJlCmJpb3BfdGFibGVzX3N2YVtbInBsb3RzIl1dW1siZmFpbF92c19jdXJlIl1dW1siZGVzZXFfbWFfcGxvdHMiXV0kcGxvdAoKIyMgREVTZXEyIFZvbGNhbm8gcGxvdCBvZiBmYWlsdXJlIC8gY3VyZQpiaW9wX3RhYmxlc19zdmFbWyJwbG90cyJdXVtbImZhaWxfdnNfY3VyZSJdXVtbImRlc2VxX3ZvbF9wbG90cyJdXSRwbG90CmBgYAoKIyBMb29rIGZvciBzaGFyZWQgZ2VuZXMgYW1vbmcgTW9ub2N5dGVzL05ldXRyb3BoaWxzL0Vvc2lub3BoaWxzCgpXZSBoYXZlIHRocmVlIHZhcmlhYmxlcyBjb250YWluaW5nIHRoZSAnc2lnbmlmaWNhbnQnIERFIGdlbmVzIGZvciB0aGUKdGhyZWUgY2VsbCB0eXBlcy4gIEZvciB0aGlzIEkgYW0gY2hvb3NpbmcgKGZvciB0aGUgbW9tZW50KSB0byB1c2UgdGhlCnN2YSBkYXRhLgoKYGBge3Igc2hhcmVkX2J5X3R5cGV9CiMjIG1vbm9fc2lnX3N2YSwgbmV1dF9zaWdfc3ZhLCBlb19zaWdfc3ZhCnNpZ192ZWN0b3JzIDwtIGxpc3QoCiAgICAibW9ub2N5dGVzIiA9IGMocm93bmFtZXMobW9ub19zaWdfc3ZhW1siZGVzZXEiXV1bWyJ1cHMiXV1bWyJmYWlsX3ZzX2N1cmUiXV0pLAogICAgICAgICAgICAgICAgICAgIHJvd25hbWVzKG1vbm9fc2lnX3N2YVtbImRlc2VxIl1dW1siZG93bnMiXV1bWyJmYWlsX3ZzX2N1cmUiXV0pKSwKICAgICJuZXV0cm9waGlscyIgPSBjKHJvd25hbWVzKG5ldXRfc2lnX3N2YVtbImRlc2VxIl1dW1sidXBzIl1dW1siZmFpbF92c19jdXJlIl1dKSwKICAgICAgICAgICAgICAgICAgICAgIHJvd25hbWVzKG5ldXRfc2lnX3N2YVtbImRlc2VxIl1dW1siZG93bnMiXV1bWyJmYWlsX3ZzX2N1cmUiXV0pKSwKICAgICJlb3Npbm9waGlscyIgPSAgYyhyb3duYW1lcyhlb19zaWdfc3ZhW1siZGVzZXEiXV1bWyJ1cHMiXV1bWyJmYWlsX3ZzX2N1cmUiXV0pLAogICAgICAgICAgICAgICAgICAgICAgIHJvd25hbWVzKGVvX3NpZ19zdmFbWyJkZXNlcSJdXVtbImRvd25zIl1dW1siZmFpbF92c19jdXJlIl1dKSkpCgpzaGFyZWRfdmVjdG9yIDwtIFZlbm5lcmFibGU6OlZlbm4oU2V0cyA9IHNpZ192ZWN0b3JzKQpWZW5uZXJhYmxlOjpwbG90KHNoYXJlZF92ZWN0b3IsIGRvV2VpZ2h0cyA9IEZBTFNFKQoKc2hhcmVkX2lkcyA8LSBzaGFyZWRfdmVjdG9yQEludGVyc2VjdGlvblNldHNbWyIxMTEiXV0Kc2hhcmVkX2V4cHQgPC0gZXhjbHVkZV9nZW5lc19leHB0KGhzX2NsaW5pY2FsLCBpZHMgPSBzaGFyZWRfaWRzLCBtZXRob2QgPSAia2VlcCIpCnNoYXJlZF93cml0dGVuIDwtIHdyaXRlX2V4cHQoCiAgICBzaGFyZWRfZXhwdCwgbm9ybSA9ICJyYXciLApgYGAKCiMgTW9ub2N5dGVzIGJ5IHZpc2l0CgogMS4gQ2FuIHlvdSBwbGVhc2Ugc2hhcmUgd2l0aCB1cyBhIFBDQSAoU1ZBIGFuZCBub24tU1ZBKSBvZiB0aGUKICAgIG1vbm9jeXRlcyBvZiB0aGUgVE1SQzMgcHJvamVjdCwgYnV0IGxhYmVsaW5nIHRoZW0gYmFzZWQgb24gdGhlIHZpc2l0CiAgICAoVjEsIFYyLCBWMyk/CiAyLiBDYW4geW91IHBsZWFzZSBzaGFyZSBERSBsaXN0cyBvZiBWMSB2cyBWMiwgVjEgdnMgVjMsIFYxIHZzLiBWMitWMwogICAgYW5kIFYyIHZzIFYzPwoKYGBge3IgbW9ub2N5dGVzX2J5X3Zpc2l0fQp2aXNpdF9jb2xvcnMgPC0gY2hvc2VuX2NvbG9ycyA8LSBjKCIjRDk1RjAyIiwgIiM3NTcwQjMiLCAiIzFCOUU3NyIpCm5hbWVzKHZpc2l0X2NvbG9ycykgPC0gYygxLCAyLCAzKQptb25vX3Zpc2l0IDwtIHN1YnNldF9leHB0KGhzX3ZhbGlkLCBzdWJzZXQgPSAidHlwZW9mY2VsbHM9PSdtb25vY3l0ZXMnIikgJT4lCiAgc2V0X2V4cHRfY29uZGl0aW9ucyhmYWN0ID0gInZpc2l0bnVtYmVyIikgJT4lCiAgc2V0X2V4cHRfYmF0Y2hlcyhmYWN0ID0gImNsaW5pY2Fsb3V0Y29tZSIpICU+JQogIHNldF9leHB0X2NvbG9ycyhjb2xvcnMgPSBjaG9zZW5fY29sb3JzKQoKbW9ub192aXNpdF9ub3JtIDwtIG5vcm1hbGl6ZV9leHB0KG1vbm9fdmlzaXQsIGZpbHRlciA9IFRSVUUsIG5vcm0gPSAicXVhbnQiLCBjb252ZXJ0ID0gImNwbSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB0cmFuc2Zvcm0gPSAibG9nMiIpCm1vbm9fdmlzaXRfcGNhIDwtIHBsb3RfcGNhKG1vbm9fdmlzaXRfbm9ybSkKcHAoZmlsZSA9ICJpbWFnZXMvbW9ub2N5dGVfYnlfdmlzaXQucG5nIiwgaW1hZ2UgPSBtb25vX3Zpc2l0X3BjYSRwbG90KQoKbW9ub192aXNpdF9uYiA8LSBub3JtYWxpemVfZXhwdChtb25vX3Zpc2l0LCBmaWx0ZXIgPSBUUlVFLCBjb252ZXJ0ID0gImNwbSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYmF0Y2ggPSAic3Zhc2VxIiwgdHJhbnNmb3JtID0gImxvZzIiKQptb25vX3Zpc2l0X25iX3BjYSA8LSBwbG90X3BjYShtb25vX3Zpc2l0X25iKQpwcChmaWxlID0gImltYWdlcy9tb25vY3l0ZV9ieV92aXNpdF9uYi5wbmciLCBpbWFnZSA9IG1vbm9fdmlzaXRfbmJfcGNhJHBsb3QpCgp0YWJsZShwRGF0YShtb25vX3Zpc2l0X25vcm0pJGJhdGNoKQpgYGAKCmBgYHtyIG1vbm9fdmlzaXRfZGUsIGZpZy5zaG93ID0gImhpZGUifQprZWVwZXJzIDwtIGxpc3QoCiAgICAic2Vjb25kX3ZzX2ZpcnN0IiA9IGMoImMyIiwgImMxIiksCiAgICAidGhpcmRfdnNfc2Vjb25kIiA9IGMoImMzIiwgImMyIiksCiAgICAidGhpcmRfdnNfZmlyc3QiID0gYygiYzMiLCAiYzEiKSkKbW9ub192aXNpdF9kZSA8LSBhbGxfcGFpcndpc2UobW9ub192aXNpdCwgbW9kZWxfYmF0Y2ggPSAic3Zhc2VxIiwgZmlsdGVyID0gVFJVRSkKCm1vbm9fdmlzaXRfdGFibGVzIDwtIGNvbWJpbmVfZGVfdGFibGVzKAogICAgbW9ub192aXNpdF9kZSwKICAgIGtlZXBlcnMgPSBrZWVwZXJzLAogICAgZXhjZWwgPSBnbHVlOjpnbHVlKCJleGNlbC9tb25vX3Zpc2l0X3RhYmxlcy12e3Zlcn0ueGxzeCIpKQpgYGAKCmBgYHtyIHYxX3ZzX2FsbH0KbmV3X2ZhY3RvciA8LSBhcy5jaGFyYWN0ZXIocERhdGEobW9ub192aXNpdClbWyJ2aXNpdG51bWJlciJdXSkKbm90X29uZV9pZHggPC0gbmV3X2ZhY3RvciAhPSAxCm5ld19mYWN0b3Jbbm90X29uZV9pZHhdIDwtICJub3RfMSIKbW9ub19vbmVfdnMgPC0gc2V0X2V4cHRfY29uZGl0aW9ucyhtb25vX3Zpc2l0LCBuZXdfZmFjdG9yKQoKbW9ub19vbmVfdnNfZGUgPC0gYWxsX3BhaXJ3aXNlKG1vbm9fb25lX3ZzLCBtb2RlbF9iYXRjaCA9ICJzdmFzZXEiLCBmaWx0ZXIgPSBUUlVFKQoKbW9ub19vbmVfdnNfdGFibGVzIDwtIGNvbWJpbmVfZGVfdGFibGVzKAogICAgbW9ub19vbmVfdnNfZGUsCiAgICBleGNlbCA9IGdsdWU6OmdsdWUoImV4Y2VsL21vbm9fb25lX3ZzX3RhYmxlcy12e3Zlcn0ueGxzeCIpKQpgYGAKCiMgVGVzdCBUU1AKCkluIHdyaXRpbmcgdGhlIGZvbGxvd2luZywgSSBxdWlja2x5IHJlYWxpemVkIHRoYXQgdHNwYWlyIHdhcyBub3QKam9raW5nIHdoZW4gaXQgc2FpZCBpdCBpcyBpbnRlbmRlZCBmb3Igc21hbGwgbnVtYmVycyBvZiBnZW5lcy4gIEZvciBhCmZ1bGwgZXhwcmVzc2lvbnNldCBvZiBodW1hbiBkYXRhIGl0IGlzIHN0cnVnZ2xpbmcuICBJIGxpa2UgdGhlIGlkZWEsCml0IG1heSBwcm92ZSB3b3J0aCB3aGlsZSB0byBzcGVuZCBzb21lIHRpbWUgb3B0aW1pemluZyB0aGUgcGFja2FnZSBzbwp0aGF0IGl0IGlzIG1vcmUgdXNhYmxlLgoKYGBge3IgdHNwLCBldmFsID0gRkFMU0V9CmV4cHQgPC0gaHNfY2xpbmljYWxfbm9iaW9wCgpzaW1wbGVfdHNwIDwtIGZ1bmN0aW9uKGV4cHQsIGNvbHVtbiA9ICJjb25kaXRpb24iKSB7CiAgZmFjdHMgPC0gbGV2ZWxzKGFzLmZhY3RvcihwRGF0YShleHB0KVtbY29sdW1uXV0pKQogIHJldGxpc3QgPC0gbGlzdCgpCiAgaWYgKGxlbmd0aChmYWN0cykgPCAyKSB7CiAgICBzdG9wKCJUaGlzIHJlcXVpcmVzIGZhY3RvcnMgd2l0aCBhdCBsZWFzdCAyIGxldmVscy4iKQogIH0gZWxzZSBpZiAobGVuZ3RoKGZhY3RzKSA9PSAyKSB7CiAgICByZXRsaXN0IDwtIHNpbXBsZV90c3BfcGFpcihleHB0LCBjb2x1bW4gPSBjb2x1bW4pCiAgfSBlbHNlIHsKICAgIGZvciAoZmlyc3QgaW4gMToobGVuZ3RoKGZhY3RzKSAtIDEpKSB7CiAgICAgIGZvciAoc2Vjb25kIGluIDI6KGxlbmd0aChmYWN0cykpKSB7CiAgICAgICAgaWYgKGZpcnN0IDwgc2Vjb25kKSB7CiAgICAgICAgICBuYW1lIDwtIGdsdWU6OmdsdWUoIntmYWN0c1tmaXJzdF19X3ZzX3tmYWN0c1tzZWNvbmRdfSIpCiAgICAgICAgICBtZXNzYWdlKCJTdGFydGluZyAiLCBuYW1lLCAiLiIpCiAgICAgICAgICBzdWJzdHJpbmcgPC0gZ2x1ZTo6Z2x1ZSgie2NvbHVtbn09PSd7ZmFjdHNbZmlyc3RdfSd8e2NvbHVtbn09PSd7ZmFjdHNbc2Vjb25kXX0nIikKICAgICAgICAgIHN1YmJ5IDwtIHN1YnNldF9leHB0KGV4cHQsIHN1YnNldD1hcy5jaGFyYWN0ZXIoc3Vic3RyaW5nKSkKICAgICAgICAgIHJldGxpc3RbW25hbWVdXSA8LSBzaW1wbGVfdHNwX3BhaXIoc3ViYnksIGNvbHVtbiA9IGNvbHVtbikKICAgICAgICB9CiAgICAgIH0KICAgIH0KICB9Cn0KCnNpbXBsZV90c3BfcGFpciA8LSBmdW5jdGlvbihzdWJieSwgY29sdW1uID0gImNvbmRpdGlvbiIsIHJlcGV0aXRpb25zID0gNTApIHsKICB0c3BfaW5wdXQgPC0gc3ViYnlbWyJleHByZXNzaW9uc2V0Il1dCiAgdHNwX291dHB1dCA8LSB0c3BjYWxjKHRzcF9pbnB1dCwgY29sdW1uKQogIHRzcF9zY29yZXMgPC0gdHNwc2lnKHRzcF9pbnB1dCwgY29sdW1uLCBCID0gcmVwZXRpdGlvbnMpCn0KCnRzcDEgPC0gdHNwY2FsYyh0c3BfaW5wdXQsICJjb25kaXRpb24iKQoKYGBgCgpgYGB7ciBzYXZlbWV9CmlmICghaXNUUlVFKGdldDAoInNraXBfbG9hZCIpKSkgewogIHBhbmRlcjo6cGFuZGVyKHNlc3Npb25JbmZvKCkpCiAgbWVzc2FnZShwYXN0ZTAoIlRoaXMgaXMgaHBnbHRvb2xzIGNvbW1pdDogIiwgZ2V0X2dpdF9jb21taXQoKSkpCiAgbWVzc2FnZShwYXN0ZTAoIlNhdmluZyB0byAiLCBzYXZlZmlsZSkpCiAgdG1wIDwtIHNtKHNhdmVtZShmaWxlbmFtZSA9IHNhdmVmaWxlKSkKfQpgYGAKCmBgYHtyIGxvYWRtZV9hZnRlciwgZXZhbCA9IEZBTFNFfQp0bXAgPC0gbG9hZG1lKGZpbGVuYW1lID0gc2F2ZWZpbGUpCmBgYAo=