Back to Multiple platform build/check report for BioC 3.20: simplified long |
|
This page was generated on 2024-07-16 11:40 -0400 (Tue, 16 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" | 4677 |
palomino6 | Windows Server 2022 Datacenter | x64 | 4.4.1 (2024-06-14 ucrt) -- "Race for Your Life" | 4416 |
lconway | macOS 12.7.1 Monterey | x86_64 | 4.4.1 (2024-06-14) -- "Race for Your Life" | 4444 |
kjohnson3 | macOS 13.6.5 Ventura | arm64 | 4.4.1 (2024-06-14) -- "Race for Your Life" | 4393 |
palomino8 | Windows Server 2022 Datacenter | x64 | 4.4.1 (2024-06-14 ucrt) -- "Race for Your Life" | 4373 |
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 436/2243 | 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 | ERROR | |||||||||
palomino6 | 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 | |||||||||
palomino8 | Windows Server 2022 Datacenter / x64 | 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: C:\Users\biocbuild\bbs-3.20-bioc\R\bin\R.exe CMD check --no-multiarch --install=check:COTAN.install-out.txt --library=C:\Users\biocbuild\bbs-3.20-bioc\R\library --no-vignettes --timings COTAN_2.5.6.tar.gz |
StartedAt: 2024-07-15 23:00:47 -0400 (Mon, 15 Jul 2024) |
EndedAt: 2024-07-15 23:26:41 -0400 (Mon, 15 Jul 2024) |
EllapsedTime: 1554.3 seconds |
RetCode: 0 |
Status: OK |
CheckDir: COTAN.Rcheck |
Warnings: 0 |
############################################################################## ############################################################################## ### ### Running command: ### ### C:\Users\biocbuild\bbs-3.20-bioc\R\bin\R.exe CMD check --no-multiarch --install=check:COTAN.install-out.txt --library=C:\Users\biocbuild\bbs-3.20-bioc\R\library --no-vignettes --timings COTAN_2.5.6.tar.gz ### ############################################################################## ############################################################################## * using log directory 'C:/Users/biocbuild/bbs-3.20-bioc/meat/COTAN.Rcheck' * using R version 4.4.1 (2024-06-14 ucrt) * using platform: x86_64-w64-mingw32 * R was compiled by gcc.exe (GCC) 13.2.0 GNU Fortran (GCC) 13.2.0 * running under: Windows Server 2022 x64 (build 20348) * 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 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 239.80 1.94 242.30 CalculatingCOEX 97.22 1.59 97.49 HeatmapPlots 66.21 2.04 70.90 ParametersEstimations 64.47 0.84 64.72 HandlingClusterizations 40.59 1.70 41.99 COTANObjectCreation 16.36 2.48 16.38 GenesStatistics 14.84 0.25 14.92 * 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 'C:/Users/biocbuild/bbs-3.20-bioc/meat/COTAN.Rcheck/00check.log' for details.
COTAN.Rcheck/00install.out
############################################################################## ############################################################################## ### ### Running command: ### ### C:\Users\biocbuild\bbs-3.20-bioc\R\bin\R.exe CMD INSTALL COTAN ### ############################################################################## ############################################################################## * installing to library 'C:/Users/biocbuild/bbs-3.20-bioc/R/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 ucrt) -- "Race for Your Life" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 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.18 0.03 0.18
COTAN.Rcheck/tests/spelling.Rout
R version 4.4.1 (2024-06-14 ucrt) -- "Race for Your Life" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 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.14 0.09 0.18
COTAN.Rcheck/tests/testthat.Rout
R version 4.4.1 (2024-06-14 ucrt) -- "Race for Your Life" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 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: 2.24153900146484 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00124788284301758 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00718903541564941 Estimating genes' coex Calculate genes' coex elapsed time: 0.737301111221313 Total calculations elapsed time: 2.98727703094482 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: 3.29573392868042 Calculating genes' coex normalization factor Fraction of genes' with very low expected contingency tables: 0.181818181818182 Calculate genes' normalization factor elapsed time: 0.000890970230102539 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00357508659362793 Estimating genes' coex Calculate genes' coex elapsed time: 0.815598011016846 Total calculations elapsed time: 4.115797996521 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: 2.57664179801941 Calculating cells' coex normalization factor Fraction of cells' with very low expected contingency tables: 0 Calculate cells' normalization factor elapsed time: 0.00110816955566406 Retrieving observed cells' yes/yes contingency table calculating YY.. done Observed cells' contingency table elapsed time: 0.00313282012939453 Estimating cells' coex Calculate cells' coex elapsed time: 0.863232135772705 Total calculations elapsed time: 3.44411492347717 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.619375430892821 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.0262702142278251 | 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: 2.5265040397644 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.000803947448730469 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00197792053222656 Estimating genes' coex Calculate genes' coex elapsed time: 0.84185004234314 Total calculations elapsed time: 3.3711359500885 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: 2.46900296211243 Calculating cells' coex normalization factor Calculate cells' normalization factor elapsed time: 0.000830888748168945 Retrieving observed cells' yes/yes contingency table calculating YY.. done Observed cells' contingency table elapsed time: 0.00194001197814941 Estimating cells' coex Calculate cells' coex elapsed time: 0.80699610710144 Total calculations elapsed time: 3.27876996994019 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: 2.54581904411316 Calculating genes' coex normalization factor Fraction of genes' with very low expected contingency tables: 0.181818181818182 Calculate genes' normalization factor elapsed time: 0.000873088836669922 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0017249584197998 Estimating genes' coex Calculate genes' coex elapsed time: 0.785531044006348 Total calculations elapsed time: 3.33394813537598 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: 2.38498115539551 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.000697851181030273 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00174903869628906 Estimating genes' coex Calculate genes' coex elapsed time: 0.777189970016479 Total calculations elapsed time: 3.16461801528931 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: 2.36004018783569 Calculating cells' coex normalization factor Fraction of cells' with very low expected contingency tables: 0 Calculate cells' normalization factor elapsed time: 0.000901937484741211 Retrieving observed cells' yes/yes contingency table calculating YY.. done Observed cells' contingency table elapsed time: 0.00177502632141113 Estimating cells' coex Calculate cells' coex elapsed time: 0.77904200553894 Total calculations elapsed time: 3.14175915718079 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: 2.90378189086914 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.000733137130737305 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00178790092468262 Estimating genes' coex Calculate genes' coex elapsed time: 0.774948120117188 Total calculations elapsed time: 3.68125104904175 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: 6.47106909751892 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: 2.88798403739929 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00615310668945312 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0456068515777588 Estimating genes' coex Calculate genes' coex elapsed time: 0.748746156692505 Total calculations elapsed time: 3.68849015235901 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.76021003723145 Total elapsed time: 12.0154001712799 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: 6.40040707588196 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: 2.91062808036804 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00607609748840332 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0455708503723145 Estimating genes' coex Calculate genes' coex elapsed time: 0.809651136398315 Total calculations elapsed time: 3.77192616462708 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.52281999588013 Total elapsed time: 11.7322289943695 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: 6.43912601470947 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: 3.00715088844299 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00593209266662598 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.05023193359375 Estimating genes' coex Calculate genes' coex elapsed time: 0.782701015472412 Total calculations elapsed time: 3.84601593017578 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.74821996688843 Total elapsed time: 11.9683330059052 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: 8.91274499893188 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: 3.00428509712219 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00642895698547363 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0452489852905273 Estimating genes' coex Calculate genes' coex elapsed time: 2.21301007270813 Total calculations elapsed time: 5.26897311210632 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 6.08187198638916 Total elapsed time: 15.7800438404083 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: 4.38916897773743 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: 2.58492493629456 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00753593444824219 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00754714012145996 Estimating genes' coex Calculate genes' coex elapsed time: 0.845319032669067 Total calculations elapsed time: 3.44532704353333 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: 2.72854995727539 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00680994987487793 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00918006896972656 Estimating genes' coex Calculate genes' coex elapsed time: 0.806274890899658 Total calculations elapsed time: 3.55081486701965 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: 2.4420599937439 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00818896293640137 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.00728988647460938 Estimating genes' coex Calculate genes' coex elapsed time: 0.778784036636353 Total calculations elapsed time: 3.23632287979126 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: 11.3320739269257 Saving elaborated data locally at: C:\Users\biocbuild\bbs-3.20-bioc\tmpdir\Rtmpa46qY2/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-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 Negative: 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 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 23:20:02 UMAP embedding parameters a = 0.9922 b = 1.112 23:20:02 Read 1200 rows and found 50 numeric columns 23:20:02 Using Annoy for neighbor search, n_neighbors = 30 23:20:02 Building Annoy index with metric = cosine, n_trees = 50 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| 23:20:02 Writing NN index file to temp file C:\Users\biocbuild\bbs-3.20-bioc\tmpdir\Rtmpa46qY2\file1921850b57ad3 23:20:02 Searching Annoy index using 1 thread, search_k = 3000 23:20:03 Annoy recall = 100% 23:20:03 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30 23:20:05 Initializing from normalized Laplacian + noise (using RSpectra) 23:20:05 Commencing optimization for 500 epochs, with 42270 positive edges Using method 'umap' 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| 23:20:10 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: 3.67278480529785 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: 3.0281400680542 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0113208293914795 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0403671264648438 Estimating genes' coex Calculate genes' coex elapsed time: 0.821187973022461 Total calculations elapsed time: 3.90101599693298 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.73664903640747 Total elapsed time: 9.25002098083496 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: 3.3020179271698 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: 3.03884696960449 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00733709335327148 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0327038764953613 Estimating genes' coex Calculate genes' coex elapsed time: 0.836619138717651 Total calculations elapsed time: 3.91550707817078 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.75269913673401 Total elapsed time: 8.90082001686096 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: 3.18818402290344 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: 2.8837468624115 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00845003128051758 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.016510009765625 Estimating genes' coex Calculate genes' coex elapsed time: 0.814328193664551 Total calculations elapsed time: 3.72303509712219 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.56463193893433 Total elapsed time: 8.58449697494507 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: 3.20858693122864 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: 3.00737380981445 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.013375997543335 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0217511653900146 Estimating genes' coex Calculate genes' coex elapsed time: 1.43826699256897 Total calculations elapsed time: 4.48076796531677 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 5.33369493484497 Total elapsed time: 10.2655048370361 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: 3.9691309928894 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: 2.92781591415405 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0127019882202148 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.02081298828125 Estimating genes' coex Calculate genes' coex elapsed time: 1.01730108261108 Total calculations elapsed time: 3.9786319732666 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.85404109954834 Total elapsed time: 9.70939803123474 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-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 Negative: 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 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 23:21:36 UMAP embedding parameters a = 0.9922 b = 1.112 23:21:36 Read 1200 rows and found 50 numeric columns 23:21:36 Using Annoy for neighbor search, n_neighbors = 30 23:21:36 Building Annoy index with metric = cosine, n_trees = 50 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| 23:21:36 Writing NN index file to temp file C:\Users\biocbuild\bbs-3.20-bioc\tmpdir\Rtmpa46qY2\file192185a0d5fb9 23:21:36 Searching Annoy index using 1 thread, search_k = 3000 23:21:36 Annoy recall = 100% 23:21:37 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30 23:21:38 Initializing from normalized Laplacian + noise (using RSpectra) 23:21:38 Commencing optimization for 500 epochs, with 42270 positive edges Using method 'umap' 0% 10 20 30 40 50 60 70 80 90 100% [----|----|----|----|----|----|----|----|----|----| **************************************************| 23:21:43 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: 4.82938098907471 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: 3.05662798881531 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0110421180725098 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0378220081329346 Estimating genes' coex Calculate genes' coex elapsed time: 0.938365936279297 Total calculations elapsed time: 4.04385805130005 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.89435815811157 Total elapsed time: 10.5452401638031 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: 4.14039897918701 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: 2.89729595184326 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0105090141296387 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0431411266326904 Estimating genes' coex Calculate genes' coex elapsed time: 0.868382930755615 Total calculations elapsed time: 3.81932902336121 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.67956113815308 Total elapsed time: 9.71559500694275 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.538244816501366 Handling cluster '3' with mean UDE 1.43530796540674 Handling cluster '4' with mean UDE 0.632684489354434 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.538244816501366 Handling cluster '3' with mean UDE 1.43530796540674 Handling cluster '4' with mean UDE 0.632684489354434 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.538244816501366 Handling cluster '3' with mean UDE 1.43530796540674 Handling cluster '4' with mean UDE 0.632684489354434 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: 3.35174202919006 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: 2.78860807418823 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00755190849304199 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0348360538482666 Estimating genes' coex Calculate genes' coex elapsed time: 0.830036878585815 Total calculations elapsed time: 3.66103291511536 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.54759216308594 Total elapsed time: 8.75048017501831 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: 3.46594309806824 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: 2.80252695083618 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00744318962097168 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0364468097686768 Estimating genes' coex Calculate genes' coex elapsed time: 0.83777117729187 Total calculations elapsed time: 3.6841881275177 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.57472109794617 Total elapsed time: 8.89262008666992 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: 6.12854099273682 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: 3.06846189498901 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00612092018127441 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0449461936950684 Estimating genes' coex Calculate genes' coex elapsed time: 0.822009801864624 Total calculations elapsed time: 3.94153881072998 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.78909516334534 Total elapsed time: 11.765328168869 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: 7.77938485145569 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.43530796540674 Handling cluster '2' with mean UDE 0.640931107489519 Handling cluster '3' with mean UDE 0.546546914955575 Handling cluster '4' with mean UDE 1.22034802329657 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: 4.05067610740662 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: 2.82808685302734 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0133261680603027 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0458419322967529 Estimating genes' coex Calculate genes' coex elapsed time: 0.841284036636353 Total calculations elapsed time: 3.72853899002075 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.58864593505859 Total elapsed time: 9.48018598556519 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: 4.15282797813416 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: 2.85358881950378 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0109570026397705 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0380070209503174 Estimating genes' coex Calculate genes' coex elapsed time: 0.834509134292603 Total calculations elapsed time: 3.73706197738647 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.58983087539673 Total elapsed time: 9.57771182060242 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: 3.82480907440186 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: 2.87442398071289 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0109419822692871 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0306451320648193 Estimating genes' coex Calculate genes' coex elapsed time: 0.837031841278076 Total calculations elapsed time: 3.75304293632507 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.69163799285889 Total elapsed time: 9.36018419265747 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: 3.91703915596008 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: 2.84470415115356 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0107719898223877 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0386388301849365 Estimating genes' coex Calculate genes' coex elapsed time: 0.811388969421387 Total calculations elapsed time: 3.70550394058228 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.51492595672607 Total elapsed time: 9.28242516517639 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: 4.40789103507996 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: 2.86508512496948 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0111038684844971 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0430841445922852 Estimating genes' coex Calculate genes' coex elapsed time: 0.820337057113647 Total calculations elapsed time: 3.73961019515991 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.55813789367676 Total elapsed time: 9.83441686630249 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: 4.49824595451355 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: 2.89403915405273 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0110199451446533 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.045583963394165 Estimating genes' coex Calculate genes' coex elapsed time: 0.845510959625244 Total calculations elapsed time: 3.7961540222168 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.76588916778564 Total elapsed time: 10.2023751735687 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: 5.09750199317932 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: 3.43088316917419 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0108070373535156 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0415289402008057 Estimating genes' coex Calculate genes' coex elapsed time: 0.83339786529541 Total calculations elapsed time: 4.31661701202393 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 5.19287514686584 Total elapsed time: 11.1462140083313 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: 5.17461395263672 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: 2.98239803314209 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00604891777038574 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0429868698120117 Estimating genes' coex Calculate genes' coex elapsed time: 0.80011510848999 Total calculations elapsed time: 3.83154892921448 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.68909907341003 Total elapsed time: 10.6970009803772 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: 6.1833221912384 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: 3.15712404251099 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.00659704208374023 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.045166015625 Estimating genes' coex Calculate genes' coex elapsed time: 0.820941925048828 Total calculations elapsed time: 4.02982902526855 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.93459701538086 Total elapsed time: 11.9361910820007 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: 4.23837804794312 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: 2.86683893203735 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0110781192779541 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0403409004211426 Estimating genes' coex Calculate genes' coex elapsed time: 0.813841104507446 Total calculations elapsed time: 3.7320990562439 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.56960988044739 Total elapsed time: 9.72545790672302 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: 4.3111310005188 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: 2.94214487075806 Calculating genes' coex normalization factor Calculate genes' normalization factor elapsed time: 0.0110390186309814 Retrieving observed genes' yes/yes contingency table calculating YY.. done Observed genes' contingency table elapsed time: 0.0381290912628174 Estimating genes' coex Calculate genes' coex elapsed time: 0.882302045822144 Total calculations elapsed time: 3.873615026474 Calculate genes' coex (legacy): DONE Only genes' coex elapsed time: 4.75784707069397 Total elapsed time: 9.91840410232544 Calculate GDI dataframe: START Using S Calculating S: START Calculating S: DONE S matrix sorted Calculate GDI dataframe: DONE [ FAIL 0 | WARN 41 | SKIP 0 | PASS 452 ] [ FAIL 0 | WARN 41 | SKIP 0 | PASS 452 ] > > proc.time() user system elapsed 759.53 6.31 766.73
COTAN.Rcheck/COTAN-Ex.timings
name | user | system | elapsed | |
COTAN-Legacy | 0 | 0 | 0 | |
COTANObjectCreation | 16.36 | 2.48 | 16.38 | |
CalculatingCOEX | 97.22 | 1.59 | 97.49 | |
ClustersList | 0 | 0 | 0 | |
GenesStatistics | 14.84 | 0.25 | 14.92 | |
HandleMetaData | 0.08 | 0.00 | 0.08 | |
HandlingClusterizations | 40.59 | 1.70 | 41.99 | |
HandlingConditions | 0.11 | 0.04 | 0.14 | |
HeatmapPlots | 66.21 | 2.04 | 70.90 | |
LoggingFunctions | 0.00 | 0.02 | 0.02 | |
ParametersEstimations | 64.47 | 0.84 | 64.72 | |
RawDataCleaning | 4.12 | 0.08 | 4.19 | |
RawDataGetters | 0.08 | 0.00 | 0.07 | |
UniformClusters | 239.80 | 1.94 | 242.30 | |
getColorsVector | 0 | 0 | 0 | |