Thesis SIGNAL Mutational Signatures GLMM Comparison Script
ALL IMPORTED FILES MUST BE IN THE SAME DIRECTORY AS THIS SCRIPT
High-performance computing cluster (e.g. using RStudio Server) strongly recommended, due to large memory requirements
Load required packages
library("tidyverse")
library("readbulk")
library("glmmTMB")
library("betareg")
library("DHARMa")Import data; each file comprises 1000 SIGNAL bootstrap solutions (generated using https://signal.mutationalsignatures.com/analyse2) using the OvCa-specific signature set (OVARY_A/B/C/D/E/F/G- see https://signal.mutationalsignatures.com/explore/studyTissueType/1-15) for the 12 samples comprising the pooled reference group (‘HGSrefMutSig’), along with either the 2-3 samples comprising each pooled test group (‘ATM/PALB2a/PALB2b/LLGL2/SCYL3’) or an individual sample (‘sample’); refer to Chapters 4.2.2.3-4 for further details
mut_sig_SIGNAL_exposures_ATM <- read_tsv("masterfile_tumour_samples_SIGNAL_HGSrefMutSig_ATM.tsv")
mut_sig_SIGNAL_exposures_PALB2a <- read_tsv("masterfile_tumour_samples_SIGNAL_HGSrefMutSig_PALB2a.tsv")
mut_sig_SIGNAL_exposures_PALB2b <- read_tsv("masterfile_tumour_samples_SIGNAL_HGSrefMutSig_PALB2b.tsv")
mut_sig_SIGNAL_exposures_LLGL2 <- read_tsv("masterfile_tumour_samples_SIGNAL_HGSrefMutSig_LLGL2.tsv")
mut_sig_SIGNAL_exposures_SCYL3 <- read_tsv("masterfile_tumour_samples_SIGNAL_HGSrefMutSig_SCYL3.tsv")
mut_sig_SIGNAL_exposures_sample <- read_tsv("masterfile_tumour_samples_SIGNAL_HGSrefMutSig_sample.tsv")Round down 1s and convert data into factors
mut_sig_SIGNAL_exposures_ATM$Exposures <- replace(mut_sig_SIGNAL_exposures_ATM$Exposures, mut_sig_SIGNAL_exposures_ATM$Exposures==1, 0.9999999)
mut_sig_SIGNAL_exposures_ATM$Group <- as.factor(mut_sig_SIGNAL_exposures_ATM$Group)
mut_sig_SIGNAL_exposures_ATM$Tumour_Sample <- as.factor(mut_sig_SIGNAL_exposures_ATM$Tumour_Sample)
mut_sig_SIGNAL_exposures_ATM$Solution <- as.factor(mut_sig_SIGNAL_exposures_ATM$Solution)
mut_sig_SIGNAL_exposures_ATM$Solu_Sample <- paste(mut_sig_SIGNAL_exposures_ATM$Solution,
mut_sig_SIGNAL_exposures_ATM$Tumour_Sample,
sep=":") %>%
as.factor()
mut_sig_SIGNAL_exposures_ATM$Signature <- as.factor(mut_sig_SIGNAL_exposures_ATM$Signature)
mut_sig_SIGNAL_exposures_PALB2a$Exposures <- replace(mut_sig_SIGNAL_exposures_PALB2a$Exposures, mut_sig_SIGNAL_exposures_PALB2a$Exposures==1, 0.9999999)
mut_sig_SIGNAL_exposures_PALB2a$Group <- as.factor(mut_sig_SIGNAL_exposures_PALB2a$Group)
mut_sig_SIGNAL_exposures_PALB2a$Tumour_Sample <- as.factor(mut_sig_SIGNAL_exposures_PALB2a$Tumour_Sample)
mut_sig_SIGNAL_exposures_PALB2a$Solution <- as.factor(mut_sig_SIGNAL_exposures_PALB2a$Solution)
mut_sig_SIGNAL_exposures_PALB2a$Solu_Sample <- paste(mut_sig_SIGNAL_exposures_PALB2a$Solution,
mut_sig_SIGNAL_exposures_PALB2a$Tumour_Sample,
sep=":") %>%
as.factor()
mut_sig_SIGNAL_exposures_PALB2a$Signature <- as.factor(mut_sig_SIGNAL_exposures_PALB2a$Signature)
mut_sig_SIGNAL_exposures_PALB2b$Exposures <- replace(mut_sig_SIGNAL_exposures_PALB2b$Exposures, mut_sig_SIGNAL_exposures_PALB2b$Exposures==1, 0.9999999)
mut_sig_SIGNAL_exposures_PALB2b$Group <- as.factor(mut_sig_SIGNAL_exposures_PALB2b$Group)
mut_sig_SIGNAL_exposures_PALB2b$Tumour_Sample <- as.factor(mut_sig_SIGNAL_exposures_PALB2b$Tumour_Sample)
mut_sig_SIGNAL_exposures_PALB2b$Solution <- as.factor(mut_sig_SIGNAL_exposures_PALB2b$Solution)
mut_sig_SIGNAL_exposures_PALB2b$Solu_Sample <- paste(mut_sig_SIGNAL_exposures_PALB2b$Solution,
mut_sig_SIGNAL_exposures_PALB2b$Tumour_Sample,
sep=":") %>%
as.factor()
mut_sig_SIGNAL_exposures_PALB2b$Signature <- as.factor(mut_sig_SIGNAL_exposures_PALB2b$Signature)
mut_sig_SIGNAL_exposures_LLGL2$Exposures <- replace(mut_sig_SIGNAL_exposures_LLGL2$Exposures, mut_sig_SIGNAL_exposures_LLGL2$Exposures==1, 0.9999999)
mut_sig_SIGNAL_exposures_LLGL2$Group <- as.factor(mut_sig_SIGNAL_exposures_LLGL2$Group)
mut_sig_SIGNAL_exposures_LLGL2$Tumour_Sample <- as.factor(mut_sig_SIGNAL_exposures_LLGL2$Tumour_Sample)
mut_sig_SIGNAL_exposures_LLGL2$Solution <- as.factor(mut_sig_SIGNAL_exposures_LLGL2$Solution)
mut_sig_SIGNAL_exposures_LLGL2$Solu_Sample <- paste(mut_sig_SIGNAL_exposures_LLGL2$Solution,
mut_sig_SIGNAL_exposures_LLGL2$Tumour_Sample,
sep=":") %>%
as.factor()
mut_sig_SIGNAL_exposures_LLGL2$Signature <- as.factor(mut_sig_SIGNAL_exposures_LLGL2$Signature)
mut_sig_SIGNAL_exposures_SCYL3$Exposures <- replace(mut_sig_SIGNAL_exposures_SCYL3$Exposures, mut_sig_SIGNAL_exposures_SCYL3$Exposures==1, 0.9999999)
mut_sig_SIGNAL_exposures_SCYL3$Group <- as.factor(mut_sig_SIGNAL_exposures_SCYL3$Group)
mut_sig_SIGNAL_exposures_SCYL3$Tumour_Sample <- as.factor(mut_sig_SIGNAL_exposures_SCYL3$Tumour_Sample)
mut_sig_SIGNAL_exposures_SCYL3$Solution <- as.factor(mut_sig_SIGNAL_exposures_SCYL3$Solution)
mut_sig_SIGNAL_exposures_SCYL3$Solu_Sample <- paste(mut_sig_SIGNAL_exposures_SCYL3$Solution,
mut_sig_SIGNAL_exposures_SCYL3$Tumour_Sample,
sep=":") %>%
as.factor()
mut_sig_SIGNAL_exposures_SCYL3$Signature <- as.factor(mut_sig_SIGNAL_exposures_SCYL3$Signature)
mut_sig_SIGNAL_exposures_SCYL3$Exposures <- replace(mut_sig_SIGNAL_exposures_SCYL3$Exposures, mut_sig_SIGNAL_exposures_SCYL3$Exposures==1, 0.9999999)
mut_sig_SIGNAL_exposures_SCYL3$Group <- as.factor(mut_sig_SIGNAL_exposures_SCYL3$Group)
mut_sig_SIGNAL_exposures_SCYL3$Tumour_Sample <- as.factor(mut_sig_SIGNAL_exposures_SCYL3$Tumour_Sample)
mut_sig_SIGNAL_exposures_SCYL3$Solution <- as.factor(mut_sig_SIGNAL_exposures_SCYL3$Solution)
mut_sig_SIGNAL_exposures_SCYL3$Solu_Sample <- paste(mut_sig_SIGNAL_exposures_SCYL3$Solution,
mut_sig_SIGNAL_exposures_SCYL3$Tumour_Sample,
sep=":") %>%
as.factor()
mut_sig_SIGNAL_exposures_SCYL3$Signature <- as.factor(mut_sig_SIGNAL_exposures_SCYL3$Signature)
mut_sig_SIGNAL_exposures_sample$Exposures <- replace(mut_sig_SIGNAL_exposures_sample$Exposures, mut_sig_SIGNAL_exposures_sample$Exposures==1, 0.9999999)
mut_sig_SIGNAL_exposures_sample$Group <- as.factor(mut_sig_SIGNAL_exposures_sample$Group)
mut_sig_SIGNAL_exposures_sample$Tumour_Sample <- as.factor(mut_sig_SIGNAL_exposures_sample$Tumour_Sample)
mut_sig_SIGNAL_exposures_sample$Solution <- as.factor(mut_sig_SIGNAL_exposures_sample$Solution)
mut_sig_SIGNAL_exposures_sample$Solu_Sample <- paste(mut_sig_SIGNAL_exposures_sample$Solution,
mut_sig_SIGNAL_exposures_sample$Tumour_Sample,
sep=":") %>%
as.factor()
mut_sig_SIGNAL_exposures_sample$Signature <- as.factor(mut_sig_SIGNAL_exposures_sample$Signature)Run GLMM with beta-regression for combined HGSrefMutSig vs combined and individual samples of interest, and simulate residuals for model using DHARMa
mut_sig_mixed_lmer_bin_ATM <- glmmTMB(Exposures ~ Age + Group +
(1|Solu_Sample) + (1|Signature),
data = mut_sig_SIGNAL_exposures_ATM,
family=beta_family(),
ziformula = ~1,
control = glmmTMBControl(parallel = 48) #parallel = no. of available processor cores/threads
)
simres_mut_sig_mixed_lmer_bin_ATM <- simulateResiduals(mut_sig_mixed_lmer_bin_ATM)
mut_sig_mixed_lmer_bin_PALB2a <- glmmTMB(Exposures ~ Age + Group +
(1|Solu_Sample) + (1|Signature),
data = mut_sig_SIGNAL_exposures_PALB2a,
family=beta_family(),
ziformula = ~1,
control = glmmTMBControl(parallel = 48)
)
simres_mut_sig_mixed_lmer_bin_PALB2a <- simulateResiduals(mut_sig_mixed_lmer_bin_PALB2a)
mut_sig_mixed_lmer_bin_PALB2b <- glmmTMB(Exposures ~ Age + Group +
(1|Solu_Sample) + (1|Signature),
data = mut_sig_SIGNAL_exposures_PALB2b,
family=beta_family(),
ziformula = ~1,
control = glmmTMBControl(parallel = 48)
)
simres_mut_sig_mixed_lmer_bin_PALB2b <- simulateResiduals(mut_sig_mixed_lmer_bin_PALB2b)
mut_sig_mixed_lmer_bin_LLGL2 <- glmmTMB(Exposures ~ Age + Group +
(1|Solu_Sample) + (1|Signature),
data = mut_sig_SIGNAL_exposures_LLGL2,
family=beta_family(),
ziformula = ~1,
control = glmmTMBControl(parallel = 48)
)
simres_mut_sig_mixed_lmer_bin_LLGL2 <- simulateResiduals(mut_sig_mixed_lmer_bin_LLGL2)
mut_sig_mixed_lmer_bin_SCYL3 <- glmmTMB(Exposures ~ Age + Group +
(1|Solu_Sample) + (1|Signature),
data = mut_sig_SIGNAL_exposures_SCYL3,
family=beta_family(),
ziformula = ~1,
control = glmmTMBControl(parallel = 48)
)
simres_mut_sig_mixed_lmer_bin_SCYL3 <- simulateResiduals(mut_sig_mixed_lmer_bin_SCYL3)
mut_sig_mixed_lmer_bin_sample <- glmmTMB(Exposures ~ Age + Group +
(1|Solu_Sample) + (1|Signature),
data = mut_sig_SIGNAL_exposures_sample,
family=beta_family(),
ziformula = ~1,
control = glmmTMBControl(parallel = 48)
)
simres_mut_sig_mixed_lmer_bin_sample <- simulateResiduals(mut_sig_mixed_lmer_bin_sample)Summary of GLMM results
summary(mut_sig_mixed_lmer_bin_ATM) Family: beta ( logit )
Formula: Exposures ~ Age + Group + (1 | Solu_Sample) + (1 | Signature)
Zero inflation: ~1
Data: mut_sig_SIGNAL_exposures_ATM
AIC BIC logLik deviance df.resid
44048.1 44115.0 -22017.0 44034.1 104993
Random effects:
Conditional model:
Groups Name Variance Std.Dev.
Solu_Sample (Intercept) 0.05436 0.2332
Signature (Intercept) 0.17581 0.4193
Number of obs: 105000, groups: Solu_Sample, 15000; Signature, 7
Dispersion parameter for beta family (): 5.39
Conditional model:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.1272493 0.1632022 -13.03 <2e-16 ***
Age 0.0129747 0.0005831 22.25 <2e-16 ***
GroupHGSrefMutSig -0.1288221 0.0101442 -12.70 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Zero-inflation model:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.633850 0.006485 -97.75 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1summary(mut_sig_mixed_lmer_bin_PALB2a) Family: beta ( logit )
Formula: Exposures ~ Age + Group + (1 | Solu_Sample) + (1 | Signature)
Zero inflation: ~1
Data: mut_sig_SIGNAL_exposures_PALB2a
AIC BIC logLik deviance df.resid
30004.1 30070.5 -14995.0 29990.1 97993
Random effects:
Conditional model:
Groups Name Variance Std.Dev.
Solu_Sample (Intercept) 0.04277 0.2068
Signature (Intercept) 0.19685 0.4437
Number of obs: 98000, groups: Solu_Sample, 14000; Signature, 7
Dispersion parameter for beta family (): 5.82
Conditional model:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.2596816 0.1711655 -13.202 < 2e-16 ***
Age 0.0124926 0.0005535 22.571 < 2e-16 ***
GroupPALB2a 0.0528082 0.0120234 4.392 1.12e-05 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Zero-inflation model:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.732795 0.006822 -107.4 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1summary(mut_sig_mixed_lmer_bin_PALB2b) Family: beta ( logit )
Formula: Exposures ~ Age + Group + (1 | Solu_Sample) + (1 | Signature)
Zero inflation: ~1
Data: mut_sig_SIGNAL_exposures_PALB2b
AIC BIC logLik deviance df.resid
31932.4 31999.3 -15959.2 31918.4 104993
Random effects:
Conditional model:
Groups Name Variance Std.Dev.
Solu_Sample (Intercept) 0.04491 0.2119
Signature (Intercept) 0.17691 0.4206
Number of obs: 105000, groups: Solu_Sample, 15000; Signature, 7
Dispersion parameter for beta family (): 5.69
Conditional model:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.2644555 0.1626662 -13.921 < 2e-16 ***
Age 0.0128562 0.0005562 23.116 < 2e-16 ***
GroupPALB2b 0.0764734 0.0110615 6.913 4.73e-12 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Zero-inflation model:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.750185 0.006611 -113.5 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1summary(mut_sig_mixed_lmer_bin_LLGL2) Family: beta ( logit )
Formula: Exposures ~ Age + Group + (1 | Solu_Sample) + (1 | Signature)
Zero inflation: ~1
Data: mut_sig_SIGNAL_exposures_LLGL2
AIC BIC logLik deviance df.resid
60057.5 60124.4 -30021.7 60043.5 104993
Random effects:
Conditional model:
Groups Name Variance Std.Dev.
Solu_Sample (Intercept) 0.05634 0.2374
Signature (Intercept) 0.21470 0.4634
Number of obs: 105000, groups: Solu_Sample, 15000; Signature, 7
Dispersion parameter for beta family (): 4.08
Conditional model:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.0859905 0.1788675 -11.66 <2e-16 ***
Age 0.0112367 0.0005847 19.22 <2e-16 ***
GroupLLGL2b 0.1460556 0.0115641 12.63 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Zero-inflation model:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.531563 0.006391 -83.17 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1summary(mut_sig_mixed_lmer_bin_SCYL3) Family: beta ( logit )
Formula: Exposures ~ Age + Group + (1 | Solu_Sample) + (1 | Signature)
Zero inflation: ~1
Data: mut_sig_SIGNAL_exposures_SCYL3
AIC BIC logLik deviance df.resid
35547.5 35613.9 -17766.7 35533.5 97993
Random effects:
Conditional model:
Groups Name Variance Std.Dev.
Solu_Sample (Intercept) 0.04149 0.2037
Signature (Intercept) 0.21897 0.4679
Number of obs: 98000, groups: Solu_Sample, 14000; Signature, 7
Dispersion parameter for beta family (): 5.67
Conditional model:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.1432706 0.1797911 -11.921 <2e-16 ***
Age 0.0105856 0.0005209 20.322 <2e-16 ***
GroupSCYL3 0.0187499 0.0110561 1.696 0.0899 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Zero-inflation model:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.664613 0.006745 -98.54 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1# summary(mut_sig_mixed_lmer_bin_sample)Plots of simulated residuals
plot(simres_mut_sig_mixed_lmer_bin_ATM)
plot(simres_mut_sig_mixed_lmer_bin_PALB2a)
plot(simres_mut_sig_mixed_lmer_bin_PALB2b)
plot(simres_mut_sig_mixed_lmer_bin_LLGL2)
plot(simres_mut_sig_mixed_lmer_bin_SCYL3)
Summary and plots of dispersion
testDispersion(simres_mut_sig_mixed_lmer_bin_ATM)
DHARMa nonparametric dispersion test via sd of residuals fitted vs. simulated
data: simulationOutput
dispersion = 1.0582, p-value = 0.616
alternative hypothesis: two.sidedtestDispersion(simres_mut_sig_mixed_lmer_bin_PALB2a)
DHARMa nonparametric dispersion test via sd of residuals fitted vs. simulated
data: simulationOutput
dispersion = 1.0424, p-value = 0.744
alternative hypothesis: two.sided
testDispersion(simres_mut_sig_mixed_lmer_bin_PALB2b)
DHARMa nonparametric dispersion test via sd of residuals fitted vs. simulated
data: simulationOutput
dispersion = 1.0152, p-value = 0.84
alternative hypothesis: two.sided
testDispersion(simres_mut_sig_mixed_lmer_bin_LLGL2)
DHARMa nonparametric dispersion test via sd of residuals fitted vs. simulated
data: simulationOutput
dispersion = 1.0324, p-value = 0.816
alternative hypothesis: two.sided
testDispersion(simres_mut_sig_mixed_lmer_bin_SCYL3)
DHARMa nonparametric dispersion test via sd of residuals fitted vs. simulated
data: simulationOutput
dispersion = 1.0557, p-value = 0.792
alternative hypothesis: two.sided
# testDispersion(simres_mut_sig_mixed_lmer_bin_sample)Summary and plots of uniformity
testUniformity(simres_mut_sig_mixed_lmer_bin_ATM)
One-sample Kolmogorov-Smirnov test
data: simulationOutput$scaledResiduals
D = 0.048133, p-value < 2.2e-16
alternative hypothesis: two-sided
testUniformity(simres_mut_sig_mixed_lmer_bin_PALB2a)
One-sample Kolmogorov-Smirnov test
data: simulationOutput$scaledResiduals
D = 0.056776, p-value < 2.2e-16
alternative hypothesis: two-sided
testUniformity(simres_mut_sig_mixed_lmer_bin_PALB2b)
One-sample Kolmogorov-Smirnov test
data: simulationOutput$scaledResiduals
D = 0.0542, p-value < 2.2e-16
alternative hypothesis: two-sided
testUniformity(simres_mut_sig_mixed_lmer_bin_LLGL2)
One-sample Kolmogorov-Smirnov test
data: simulationOutput$scaledResiduals
D = 0.052019, p-value < 2.2e-16
alternative hypothesis: two-sided
testUniformity(simres_mut_sig_mixed_lmer_bin_SCYL3)
One-sample Kolmogorov-Smirnov test
data: simulationOutput$scaledResiduals
D = 0.054255, p-value < 2.2e-16
alternative hypothesis: two-sided
# testUniformity(simres_mut_sig_mixed_lmer_bin_sample)Summary and plots of outliers
testOutliers(simres_mut_sig_mixed_lmer_bin_ATM,type="bootstrap")
DHARMa bootstrapped outlier test
data: simres_mut_sig_mixed_lmer_bin_ATM
outliers at both margin(s) = 543, observations = 105000, p-value = 0.76
alternative hypothesis: two.sided
percent confidence interval:
0.002842143 0.014148333
sample estimates:
outlier frequency (expected: 0.00645638095238095 )
0.005171429 
testOutliers(simres_mut_sig_mixed_lmer_bin_PALB2a,type="bootstrap")
DHARMa bootstrapped outlier test
data: simres_mut_sig_mixed_lmer_bin_PALB2a
outliers at both margin(s) = 427, observations = 98000, p-value = 0.6
alternative hypothesis: two.sided
percent confidence interval:
0.002389541 0.016099745
sample estimates:
outlier frequency (expected: 0.00667224489795918 )
0.004357143 
testOutliers(simres_mut_sig_mixed_lmer_bin_PALB2b,type="bootstrap")
DHARMa bootstrapped outlier test
data: simres_mut_sig_mixed_lmer_bin_PALB2b
outliers at both margin(s) = 441, observations = 105000, p-value = 0.54
alternative hypothesis: two.sided
percent confidence interval:
0.002544524 0.015225476
sample estimates:
outlier frequency (expected: 0.00646295238095238 )
0.0042 
testOutliers(simres_mut_sig_mixed_lmer_bin_LLGL2,type="bootstrap")
DHARMa bootstrapped outlier test
data: simres_mut_sig_mixed_lmer_bin_LLGL2
outliers at both margin(s) = 1094, observations = 105000, p-value = 0.2
alternative hypothesis: two.sided
percent confidence interval:
0.002551429 0.012666667
sample estimates:
outlier frequency (expected: 0.00629619047619048 )
0.01041905 
testOutliers(simres_mut_sig_mixed_lmer_bin_SCYL3,type="bootstrap")
DHARMa bootstrapped outlier test
data: simres_mut_sig_mixed_lmer_bin_SCYL3
outliers at both margin(s) = 463, observations = 98000, p-value = 0.92
alternative hypothesis: two.sided
percent confidence interval:
0.002066071 0.012079082
sample estimates:
outlier frequency (expected: 0.00571051020408163 )
0.00472449 
# testOutliers(simres_mut_sig_mixed_lmer_bin_sample,type="bootstrap")Summary and plots of zero-inflation
testZeroInflation(simres_mut_sig_mixed_lmer_bin_ATM)
DHARMa zero-inflation test via comparison to expected zeros with simulation under H0 = fitted model
data: simulationOutput
ratioObsSim = 1.0004, p-value = 0.92
alternative hypothesis: two.sided
testZeroInflation(simres_mut_sig_mixed_lmer_bin_PALB2a)
DHARMa zero-inflation test via comparison to expected zeros with simulation under H0 = fitted model
data: simulationOutput
ratioObsSim = 1.0002, p-value = 0.952
alternative hypothesis: two.sided
testZeroInflation(simres_mut_sig_mixed_lmer_bin_PALB2b)
DHARMa zero-inflation test via comparison to expected zeros with simulation under H0 = fitted model
data: simulationOutput
ratioObsSim = 1.0002, p-value = 0.984
alternative hypothesis: two.sided
testZeroInflation(simres_mut_sig_mixed_lmer_bin_LLGL2)
DHARMa zero-inflation test via comparison to expected zeros with simulation under H0 = fitted model
data: simulationOutput
ratioObsSim = 1.0005, p-value = 0.888
alternative hypothesis: two.sided
testZeroInflation(simres_mut_sig_mixed_lmer_bin_SCYL3)
DHARMa zero-inflation test via comparison to expected zeros with simulation under H0 = fitted model
data: simulationOutput
ratioObsSim = 0.99979, p-value = 0.904
alternative hypothesis: two.sided
# testZeroInflation(simres_mut_sig_mixed_lmer_bin_sample)Plots of residuals against sample/group predictors
plotResiduals(simres_mut_sig_mixed_lmer_bin_ATM, form=mut_sig_SIGNAL_exposures_ATM$Group)
plotResiduals(simres_mut_sig_mixed_lmer_bin_PALB2a, form=mut_sig_SIGNAL_exposures_PALB2a$Group)
plotResiduals(simres_mut_sig_mixed_lmer_bin_PALB2b, form=mut_sig_SIGNAL_exposures_PALB2b$Group)
plotResiduals(simres_mut_sig_mixed_lmer_bin_LLGL2, form=mut_sig_SIGNAL_exposures_LLGL2$Group)
plotResiduals(simres_mut_sig_mixed_lmer_bin_SCYL3, form=mut_sig_SIGNAL_exposures_SCYL3$Group)
# plotResiduals(simres_mut_sig_mixed_lmer_bin_sample, form=mut_sig_SIGNAL_exposures_sample$Group)Save results
sink(file="pooled_vs_HGSrefMutSig_SIGNAL_lmer_bin_all_p_values.txt",append = FALSE)
print("HGS_refMutSig_vs_ATM_lmer_bin_results")
print(coef(summary(mut_sig_mixed_lmer_bin_ATM))$cond)
print("HGS_refMutSig_vs_PALB2a_lmer_bin_results")
print(coef(summary(mut_sig_mixed_lmer_bin_PALB2a))$cond)
print("HGS_refMutSig_vs_PALB2b_lmer_bin_results")
print(coef(summary(mut_sig_mixed_lmer_bin_PALB2b))$cond)
print("HGS_refMutSig_vs_LLGL2_lmer_bin_results")
print(coef(summary(mut_sig_mixed_lmer_bin_LLGL2))$cond)
print("HGS_refMutSig_vs_SCYL3_lmer_bin_results")
print(coef(summary(mut_sig_mixed_lmer_bin_SCYL3))$cond)
# print("HGS_refMutSig_vs_sample_lmer_bin_results")
# print(coef(summary(mut_sig_mixed_lmer_bin_sample))$cond)
sink()
pdf("HGS_refMutSig_vs_ATM_SIGNAL_lmer_bin_simRes_results.pdf")
sink(file="HGS_refMutSig_vs_ATM_SIGNAL_lmer_bin_simRes_results.txt",append = FALSE)
summary(mut_sig_mixed_lmer_bin_ATM)
plot(simres_mut_sig_mixed_lmer_bin_ATM)
testDispersion(simres_mut_sig_mixed_lmer_bin_ATM)testUniformity(simres_mut_sig_mixed_lmer_bin_ATM)
testOutliers(simres_mut_sig_mixed_lmer_bin_ATM,type="bootstrap")testZeroInflation(simres_mut_sig_mixed_lmer_bin_ATM)
plotResiduals(simres_mut_sig_mixed_lmer_bin_ATM, mut_sig_SIGNAL_exposures_ATM$Group)dev.off()
sink()
pdf("HGS_refMutSig_vs_PALB2a_SIGNAL_lmer_bin_simRes_results.pdf")
sink(file="HGS_refMutSig_vs_PALB2a_SIGNAL_lmer_bin_simRes_results.txt",append = FALSE)
summary(mut_sig_mixed_lmer_bin_PALB2a)
plot(simres_mut_sig_mixed_lmer_bin_PALB2a)testDispersion(simres_mut_sig_mixed_lmer_bin_PALB2a)
testUniformity(simres_mut_sig_mixed_lmer_bin_PALB2a)testOutliers(simres_mut_sig_mixed_lmer_bin_PALB2a,type="bootstrap")
testZeroInflation(simres_mut_sig_mixed_lmer_bin_PALB2a)plotResiduals(simres_mut_sig_mixed_lmer_bin_PALB2a, mut_sig_SIGNAL_exposures_PALB2a$Group)
dev.off()
sink()
pdf("HGS_refMutSig_vs_PALB2b_SIGNAL_lmer_bin_simRes_results.pdf")sink(file="HGS_refMutSig_vs_PALB2b_SIGNAL_lmer_bin_simRes_results.txt",append = FALSE)
summary(mut_sig_mixed_lmer_bin_PALB2b)
plot(simres_mut_sig_mixed_lmer_bin_PALB2b)
testDispersion(simres_mut_sig_mixed_lmer_bin_PALB2b)testUniformity(simres_mut_sig_mixed_lmer_bin_PALB2b)
testOutliers(simres_mut_sig_mixed_lmer_bin_PALB2b,type="bootstrap")testZeroInflation(simres_mut_sig_mixed_lmer_bin_PALB2b)
plotResiduals(simres_mut_sig_mixed_lmer_bin_PALB2b, mut_sig_SIGNAL_exposures_PALB2b$Group)dev.off()
sink()
pdf("HGS_refMutSig_vs_LLGL2_SIGNAL_lmer_bin_simRes_results.pdf")
sink(file="HGS_refMutSig_vs_LLGL2_SIGNAL_lmer_bin_simRes_results.txt",append = FALSE)
summary(mut_sig_mixed_lmer_bin_LLGL2)
plot(simres_mut_sig_mixed_lmer_bin_LLGL2)testDispersion(simres_mut_sig_mixed_lmer_bin_LLGL2)
testUniformity(simres_mut_sig_mixed_lmer_bin_LLGL2)testOutliers(simres_mut_sig_mixed_lmer_bin_LLGL2,type="bootstrap")
testZeroInflation(simres_mut_sig_mixed_lmer_bin_LLGL2)plotResiduals(simres_mut_sig_mixed_lmer_bin_LLGL2, mut_sig_SIGNAL_exposures_LLGL2$Group)
dev.off()
sink()
pdf("HGS_refMutSig_vs_SCYL3_SIGNAL_lmer_bin_simRes_results.pdf")sink(file="HGS_refMutSig_vs_SCYL3_SIGNAL_lmer_bin_simRes_results.txt",append = FALSE)
summary(mut_sig_mixed_lmer_bin_SCYL3)
plot(simres_mut_sig_mixed_lmer_bin_SCYL3)
testDispersion(simres_mut_sig_mixed_lmer_bin_SCYL3)testUniformity(simres_mut_sig_mixed_lmer_bin_SCYL3)
testOutliers(simres_mut_sig_mixed_lmer_bin_SCYL3,type="bootstrap")testZeroInflation(simres_mut_sig_mixed_lmer_bin_SCYL3)
plotResiduals(simres_mut_sig_mixed_lmer_bin_SCYL3, mut_sig_SIGNAL_exposures_SCYL3$Group)dev.off()
sink()
# pdf("HGS_refMutSig_vs_sample_SIGNAL_lmer_bin_simRes_results.pdf")
# sink(file="HGS_refMutSig_vs_sample_SIGNAL_lmer_bin_simRes_results.txt",append = FALSE)
# summary(mut_sig_mixed_lmer_bin_sample)
# plot(simres_mut_sig_mixed_lmer_bin_sample)
# testDispersion(simres_mut_sig_mixed_lmer_bin_sample)
# testUniformity(simres_mut_sig_mixed_lmer_bin_sample)
# testOutliers(simres_mut_sig_mixed_lmer_bin_sample,type="bootstrap")
# testZeroInflation(simres_mut_sig_mixed_lmer_bin_sample)
# plotResiduals(simres_mut_sig_mixed_lmer_bin_sample, mut_sig_SIGNAL_exposures_sample$Group)
# dev.off()
# sink()
closeAllConnections()LS0tCnRpdGxlOiAiVGhlc2lzIFNJR05BTCBNdXRhdGlvbmFsIFNpZ25hdHVyZXMgR0xNTSBDb21wYXJpc29uIFNjcmlwdCIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKQUxMIElNUE9SVEVEIEZJTEVTIE1VU1QgQkUgSU4gVEhFIFNBTUUgRElSRUNUT1JZIEFTIFRISVMgU0NSSVBUCgpIaWdoLXBlcmZvcm1hbmNlIGNvbXB1dGluZyBjbHVzdGVyIChlLmcuIHVzaW5nIFJTdHVkaW8gU2VydmVyKSBzdHJvbmdseSByZWNvbW1lbmRlZCwgZHVlIHRvIGxhcmdlIG1lbW9yeSByZXF1aXJlbWVudHMKCkxvYWQgcmVxdWlyZWQgcGFja2FnZXMKYGBge3J9CmxpYnJhcnkoInRpZHl2ZXJzZSIpCmxpYnJhcnkoInJlYWRidWxrIikKbGlicmFyeSgiZ2xtbVRNQiIpCmxpYnJhcnkoImJldGFyZWciKQpsaWJyYXJ5KCJESEFSTWEiKQpgYGAKCkltcG9ydCBkYXRhOyBlYWNoIGZpbGUgY29tcHJpc2VzIDEwMDAgU0lHTkFMIGJvb3RzdHJhcCBzb2x1dGlvbnMgKGdlbmVyYXRlZCB1c2luZyBodHRwczovL3NpZ25hbC5tdXRhdGlvbmFsc2lnbmF0dXJlcy5jb20vYW5hbHlzZTIpIHVzaW5nIHRoZSBPdkNhLXNwZWNpZmljIHNpZ25hdHVyZSBzZXQgKE9WQVJZX0EvQi9DL0QvRS9GL0ctIHNlZSBodHRwczovL3NpZ25hbC5tdXRhdGlvbmFsc2lnbmF0dXJlcy5jb20vZXhwbG9yZS9zdHVkeVRpc3N1ZVR5cGUvMS0xNSkgZm9yIHRoZSAxMiBzYW1wbGVzIGNvbXByaXNpbmcgdGhlIHBvb2xlZCByZWZlcmVuY2UgZ3JvdXAgKCdIR1NyZWZNdXRTaWcnKSwgYWxvbmcgd2l0aCBlaXRoZXIgdGhlIDItMyBzYW1wbGVzIGNvbXByaXNpbmcgZWFjaCBwb29sZWQgdGVzdCBncm91cCAoJ0FUTS9QQUxCMmEvUEFMQjJiL0xMR0wyL1NDWUwzJykgb3IgYW4gaW5kaXZpZHVhbCBzYW1wbGUgKCdzYW1wbGUnKTsgcmVmZXIgdG8gQ2hhcHRlcnMgNC4yLjIuMy00IGZvciBmdXJ0aGVyIGRldGFpbHMKYGBge3J9Cm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19BVE0gPC0gcmVhZF90c3YoIm1hc3RlcmZpbGVfdHVtb3VyX3NhbXBsZXNfU0lHTkFMX0hHU3JlZk11dFNpZ19BVE0udHN2IikKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYSA8LSByZWFkX3RzdigibWFzdGVyZmlsZV90dW1vdXJfc2FtcGxlc19TSUdOQUxfSEdTcmVmTXV0U2lnX1BBTEIyYS50c3YiKQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJiIDwtIHJlYWRfdHN2KCJtYXN0ZXJmaWxlX3R1bW91cl9zYW1wbGVzX1NJR05BTF9IR1NyZWZNdXRTaWdfUEFMQjJiLnRzdiIpCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19MTEdMMiA8LSByZWFkX3RzdigibWFzdGVyZmlsZV90dW1vdXJfc2FtcGxlc19TSUdOQUxfSEdTcmVmTXV0U2lnX0xMR0wyLnRzdiIpCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19TQ1lMMyA8LSByZWFkX3RzdigibWFzdGVyZmlsZV90dW1vdXJfc2FtcGxlc19TSUdOQUxfSEdTcmVmTXV0U2lnX1NDWUwzLnRzdiIpCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19zYW1wbGUgPC0gcmVhZF90c3YoIm1hc3RlcmZpbGVfdHVtb3VyX3NhbXBsZXNfU0lHTkFMX0hHU3JlZk11dFNpZ19zYW1wbGUudHN2IikKYGBgCgpSb3VuZCBkb3duIDFzIGFuZCBjb252ZXJ0IGRhdGEgaW50byBmYWN0b3JzCmBgYHtyfQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfQVRNJEV4cG9zdXJlcyA8LSByZXBsYWNlKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19BVE0kRXhwb3N1cmVzLCBtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfQVRNJEV4cG9zdXJlcz09MSwgMC45OTk5OTk5KQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfQVRNJEdyb3VwIDwtIGFzLmZhY3RvcihtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfQVRNJEdyb3VwKQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfQVRNJFR1bW91cl9TYW1wbGUgPC0gYXMuZmFjdG9yKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19BVE0kVHVtb3VyX1NhbXBsZSkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX0FUTSRTb2x1dGlvbiA8LSBhcy5mYWN0b3IobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX0FUTSRTb2x1dGlvbikKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX0FUTSRTb2x1X1NhbXBsZSA8LSBwYXN0ZShtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfQVRNJFNvbHV0aW9uLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19BVE0kVHVtb3VyX1NhbXBsZSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzZXA9IjoiKSAlPiUgCiAgYXMuZmFjdG9yKCkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX0FUTSRTaWduYXR1cmUgPC0gYXMuZmFjdG9yKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19BVE0kU2lnbmF0dXJlKQoKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYSRFeHBvc3VyZXMgPC0gcmVwbGFjZShtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJhJEV4cG9zdXJlcywgbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYSRFeHBvc3VyZXM9PTEsIDAuOTk5OTk5OSkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYSRHcm91cCA8LSBhcy5mYWN0b3IobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYSRHcm91cCkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYSRUdW1vdXJfU2FtcGxlIDwtIGFzLmZhY3RvcihtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJhJFR1bW91cl9TYW1wbGUpCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19QQUxCMmEkU29sdXRpb24gPC0gYXMuZmFjdG9yKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19QQUxCMmEkU29sdXRpb24pCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19QQUxCMmEkU29sdV9TYW1wbGUgPC0gcGFzdGUobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYSRTb2x1dGlvbiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19QQUxCMmEkVHVtb3VyX1NhbXBsZSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHNlcD0iOiIpICU+JSAKICBhcy5mYWN0b3IoKQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJhJFNpZ25hdHVyZSA8LSBhcy5mYWN0b3IobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYSRTaWduYXR1cmUpCgptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJiJEV4cG9zdXJlcyA8LSByZXBsYWNlKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19QQUxCMmIkRXhwb3N1cmVzLCBtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJiJEV4cG9zdXJlcz09MSwgMC45OTk5OTk5KQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJiJEdyb3VwIDwtIGFzLmZhY3RvcihtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJiJEdyb3VwKQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJiJFR1bW91cl9TYW1wbGUgPC0gYXMuZmFjdG9yKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19QQUxCMmIkVHVtb3VyX1NhbXBsZSkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYiRTb2x1dGlvbiA8LSBhcy5mYWN0b3IobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYiRTb2x1dGlvbikKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYiRTb2x1X1NhbXBsZSA8LSBwYXN0ZShtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJiJFNvbHV0aW9uLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19QQUxCMmIkVHVtb3VyX1NhbXBsZSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzZXA9IjoiKSAlPiUgCiAgYXMuZmFjdG9yKCkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYiRTaWduYXR1cmUgPC0gYXMuZmFjdG9yKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19QQUxCMmIkU2lnbmF0dXJlKQoKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX0xMR0wyJEV4cG9zdXJlcyA8LSByZXBsYWNlKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19MTEdMMiRFeHBvc3VyZXMsIG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19MTEdMMiRFeHBvc3VyZXM9PTEsIDAuOTk5OTk5OSkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX0xMR0wyJEdyb3VwIDwtIGFzLmZhY3RvcihtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfTExHTDIkR3JvdXApCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19MTEdMMiRUdW1vdXJfU2FtcGxlIDwtIGFzLmZhY3RvcihtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfTExHTDIkVHVtb3VyX1NhbXBsZSkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX0xMR0wyJFNvbHV0aW9uIDwtIGFzLmZhY3RvcihtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfTExHTDIkU29sdXRpb24pCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19MTEdMMiRTb2x1X1NhbXBsZSA8LSBwYXN0ZShtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfTExHTDIkU29sdXRpb24sCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX0xMR0wyJFR1bW91cl9TYW1wbGUsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc2VwPSI6IikgJT4lIAogIGFzLmZhY3RvcigpCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19MTEdMMiRTaWduYXR1cmUgPC0gYXMuZmFjdG9yKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19MTEdMMiRTaWduYXR1cmUpCgptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkRXhwb3N1cmVzIDwtIHJlcGxhY2UobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1NDWUwzJEV4cG9zdXJlcywgbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1NDWUwzJEV4cG9zdXJlcz09MSwgMC45OTk5OTk5KQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkR3JvdXAgPC0gYXMuZmFjdG9yKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19TQ1lMMyRHcm91cCkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1NDWUwzJFR1bW91cl9TYW1wbGUgPC0gYXMuZmFjdG9yKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19TQ1lMMyRUdW1vdXJfU2FtcGxlKQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkU29sdXRpb24gPC0gYXMuZmFjdG9yKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19TQ1lMMyRTb2x1dGlvbikKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1NDWUwzJFNvbHVfU2FtcGxlIDwtIHBhc3RlKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19TQ1lMMyRTb2x1dGlvbiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkVHVtb3VyX1NhbXBsZSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzZXA9IjoiKSAlPiUgCiAgYXMuZmFjdG9yKCkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1NDWUwzJFNpZ25hdHVyZSA8LSBhcy5mYWN0b3IobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1NDWUwzJFNpZ25hdHVyZSkKCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19TQ1lMMyRFeHBvc3VyZXMgPC0gcmVwbGFjZShtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkRXhwb3N1cmVzLCBtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkRXhwb3N1cmVzPT0xLCAwLjk5OTk5OTkpCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19TQ1lMMyRHcm91cCA8LSBhcy5mYWN0b3IobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1NDWUwzJEdyb3VwKQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkVHVtb3VyX1NhbXBsZSA8LSBhcy5mYWN0b3IobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1NDWUwzJFR1bW91cl9TYW1wbGUpCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19TQ1lMMyRTb2x1dGlvbiA8LSBhcy5mYWN0b3IobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1NDWUwzJFNvbHV0aW9uKQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkU29sdV9TYW1wbGUgPC0gcGFzdGUobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1NDWUwzJFNvbHV0aW9uLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19TQ1lMMyRUdW1vdXJfU2FtcGxlLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHNlcD0iOiIpICU+JSAKICBhcy5mYWN0b3IoKQptdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkU2lnbmF0dXJlIDwtIGFzLmZhY3RvcihtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkU2lnbmF0dXJlKQoKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX3NhbXBsZSRFeHBvc3VyZXMgPC0gcmVwbGFjZShtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfc2FtcGxlJEV4cG9zdXJlcywgbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX3NhbXBsZSRFeHBvc3VyZXM9PTEsIDAuOTk5OTk5OSkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX3NhbXBsZSRHcm91cCA8LSBhcy5mYWN0b3IobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX3NhbXBsZSRHcm91cCkKbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX3NhbXBsZSRUdW1vdXJfU2FtcGxlIDwtIGFzLmZhY3RvcihtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfc2FtcGxlJFR1bW91cl9TYW1wbGUpCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19zYW1wbGUkU29sdXRpb24gPC0gYXMuZmFjdG9yKG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19zYW1wbGUkU29sdXRpb24pCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19zYW1wbGUkU29sdV9TYW1wbGUgPC0gcGFzdGUobXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX3NhbXBsZSRTb2x1dGlvbiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfc2FtcGxlJFR1bW91cl9TYW1wbGUsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgc2VwPSI6IikgJT4lIAogIGFzLmZhY3RvcigpCm11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19zYW1wbGUkU2lnbmF0dXJlIDwtIGFzLmZhY3RvcihtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfc2FtcGxlJFNpZ25hdHVyZSkKYGBgCgpSdW4gR0xNTSB3aXRoIGJldGEtcmVncmVzc2lvbiBmb3IgY29tYmluZWQgSEdTcmVmTXV0U2lnIHZzIGNvbWJpbmVkIGFuZCBpbmRpdmlkdWFsIHNhbXBsZXMgb2YgaW50ZXJlc3QsIGFuZCBzaW11bGF0ZSByZXNpZHVhbHMgZm9yIG1vZGVsIHVzaW5nIERIQVJNYQpgYGB7cn0KbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9BVE0gPC0gZ2xtbVRNQihFeHBvc3VyZXMgfiBBZ2UgKyBHcm91cCArIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAoMXxTb2x1X1NhbXBsZSkgKyAoMXxTaWduYXR1cmUpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBkYXRhID0gbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX0FUTSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZmFtaWx5PWJldGFfZmFtaWx5KCksIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB6aWZvcm11bGEgPSB+MSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29udHJvbCA9IGdsbW1UTUJDb250cm9sKHBhcmFsbGVsID0gNDgpICNwYXJhbGxlbCA9IG5vLiBvZiBhdmFpbGFibGUgcHJvY2Vzc29yIGNvcmVzL3RocmVhZHMKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgKQpzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9BVE0gPC0gc2ltdWxhdGVSZXNpZHVhbHMobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9BVE0pCgptdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYSA8LSBnbG1tVE1CKEV4cG9zdXJlcyB+IEFnZSArIEdyb3VwICsgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICgxfFNvbHVfU2FtcGxlKSArICgxfFNpZ25hdHVyZSksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGRhdGEgPSBtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJhLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBmYW1pbHk9YmV0YV9mYW1pbHkoKSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHppZm9ybXVsYSA9IH4xLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb250cm9sID0gZ2xtbVRNQkNvbnRyb2wocGFyYWxsZWwgPSA0OCkKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgKQpzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmEgPC0gc2ltdWxhdGVSZXNpZHVhbHMobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmEpCgptdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYiA8LSBnbG1tVE1CKEV4cG9zdXJlcyB+IEFnZSArIEdyb3VwICsgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICgxfFNvbHVfU2FtcGxlKSArICgxfFNpZ25hdHVyZSksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGRhdGEgPSBtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJiLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZmFtaWx5PWJldGFfZmFtaWx5KCksIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB6aWZvcm11bGEgPSB+MSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29udHJvbCA9IGdsbW1UTUJDb250cm9sKHBhcmFsbGVsID0gNDgpCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICkKc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJiIDwtIHNpbXVsYXRlUmVzaWR1YWxzKG11dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJiKQoKbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9MTEdMMiA8LSBnbG1tVE1CKEV4cG9zdXJlcyB+IEFnZSArIEdyb3VwICsKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAoMXxTb2x1X1NhbXBsZSkgKyAoMXxTaWduYXR1cmUpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGRhdGEgPSBtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfTExHTDIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGZhbWlseT1iZXRhX2ZhbWlseSgpLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB6aWZvcm11bGEgPSB+MSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb250cm9sID0gZ2xtbVRNQkNvbnRyb2wocGFyYWxsZWwgPSA0OCkKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICApCnNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX0xMR0wyIDwtIHNpbXVsYXRlUmVzaWR1YWxzKG11dF9zaWdfbWl4ZWRfbG1lcl9iaW5fTExHTDIpCgptdXRfc2lnX21peGVkX2xtZXJfYmluX1NDWUwzIDwtIGdsbW1UTUIoRXhwb3N1cmVzIH4gQWdlICsgR3JvdXAgKyAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAoMXxTb2x1X1NhbXBsZSkgKyAoMXxTaWduYXR1cmUpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGRhdGEgPSBtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZmFtaWx5PWJldGFfZmFtaWx5KCksIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHppZm9ybXVsYSA9IH4xLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbnRyb2wgPSBnbG1tVE1CQ29udHJvbChwYXJhbGxlbCA9IDQ4KQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICkKc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fU0NZTDMgPC0gc2ltdWxhdGVSZXNpZHVhbHMobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9TQ1lMMykKCm11dF9zaWdfbWl4ZWRfbG1lcl9iaW5fc2FtcGxlIDwtIGdsbW1UTUIoRXhwb3N1cmVzIH4gQWdlICsgR3JvdXAgKyAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgKDF8U29sdV9TYW1wbGUpICsgKDF8U2lnbmF0dXJlKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZGF0YSA9IG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19zYW1wbGUsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGZhbWlseT1iZXRhX2ZhbWlseSgpLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgemlmb3JtdWxhID0gfjEsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbnRyb2wgPSBnbG1tVE1CQ29udHJvbChwYXJhbGxlbCA9IDQ4KQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICApCnNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX3NhbXBsZSA8LSBzaW11bGF0ZVJlc2lkdWFscyhtdXRfc2lnX21peGVkX2xtZXJfYmluX3NhbXBsZSkKYGBgCgpTdW1tYXJ5IG9mIEdMTU0gcmVzdWx0cwpgYGB7cn0Kc3VtbWFyeShtdXRfc2lnX21peGVkX2xtZXJfYmluX0FUTSkKYGBgCmBgYHtyfQpzdW1tYXJ5KG11dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJhKQpgYGAKYGBge3J9CnN1bW1hcnkobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmIpCmBgYApgYGB7cn0Kc3VtbWFyeShtdXRfc2lnX21peGVkX2xtZXJfYmluX0xMR0wyKQpgYGAKYGBge3J9CnN1bW1hcnkobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9TQ1lMMykKYGBgCmBgYHtyfQojIHN1bW1hcnkobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9zYW1wbGUpCmBgYAoKUGxvdHMgb2Ygc2ltdWxhdGVkIHJlc2lkdWFscwpgYGB7cn0KcGxvdChzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9BVE0pCnBsb3Qoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJhKQpwbG90KHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYikKcGxvdChzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9MTEdMMikKcGxvdChzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9TQ1lMMykKIyBwbG90KHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX3NhbXBsZSkKYGBgCgpTdW1tYXJ5IGFuZCBwbG90cyBvZiBkaXNwZXJzaW9uCmBgYHtyfQp0ZXN0RGlzcGVyc2lvbihzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9BVE0pCmBgYApgYGB7cn0KdGVzdERpc3BlcnNpb24oc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJhKQpgYGAKYGBge3J9CnRlc3REaXNwZXJzaW9uKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYikKYGBgCmBgYHtyfQp0ZXN0RGlzcGVyc2lvbihzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9MTEdMMikKYGBgCmBgYHtyfQp0ZXN0RGlzcGVyc2lvbihzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9TQ1lMMykKYGBgCmBgYHtyfQojIHRlc3REaXNwZXJzaW9uKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX3NhbXBsZSkKYGBgCgpTdW1tYXJ5IGFuZCBwbG90cyBvZiB1bmlmb3JtaXR5CmBgYHtyfQp0ZXN0VW5pZm9ybWl0eShzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9BVE0pCmBgYApgYGB7cn0KdGVzdFVuaWZvcm1pdHkoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJhKQpgYGAKYGBge3J9CnRlc3RVbmlmb3JtaXR5KHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYikKYGBgCmBgYHtyfQp0ZXN0VW5pZm9ybWl0eShzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9MTEdMMikKYGBgCmBgYHtyfQp0ZXN0VW5pZm9ybWl0eShzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9TQ1lMMykKYGBgCmBgYHtyfQojIHRlc3RVbmlmb3JtaXR5KHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX3NhbXBsZSkKYGBgCgpTdW1tYXJ5IGFuZCBwbG90cyBvZiBvdXRsaWVycwpgYGB7cn0KdGVzdE91dGxpZXJzKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX0FUTSx0eXBlPSJib290c3RyYXAiKQpgYGAKYGBge3J9CnRlc3RPdXRsaWVycyhzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmEsdHlwZT0iYm9vdHN0cmFwIikKYGBgCmBgYHtyfQp0ZXN0T3V0bGllcnMoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJiLHR5cGU9ImJvb3RzdHJhcCIpCmBgYApgYGB7cn0KdGVzdE91dGxpZXJzKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX0xMR0wyLHR5cGU9ImJvb3RzdHJhcCIpCmBgYApgYGB7cn0KdGVzdE91dGxpZXJzKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1NDWUwzLHR5cGU9ImJvb3RzdHJhcCIpCmBgYApgYGB7cn0KIyB0ZXN0T3V0bGllcnMoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fc2FtcGxlLHR5cGU9ImJvb3RzdHJhcCIpCmBgYAoKU3VtbWFyeSBhbmQgcGxvdHMgb2YgemVyby1pbmZsYXRpb24KYGBge3J9CnRlc3RaZXJvSW5mbGF0aW9uKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX0FUTSkKYGBgCmBgYHtyfQp0ZXN0WmVyb0luZmxhdGlvbihzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmEpCmBgYApgYGB7cn0KdGVzdFplcm9JbmZsYXRpb24oc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJiKQpgYGAKYGBge3J9CnRlc3RaZXJvSW5mbGF0aW9uKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX0xMR0wyKQpgYGAKYGBge3J9CnRlc3RaZXJvSW5mbGF0aW9uKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1NDWUwzKQpgYGAKYGBge3J9CiMgdGVzdFplcm9JbmZsYXRpb24oc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fc2FtcGxlKQpgYGAKClBsb3RzIG9mIHJlc2lkdWFscyBhZ2FpbnN0IHNhbXBsZS9ncm91cCBwcmVkaWN0b3JzCmBgYHtyfQpwbG90UmVzaWR1YWxzKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX0FUTSwgZm9ybT1tdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfQVRNJEdyb3VwKQpwbG90UmVzaWR1YWxzKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYSwgZm9ybT1tdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJhJEdyb3VwKQpwbG90UmVzaWR1YWxzKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYiwgZm9ybT1tdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJiJEdyb3VwKQpwbG90UmVzaWR1YWxzKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX0xMR0wyLCBmb3JtPW11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19MTEdMMiRHcm91cCkKcGxvdFJlc2lkdWFscyhzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9TQ1lMMywgZm9ybT1tdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfU0NZTDMkR3JvdXApCiMgcGxvdFJlc2lkdWFscyhzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9zYW1wbGUsIGZvcm09bXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX3NhbXBsZSRHcm91cCkKYGBgCgpTYXZlIHJlc3VsdHMKYGBge3J9CnNpbmsoZmlsZT0icG9vbGVkX3ZzX0hHU3JlZk11dFNpZ19TSUdOQUxfbG1lcl9iaW5fYWxsX3BfdmFsdWVzLnR4dCIsYXBwZW5kID0gRkFMU0UpCnByaW50KCJIR1NfcmVmTXV0U2lnX3ZzX0FUTV9sbWVyX2Jpbl9yZXN1bHRzIikKcHJpbnQoY29lZihzdW1tYXJ5KG11dF9zaWdfbWl4ZWRfbG1lcl9iaW5fQVRNKSkkY29uZCkKcHJpbnQoIkhHU19yZWZNdXRTaWdfdnNfUEFMQjJhX2xtZXJfYmluX3Jlc3VsdHMiKQpwcmludChjb2VmKHN1bW1hcnkobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmEpKSRjb25kKQpwcmludCgiSEdTX3JlZk11dFNpZ192c19QQUxCMmJfbG1lcl9iaW5fcmVzdWx0cyIpCnByaW50KGNvZWYoc3VtbWFyeShtdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYikpJGNvbmQpCnByaW50KCJIR1NfcmVmTXV0U2lnX3ZzX0xMR0wyX2xtZXJfYmluX3Jlc3VsdHMiKQpwcmludChjb2VmKHN1bW1hcnkobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9MTEdMMikpJGNvbmQpCnByaW50KCJIR1NfcmVmTXV0U2lnX3ZzX1NDWUwzX2xtZXJfYmluX3Jlc3VsdHMiKQpwcmludChjb2VmKHN1bW1hcnkobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9TQ1lMMykpJGNvbmQpCiMgcHJpbnQoIkhHU19yZWZNdXRTaWdfdnNfc2FtcGxlX2xtZXJfYmluX3Jlc3VsdHMiKQojIHByaW50KGNvZWYoc3VtbWFyeShtdXRfc2lnX21peGVkX2xtZXJfYmluX3NhbXBsZSkpJGNvbmQpCnNpbmsoKQoKcGRmKCJIR1NfcmVmTXV0U2lnX3ZzX0FUTV9TSUdOQUxfbG1lcl9iaW5fc2ltUmVzX3Jlc3VsdHMucGRmIikKc2luayhmaWxlPSJIR1NfcmVmTXV0U2lnX3ZzX0FUTV9TSUdOQUxfbG1lcl9iaW5fc2ltUmVzX3Jlc3VsdHMudHh0IixhcHBlbmQgPSBGQUxTRSkKc3VtbWFyeShtdXRfc2lnX21peGVkX2xtZXJfYmluX0FUTSkKcGxvdChzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9BVE0pCnRlc3REaXNwZXJzaW9uKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX0FUTSkKdGVzdFVuaWZvcm1pdHkoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fQVRNKQp0ZXN0T3V0bGllcnMoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fQVRNLHR5cGU9ImJvb3RzdHJhcCIpCnRlc3RaZXJvSW5mbGF0aW9uKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX0FUTSkKcGxvdFJlc2lkdWFscyhzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9BVE0sIG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19BVE0kR3JvdXApCmRldi5vZmYoKQpzaW5rKCkKCnBkZigiSEdTX3JlZk11dFNpZ192c19QQUxCMmFfU0lHTkFMX2xtZXJfYmluX3NpbVJlc19yZXN1bHRzLnBkZiIpCnNpbmsoZmlsZT0iSEdTX3JlZk11dFNpZ192c19QQUxCMmFfU0lHTkFMX2xtZXJfYmluX3NpbVJlc19yZXN1bHRzLnR4dCIsYXBwZW5kID0gRkFMU0UpCnN1bW1hcnkobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmEpCnBsb3Qoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJhKQp0ZXN0RGlzcGVyc2lvbihzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmEpCnRlc3RVbmlmb3JtaXR5KHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYSkKdGVzdE91dGxpZXJzKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYSx0eXBlPSJib290c3RyYXAiKQp0ZXN0WmVyb0luZmxhdGlvbihzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmEpCnBsb3RSZXNpZHVhbHMoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJhLCBtdXRfc2lnX1NJR05BTF9leHBvc3VyZXNfUEFMQjJhJEdyb3VwKQpkZXYub2ZmKCkKc2luaygpCgpwZGYoIkhHU19yZWZNdXRTaWdfdnNfUEFMQjJiX1NJR05BTF9sbWVyX2Jpbl9zaW1SZXNfcmVzdWx0cy5wZGYiKQpzaW5rKGZpbGU9IkhHU19yZWZNdXRTaWdfdnNfUEFMQjJiX1NJR05BTF9sbWVyX2Jpbl9zaW1SZXNfcmVzdWx0cy50eHQiLGFwcGVuZCA9IEZBTFNFKQpzdW1tYXJ5KG11dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJiKQpwbG90KHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYikKdGVzdERpc3BlcnNpb24oc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJiKQp0ZXN0VW5pZm9ybWl0eShzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmIpCnRlc3RPdXRsaWVycyhzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9QQUxCMmIsdHlwZT0iYm9vdHN0cmFwIikKdGVzdFplcm9JbmZsYXRpb24oc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fUEFMQjJiKQpwbG90UmVzaWR1YWxzKHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1BBTEIyYiwgbXV0X3NpZ19TSUdOQUxfZXhwb3N1cmVzX1BBTEIyYiRHcm91cCkKZGV2Lm9mZigpCnNpbmsoKQoKcGRmKCJIR1NfcmVmTXV0U2lnX3ZzX0xMR0wyX1NJR05BTF9sbWVyX2Jpbl9zaW1SZXNfcmVzdWx0cy5wZGYiKQpzaW5rKGZpbGU9IkhHU19yZWZNdXRTaWdfdnNfTExHTDJfU0lHTkFMX2xtZXJfYmluX3NpbVJlc19yZXN1bHRzLnR4dCIsYXBwZW5kID0gRkFMU0UpCnN1bW1hcnkobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9MTEdMMikKcGxvdChzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9MTEdMMikKdGVzdERpc3BlcnNpb24oc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fTExHTDIpCnRlc3RVbmlmb3JtaXR5KHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX0xMR0wyKQp0ZXN0T3V0bGllcnMoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fTExHTDIsdHlwZT0iYm9vdHN0cmFwIikKdGVzdFplcm9JbmZsYXRpb24oc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fTExHTDIpCnBsb3RSZXNpZHVhbHMoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fTExHTDIsIG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19MTEdMMiRHcm91cCkKZGV2Lm9mZigpCnNpbmsoKQoKcGRmKCJIR1NfcmVmTXV0U2lnX3ZzX1NDWUwzX1NJR05BTF9sbWVyX2Jpbl9zaW1SZXNfcmVzdWx0cy5wZGYiKQpzaW5rKGZpbGU9IkhHU19yZWZNdXRTaWdfdnNfU0NZTDNfU0lHTkFMX2xtZXJfYmluX3NpbVJlc19yZXN1bHRzLnR4dCIsYXBwZW5kID0gRkFMU0UpCnN1bW1hcnkobXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9TQ1lMMykKcGxvdChzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9TQ1lMMykKdGVzdERpc3BlcnNpb24oc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fU0NZTDMpCnRlc3RVbmlmb3JtaXR5KHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX1NDWUwzKQp0ZXN0T3V0bGllcnMoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fU0NZTDMsdHlwZT0iYm9vdHN0cmFwIikKdGVzdFplcm9JbmZsYXRpb24oc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fU0NZTDMpCnBsb3RSZXNpZHVhbHMoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fU0NZTDMsIG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19TQ1lMMyRHcm91cCkKZGV2Lm9mZigpCnNpbmsoKQoKIyBwZGYoIkhHU19yZWZNdXRTaWdfdnNfc2FtcGxlX1NJR05BTF9sbWVyX2Jpbl9zaW1SZXNfcmVzdWx0cy5wZGYiKQojIHNpbmsoZmlsZT0iSEdTX3JlZk11dFNpZ192c19zYW1wbGVfU0lHTkFMX2xtZXJfYmluX3NpbVJlc19yZXN1bHRzLnR4dCIsYXBwZW5kID0gRkFMU0UpCiMgc3VtbWFyeShtdXRfc2lnX21peGVkX2xtZXJfYmluX3NhbXBsZSkKIyBwbG90KHNpbXJlc19tdXRfc2lnX21peGVkX2xtZXJfYmluX3NhbXBsZSkKIyB0ZXN0RGlzcGVyc2lvbihzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9zYW1wbGUpCiMgdGVzdFVuaWZvcm1pdHkoc2ltcmVzX211dF9zaWdfbWl4ZWRfbG1lcl9iaW5fc2FtcGxlKQojIHRlc3RPdXRsaWVycyhzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9zYW1wbGUsdHlwZT0iYm9vdHN0cmFwIikKIyB0ZXN0WmVyb0luZmxhdGlvbihzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9zYW1wbGUpCiMgcGxvdFJlc2lkdWFscyhzaW1yZXNfbXV0X3NpZ19taXhlZF9sbWVyX2Jpbl9zYW1wbGUsIG11dF9zaWdfU0lHTkFMX2V4cG9zdXJlc19zYW1wbGUkR3JvdXApCiMgZGV2Lm9mZigpCiMgc2luaygpCgpjbG9zZUFsbENvbm5lY3Rpb25zKCkKYGBgCg==