Warning message:
In readChar(file, size, TRUE) : truncating string with embedded nuls
    p_values[lower.tri(p_values,diag=TRUE)] <- NA

    p_values2 = as.data.frame(as.table(as.matrix(p_values))) 
    
    if(nrow(p_values2) > 10000){
      p_values2 = p_values2[sample(nrow(p_values2), (nrow(p_values2)/20)) , ]
    }
    p_values2 = p_values2[complete.cases(p_values2),]
    p_values2 = p_values2[order(p_values2$Freq, decreasing = F),]
    p_values2$n = c(1:nrow(p_values2))/nrow(p_values2)
    
    p_values = get.pval(obj,type_stat = "G")
    p_values[lower.tri(p_values,diag=TRUE)] <- NA
    
    p_values3 = as.data.frame(as.table(as.matrix(p_values))) 
    if(nrow(p_values3) > 10000){
      p_values3 = p_values3[sample(nrow(p_values3), (nrow(p_values3)/20)),]
    }
    
    p_values3 = p_values3[complete.cases(p_values3),]
    p_values3 = p_values3[order(p_values3$Freq, decreasing = F),]
    p_values3$n = c(1:nrow(p_values3))/nrow(p_values3)
    print(dim(p_values2))
    print(dim(p_values3))
    
    
    p_values2$Type = "Chi-squared test"
    p_values3$Type = "G-test"
    
    p_values = rbind(p_values2,p_values3)
    
    
    #---------for p-values form pearson correlation usign Seurat normalization
    seur.obj <- CreateSeuratObject(counts = as.matrix(obj@raw), project = "neg", 
                                   min.cells = 0, min.features = 2)
    seur.obj <- NormalizeData(seur.obj)

    #seur.obj[["RNA"]]@data[1:10,1:10]

    seurat.data = as.matrix(seur.obj[["RNA"]]@data)
    p_val.pearson = rcorr(t(seurat.data), type="pearson")
    
    p_values4 = as.data.frame(as.table(as.matrix(p_val.pearson$P))) 
    
    if(nrow(p_values4) > 10000){
        p_values4 = p_values4[sample(nrow(p_values4), (nrow(p_values4)/20)) , ]
    }
    p_values4 = p_values4[complete.cases(p_values4),]
    p_values4 = p_values4[order(p_values4$Freq, decreasing = F),]
    p_values4$n = c(1:nrow(p_values4))/nrow(p_values4)
    p_values4$Type = "Pearson"
    
    p_values = rbind(p_values,p_values4)
    #-----------------------------
    
    plot_p = ggplot(p_values, aes(x = Freq, y = n,colour = Type)) +
      theme(axis.text.x = element_text(size = 12, angle = 0, hjust = .5, vjust = .5, face = "plain"),
            axis.text.y = element_text( size = 12, angle = 0, hjust = 0, vjust = .5, face = "plain"),  
            axis.title.x = element_text( size = 12, angle = 0, hjust = .5, vjust = 0, face = "plain"),
            axis.title.y = element_text( size = 12, angle = 90, hjust = .5, vjust = .5, face = "plain"))       + 
      labs(x = "p-value", y = "percentile" ) +
      geom_line(size = 1.5) + 
      scale_x_log10(breaks = trans_breaks("log10", function(x) 10^x),
              labels = trans_format("log10", math_format(10^.x))) +
      scale_y_log10(breaks = trans_breaks("log10", function(x) 10^x),
              labels = trans_format("log10", math_format(10^.x))) +
      annotation_logticks()
    return(plot_p)
}

1 Technical negative dataset: ERCC 10x

NULL
[1] "Get p-values genome wide on columns genome wide on rows"
[1] "Using function S"
[1] "function to generate S "
NULL
[1] "Get p-values genome wide on columns genome wide on rows"
[1] "Using function G"
[1] "function to generate G "
[1] "Generating contingency tables for observed data"
[1] "mu estimator creation"
[1] "expected contingency tables creation"
[1] "The distance between estimated n of zeros and observed number of zero is 0.00446406311693362 over 66"
[1] "Done"
[1] "G estimation"
[1] 2145    4
[1] 2145    4
Performing log-normalization
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|

