Back to Multiple platform build/check report for BioC 3.20: simplified long |
|
This page was generated on 2024-07-23 11:42 -0400 (Tue, 23 Jul 2024).
Hostname | OS | Arch (*) | R version | Installed pkgs |
---|---|---|---|---|
nebbiolo2 | Linux (Ubuntu 22.04.3 LTS) | x86_64 | 4.4.1 (2024-06-14) -- "Race for Your Life" | 4688 |
palomino8 | Windows Server 2022 Datacenter | x64 | 4.4.1 (2024-06-14 ucrt) -- "Race for Your Life" | 4280 |
lconway | macOS 12.7.1 Monterey | x86_64 | 4.4.1 (2024-06-14) -- "Race for Your Life" | 4455 |
kjohnson3 | macOS 13.6.5 Ventura | arm64 | 4.4.1 (2024-06-14) -- "Race for Your Life" | 4404 |
Click on any hostname to see more info about the system (e.g. compilers) (*) as reported by 'uname -p', except on Windows and Mac OS X |
Package 439/2248 | Hostname | OS / Arch | INSTALL | BUILD | CHECK | BUILD BIN | ||||||||
COTAN 2.5.6 (landing page) Galfrè Silvia Giulia
| nebbiolo2 | Linux (Ubuntu 22.04.3 LTS) / x86_64 | OK | OK | OK | ![]() | ||||||||
palomino8 | Windows Server 2022 Datacenter / x64 | OK | OK | OK | OK | ![]() | ||||||||
lconway | macOS 12.7.1 Monterey / x86_64 | OK | OK | OK | OK | ![]() | ||||||||
kjohnson3 | macOS 13.6.5 Ventura / arm64 | OK | OK | OK | OK | ![]() | ||||||||
To the developers/maintainers of the COTAN package: - Allow up to 24 hours (and sometimes 48 hours) for your latest push to git@git.bioconductor.org:packages/COTAN.git to reflect on this report. See Troubleshooting Build Report for more information. - Use the following Renviron settings to reproduce errors and warnings. - If 'R CMD check' started to fail recently on the Linux builder(s) over a missing dependency, add the missing dependency to 'Suggests:' in your DESCRIPTION file. See Renviron.bioc for more information. |
Package: COTAN |
Version: 2.5.6 |
Command: /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:COTAN.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings COTAN_2.5.6.tar.gz |
StartedAt: 2024-07-22 20:53:26 -0400 (Mon, 22 Jul 2024) |
EndedAt: 2024-07-22 21:02:29 -0400 (Mon, 22 Jul 2024) |
EllapsedTime: 542.9 seconds |
RetCode: 0 |
Status: OK |
CheckDir: COTAN.Rcheck |
Warnings: 0 |
############################################################################## ############################################################################## ### ### Running command: ### ### /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:COTAN.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings COTAN_2.5.6.tar.gz ### ############################################################################## ############################################################################## * using log directory ‘/Users/biocbuild/bbs-3.20-bioc/meat/COTAN.Rcheck’ * using R version 4.4.1 (2024-06-14) * using platform: aarch64-apple-darwin20 * R was compiled by Apple clang version 14.0.0 (clang-1400.0.29.202) GNU Fortran (GCC) 12.2.0 * running under: macOS Ventura 13.6.7 * using session charset: UTF-8 * using option ‘--no-vignettes’ * checking for file ‘COTAN/DESCRIPTION’ ... OK * checking extension type ... Package * this is package ‘COTAN’ version ‘2.5.6’ * package encoding: UTF-8 * checking package namespace information ... OK * checking package dependencies ... OK * checking if this is a source package ... OK * checking if there is a namespace ... OK * checking for hidden files and directories ... OK * checking for portable file names ... OK * checking for sufficient/correct file permissions ... OK * checking whether package ‘COTAN’ can be installed ... OK * checking installed package size ... OK * checking package directory ... OK * checking ‘build’ directory ... OK * checking DESCRIPTION meta-information ... OK * checking top-level files ... OK * checking for left-over files ... OK * checking index information ... OK * checking package subdirectories ... OK * checking code files for non-ASCII characters ... OK * checking R files for syntax errors ... OK * checking whether the package can be loaded ... OK * checking whether the package can be loaded with stated dependencies ... OK * checking whether the package can be unloaded cleanly ... OK * checking whether the namespace can be loaded with stated dependencies ... OK * checking whether the namespace can be unloaded cleanly ... OK * checking dependencies in R code ... NOTE 'library' or 'require' call to ‘torch’ in package code. Please use :: or requireNamespace() instead. See section 'Suggested packages' in the 'Writing R Extensions' manual. Unexported object imported by a ':::' call: ‘ggplot2:::ggname’ See the note in ?`:::` about the use of this operator. * checking S3 generic/method consistency ... OK * checking replacement functions ... OK * checking foreign function calls ... OK * checking R code for possible problems ... NOTE GDIPlot: no visible binding for global variable ‘sum.raw.norm’ GDIPlot: no visible binding for global variable ‘GDI’ UMAPPlot: no visible binding for global variable ‘x’ UMAPPlot: no visible binding for global variable ‘y’ UMAPPlot: no visible binding for global variable ‘types’ calculateCoex_Legacy: no visible binding for global variable ‘expectedNN’ calculateCoex_Legacy: no visible binding for global variable ‘expectedNY’ calculateCoex_Legacy: no visible binding for global variable ‘expectedYN’ calculateCoex_Legacy: no visible binding for global variable ‘expectedYY’ calculateCoex_Legacy: no visible binding for global variable ‘observedYY’ calculateCoex_Legacy: no visible binding for global variable ‘.’ calculateG: no visible binding for global variable ‘observedNN’ calculateG: no visible binding for global variable ‘observedNY’ calculateG: no visible binding for global variable ‘observedYN’ calculateG: no visible binding for global variable ‘observedYY’ calculateG: no visible binding for global variable ‘expectedNN’ calculateG: no visible binding for global variable ‘expectedNY’ calculateG: no visible binding for global variable ‘expectedYN’ calculateG: no visible binding for global variable ‘expectedYY’ calculatePartialCoex: no visible binding for global variable ‘expectedNN’ calculatePartialCoex: no visible binding for global variable ‘expectedNY’ calculatePartialCoex: no visible binding for global variable ‘expectedYN’ calculatePartialCoex: no visible binding for global variable ‘expectedYY’ calculatePartialCoex: no visible binding for global variable ‘observedYY’ calculatePartialCoex: no visible binding for global variable ‘.’ cellsUniformClustering: no visible binding for global variable ‘objSeurat’ cellsUniformClustering: no visible binding for global variable ‘usedMaxResolution’ cellsUniformClustering: no visible binding for global variable ‘permMap’ checkClusterUniformity: ... may be used in an incorrect context: ‘c(..., nuPlot, zoomedNuPlot)’ checkClusterUniformity: no visible binding for global variable ‘nuPlot’ checkClusterUniformity: no visible binding for global variable ‘zoomedNuPlot’ cleanPlots: no visible binding for global variable ‘PC1’ cleanPlots: no visible binding for global variable ‘PC2’ cleanPlots: no visible binding for global variable ‘n’ cleanPlots: no visible binding for global variable ‘means’ cleanPlots: no visible binding for global variable ‘nu’ clustersMarkersHeatmapPlot: no visible binding for global variable ‘condName’ clustersMarkersHeatmapPlot: no visible binding for global variable ‘conditions’ clustersSummaryPlot: no visible binding for global variable ‘keys’ clustersSummaryPlot: no visible binding for global variable ‘values’ clustersSummaryPlot: no visible binding for global variable ‘CellNumber’ clustersSummaryPlot: no visible binding for global variable ‘ExpGenes’ clustersSummaryPlot: no visible binding for global variable ‘Cluster’ clustersSummaryPlot: no visible binding for global variable ‘Condition’ establishGenesClusters: no visible binding for global variable ‘secondaryMarkers’ establishGenesClusters: no visible binding for global variable ‘GCS’ establishGenesClusters: no visible binding for global variable ‘rankGenes’ expectedContingencyTables: no visible binding for global variable ‘expectedN’ expectedPartialContingencyTables: no visible binding for global variable ‘expectedNN’ expectedPartialContingencyTables: no visible binding for global variable ‘expectedN’ geom_flat_violin : <anonymous>: no visible binding for global variable ‘group’ geom_flat_violin : <anonymous>: no visible binding for global variable ‘y’ geom_flat_violin : <anonymous>: no visible binding for global variable ‘x’ geom_flat_violin : <anonymous>: no visible binding for global variable ‘width’ geom_flat_violin : <anonymous>: no visible binding for global variable ‘violinwidth’ geom_flat_violin : <anonymous>: no visible binding for global variable ‘xmax’ geom_flat_violin : <anonymous>: no visible binding for global variable ‘xminv’ geom_flat_violin : <anonymous>: no visible binding for global variable ‘xmaxv’ heatmapPlot: no visible binding for global variable ‘obj’ heatmapPlot: no visible binding for global variable ‘type’ heatmapPlot: no visible binding for global variable ‘g2’ heatmapPlot: no visible binding for global variable ‘coex’ mergeUniformCellsClusters : testPairListMerge: no visible binding for global variable ‘cl1’ mergeUniformCellsClusters : testPairListMerge: no visible binding for global variable ‘cl2’ mergeUniformCellsClusters : testPairListMerge : <anonymous>: no visible binding for global variable ‘cl1’ mergeUniformCellsClusters : testPairListMerge : <anonymous>: no visible binding for global variable ‘cl2’ mergeUniformCellsClusters: no visible binding for global variable ‘permMap’ mitochondrialPercentagePlot: no visible binding for global variable ‘mit.percentage’ observedContingencyTables: no visible binding for global variable ‘observedY’ observedPartialContingencyTables: no visible binding for global variable ‘observedYY’ observedPartialContingencyTables: no visible binding for global variable ‘observedY’ scatterPlot: no visible binding for global variable ‘.x’ calculateCoex,COTAN: no visible binding for global variable ‘useTorch’ coerce,COTAN-scCOTAN: no visible binding for global variable ‘rawNorm’ coerce,COTAN-scCOTAN: no visible binding for global variable ‘nu’ coerce,COTAN-scCOTAN: no visible binding for global variable ‘lambda’ coerce,COTAN-scCOTAN: no visible binding for global variable ‘a’ coerce,COTAN-scCOTAN: no visible binding for global variable ‘hk’ coerce,COTAN-scCOTAN: no visible binding for global variable ‘clusters’ coerce,COTAN-scCOTAN: no visible binding for global variable ‘clusterData’ Undefined global functions or variables: . .x CellNumber Cluster Condition ExpGenes GCS GDI PC1 PC2 a cl1 cl2 clusterData clusters coex condName conditions expectedN expectedNN expectedNY expectedYN expectedYY g2 group hk keys lambda means mit.percentage n nu nuPlot obj objSeurat observedNN observedNY observedY observedYN observedYY permMap rankGenes rawNorm secondaryMarkers sum.raw.norm type types useTorch usedMaxResolution values violinwidth width x xmax xmaxv xminv y zoomedNuPlot * checking Rd files ... OK * checking Rd metadata ... OK * checking Rd cross-references ... OK * checking for missing documentation entries ... OK * checking for code/documentation mismatches ... OK * checking Rd \usage sections ... OK * checking Rd contents ... OK * checking for unstated dependencies in examples ... OK * checking contents of ‘data’ directory ... OK * checking data for non-ASCII characters ... OK * checking data for ASCII and uncompressed saves ... OK * checking files in ‘vignettes’ ... OK * checking examples ... OK Examples with CPU (user + system) or elapsed time > 5s user system elapsed UniformClusters 84.655 0.700 85.318 CalculatingCOEX 23.399 0.398 23.812 HeatmapPlots 16.718 0.378 17.098 HandlingClusterizations 12.342 0.313 12.643 ParametersEstimations 10.824 0.339 11.169 * checking for unstated dependencies in ‘tests’ ... OK * checking tests ... Running ‘outputTestDatasetCreation.R’ Running ‘spelling.R’ Running ‘testthat.R’ OK * checking for unstated dependencies in vignettes ... OK * checking package vignettes ... OK * checking running R code from vignettes ... SKIPPED * checking re-building of vignette outputs ... SKIPPED * checking PDF version of manual ... OK * DONE Status: 2 NOTEs See ‘/Users/biocbuild/bbs-3.20-bioc/meat/COTAN.Rcheck/00check.log’ for details.
COTAN.Rcheck/00install.out
############################################################################## ############################################################################## ### ### Running command: ### ### /Library/Frameworks/R.framework/Resources/bin/R CMD INSTALL COTAN ### ############################################################################## ############################################################################## * installing to library ‘/Library/Frameworks/R.framework/Versions/4.4-arm64/Resources/library’ * installing *source* package ‘COTAN’ ... ** using staged installation ** R ** data ** inst ** byte-compile and prepare package for lazy loading Note: ... may be used in an incorrect context ** help *** installing help indices ** building package indices ** installing vignettes ** testing if installed package can be loaded from temporary location ** testing if installed package can be loaded from final location ** testing if installed package keeps a record of temporary installation path * DONE (COTAN)
COTAN.Rcheck/tests/outputTestDatasetCreation.Rout
R version 4.4.1 (2024-06-14) -- "Race for Your Life" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: aarch64-apple-darwin20 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > > # Creates the files to be reloaded by the tests for comparisons > library(zeallot) > > outputTestDatasetCreation <- function(testsDir = file.path("tests", + "testthat")) { + utils::data("test.dataset", package = "COTAN") + options(parallelly.fork.enable = TRUE) + + obj <- COTAN(raw = test.dataset) + obj <- initializeMetaDataset(obj, GEO = " ", + sequencingMethod = "artificial", + sampleCondition = "test") + + obj <- proceedToCoex(obj, cores = 6L, saveObj = FALSE) + #saveRDS(obj, file = file.path(testsDir,"temp.RDS")) + + cells.names.test <- getCells(obj)[c(1L:10L, 591L:610L, 991L:1000L)] + genes.names.test <- getGenes(obj)[c(1L:10L, 291L:310L, 591L: 600L)] + saveRDS(cells.names.test, file.path(testsDir, "cells.names.test.RDS")) + saveRDS(genes.names.test, file.path(testsDir, "genes.names.test.RDS")) + + dispersion.test <- getDispersion(obj)[genes.names.test] + saveRDS(dispersion.test, file.path(testsDir, "dispersion.test.RDS")) + + raw.norm.test <- getNuNormData(obj)[genes.names.test, cells.names.test] + saveRDS(raw.norm.test, file.path(testsDir, "raw.norm.test.RDS")) + + coex.test <- getGenesCoex(obj, genes = genes.names.test, zeroDiagonal = FALSE) + saveRDS(coex.test, file.path(testsDir, "coex.test.RDS")) + + lambda.test <- getLambda(obj)[genes.names.test] + saveRDS(lambda.test, file.path(testsDir, "lambda.test.RDS")) + + GDI.test <- calculateGDI(obj) + GDI.test <- GDI.test[genes.names.test, ] + saveRDS(GDI.test, file.path(testsDir, "GDI.test.RDS")) + + nu.test <- getNu(obj)[cells.names.test] + saveRDS(nu.test, file.path(testsDir, "nu.test.RDS")) + + pval.test <- calculatePValue(obj, geneSubsetCol = genes.names.test) + saveRDS(pval.test, file.path(testsDir, "pval.test.RDS")) + + GDIThreshold <- 1.46 + ratioAboveThreshold <- 0.01 + initialResolution <- 0.8 + clusters <- cellsUniformClustering(obj, GDIThreshold = GDIThreshold, + ratioAboveThreshold = ratioAboveThreshold, + initialResolution = initialResolution, + cores = 6L, saveObj = FALSE)[["clusters"]] + saveRDS(clusters, file.path(testsDir, "clusters1.RDS")) + + coexDF <- DEAOnClusters(obj, clusters = clusters) + obj <- addClusterization(obj, clName = "clusters", + clusters = clusters, coexDF = coexDF) + + saveRDS(coexDF[genes.names.test, ], + file.path(testsDir, "coex.test.cluster1.RDS")) + + pvalDF <- pValueFromDEA(coexDF, getNumCells(obj), adjustmentMethod = "none") + + saveRDS(pvalDF[genes.names.test, ], + file.path(testsDir, "pval.test.cluster1.RDS")) + + c(mergedClusters, mCoexDF) %<-% + mergeUniformCellsClusters(objCOTAN = obj, + clusters = NULL, + GDIThreshold = GDIThreshold, + cores = 6L, + distance = "cosine", + hclustMethod = "ward.D2", + saveObj = FALSE) + + saveRDS(mergedClusters[genes.names.test], + file.path(testsDir, "cluster_data_merged.RDS")) + } > > proc.time() user system elapsed 0.078 0.018 0.092
COTAN.Rcheck/tests/spelling.Rout
R version 4.4.1 (2024-06-14) -- "Race for Your Life" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: aarch64-apple-darwin20 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > if (requireNamespace("spelling", quietly = TRUE)) + spelling::spell_check_test(vignettes = TRUE, error = FALSE, + skip_on_cran = TRUE) NULL > > proc.time() user system elapsed 0.072 0.020 0.087
COTAN.Rcheck/tests/testthat.Rout
R version 4.4.1 (2024-06-14) -- "Race for Your Life" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: aarch64-apple-darwin20 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > Sys.setenv(R_TESTS = "") > library(testthat) > library(COTAN) > test_check("COTAN") Setting new log level to 4 Initializing `COTAN` meta-data Genes/cells selection done: dropped [0] genes and [0] cells Working on [10] genes and [20] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.9033203125 | max: 4.6796875 | % negative: 10 Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.05829095840454 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.000468969345092773 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00231599807739258 Estimating genes' coex Calculate genes' coex elapsed time: 0.375002145767212 Total calculations elapsed time: 1.43607807159424 Calculate genes' coex (legacy): DONE Genes/cells selection done: dropped [0] genes and [0] cells Working on [10] genes and [20] cells Estimate dispersion: START Effective number of cores used: 1 Executing 4 genes batches from [1:2] to [7:8] Executing 1 genes batches from [9:10] to [9:10] Estimate dispersion: DONE dispersion | min: 0.9033203125 | max: 4.6796875 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 4 cells batches from [1:3] to [10:11] Executing 3 cells batches from [12:14] to [18:20] Estimate nu: DONE nu change (abs) | max: 1.75595238095238 | median: 1.07174634176587 | mean: 1.07174634176587 Estimate 'dispersion'/'nu': START Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 1.0362548828125 | max: 4.60986328125 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.0265938895089288 | median: 0.0144680038331048 | mean: 0.0144680038331048 Nu mean: 1.69633192486233 Marginal errors | max: 1.95570586131367 | median 1.32068160171502 | mean: 1.33375826507259 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.058837890625 | max: 3.528076171875 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.416683423613994 | median: 0.239880630367975 | mean: 0.239880630367975 Nu mean: 0.823197206753982 Marginal errors | max: 0.836359531101206 | median 0.703684202571891 | mean: 0.645537958989614 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.32879638671875 | max: 4.0302734375 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.164237872898673 | median: 0.0955985184389135 | mean: 0.0955985184389135 Nu mean: 1.06863935445976 Marginal errors | max: 0.259872988828242 | median 0.213703042752633 | mean: 0.197386407582083 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.2294921875 | max: 3.8720703125 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.055185575120883 | median: 0.0319991762044448 | mean: 0.0319991762044448 Nu mean: 0.976813601083562 Marginal errors | max: 0.0951586919577032 | median 0.079429703709466 | mean: 0.0724140148396648 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.2637939453125 | max: 3.929443359375 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.0196211148938294 | median: 0.01138609597457 | mean: 0.01138609597457 Nu mean: 1.00823501891926 Marginal errors | max: 0.0327747321002292 | median 0.0272104747849538 | mean: 0.0248963830312038 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.25177001953125 | max: 3.90966796875 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.00670099066960717 | median: 0.00388888266671264 | mean: 0.00388888266671264 Nu mean: 0.997187891997105 Marginal errors | max: 0.0114324509186883 | median 0.0094232649770607 | mean: 0.00863113610779571 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.25592041015625 | max: 3.91650390625 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.00230093811689414 | median: 0.00132529122122246 | mean: 0.00132529122122246 Nu mean: 1.00097564689567 Marginal errors | max: 0.00387133150664631 | median 0.0031091017608853 | mean: 0.00286071175800213 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.2545166015625 | max: 3.914306640625 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.000837011646904529 | median: 0.000470837393363011 | mean: 0.000470837393363011 Nu mean: 0.999633825746458 Marginal errors | max: 0.00122501723202006 | median 0.00102126435760308 | mean: 0.000943992659051318 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.2550048828125 | max: 3.9150390625 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.000209227351122054 | median: 0.000122070312500028 | mean: 0.000122070312500028 Nu mean: 1.00008715703862 Marginal errors | max: 0.000364602956581805 | median 0.000313956936819793 | mean: 0.000282899574318485 Estimate dispersion/nu: DONE Estimate 'dispersion'/'nu': START Initializing `COTAN` meta-data Genes/cells selection done: dropped [0] genes and [0] cells Working on [10] genes and [20] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.9033203125 | max: 4.6796875 | % negative: 10 While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.06974220275879 Calculating genes' coex normalization factor Fraction of genes' with very low expected contingency tables: 0.181818181818182 Calculate genes' normalization factor elapsed time: 0.000321865081787109 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0013129711151123 Estimating genes' coex Calculate genes' coex elapsed time: 0.357406139373779 Total calculations elapsed time: 1.42878317832947 Calculate genes' coex (legacy): DONE Calculate cells' coex (legacy): START Retrieving expected cells' contingency table calculating NN.. done calculating YN..NY..YY..t().. done Expected cells' contingency table elapsed time: 1.0862832069397 Calculating cells' coex normalization factor Fraction of cells' with very low expected contingency tables: 0 Calculate cells' normalization factor elapsed time: 0.000276803970336914 Retrieving observed cells' yes/yes contingency table calculating YY.. done Observed cells' contingency table elapsed time: 0.00083613395690918 Estimating cells' coex Calculate cells' coex elapsed time: 0.350843906402588 Total calculations elapsed time: 1.43824005126953 Calculate cells' coex (legacy): DONE Initializing `COTAN` meta-data Initializing `COTAN` meta-data Genes/cells selection done: dropped [0] genes and [0] cells Working on [10] genes and [20] cells Estimate 'dispersion'/'nu': START Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.903564453125 | max: 4.679443359375 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 1.75719246031746 | median: 1.07229953342014 | mean: 1.07229953342014 Nu mean: 1.68489292689732 Marginal errors | max: 1.73564890252257 | median 1.37996360874076 | mean: 1.32180348113228 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.0655517578125 | max: 3.5439453125 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.402649984216273 | median: 0.231868788425666 | mean: 0.231868788425666 Nu mean: 0.829218804209393 Marginal errors | max: 0.80321315986594 | median 0.677497553540581 | mean: 0.61937543089282 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.3260498046875 | max: 4.026123046875 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.158004893526231 | median: 0.0919692884670312 | mean: 0.0919692884670312 Nu mean: 1.0660356050592 Marginal errors | max: 0.250724014302326 | median 0.206232152124436 | mean: 0.190425623677198 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.23040771484375 | max: 3.8736572265625 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.0532774732102337 | median: 0.0308837890624999 | mean: 0.0308837890624999 Nu mean: 0.977606315852266 Marginal errors | max: 0.0916983669060123 | median 0.0765266929824948 | mean: 0.0697593208689684 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.26348876953125 | max: 3.928955078125 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.0189966206463044 | median: 0.0110199320575908 | mean: 0.0110199320575908 Nu mean: 1.00797668858871 Marginal errors | max: 0.0317151207459254 | median 0.026270214227825 | mean: 0.0240886952086962 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.2518310546875 | max: 3.9097900390625 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.00670088501353994 | median: 0.00388888101583662 | mean: 0.00388888101583662 Nu mean: 0.997187996002297 Marginal errors | max: 0.0113693316356223 | median 0.00939669372836249 | mean: 0.00860715734056932 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.2559814453125 | max: 3.9166259765625 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.00251007446998996 | median: 0.00144735987374958 | mean: 0.00144735987374958 Nu mean: 1.00106271459624 Marginal errors | max: 0.00406746973787442 | median 0.00343393462175534 | mean: 0.00313496757119527 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.25445556640625 | max: 3.9140625 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.000837027590858019 | median: 0.000488281249999889 | mean: 0.000488281249999889 Nu mean: 0.999651253659142 Marginal errors | max: 0.00143433714371533 | median 0.00116636244706747 | mean: 0.00109289166947839 Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:10] to [1:10] Estimate dispersion: DONE dispersion | min: 0.2550048828125 | max: 3.9150390625 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 1 cells batches from [1:20] to [1:20] Estimate nu: DONE nu change (abs) | max: 0.000209227688885871 | median: 0.0001220703125 | mean: 0.0001220703125 Nu mean: 1.00008715737639 Marginal errors | max: 0.000379524846207957 | median 0.000325685250686547 | mean: 0.000295532844332236 Estimate dispersion/nu: DONE While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.13167023658752 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.000274896621704102 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.000647068023681641 Estimating genes' coex Calculate genes' coex elapsed time: 0.386898994445801 Total calculations elapsed time: 1.51949119567871 Calculate genes' coex (legacy): DONE Calculate cells' coex (legacy): START Retrieving expected cells' contingency table calculating NN.. done calculating YN..NY..YY..t().. done Expected cells' contingency table elapsed time: 1.13314008712769 Calculating cells' coex normalization factor Calculate cells' normalization factor elapsed time: 0.000248908996582031 Retrieving observed cells' yes/yes contingency table calculating YY.. done Observed cells' contingency table elapsed time: 0.000579118728637695 Estimating cells' coex Calculate cells' coex elapsed time: 0.395025014877319 Total calculations elapsed time: 1.52899312973022 Calculate cells' coex (legacy): DONE Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE Calculate genes' partial coex: START Retrieving expected genes' partial contingency table calculating partial NN.. done calculating partial NY..YN..YY.. done Calculating genes' partial coex normalization factor Fraction of genes with very low expected contingency tables: 0.325 Retrieving observed genes' yes/yes partial contingency table calculating partial YY.. done Estimating genes' partial coex Calculate genes' partial coex: DONE Calculate cells' partial coex: START Retrieving expected cells' partial contingency table calculating partial NN.. done calculating partial YN..NY..YY.. done Calculating cells' partial coex normalization factor Fraction of genes with very low expected contingency tables: 0 Retrieving observed cells' yes/yes partial contingency table calculating partial YY.. done Estimating cells' partial coex Calculate cells' partial coex: DONE Asked to drop 2 genes and 0 cells Asked to drop 0 genes and 4 cells Asked to drop 2 genes and 2 cells Attaching package: 'rlang' The following objects are masked from 'package:testthat': is_false, is_null, is_true Genes/cells selection done: dropped [0] genes and [0] cells Working on [10] genes and [20] cells calculating YY.. done calculating YY.. done calculating YN..NY..NN..t().. done Estimate dispersion: START Effective number of cores used: 1 Executing 3 genes batches from [1:3] to [8:10] Estimate dispersion: DONE dispersion | min: 0.9033203125 | max: 4.6796875 | % negative: 10 calculating NN.. done calculating NN.. done calculating NY..YN..YY..t().. done Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.06967806816101 Calculating genes' coex normalization factor Fraction of genes' with very low expected contingency tables: 0.181818181818182 Calculate genes' normalization factor elapsed time: 0.000278949737548828 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.000528097152709961 Estimating genes' coex Calculate genes' coex elapsed time: 0.386295080184937 Total calculations elapsed time: 1.45678019523621 Calculate genes' coex (legacy): DONE While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.04586720466614 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.000208854675292969 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.000509977340698242 Estimating genes' coex Calculate genes' coex elapsed time: 0.362828969955444 Total calculations elapsed time: 1.40941500663757 Calculate genes' coex (legacy): DONE Calculate genes' partial coex: START Retrieving expected genes' partial contingency table calculating partial NN.. done calculating partial NY..YN..YY.. done Calculating genes' partial coex normalization factor Fraction of genes with very low expected contingency tables: 0.1 Retrieving observed genes' yes/yes partial contingency table calculating partial YY.. done Estimating genes' partial coex Calculate genes' partial coex: DONE Calculate genes' partial coex: START Retrieving expected genes' partial contingency table calculating partial NN.. done calculating partial NY..YN..YY.. done Calculating genes' partial coex normalization factor Fraction of genes with very low expected contingency tables: 0.4 Retrieving observed genes' yes/yes partial contingency table calculating partial YY.. done Estimating genes' partial coex Calculate genes' partial coex: DONE Genes/cells selection done: dropped [0] genes and [0] cells Working on [10] genes and [20] cells calculating YY.. done calculating YY.. done calculating NY..YN..NN..t().. done Estimate 'dispersion'/'nu': START Estimate dispersion: START Effective number of cores used: 1 Executing 3 genes batches from [1:3] to [8:10] Estimate dispersion: DONE dispersion | min: 0.903564453125 | max: 4.679443359375 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 4 cells batches from [1:4] to [13:16] Executing 1 cells batches from [17:20] to [17:20] Estimate nu: DONE nu change (abs) | max: 1.75719246031746 | median: 1.07229953342014 | mean: 1.07229953342014 Nu mean: 1.68489292689732 Marginal errors | max: 0.25535328937316 | median 0.0807577993228135 | mean: 0.101980750205762 Estimate dispersion: START Effective number of cores used: 1 Executing 3 genes batches from [1:3] to [8:10] Estimate dispersion: DONE dispersion | min: 1.037109375 | max: 4.6107177734375 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 4 cells batches from [1:4] to [13:16] Executing 1 cells batches from [17:20] to [17:20] Estimate nu: DONE nu change (abs) | max: 0.0273438105507502 | median: 0.0148852611818011 | mean: 0.0148852611818011 Nu mean: 1.69735147626627 Marginal errors | max: 0.00326864580272002 | median 0.00111524657842743 | mean: 0.00131556083122533 Estimate dispersion: START Effective number of cores used: 1 Executing 3 genes batches from [1:3] to [8:10] Estimate dispersion: DONE dispersion | min: 1.03887939453125 | max: 4.6097412109375 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 4 cells batches from [1:4] to [13:16] Executing 1 cells batches from [17:20] to [17:20] Estimate nu: DONE nu change (abs) | max: 0 | median: 0 | mean: 0 Nu mean: 1.69735147626627 Marginal errors | max: 7.56328383637594e-05 | median 1.72948087246994e-05 | mean: 2.99252342145451e-05 Estimate dispersion/nu: DONE calculating NN.. done calculating NN.. done calculating YN..NY..YY..t().. done Calculate cells' coex (legacy): START Retrieving expected cells' contingency table calculating NN.. done calculating YN..NY..YY..t().. done Expected cells' contingency table elapsed time: 1.10965800285339 Calculating cells' coex normalization factor Fraction of cells' with very low expected contingency tables: 0 Calculate cells' normalization factor elapsed time: 0.000301122665405273 Retrieving observed cells' yes/yes contingency table calculating YY.. done Observed cells' contingency table elapsed time: 0.000562906265258789 Estimating cells' coex Calculate cells' coex elapsed time: 0.369242906570435 Total calculations elapsed time: 1.47976493835449 Calculate cells' coex (legacy): DONE Calculate cells' partial coex: START Retrieving expected cells' partial contingency table calculating partial NN.. done calculating partial YN..NY..YY.. done Calculating cells' partial coex normalization factor Fraction of genes with very low expected contingency tables: 0 Retrieving observed cells' yes/yes partial contingency table calculating partial YY.. done Estimating cells' partial coex Calculate cells' partial coex: DONE Calculate cells' partial coex: START Retrieving expected cells' partial contingency table calculating partial NN.. done calculating partial YN..NY..YY.. done Calculating cells' partial coex normalization factor Fraction of genes with very low expected contingency tables: 0 Retrieving observed cells' yes/yes partial contingency table calculating partial YY.. done Estimating cells' partial coex Calculate cells' partial coex: DONE Genes/cells selection done: dropped [0] genes and [0] cells Working on [10] genes and [20] cells Estimate 'dispersion'/'nu': START Estimate dispersion: START Effective number of cores used: 1 Executing 3 genes batches from [1:3] to [8:10] Estimate dispersion: DONE dispersion | min: 0.903564453125 | max: 4.679443359375 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 4 cells batches from [1:4] to [13:16] Executing 1 cells batches from [17:20] to [17:20] Estimate nu: DONE nu change (abs) | max: 1.75719246031746 | median: 1.07229953342014 | mean: 1.07229953342014 Nu mean: 1.68489292689732 Marginal errors | max: 0.25535328937316 | median 0.0807577993228135 | mean: 0.101980750205762 Estimate dispersion: START Effective number of cores used: 1 Executing 3 genes batches from [1:3] to [8:10] Estimate dispersion: DONE dispersion | min: 1.037109375 | max: 4.6107177734375 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 4 cells batches from [1:4] to [13:16] Executing 1 cells batches from [17:20] to [17:20] Estimate nu: DONE nu change (abs) | max: 0.0273438105507502 | median: 0.0148852611818011 | mean: 0.0148852611818011 Nu mean: 1.69735147626627 Marginal errors | max: 0.00326864580272002 | median 0.00111524657842743 | mean: 0.00131556083122533 Estimate dispersion: START Effective number of cores used: 1 Executing 3 genes batches from [1:3] to [8:10] Estimate dispersion: DONE dispersion | min: 1.03887939453125 | max: 4.6097412109375 | % negative: 10 Estimate nu: START Effective number of cores used: 1 Executing 4 cells batches from [1:4] to [13:16] Executing 1 cells batches from [17:20] to [17:20] Estimate nu: DONE nu change (abs) | max: 0 | median: 0 | mean: 0 Nu mean: 1.69735147626627 Marginal errors | max: 7.56328383637594e-05 | median 1.72948087246994e-05 | mean: 2.99252342145451e-05 Estimate dispersion/nu: DONE Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.01087284088135 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.000220060348510742 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.000519037246704102 Estimating genes' coex Calculate genes' coex elapsed time: 0.355808019638062 Total calculations elapsed time: 1.36741995811462 Calculate genes' coex (legacy): DONE Calculating S: START Calculating S: DONE Calculating G: START calculating YY.. done calculating YN..NY..NN..t().. done calculating NN.. done calculating NY..YN..YY..t().. done Estimating G Calculating G: DONE Using S Calculating S: START Calculating S: DONE calculating PValues: START Get p-values genome wide on columns and genome wide on rows calculating PValues: DONE Using G Calculating G: START calculating YY.. done calculating YN..NY..NN..t().. done calculating NN.. done calculating NY..YN..YY..t().. done Estimating G Calculating G: DONE calculating PValues: START Get p-values on a set of genes on columns and on a set of genes on rows calculating PValues: DONE Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE Calculate GDI dataframe: START Using G Calculating G: START calculating YY.. done calculating YN..NY..NN..t().. done calculating NN.. done calculating NY..YN..YY..t().. done Estimating G Calculating G: DONE S matrix sorted Calculate GDI dataframe: DONE Calculating S: START Calculating S: DONE S matrix sorted S matrix sorted Initializing `COTAN` meta-data Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [1200] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0 Only analysis elapsed time: 1.32460713386536 Cotan genes' coex estimation started Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.61635899543762 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00331902503967285 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0227189064025879 Estimating genes' coex Calculate genes' coex elapsed time: 0.340723991394043 Total calculations elapsed time: 1.98312091827393 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.33641695976257 Total elapsed time: 4.01845407485962 Initializing `COTAN` meta-data Condition test n cells 1200 Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [1200] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0 Only analysis elapsed time: 1.34596490859985 Cotan genes' coex estimation started Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.61781883239746 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00303196907043457 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0222330093383789 Estimating genes' coex Calculate genes' coex elapsed time: 0.339533090591431 Total calculations elapsed time: 1.9826169013977 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.35552906990051 Total elapsed time: 4.05354714393616 Using S Calculating S: START Calculating S: DONE calculating PValues: START Get p-values on a set of genes on columns and genome wide on rows calculating PValues: DONE Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE Initializing `COTAN` meta-data Condition test n cells 1200 Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [1200] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0 Only analysis elapsed time: 1.35616898536682 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.57804417610168 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00283598899841309 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0223088264465332 Estimating genes' coex Calculate genes' coex elapsed time: 0.347666025161743 Total calculations elapsed time: 1.95085501670837 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.26392793655396 Total elapsed time: 3.95630502700806 Using S Calculating S: START Calculating S: DONE calculating PValues: START Get p-values on a set of genes on columns and genome wide on rows calculating PValues: DONE Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE Initializing `COTAN` meta-data Condition test n cells 1200 Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [1200] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0 Only analysis elapsed time: 1.30555701255798 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.59665489196777 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00460910797119141 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0240678787231445 Estimating genes' coex Calculate genes' coex elapsed time: 0.370126962661743 Total calculations elapsed time: 1.99545884132385 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.38195490837097 Total elapsed time: 4.05131793022156 Using S Calculating S: START Calculating S: DONE calculating PValues: START Get p-values on a set of genes on columns and genome wide on rows calculating PValues: DONE Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE Initializing `COTAN` meta-data Asked to drop 0 genes and 1160 cells Cotan analysis functions started Asked to drop 31 genes and 0 cells Genes/cells selection done: dropped [31] genes and [0] cells Working on [569] genes and [40] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:569] to [1:569] Estimate dispersion: DONE dispersion | min: -0.2119140625 | max: 85.5 | % negative: 62.2144112478032 Cotan genes' coex estimation not requested Total elapsed time: 1.69538593292236 Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.02721405029297 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00330185890197754 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00267815589904785 Estimating genes' coex Calculate genes' coex elapsed time: 0.334754943847656 Total calculations elapsed time: 1.36794900894165 Calculate genes' coex (legacy): DONE While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.10032200813293 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00295591354370117 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00291109085083008 Estimating genes' coex Calculate genes' coex elapsed time: 0.332762002944946 Total calculations elapsed time: 1.43895101547241 Calculate genes' coex (legacy): DONE While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.09233713150024 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00284004211425781 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00256180763244629 Estimating genes' coex Calculate genes' coex elapsed time: 0.341526985168457 Total calculations elapsed time: 1.43926596641541 Calculate genes' coex (legacy): DONE Hangling COTAN object with condition: test Using S Calculating S: START Calculating S: DONE calculating PValues: START Get p-values on a set of genes on columns and on a set of genes on rows calculating PValues: DONE Handling genes type: G2 Handling genes type: G3 Initializing `COTAN` meta-data Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [1200] cells PCA: START PCA: DONE Hierarchical clustering: START Hierarchical clustering: DONE Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0 Cotan genes' coex estimation not requested Total elapsed time: 3.57572889328003 Saving elaborated data locally at: /var/folders/r0/l4fjk6cj5xj0j3brt4bplpl40000gt/T//RtmpHONuQh/test.cotan.RDS Creating cells' uniform clustering: START In iteration 1 the number of cells to re-cluster is 1200 cells belonging to 0 clusters Creating Seurat object: START Normalizing layer: counts Performing log-normalization 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| Finding variable features for layer counts Calculating gene variances 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| Calculating feature variances of standardized and clipped values 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| Centering and scaling data matrix | | | 0% | |======================================================================| 100% PC_ 1 Positive: g-000558, g-000570, g-000499, g-000504, g-000546, g-000506, g-000503, g-000517, g-000596, g-000528 g-000527, g-000580, g-000592, g-000578, g-000509, g-000488, g-000555, g-000577, g-000534, g-000583 g-000598, g-000535, g-000512, g-000554, g-000519, g-000525, g-000548, g-000544, g-000502, g-000541 Negative: g-000133, g-000007, g-000074, g-000141, g-000057, g-000235, g-000170, g-000019, g-000195, g-000140 g-000183, g-000031, g-000046, g-000178, g-000177, g-000161, g-000157, g-000139, g-000011, g-000135 g-000125, g-000208, g-000061, g-000085, g-000204, g-000104, g-000237, g-000004, g-000038, g-000128 PC_ 2 Positive: g-000039, g-000050, g-000175, g-000078, g-000116, g-000189, g-000135, g-000047, g-000072, g-000087 g-000063, g-000235, g-000066, g-000109, g-000018, g-000074, g-000231, g-000136, g-000034, g-000207 g-000128, g-000167, g-000171, g-000049, g-000182, g-000013, g-000054, g-000062, g-000240, g-000158 Negative: g-000584, g-000583, g-000544, g-000519, g-000575, g-000516, g-000585, g-000486, g-000489, g-000539 g-000484, g-000502, g-000523, g-000595, g-000305, g-000574, g-000599, g-000589, g-000509, g-000538 g-000526, g-000551, g-000579, g-000590, g-000445, g-000556, g-000543, g-000501, g-000504, g-000570 PC_ 3 Positive: g-000015, g-000575, g-000483, g-000316, g-000025, g-000364, g-000050, g-000278, g-000443, g-000360 g-000332, g-000124, g-000212, g-000387, g-000536, g-000252, g-000251, g-000321, g-000501, g-000470 g-000582, g-000106, g-000455, g-000368, g-000081, g-000104, g-000437, g-000288, g-000386, g-000317 Negative: g-000211, g-000337, g-000129, g-000185, g-000397, g-000403, g-000253, g-000098, g-000390, g-000303 g-000052, g-000088, g-000463, g-000468, g-000236, g-000209, g-000005, g-000375, g-000342, g-000262 g-000388, g-000091, g-000413, g-000285, g-000003, g-000095, g-000142, g-000205, g-000432, g-000241 PC_ 4 Positive: g-000379, g-000193, g-000212, g-000434, g-000593, g-000513, g-000177, g-000223, g-000069, g-000131 g-000162, g-000345, g-000462, g-000484, g-000448, g-000229, g-000365, g-000302, g-000010, g-000366 g-000051, g-000535, g-000269, g-000270, g-000155, g-000529, g-000373, g-000008, g-000393, g-000306 Negative: g-000334, g-000398, g-000292, g-000095, g-000097, g-000202, g-000382, g-000195, g-000007, g-000079 g-000086, g-000240, g-000263, g-000317, g-000576, g-000557, g-000160, g-000154, g-000214, g-000228 g-000313, g-000053, g-000524, g-000374, g-000568, g-000188, g-000358, g-000528, g-000362, g-000150 PC_ 5 Positive: g-000518, g-000108, g-000186, g-000170, g-000401, g-000337, g-000047, g-000599, g-000432, g-000578 g-000042, g-000065, g-000493, g-000261, g-000533, g-000256, g-000560, g-000596, g-000368, g-000381 g-000535, g-000338, g-000215, g-000159, g-000365, g-000234, g-000173, g-000387, g-000225, g-000272 Negative: g-000451, g-000339, g-000295, g-000328, g-000544, g-000061, g-000227, g-000391, g-000556, g-000237 g-000067, g-000165, g-000449, g-000591, g-000087, g-000129, g-000197, g-000203, g-000487, g-000505 g-000333, g-000029, g-000271, g-000064, g-000583, g-000156, g-000448, g-000153, g-000526, g-000393 Computing nearest neighbor graph Computing SNN Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck Number of nodes: 1200 Number of edges: 55489 Running Louvain algorithm with multilevel refinement... 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| Maximum modularity in 10 random starts: 0.5973 Number of communities: 4 Elapsed time: 0 seconds Used resolution for Seurat clusterization is: 0.8 20:59:56 UMAP embedding parameters a = 0.9922 b = 1.112 20:59:56 Read 1200 rows and found 50 numeric columns 20:59:56 Using Annoy for neighbor search, n_neighbors = 30 20:59:56 Building Annoy index with metric = cosine, n_trees = 50 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| 20:59:56 Writing NN index file to temp file /var/folders/r0/l4fjk6cj5xj0j3brt4bplpl40000gt/T//RtmpHONuQh/file1207520d3334d 20:59:56 Searching Annoy index using 1 thread, search_k = 3000 20:59:57 Annoy recall = 100% 20:59:57 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30 20:59:58 Initializing from normalized Laplacian + noise (using RSpectra) 20:59:58 Commencing optimization for 500 epochs, with 42274 positive edges Using method 'umap' 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| 20:59:59 Optimization finished Creating Seurat object: DONE * checking uniformity of cluster '0' of 4 clusters Asked to drop 0 genes and 847 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [353] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.038818359375 | max: 11.109375 | % negative: 5 Only analysis elapsed time: 1.10217094421387 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.17173099517822 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0027930736541748 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.016963005065918 Estimating genes' coex Calculate genes' coex elapsed time: 0.345848083496094 Total calculations elapsed time: 1.53733515739441 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 1.89135909080505 Total elapsed time: 3.33467602729797 Checking uniformity for the cluster '01_0000' with 353 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 01_0000, with size 353, is uniform 0.83% of the genes is above the given GDI threshold 1.46 highest 1% GDI quantile is at 1.4439 cluster 01_0000 is uniform * checking uniformity of cluster '1' of 4 clusters Asked to drop 0 genes and 879 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [321] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.0501708984375 | max: 14.515625 | % negative: 8.5 Only analysis elapsed time: 1.10461091995239 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.13993787765503 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00388312339782715 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.016010046005249 Estimating genes' coex Calculate genes' coex elapsed time: 0.349992990493774 Total calculations elapsed time: 1.50982403755188 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 1.90003800392151 Total elapsed time: 3.35067296028137 Checking uniformity for the cluster '01_0001' with 321 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 01_0001, with size 321, is uniform 0.5% of the genes is above the given GDI threshold 1.46 highest 1% GDI quantile is at 1.4305 cluster 01_0001 is uniform * checking uniformity of cluster '2' of 4 clusters Asked to drop 0 genes and 905 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [295] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.054443359375 | max: 104 | % negative: 36.5 Only analysis elapsed time: 1.18277406692505 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.2152099609375 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0036320686340332 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00765299797058105 Estimating genes' coex Calculate genes' coex elapsed time: 0.396613836288452 Total calculations elapsed time: 1.62310886383057 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.00606203079224 Total elapsed time: 3.55930089950562 Checking uniformity for the cluster '01_0002' with 295 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 01_0002, with size 295, is uniform 0% of the genes is above the given GDI threshold 1.46 highest 1% GDI quantile is at 1.3925 cluster 01_0002 is uniform * checking uniformity of cluster '3' of 4 clusters Asked to drop 0 genes and 969 cells Cotan analysis functions started Asked to drop 1 genes and 0 cells Genes/cells selection done: dropped [1] genes and [0] cells Working on [599] genes and [231] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:599] to [1:599] Estimate dispersion: DONE dispersion | min: -0.0662841796875 | max: 82.5 | % negative: 33.889816360601 Only analysis elapsed time: 1.16943502426147 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.30145001411438 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00429296493530273 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00952506065368652 Estimating genes' coex Calculate genes' coex elapsed time: 0.415613889694214 Total calculations elapsed time: 1.73088192939758 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.1038830280304 Total elapsed time: 3.63622903823853 Checking uniformity for the cluster '01_0003' with 231 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 01_0003, with size 231, is uniform 0% of the genes is above the given GDI threshold 1.46 highest 1% GDI quantile is at 1.373 cluster 01_0003 is uniform Found 4 uniform and 0 non-uniform clusters NO new possible uniform clusters! Unclustered cell left: 0 The final raw clusterization contains [ 4 ] different clusters: 01_0000, 01_0001, 01_0002, 01_0003 Differential Expression Analysis - START * DEA on cluster '1' with 353 cells * DEA on cluster '2' with 321 cells * DEA on cluster '3' with 295 cells * DEA on cluster '4' with 231 cells Differential Expression Analysis - DONE Applied reordering to clusterization is: 1 -> 3, 2 -> 1, 3 -> 2, 4 -> 4 Creating cells' uniform clustering: DONE Applied reordering to clusterization is: 1 -> 1, 2 -> 2, 3 -> 3, 4 -> 4 Asked to drop 0 genes and 905 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [295] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.054443359375 | max: 104 | % negative: 36.5 Only analysis elapsed time: 1.20586800575256 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.22621393203735 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0038149356842041 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00797510147094727 Estimating genes' coex Calculate genes' coex elapsed time: 0.355824947357178 Total calculations elapsed time: 1.59382891654968 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 1.95308804512024 Total elapsed time: 3.53132200241089 Checking uniformity for the cluster 'Cluster_2' with 295 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster Cluster_2, with size 295, is uniform 0% of the genes is above the given GDI threshold 1.46 highest 1% GDI quantile is at 1.3925 Differential Expression Analysis - START * DEA on cluster '-1' with 50 cells * DEA on cluster '1' with 321 cells * DEA on cluster '2' with 245 cells * DEA on cluster '3' with 353 cells * DEA on cluster '4' with 231 cells Differential Expression Analysis - DONE Applied reordering to clusterization is: 1 -> 2, 2 -> 1, 3 -> 4, 4 -> 3, -1 -> -1 Applied reordering to clusterization is: 1 -> 2, 2 -> 3, 3 -> 1, 4 -> 4, -1 -> -1 Creating cells' uniform clustering: START In iteration 1 the number of cells to re-cluster is 1200 cells belonging to 0 clusters Creating Seurat object: START Normalizing layer: counts Performing log-normalization 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| Finding variable features for layer counts Calculating gene variances 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| Calculating feature variances of standardized and clipped values 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| Centering and scaling data matrix | | | 0% | |======================================================================| 100% PC_ 1 Positive: g-000558, g-000570, g-000499, g-000504, g-000546, g-000506, g-000503, g-000517, g-000596, g-000528 g-000527, g-000580, g-000592, g-000578, g-000509, g-000488, g-000555, g-000577, g-000534, g-000583 g-000598, g-000535, g-000512, g-000554, g-000519, g-000525, g-000548, g-000544, g-000502, g-000541 Negative: g-000133, g-000007, g-000074, g-000141, g-000057, g-000235, g-000170, g-000019, g-000195, g-000140 g-000183, g-000031, g-000046, g-000178, g-000177, g-000161, g-000157, g-000139, g-000011, g-000135 g-000125, g-000208, g-000061, g-000085, g-000204, g-000104, g-000237, g-000004, g-000038, g-000128 PC_ 2 Positive: g-000039, g-000050, g-000175, g-000078, g-000116, g-000189, g-000135, g-000047, g-000072, g-000087 g-000063, g-000235, g-000066, g-000109, g-000018, g-000074, g-000231, g-000136, g-000034, g-000207 g-000128, g-000167, g-000171, g-000049, g-000182, g-000013, g-000054, g-000062, g-000240, g-000158 Negative: g-000584, g-000583, g-000544, g-000519, g-000575, g-000516, g-000585, g-000486, g-000489, g-000539 g-000484, g-000502, g-000523, g-000595, g-000305, g-000574, g-000599, g-000589, g-000509, g-000538 g-000526, g-000551, g-000579, g-000590, g-000445, g-000556, g-000543, g-000501, g-000504, g-000570 PC_ 3 Positive: g-000015, g-000575, g-000483, g-000316, g-000025, g-000364, g-000050, g-000278, g-000443, g-000360 g-000332, g-000124, g-000212, g-000387, g-000536, g-000252, g-000251, g-000321, g-000501, g-000470 g-000582, g-000106, g-000455, g-000368, g-000081, g-000104, g-000437, g-000288, g-000386, g-000317 Negative: g-000211, g-000337, g-000129, g-000185, g-000397, g-000403, g-000253, g-000098, g-000390, g-000303 g-000052, g-000088, g-000463, g-000468, g-000236, g-000209, g-000005, g-000375, g-000342, g-000262 g-000388, g-000091, g-000413, g-000285, g-000003, g-000095, g-000142, g-000205, g-000432, g-000241 PC_ 4 Positive: g-000379, g-000193, g-000212, g-000434, g-000593, g-000513, g-000177, g-000223, g-000069, g-000131 g-000162, g-000345, g-000462, g-000484, g-000448, g-000229, g-000365, g-000302, g-000010, g-000366 g-000051, g-000535, g-000269, g-000270, g-000155, g-000529, g-000373, g-000008, g-000393, g-000306 Negative: g-000334, g-000398, g-000292, g-000095, g-000097, g-000202, g-000382, g-000195, g-000007, g-000079 g-000086, g-000240, g-000263, g-000317, g-000576, g-000557, g-000160, g-000154, g-000214, g-000228 g-000313, g-000053, g-000524, g-000374, g-000568, g-000188, g-000358, g-000528, g-000362, g-000150 PC_ 5 Positive: g-000518, g-000108, g-000186, g-000170, g-000401, g-000337, g-000047, g-000599, g-000432, g-000578 g-000042, g-000065, g-000493, g-000261, g-000533, g-000256, g-000560, g-000596, g-000368, g-000381 g-000535, g-000338, g-000215, g-000159, g-000365, g-000234, g-000173, g-000387, g-000225, g-000272 Negative: g-000451, g-000339, g-000295, g-000328, g-000544, g-000061, g-000227, g-000391, g-000556, g-000237 g-000067, g-000165, g-000449, g-000591, g-000087, g-000129, g-000197, g-000203, g-000487, g-000505 g-000333, g-000029, g-000271, g-000064, g-000583, g-000156, g-000448, g-000153, g-000526, g-000393 Computing nearest neighbor graph Computing SNN Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck Number of nodes: 1200 Number of edges: 55489 Running Louvain algorithm with multilevel refinement... 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| Maximum modularity in 10 random starts: 0.5973 Number of communities: 4 Elapsed time: 0 seconds Used resolution for Seurat clusterization is: 0.8 21:00:31 UMAP embedding parameters a = 0.9922 b = 1.112 21:00:31 Read 1200 rows and found 50 numeric columns 21:00:31 Using Annoy for neighbor search, n_neighbors = 30 21:00:31 Building Annoy index with metric = cosine, n_trees = 50 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| 21:00:31 Writing NN index file to temp file /var/folders/r0/l4fjk6cj5xj0j3brt4bplpl40000gt/T//RtmpHONuQh/file1207510916398 21:00:31 Searching Annoy index using 1 thread, search_k = 3000 21:00:31 Annoy recall = 100% 21:00:32 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30 21:00:33 Initializing from normalized Laplacian + noise (using RSpectra) 21:00:33 Commencing optimization for 500 epochs, with 42274 positive edges Using method 'umap' 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| 21:00:34 Optimization finished Creating Seurat object: DONE Using passed in clusterization * checking uniformity of cluster '1' of 2 clusters Asked to drop 0 genes and 600 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [600] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.0386962890625 | max: 19.40625 | % negative: 6.5 Only analysis elapsed time: 1.24940013885498 Cotan genes' coex estimation started Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.38002395629883 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00346088409423828 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0181970596313477 Estimating genes' coex Calculate genes' coex elapsed time: 0.401204109191895 Total calculations elapsed time: 1.80288600921631 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.19585585594177 Total elapsed time: 3.82551884651184 Checking uniformity for the cluster '01_0001' with 600 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 01_0001, with size 600, is uniform 1% of the genes is above the given GDI threshold 1.46 highest 1% GDI quantile is at 1.4534 cluster 01_0001 is uniform * checking uniformity of cluster '2' of 2 clusters Asked to drop 0 genes and 600 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [600] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.032958984375 | max: 10.0859375 | % negative: 3.66666666666667 Only analysis elapsed time: 1.24108004570007 Cotan genes' coex estimation started Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.43032383918762 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00355815887451172 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0209739208221436 Estimating genes' coex Calculate genes' coex elapsed time: 0.370419025421143 Total calculations elapsed time: 1.82527494430542 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.21545815467834 Total elapsed time: 3.86642813682556 Checking uniformity for the cluster '01_0002' with 600 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 01_0002, with size 600, is uniform 0.5% of the genes is above the given GDI threshold 1.46 highest 1% GDI quantile is at 1.4378 cluster 01_0002 is uniform Found 2 uniform and 0 non-uniform clusters NO new possible uniform clusters! Unclustered cell left: 0 The final raw clusterization contains [ 2 ] different clusters: 01_0001, 01_0002 Differential Expression Analysis - START * DEA on cluster '1' with 600 cells * DEA on cluster '2' with 600 cells Differential Expression Analysis - DONE Applied reordering to clusterization is: 1 -> 1, 2 -> 2 Creating cells' uniform clustering: DONE findClustersMarkers - START clustersDeltaExpression - START Handling cluster '1' with mean UDE 1.20998115359079 Handling cluster '2' with mean UDE 0.538244816501365 Handling cluster '3' with mean UDE 1.43530796540673 Handling cluster '4' with mean UDE 0.632684489354433 clustersDeltaExpression - DONE Log Fold Change Analysis - START * Analysis of cluster: '1' * Analysis of cluster: '2' * Analysis of cluster: '3' * Analysis of cluster: '4' Log Fold Change Analysis - DONE findClustersMarkers - DONE findClustersMarkers - START clustersDeltaExpression - START Handling cluster '1' with mean UDE 1.20998115359079 Handling cluster '2' with mean UDE 0.538244816501365 Handling cluster '3' with mean UDE 1.43530796540673 Handling cluster '4' with mean UDE 0.632684489354433 clustersDeltaExpression - DONE Log Fold Change Analysis - START * Analysis of cluster: '1' * Analysis of cluster: '2' * Analysis of cluster: '3' * Analysis of cluster: '4' Log Fold Change Analysis - DONE findClustersMarkers - DONE findClustersMarkers - START Differential Expression Analysis - START * DEA on cluster '1' with 321 cells * DEA on cluster '2' with 295 cells * DEA on cluster '3' with 353 cells * DEA on cluster '4' with 231 cells Differential Expression Analysis - DONE clustersDeltaExpression - START Handling cluster '1' with mean UDE 1.20998115359079 Handling cluster '2' with mean UDE 0.538244816501365 Handling cluster '3' with mean UDE 1.43530796540673 Handling cluster '4' with mean UDE 0.632684489354433 clustersDeltaExpression - DONE Log Fold Change Analysis - START * Analysis of cluster: '1' * Analysis of cluster: '2' * Analysis of cluster: '3' * Analysis of cluster: '4' Log Fold Change Analysis - DONE findClustersMarkers - DONE [1] "1" Asked to drop 0 genes and 879 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [321] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.0501708984375 | max: 14.515625 | % negative: 8.5 Only analysis elapsed time: 1.14441108703613 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.25949120521545 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00286293029785156 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0161068439483643 Estimating genes' coex Calculate genes' coex elapsed time: 0.388387203216553 Total calculations elapsed time: 1.66684818267822 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.07299900054932 Total elapsed time: 3.60918712615967 Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE [1] "3" Asked to drop 0 genes and 847 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [353] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.038818359375 | max: 11.109375 | % negative: 5 Only analysis elapsed time: 1.15467190742493 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.25313806533813 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00360202789306641 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0190219879150391 Estimating genes' coex Calculate genes' coex elapsed time: 0.382344007492065 Total calculations elapsed time: 1.65810608863831 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.06196808815002 Total elapsed time: 3.59350109100342 Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE Genes/cells selection done: dropped [0] genes and [0] cells Working on [10] genes and [20] cells Initializing `COTAN` meta-data Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [1200] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0 Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [1200] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0 Only analysis elapsed time: 1.39376902580261 Cotan genes' coex estimation started Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.61140704154968 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00397300720214844 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0225419998168945 Estimating genes' coex Calculate genes' coex elapsed time: 0.354851007461548 Total calculations elapsed time: 1.99277305603027 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.36424112319946 Total elapsed time: 4.10208296775818 Calculating gene co-expression space - START Using S Calculating S: START Calculating S: DONE calculating PValues: START Get p-values on a set of genes on columns and genome wide on rows calculating PValues: DONE Number of selected secondary markers: 109 Calculating S: START Calculating S: DONE S matrix sorted Number of columns (V set - secondary markers): 109 Number of rows (U set): 60 Calculating gene co-expression space - DONE Establishing gene clusters - START Calculating gene co-expression space - START Using S Calculating S: START Calculating S: DONE calculating PValues: START Get p-values on a set of genes on columns and genome wide on rows calculating PValues: DONE Number of selected secondary markers: 109 Calculating S: START Calculating S: DONE S matrix sorted Number of columns (V set - secondary markers): 109 Number of rows (U set): 60 Calculating gene co-expression space - DONE Establishing gene clusters - DONE Initializing `COTAN` meta-data Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [1200] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0 Cotan genes' coex estimation not requested Total elapsed time: 2.1904399394989 Differential Expression Analysis - START * DEA on cluster '1' with 353 cells * DEA on cluster '2' with 221 cells * DEA on cluster '3' with 315 cells * DEA on cluster '4' with 311 cells Differential Expression Analysis - DONE Log Fold Change Analysis - START * Analysis of cluster: '1' * Analysis of cluster: '2' * Analysis of cluster: '3' * Analysis of cluster: '4' Log Fold Change Analysis - DONE clustersDeltaExpression - START Handling cluster '1' with mean UDE 1.43530796540673 Handling cluster '2' with mean UDE 0.640931107489518 Handling cluster '3' with mean UDE 0.546546914955574 Handling cluster '4' with mean UDE 1.22034802329656 clustersDeltaExpression - DONE In group G1 there are 3 detected over 3 genes In group G2 there are 2 detected over 2 genes In group G3 there are 5 detected over 5 genes Merging cells' uniform clustering: START Start merging nearest clusters: threshold 1.37 Start merging nearest clusters: iteration 1 Differential Expression Analysis - START * DEA on cluster '1' with 353 cells * DEA on cluster '2' with 221 cells * DEA on cluster '3' with 315 cells * DEA on cluster '4' with 311 cells Differential Expression Analysis - DONE New clusters pairs to be tested for merging: c("1", "2") c("3", "4") *1_2-merge Asked to drop 0 genes and 626 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [574] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.0340576171875 | max: 10.4375 | % negative: 4.33333333333333 Only analysis elapsed time: 1.20503401756287 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.3079400062561 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00294899940490723 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0190050601959229 Estimating genes' coex Calculate genes' coex elapsed time: 0.364534854888916 Total calculations elapsed time: 1.69442892074585 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.06427192687988 Total elapsed time: 3.6348819732666 Checking uniformity for the cluster '1_2-merge' with 574 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 1_2-merge, with size 574, is not uniform 6.83% of the genes is above the given GDI threshold 1.37 highest 1% GDI quantile is at 1.4297 Clusters 1 and 2 cannot be merged *3_4-merge Asked to drop 0 genes and 574 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [626] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.03839111328125 | max: 17.40625 | % negative: 6.66666666666667 Only analysis elapsed time: 1.20072603225708 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.32352519035339 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00325298309326172 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0181488990783691 Estimating genes' coex Calculate genes' coex elapsed time: 0.333503007888794 Total calculations elapsed time: 1.67843008041382 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.04006099700928 Total elapsed time: 3.6287841796875 Checking uniformity for the cluster '3_4-merge' with 626 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 3_4-merge, with size 626, is not uniform 7.17% of the genes is above the given GDI threshold 1.37 highest 1% GDI quantile is at 1.4544 Clusters 3 and 4 cannot be merged Finished the first batch - no merges were executed Start merging nearest clusters: iteration 2 Differential Expression Analysis - START * DEA on cluster '1' with 353 cells * DEA on cluster '2' with 221 cells * DEA on cluster '3' with 315 cells * DEA on cluster '4' with 311 cells Differential Expression Analysis - DONE New clusters pairs to be tested for merging: c("2", "3") c("2", "4") *2_3-merge Asked to drop 0 genes and 664 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [536] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.0390625 | max: 12.59375 | % negative: 9.16666666666667 Only analysis elapsed time: 1.16221189498901 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.29215312004089 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0030820369720459 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0137438774108887 Estimating genes' coex Calculate genes' coex elapsed time: 0.354797124862671 Total calculations elapsed time: 1.6637761592865 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.01688003540039 Total elapsed time: 3.54167699813843 Checking uniformity for the cluster '2_3-merge' with 536 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 2_3-merge, with size 536, is not uniform 10.17% of the genes is above the given GDI threshold 1.37 highest 1% GDI quantile is at 1.4902 Clusters 2 and 3 cannot be merged *2_4-merge Asked to drop 0 genes and 668 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [532] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.04754638671875 | max: 10.234375 | % negative: 5 Only analysis elapsed time: 1.24563217163086 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.4274799823761 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00409483909606934 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0183320045471191 Estimating genes' coex Calculate genes' coex elapsed time: 0.42494797706604 Total calculations elapsed time: 1.87485480308533 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.30237698554993 Total elapsed time: 3.97582316398621 Checking uniformity for the cluster '2_4-merge' with 532 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 2_4-merge, with size 532, is not uniform 22.5% of the genes is above the given GDI threshold 1.37 highest 1% GDI quantile is at 1.5909 Clusters 2 and 4 cannot be merged None of the 4 tested cluster pairs could be merged Executed all merges for threshold 1.37 Start merging nearest clusters: threshold 1.4 Start merging nearest clusters: iteration 3 Differential Expression Analysis - START * DEA on cluster '1' with 353 cells * DEA on cluster '2' with 221 cells * DEA on cluster '3' with 315 cells * DEA on cluster '4' with 311 cells Differential Expression Analysis - DONE New clusters pairs to be tested for merging: c("1", "3") c("1", "4") *1_3-merge Asked to drop 0 genes and 532 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [668] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.046630859375 | max: 9.03515625 | % negative: 1.33333333333333 Only analysis elapsed time: 1.27050805091858 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.41300582885742 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00327014923095703 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0206840038299561 Estimating genes' coex Calculate genes' coex elapsed time: 0.351626873016357 Total calculations elapsed time: 1.78858685493469 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.10658001899719 Total elapsed time: 3.77152585983276 Checking uniformity for the cluster '1_3-merge' with 668 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 1_3-merge, with size 668, is not uniform 32.67% of the genes is above the given GDI threshold 1.4 highest 1% GDI quantile is at 1.657 Clusters 1 and 3 cannot be merged *1_4-merge Asked to drop 0 genes and 536 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [664] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.0164794921875 | max: 5.5595703125 | % negative: 0.166666666666667 Only analysis elapsed time: 1.25474500656128 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.37300491333008 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00282597541809082 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0193600654602051 Estimating genes' coex Calculate genes' coex elapsed time: 0.364516973495483 Total calculations elapsed time: 1.75970792770386 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.09194684028625 Total elapsed time: 3.7419159412384 Checking uniformity for the cluster '1_4-merge' with 664 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 1_4-merge, with size 664, is not uniform 34.17% of the genes is above the given GDI threshold 1.4 highest 1% GDI quantile is at 1.5911 Clusters 1 and 4 cannot be merged None of the 6 tested cluster pairs could be merged Executed all merges for threshold 1.4 Start merging nearest clusters: threshold 1.43 Start merging nearest clusters: iteration 4 Differential Expression Analysis - START * DEA on cluster '1' with 353 cells * DEA on cluster '2' with 221 cells * DEA on cluster '3' with 315 cells * DEA on cluster '4' with 311 cells Differential Expression Analysis - DONE New clusters pairs to be tested for merging: Clusters 1 and 2 will be merged Executed 1 merges out of potentially 6 Start merging nearest clusters: iteration 5 Differential Expression Analysis - START * DEA on cluster '1_2-merge' with 574 cells * DEA on cluster '3' with 315 cells * DEA on cluster '4' with 311 cells Differential Expression Analysis - DONE New clusters pairs to be tested for merging: c("1_2-merge", "3") c("1_2-merge", "4") *1_2-merge_3-merge Asked to drop 0 genes and 311 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [889] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.04052734375 | max: 9.234375 | % negative: 1.16666666666667 Only analysis elapsed time: 1.35753583908081 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.53120493888855 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00438094139099121 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0199089050292969 Estimating genes' coex Calculate genes' coex elapsed time: 0.388054132461548 Total calculations elapsed time: 1.94354891777039 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.31879997253418 Total elapsed time: 4.06150698661804 Checking uniformity for the cluster '1_2-merge_3-merge' with 889 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 1_2-merge_3-merge, with size 889, is not uniform 15.33% of the genes is above the given GDI threshold 1.43 highest 1% GDI quantile is at 1.5775 Clusters 1_2-merge and 3 cannot be merged *1_2-merge_4-merge Asked to drop 0 genes and 315 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [885] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.01153564453125 | max: 6.0224609375 | % negative: 0.166666666666667 Only analysis elapsed time: 1.34677600860596 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.52490711212158 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00350594520568848 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0224239826202393 Estimating genes' coex Calculate genes' coex elapsed time: 0.398932933807373 Total calculations elapsed time: 1.94976997375488 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.34554386138916 Total elapsed time: 4.04590106010437 Checking uniformity for the cluster '1_2-merge_4-merge' with 885 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 1_2-merge_4-merge, with size 885, is not uniform 22.67% of the genes is above the given GDI threshold 1.43 highest 1% GDI quantile is at 1.6496 Clusters 1_2-merge and 4 cannot be merged None of the 8 tested cluster pairs could be merged Executed all merges for threshold 1.43 Start merging nearest clusters: threshold 1.46 Start merging nearest clusters: iteration 6 Differential Expression Analysis - START * DEA on cluster '1_2-merge' with 574 cells * DEA on cluster '3' with 315 cells * DEA on cluster '4' with 311 cells Differential Expression Analysis - DONE New clusters pairs to be tested for merging: Clusters 3 and 4 will be merged Executed 1 merges out of potentially 8 Start merging nearest clusters: iteration 7 Differential Expression Analysis - START * DEA on cluster '1_2-merge' with 574 cells * DEA on cluster '3_4-merge' with 626 cells Differential Expression Analysis - DONE New clusters pairs to be tested for merging: c("1_2-merge", "3_4-merge") *1_2-merge_3_4-merge-merge Asked to drop no genes or cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [1200] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0 Only analysis elapsed time: 1.34140086174011 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.64638996124268 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00340414047241211 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0226988792419434 Estimating genes' coex Calculate genes' coex elapsed time: 0.393362998962402 Total calculations elapsed time: 2.06585597991943 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.45708608627319 Total elapsed time: 4.14800500869751 Checking uniformity for the cluster '1_2-merge_3_4-merge-merge' with 1200 cells Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE GDI plot Removed 0 low GDI genes (such as the fully-expressed) in GDI plot Cluster 1_2-merge_3_4-merge-merge, with size 1200, is not uniform 22.17% of the genes is above the given GDI threshold 1.46 highest 1% GDI quantile is at 1.6405 Clusters 1_2-merge and 3_4-merge cannot be merged None of the 9 tested cluster pairs could be merged Executed all merges for threshold 1.46 The final merged clusterization contains [2] different clusters: 1_2-merge, 3_4-merge Differential Expression Analysis - START * DEA on cluster '1' with 574 cells * DEA on cluster '2' with 626 cells Differential Expression Analysis - DONE Applied reordering to clusterization is: 1 -> 1, 2 -> 2 Merging cells' uniform clustering: DONE Applied reordering to clusterization is: 1 -> 1, 2 -> 2 Log Fold Change Analysis - START * Analysis of cluster: '1' * Analysis of cluster: '2' Log Fold Change Analysis - DONE Asked to drop 0 genes and 626 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [574] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.0340576171875 | max: 10.4375 | % negative: 4.33333333333333 Only analysis elapsed time: 1.17637515068054 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.29287505149841 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00297880172729492 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0188381671905518 Estimating genes' coex Calculate genes' coex elapsed time: 0.370461940765381 Total calculations elapsed time: 1.68515396118164 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 2.0523829460144 Total elapsed time: 3.55089497566223 Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE Asked to drop 0 genes and 574 cells Cotan analysis functions started Genes/cells selection done: dropped [0] genes and [0] cells Working on [600] genes and [626] cells Estimate dispersion: START Effective number of cores used: 1 Executing 1 genes batches from [1:600] to [1:600] Estimate dispersion: DONE dispersion | min: -0.03839111328125 | max: 17.40625 | % negative: 6.66666666666667 Only analysis elapsed time: 1.15474200248718 Cotan genes' coex estimation started While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet Calculate genes' coex (legacy): START Retrieving expected genes' contingency table calculating NN.. done calculating NY..YN..YY..t().. done Expected genes' contingency table elapsed time: 1.30393409729004 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00330495834350586 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0183851718902588 Estimating genes' coex Calculate genes' coex elapsed time: 0.32604193687439 Total calculations elapsed time: 1.65166616439819 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 1.99803614616394 Total elapsed time: 3.51447606086731 Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE [ FAIL 0 | WARN 0 | SKIP 0 | PASS 452 ] > > proc.time() user system elapsed 293.859 2.024 295.875
COTAN.Rcheck/COTAN-Ex.timings
name | user | system | elapsed | |
COTAN-Legacy | 0.001 | 0.001 | 0.001 | |
COTANObjectCreation | 4.367 | 0.077 | 4.450 | |
CalculatingCOEX | 23.399 | 0.398 | 23.812 | |
ClustersList | 0.002 | 0.001 | 0.002 | |
GenesStatistics | 4.374 | 0.078 | 4.442 | |
HandleMetaData | 0.024 | 0.008 | 0.032 | |
HandlingClusterizations | 12.342 | 0.313 | 12.643 | |
HandlingConditions | 0.035 | 0.006 | 0.041 | |
HeatmapPlots | 16.718 | 0.378 | 17.098 | |
LoggingFunctions | 0.001 | 0.000 | 0.001 | |
ParametersEstimations | 10.824 | 0.339 | 11.169 | |
RawDataCleaning | 1.760 | 0.057 | 1.802 | |
RawDataGetters | 0.025 | 0.005 | 0.030 | |
UniformClusters | 84.655 | 0.700 | 85.318 | |
getColorsVector | 0.001 | 0.000 | 0.000 | |