null device 
          1

GDI plot with density

[1] "function to generate GDI dataframe"
[1] "Using S"
[1] "function to generate S "

2 Biological negative dataset: CD14+

NULL
[1] "Get p-values genome wide on columns genome wide on rows"
[1] "Using function S"
[1] "function to generate S "
NULL
[1] "Get p-values genome wide on columns genome wide on rows"
[1] "Using function G"
[1] "function to generate G "
[1] "Generating contingency tables for observed data"
[1] "mu estimator creation"
[1] "expected contingency tables creation"
[1] "The distance between estimated n of zeros and observed number of zero is 0.0486027452293039 over 7850"
[1] "Done"
[1] "G estimation"
[1] 1540172       4
[1] 1542243       4
Performing log-normalization
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|

null device 
          1

GDI plot with density

LS0tCnRpdGxlOiAiTmVnYXRpdmUgZGF0YXNldCBhbmFseXNpcyIKb3V0cHV0OgogIGh0bWxfZG9jdW1lbnQ6CiAgICBjb2xsYXBzZWQ6IG5vCiAgICBjc3M6IGh0bWwtbWQtMDEuY3NzCiAgICBmaWdfY2FwdGlvbjogeWVzCiAgICBoaWdobGlnaHQ6IGhhZGRvY2sKICAgIG51bWJlcl9zZWN0aW9uczogeWVzCiAgICB0aGVtZTogc3BhY2VsYWIKICAgIHRvYzogeWVzCiAgICB0b2NfZmxvYXQ6IHllcwogIGh0bWxfbm90ZWJvb2s6CiAgICBjb2xsYXBzZWQ6IG5vCiAgICBjc3M6IGh0bWwtbWQtMDEuY3NzCiAgICBmaWdfY2FwdGlvbjogeWVzCiAgICBoaWdobGlnaHQ6IGhhZGRvY2sKICAgIG51bWJlcl9zZWN0aW9uczogeWVzCiAgICB0aGVtZTogc3BhY2VsYWIKICAgIHRvYzogeWVzCiAgICB0b2NfZmxvYXQ6IHllcwotLS0KYGBge3IsIGluY2x1ZGUgPSBGQUxTRX0Ka25pdHI6Om9wdHNfY2h1bmskc2V0KAogIGNvbGxhcHNlID0gVFJVRSwKICBjb21tZW50ID0gIiM+IiwKICBlY2hvID0gVFJVRSwKICB3YXJuaW5nID0gRkFMU0UsCiAgdGlkeSA9IFRSVUUsIHRpZHkub3B0cz1saXN0KHdpZHRoLmN1dG9mZj00MCkKKQpvcHRpb25zKHJtYXJrZG93bi5odG1sX3ZpZ25ldHRlLmNoZWNrX3RpdGxlID0gRkFMU0UpCmBgYAoKCmBgYHtyIH0KbGlicmFyeShnZ3Bsb3QyKQpsaWJyYXJ5KGdncmVwZWwpCmxpYnJhcnkoQ09UQU4pCmxpYnJhcnkoSG1pc2MpCmxpYnJhcnkoU2V1cmF0KQpsaWJyYXJ5KHBhdGNod29yaykKbGlicmFyeShzY2FsZXMpCgpwX3ZhbHVlX3Bsb3QgPC1mdW5jdGlvbihwX3ZhbHVlcywgb2JqKXsKICAgIHBfdmFsdWVzW2xvd2VyLnRyaShwX3ZhbHVlcyxkaWFnPVRSVUUpXSA8LSBOQQoKICAgIHBfdmFsdWVzMiA9IGFzLmRhdGEuZnJhbWUoYXMudGFibGUoYXMubWF0cml4KHBfdmFsdWVzKSkpIAogICAgCiAgICBpZihucm93KHBfdmFsdWVzMikgPiAxMDAwMCl7CiAgICAgIHBfdmFsdWVzMiA9IHBfdmFsdWVzMltzYW1wbGUobnJvdyhwX3ZhbHVlczIpLCAobnJvdyhwX3ZhbHVlczIpLzIwKSkgLCBdCiAgICB9CiAgICBwX3ZhbHVlczIgPSBwX3ZhbHVlczJbY29tcGxldGUuY2FzZXMocF92YWx1ZXMyKSxdCiAgICBwX3ZhbHVlczIgPSBwX3ZhbHVlczJbb3JkZXIocF92YWx1ZXMyJEZyZXEsIGRlY3JlYXNpbmcgPSBGKSxdCiAgICBwX3ZhbHVlczIkbiA9IGMoMTpucm93KHBfdmFsdWVzMikpL25yb3cocF92YWx1ZXMyKQogICAgCiAgICBwX3ZhbHVlcyA9IGdldC5wdmFsKG9iaix0eXBlX3N0YXQgPSAiRyIpCiAgICBwX3ZhbHVlc1tsb3dlci50cmkocF92YWx1ZXMsZGlhZz1UUlVFKV0gPC0gTkEKICAgIAogICAgcF92YWx1ZXMzID0gYXMuZGF0YS5mcmFtZShhcy50YWJsZShhcy5tYXRyaXgocF92YWx1ZXMpKSkgCiAgICBpZihucm93KHBfdmFsdWVzMykgPiAxMDAwMCl7CiAgICAgIHBfdmFsdWVzMyA9IHBfdmFsdWVzM1tzYW1wbGUobnJvdyhwX3ZhbHVlczMpLCAobnJvdyhwX3ZhbHVlczMpLzIwKSksXQogICAgfQogICAgCiAgICBwX3ZhbHVlczMgPSBwX3ZhbHVlczNbY29tcGxldGUuY2FzZXMocF92YWx1ZXMzKSxdCiAgICBwX3ZhbHVlczMgPSBwX3ZhbHVlczNbb3JkZXIocF92YWx1ZXMzJEZyZXEsIGRlY3JlYXNpbmcgPSBGKSxdCiAgICBwX3ZhbHVlczMkbiA9IGMoMTpucm93KHBfdmFsdWVzMykpL25yb3cocF92YWx1ZXMzKQogICAgcHJpbnQoZGltKHBfdmFsdWVzMikpCiAgICBwcmludChkaW0ocF92YWx1ZXMzKSkKICAgIAogICAgCiAgICBwX3ZhbHVlczIkVHlwZSA9ICJDaGktc3F1YXJlZCB0ZXN0IgogICAgcF92YWx1ZXMzJFR5cGUgPSAiRy10ZXN0IgogICAgCiAgICBwX3ZhbHVlcyA9IHJiaW5kKHBfdmFsdWVzMixwX3ZhbHVlczMpCiAgICAKICAgIAogICAgIy0tLS0tLS0tLWZvciBwLXZhbHVlcyBmb3JtIHBlYXJzb24gY29ycmVsYXRpb24gdXNpZ24gU2V1cmF0IG5vcm1hbGl6YXRpb24KICAgIHNldXIub2JqIDwtIENyZWF0ZVNldXJhdE9iamVjdChjb3VudHMgPSBhcy5tYXRyaXgob2JqQHJhdyksIHByb2plY3QgPSAibmVnIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbWluLmNlbGxzID0gMCwgbWluLmZlYXR1cmVzID0gMikKICAgIHNldXIub2JqIDwtIE5vcm1hbGl6ZURhdGEoc2V1ci5vYmopCgogICAgI3NldXIub2JqW1siUk5BIl1dQGRhdGFbMToxMCwxOjEwXQoKICAgIHNldXJhdC5kYXRhID0gYXMubWF0cml4KHNldXIub2JqW1siUk5BIl1dQGRhdGEpCiAgICBwX3ZhbC5wZWFyc29uID0gcmNvcnIodChzZXVyYXQuZGF0YSksIHR5cGU9InBlYXJzb24iKQogICAgCiAgICBwX3ZhbHVlczQgPSBhcy5kYXRhLmZyYW1lKGFzLnRhYmxlKGFzLm1hdHJpeChwX3ZhbC5wZWFyc29uJFApKSkgCiAgICAKICAgIGlmKG5yb3cocF92YWx1ZXM0KSA+IDEwMDAwKXsKICAgICAgICBwX3ZhbHVlczQgPSBwX3ZhbHVlczRbc2FtcGxlKG5yb3cocF92YWx1ZXM0KSwgKG5yb3cocF92YWx1ZXM0KS8yMCkpICwgXQogICAgfQogICAgcF92YWx1ZXM0ID0gcF92YWx1ZXM0W2NvbXBsZXRlLmNhc2VzKHBfdmFsdWVzNCksXQogICAgcF92YWx1ZXM0ID0gcF92YWx1ZXM0W29yZGVyKHBfdmFsdWVzNCRGcmVxLCBkZWNyZWFzaW5nID0gRiksXQogICAgcF92YWx1ZXM0JG4gPSBjKDE6bnJvdyhwX3ZhbHVlczQpKS9ucm93KHBfdmFsdWVzNCkKICAgIHBfdmFsdWVzNCRUeXBlID0gIlBlYXJzb24gb24gU2V1cmF0IG5vcm1hbGl6ZWQgZGF0YSIKICAgIAogICAgcF92YWx1ZXMgPSByYmluZChwX3ZhbHVlcyxwX3ZhbHVlczQpCiAgICAjLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KICAgIAogICAgcGxvdF9wID0gZ2dwbG90KHBfdmFsdWVzLCBhZXMoeCA9IEZyZXEsIHkgPSBuLGNvbG91ciA9IFR5cGUpKSArCiAgICAgIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxMiwgYW5nbGUgPSAwLCBoanVzdCA9IC41LCB2anVzdCA9IC41LCBmYWNlID0gInBsYWluIiksCiAgICAgICAgICAgIGF4aXMudGV4dC55ID0gZWxlbWVudF90ZXh0KCBzaXplID0gMTIsIGFuZ2xlID0gMCwgaGp1c3QgPSAwLCB2anVzdCA9IC41LCBmYWNlID0gInBsYWluIiksICAKICAgICAgICAgICAgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KCBzaXplID0gMTIsIGFuZ2xlID0gMCwgaGp1c3QgPSAuNSwgdmp1c3QgPSAwLCBmYWNlID0gInBsYWluIiksCiAgICAgICAgICAgIGF4aXMudGl0bGUueSA9IGVsZW1lbnRfdGV4dCggc2l6ZSA9IDEyLCBhbmdsZSA9IDkwLCBoanVzdCA9IC41LCB2anVzdCA9IC41LCBmYWNlID0gInBsYWluIikpICAgICAgICsgIGdlb21fYWJsaW5lKGxpbmV0eXBlPSJkYXNoZWQiKSsKICAgICAgbGFicyh4ID0gInAtdmFsdWUiLCB5ID0gInBlcmNlbnRpbGUiICkgKwogICAgICBnZW9tX2xpbmUoc2l6ZSA9IDEuNSkgIAogICAgICAjc2NhbGVfeF9sb2cxMChicmVha3MgPSB0cmFuc19icmVha3MoImxvZzEwIiwgZnVuY3Rpb24oeCkgMTBeeCksCiAgICAgICMgICAgICAgIGxhYmVscyA9IHRyYW5zX2Zvcm1hdCgibG9nMTAiLCBtYXRoX2Zvcm1hdCgxMF4ueCkpKSArCiAgICAgICNzY2FsZV95X2xvZzEwKGJyZWFrcyA9IHRyYW5zX2JyZWFrcygibG9nMTAiLCBmdW5jdGlvbih4KSAxMF54KSwKICAgICAgIyAgICAgICAgbGFiZWxzID0gdHJhbnNfZm9ybWF0KCJsb2cxMCIsIG1hdGhfZm9ybWF0KDEwXi54KSkpICsKICAgICAgI2Fubm90YXRpb25fbG9ndGlja3MoKQogICAgcmV0dXJuKHBsb3RfcCkKfQoKcGxvdC5HREkuZGVuc2l0eSA8LSBmdW5jdGlvbihHREkuZGYpewogIHNpPTExCiAgbXljb2xvdXJzIDwtIGMoIkNvbnN0aXR1dGl2ZSIgPSAiIzAwQTA4N0ZGIiwiZGlmIj0iI0U2NEIzNUZGIiwibm9ybWFsIiA9ICIjODQ5MUI0QjIiKQogICAgICAgIHRoZW1leD0gdGhlbWUoYXhpcy50ZXh0LnggPSAKICAgICAgICAgICAgICAgICAgICAgICAgZWxlbWVudF90ZXh0KCBzaXplID0gc2ksIGFuZ2xlID0gOTAsIGhqdXN0ID0gLjUsIHZqdXN0ID0gLjUsIGZhY2UgPSAicGxhaW4iLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb2xvdXIgPSIjM0M1NDg4RkYiKSwKICAgICAgICAgICAgICAgICAgICBheGlzLnRleHQueSA9IGVsZW1lbnRfYmxhbmsoKSwgIAogICAgICAgICAgICAgICAgICAgIGF4aXMudGl0bGUueCA9IGVsZW1lbnRfYmxhbmsoKSwKICAgICAgICAgICAgICAgICAgICBheGlzLnRpdGxlLnkgPSAKICAgICAgICAgICAgICAgICAgICAgIGVsZW1lbnRfdGV4dCggc2l6ZSA9IHNpLCBhbmdsZSA9IDkwLCBoanVzdCA9IC41LCB2anVzdCA9IC41LCBmYWNlID0gInBsYWluIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbG91ciA9IiMzQzU0ODhGRiIpLAogICAgICAgICAgICAgICAgICAgIGxlZ2VuZC5wb3NpdGlvbiA9ICJub25lIikgCiAgICAgIAogICAgICB0aGVtZXkgPSB0aGVtZShheGlzLnRleHQueSA9IAogICAgICAgICAgICAgICAgICAgICAgIGVsZW1lbnRfdGV4dCggc2l6ZSA9IHNpLCBhbmdsZSA9IDAsIGhqdXN0ID0gLjUsIHZqdXN0ID0gLjUsIGZhY2UgPSAicGxhaW4iLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbG91ciA9IiMzQzU0ODhGRiIpLAogICAgICAgICAgICAgICAgICAgICBheGlzLnRpdGxlLnggPSAKICAgICAgICAgICAgICAgICAgICAgICBlbGVtZW50X3RleHQoIHNpemUgPSBzaSwgYW5nbGUgPSAwLCBoanVzdCA9IC41LCB2anVzdCA9IC41LCBmYWNlID0gInBsYWluIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb2xvdXIgPSIjM0M1NDg4RkYiKSwKICAgICAgICAgICAgICAgICAgICAgYXhpcy50ZXh0LnguYm90dG9tID0gZWxlbWVudF9ibGFuaygpLAogICAgICAgICAgICAgICAgICAgICBheGlzLnRpdGxlLnkgPSBlbGVtZW50X2JsYW5rKCksCiAgICAgICAgICAgICAgICAgICAgIGxlZ2VuZC5wb3NpdGlvbiA9ICJub25lIikgCgogIAogICAgICBmMSA9IGdncGxvdChHREkuZGYsIGFlcyh4PXN1bS5yYXcubm9ybSwgeT1HREkpKSArICBnZW9tX3BvaW50KGFscGhhID0gMC40LCBjb2xvciA9ICIjODQ5MUI0QjIiLCBzaXplPTIuNSkKICAgICAgCiAgICAgIAogICAgICBHREkuZGZfbGluID0gZjEgKyAKICAgICAgICBnZW9tX2hsaW5lKHlpbnRlcmNlcHQ9MS41LCBsaW5ldHlwZT0iZG90dGVkIiwgY29sb3IgPSAicmVkIiwgc2l6ZT0gMSkgKwogICAgICAgIHNjYWxlX2NvbG9yX21hbnVhbCgiU3RhdHVzIiwgdmFsdWVzID0gbXljb2xvdXJzKSAgKwogICAgICAgIHNjYWxlX2ZpbGxfbWFudWFsKCJTdGF0dXMiLCB2YWx1ZXMgPSBteWNvbG91cnMpICArCiAgICAgICAgeGxhYigibG9nIG5vcm1hbGl6ZWQgY291bnRzIikreWxhYigiR0RJIikrCiAgICAgICAgdGhlbWUoYXhpcy50ZXh0LnggPSAKICAgICAgICAgICAgICAgIGVsZW1lbnRfdGV4dChzaXplID0gc2ksIGFuZ2xlID0gMCwgaGp1c3QgPSAuNSwgdmp1c3QgPSAuNSwgZmFjZSA9ICJwbGFpbiIsIGNvbG91ciA9IiMzQzU0ODhGRiIgKSwKICAgICAgICAgICAgICBheGlzLnRleHQueSA9IAogICAgICAgICAgICAgICAgZWxlbWVudF90ZXh0KCBzaXplID0gc2ksIGFuZ2xlID0gMCwgaGp1c3QgPSAwLCB2anVzdCA9IC41LCBmYWNlID0gInBsYWluIiwgY29sb3VyID0iIzNDNTQ4OEZGIiksICAKICAgICAgICAgICAgICBheGlzLnRpdGxlLnggPSAKICAgICAgICAgICAgICAgIGVsZW1lbnRfdGV4dCggc2l6ZSA9IHNpLCBhbmdsZSA9IDAsIGhqdXN0ID0gLjUsIHZqdXN0ID0gMCwgZmFjZSA9ICJwbGFpbiIsIGNvbG91ciA9IiMzQzU0ODhGRiIpLAogICAgICAgICAgICAgIGF4aXMudGl0bGUueSA9IAogICAgICAgICAgICAgICAgZWxlbWVudF90ZXh0KCBzaXplID0gc2ksIGFuZ2xlID0gOTAsIGhqdXN0ID0gLjUsIHZqdXN0ID0gLjUsIGZhY2UgPSAicGxhaW4iLCBjb2xvdXIgPSIjM0M1NDg4RkYiKSwKICAgICAgICAgICAgICBsZWdlbmQudGl0bGUgPSBlbGVtZW50X2JsYW5rKCksCiAgICAgICAgICAgICAgbGVnZW5kLnRleHQgPSBlbGVtZW50X3RleHQoY29sb3IgPSAiIzNDNTQ4OEZGIixmYWNlID0iaXRhbGljIiApLAogICAgICAgICAgICAgIGxlZ2VuZC5wb3NpdGlvbiA9ICJub25lIikgCiAgICAgIAogICAgICAgIHhkZW5zaXR5R0RJLmRmIDwtIGdncGxvdChHREkuZGYsIGFlcyhzdW0ucmF3Lm5vcm0pKSArIAogICAgICAgICAgZ2VvbV9kZW5zaXR5KGFscGhhPS41LCBmaWxsID0gIiM4NDkxQjRCMiIsIGNvbG91ciA9IiM4NDkxQjRCMiIgKSArCiAgICAgICAgICB0aGVtZXggIAogICAgICAKICAgICAgICB5ZGVuc2l0eUdESS5kZiA8LSBnZ3Bsb3QoR0RJLmRmLCBhZXMoR0RJKSkgKyAKICAgICAgICAgIGdlb21fZGVuc2l0eShhbHBoYT0uNSwgZmlsbD0iIzAwQTA4N0ZGIiwgY29sb3VyPSAiIzAwQTA4N0ZGIikgKyAKICAgICAgICAgIHRoZW1leSArICBjb29yZF9mbGlwKCkgCiAgICAgICAgICAKICAgICAgCiAgICAgICAgR0RJLmRmX2xpbiA9IHhkZW5zaXR5R0RJLmRmICsgcGxvdF9zcGFjZXIoKSArIEdESS5kZl9saW4rIHlkZW5zaXR5R0RJLmRmICsgCiAgICAgICAgICBwbG90X2xheW91dChuY29sID0gMiwgbnJvdyA9IDIsIHdpZHRocyA9IGMoNCwgMSksIGhlaWdodHMgPSBjKDEsIDQpKQogICAgICAgIHJldHVybihHREkuZGZfbGluKQp9CmBgYAoKIyBUZWNobmljYWwgbmVnYXRpdmUgZGF0YXNldDogRVJDQyAxMHgKCmBgYHtyIGVjaG89VFJVRX0KRVJDQyA9IHJlYWRSRFMoIkRhdGEvbmVnYXRpdmVfZGF0YXNldHMvRVJDQy5jb3Rhbi5SRFMiKQpgYGAKYGBge3IgZWNobz1UUlVFfQpwX3ZhbHVlc19FUkNDID0gZ2V0LnB2YWwoRVJDQykKClBsb3RfRVJDQyA9IHBfdmFsdWVfcGxvdChwX3ZhbHVlc19FUkNDLCBFUkNDKQoKUGxvdF9FUkNDCmBgYApgYGB7cn0KUGxvdF9FUkNDKyB4bGltKDAsMC4wNSkrIHlsaW0oMCwwLjA1KQpgYGAKCgpgYGB7cn0KcGRmKCJFUkNDX3BfdmFsdWVfU19HX3Bsb3QucGRmIikKUGxvdF9FUkNDClBsb3RfRVJDQysgeGxpbSgwLDAuMDUpKyB5bGltKDAsMC4wNSkKZGV2Lm9mZigpCmBgYAoKCkdESSBwbG90IHdpdGggZGVuc2l0eQoKYGBge3J9CkdESV9lcmNjID0gZ2V0LkdESShFUkNDKQoKcGxvdC5HREkuZGVuc2l0eShHRElfZXJjYykKYGBgCgoKIyBCaW9sb2dpY2FsIG5lZ2F0aXZlIGRhdGFzZXQ6IENEMTQrCmBgYHtyIGVjaG89VFJVRX0KQ0QxNCA9IHJlYWRSRFMoIkRhdGEvbmVnYXRpdmVfZGF0YXNldHMvQ0QxNC5jb3Rhbi5SRFMiKQpgYGAKYGBge3IgZWNobz1UUlVFfQpwX3ZhbHVlc19DRDE0ID0gZ2V0LnB2YWwoQ0QxNCkKClBsb3RfQ0QxNCA9IHBfdmFsdWVfcGxvdChwX3ZhbHVlc19DRDE0LG9iaiA9IENEMTQpCgpQbG90X0NEMTQKYGBgCmBgYHtyfQpQbG90X0NEMTQgKyB4bGltKDAsMC4wNSkrIHlsaW0oMCwwLjA1KQpgYGAKCgpgYGB7cn0KcGRmKCJDRDE0X3BfdmFsdWVfU19HX3Bsb3QucGRmIikKUGxvdF9DRDE0ClBsb3RfQ0QxNCArIHhsaW0oMCwwLjA1KSsgeWxpbSgwLDAuMDUpCmRldi5vZmYoKQpgYGAKCgpHREkgcGxvdCB3aXRoIGRlbnNpdHkKCmBgYHtyfQpHRElfQ0QxNCA9IGdldC5HREkoQ0QxNCkKCnBsb3QuR0RJLmRlbnNpdHkoR0RJX0NEMTQpCmBgYApgYGB7cn0KZ2VuZXMgPSBsaXN0KCJnZW5lcyI9cm93bmFtZXMoR0RJX0NEMTRbR0RJX0NEMTQkR0RJPjEuNSAmIEdESV9DRDE0JHN1bS5yYXcubm9ybSA+IDUgJiBHRElfQ0QxNCRleHAuY2VsbHMgPiAyLjUsXSkpCnBsb3RfaGVhdG1hcChkZl9nZW5lcyA9IGdlbmVzLHNldHMgPSAxLGNvbmRpdGlvbnMgPSAiQ0QxNCIsZGlyID0gIkRhdGEvbmVnYXRpdmVfZGF0YXNldHMvIikKYGBgCgoKCgoKIyBTeW50ZXRpYyBuZWdhdGl2ZSBkYXRhc2V0OiA0MDAwIGNlbGxzCmBgYHtyIH0KQ0U0MDAwID0gcmVhZFJEUygiRGF0YS9uZWdhdGl2ZV9kYXRhc2V0cy9TeW1FMTdfMWNsXzQwMDBjcy5jb3Rhbi5SRFMiKQpgYGAKYGBge3IgfQpwX3ZhbHVlc19DRTQwMDAgPSBnZXQucHZhbChDRTQwMDApCgpQbG90X0NFNDAwMCA9IHBfdmFsdWVfcGxvdChwX3ZhbHVlc19DRTQwMDAsIG9iaiA9IENFNDAwMCkKClBsb3RfQ0U0MDAwCgpgYGAKCmBgYHtyfQpQbG90X0NFNDAwMCArIHhsaW0oMCwwLjA1KSsgeWxpbSgwLDAuMDUpCmBgYAoKCmBgYHtyfQpwZGYoIkNFNDAwMF9wX3ZhbHVlX1NfR19wbG90LnBkZiIpClBsb3RfQ0U0MDAwClBsb3RfQ0U0MDAwICsgeGxpbSgwLDAuMDUpKyB5bGltKDAsMC4wNSkKZGV2Lm9mZigpCmBgYAoKYGBge3J9CkdESV9DRTQwMDAgPSBnZXQuR0RJKENFNDAwMCkKCnBsb3QuR0RJLmRlbnNpdHkoR0RJX0NFNDAwMCkKYGBgCgoKIyBTeW50ZXRpYyBuZWdhdGl2ZSBkYXRhc2V0OiA4MDAgY2VsbHMKYGBge3IgZWNobz1UUlVFfQpDRTgwMCA9IHJlYWRSRFMoIkRhdGEvbmVnYXRpdmVfZGF0YXNldHMvU3ltRTE3XzFjbF84MDBjcy5jb3Rhbi5SRFMiKQpgYGAKYGBge3IgZWNobz1UUlVFfQpwX3ZhbHVlc19DRTgwMCA9IGdldC5wdmFsKENFODAwKQoKUGxvdF9DRTgwMCA9IHBfdmFsdWVfcGxvdChwX3ZhbHVlc19DRTgwMCwgb2JqID0gQ0U4MDApCgpQbG90X0NFODAwCmBgYApgYGB7cn0KUGxvdF9DRTgwMCsgeGxpbSgwLDAuMDUpKyB5bGltKDAsMC4wNSkKYGBgCgpgYGB7cn0KcGRmKCJDRTgwMF9wX3ZhbHVlX1NfR19wbG90LnBkZiIpClBsb3RfQ0U4MDAKUGxvdF9DRTgwMCsgeGxpbSgwLDAuMDUpKyB5bGltKDAsMC4wNSkKZGV2Lm9mZigpCmBgYAoKCmBgYHtyfQpzZXNzaW9uSW5mbygpCmBgYAoKCg==