Multiple platform build/check report for BioC 3.22 - CHECK results for COTAN on lconway

Hostname	OS	Arch (*)	R version	Installed pkgs
nebbiolo2	Linux (Ubuntu 24.04.3 LTS)	x86_64	4.5.1 Patched (2025-08-23 r88802) -- "Great Square Root"	4853
lconway	macOS 12.7.1 Monterey	x86_64	4.5.1 Patched (2025-09-10 r88807) -- "Great Square Root"	4640
kjohnson3	macOS 13.7.7 Ventura	arm64	4.5.1 Patched (2025-09-10 r88807) -- "Great Square Root"	4585
taishan	Linux (openEuler 24.03 LTS)	aarch64	4.5.0 (2025-04-11) -- "How About a Twenty-Six"	4576
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*

CHECK results for COTAN on lconway

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.

raw results

Summary

Package: COTAN

Version: 2.9.6

Command: /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:COTAN.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings COTAN_2.9.6.tar.gz

StartedAt: 2025-10-03 20:44:12 -0400 (Fri, 03 Oct 2025)

EndedAt: 2025-10-03 21:06:26 -0400 (Fri, 03 Oct 2025)

EllapsedTime: 1334.1 seconds

RetCode: 0

Status: OK

CheckDir: COTAN.Rcheck

Warnings: 0

Command output

##############################################################################
##############################################################################
###
### Running command:
###
###   /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:COTAN.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings COTAN_2.9.6.tar.gz
###
##############################################################################
##############################################################################


* using log directory ‘/Users/biocbuild/bbs-3.22-bioc/meat/COTAN.Rcheck’
* using R version 4.5.1 Patched (2025-09-10 r88807)
* using platform: x86_64-apple-darwin20
* R was compiled by
    Apple clang version 14.0.0 (clang-1400.0.29.202)
    GNU Fortran (GCC) 14.2.0
* running under: macOS Monterey 12.7.6
* 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.9.6’
* package encoding: UTF-8
* checking package namespace information ... OK
* checking package dependencies ... INFO
Imports includes 26 non-default packages.
Importing from so many packages makes the package vulnerable to any of
them becoming unavailable.  Move as many as possible to Suggests and
use conditionally.
* checking if this is a source package ... OK
* checking if there is a namespace ... OK
* checking for hidden files and directories ... OK
* checking for portable file names ... OK
* checking for sufficient/correct file permissions ... OK
* checking whether package ‘COTAN’ can be installed ... OK
* checking installed package size ... OK
* checking package directory ... OK
* checking ‘build’ directory ... OK
* checking DESCRIPTION meta-information ... OK
* checking top-level files ... OK
* checking for left-over files ... OK
* checking index information ... OK
* checking package subdirectories ... OK
* checking code files for non-ASCII characters ... OK
* checking R files for syntax errors ... OK
* checking whether the package can be loaded ... OK
* checking whether the package can be loaded with stated dependencies ... OK
* checking whether the package can be unloaded cleanly ... OK
* checking whether the namespace can be loaded with stated dependencies ... OK
* checking whether the namespace can be unloaded cleanly ... OK
* checking whether startup messages can be suppressed ... OK
* checking dependencies in R code ... OK
* checking S3 generic/method consistency ... OK
* checking replacement functions ... OK
* checking foreign function calls ... OK
* checking R code for possible problems ... OK
* 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         117.537  2.626 119.597
CalculatingCOEX          55.689  1.145  57.105
HandlingClusterizations  43.884  1.174  44.691
HeatmapPlots             34.730  0.941  35.883
ParametersEstimations    23.568  0.737  24.414
GenesStatistics          11.139  0.309  11.360
COTAN_ObjectCreation     10.690  0.304  11.070
RawDataCleaning           5.099  0.332   5.268
Conversions               5.021  0.128   5.175
* 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: OK

Installation output

COTAN.Rcheck/00install.out

##############################################################################
##############################################################################
###
### Running command:
###
###   /Library/Frameworks/R.framework/Resources/bin/R CMD INSTALL COTAN
###
##############################################################################
##############################################################################


* installing to library ‘/Library/Frameworks/R.framework/Versions/4.5-x86_64/Resources/library’
* installing *source* package ‘COTAN’ ...
** this is package ‘COTAN’ version ‘2.9.6’
** using staged installation
** R
** data
** inst
** byte-compile and prepare package for lazy loading
** 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)

Tests output

COTAN.Rcheck/tests/outputTestDatasetCreation.Rout


R version 4.5.1 Patched (2025-09-10 r88807) -- "Great Square Root"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin20

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> 
> # Creates the files to be reloaded by the tests for comparisons
> 
> outputTestDatasetCreation <-
+   function(testsDir = file.path("tests", "testthat")) {
+   utils::data("test.dataset", package = "COTAN")
+   options(parallelly.fork.enable = TRUE)
+   setLoggingLevel(3L)
+ 
+   obj <- COTAN(raw = test.dataset)
+   obj <- initializeMetaDataset(objCOTAN = obj, GEO = " ",
+                                sequencingMethod = "artificial",
+                                sampleCondition = "test")
+ 
+   obj <- proceedToCoex(objCOTAN = obj, cores = 6L, saveObj = FALSE)
+ 
+   if (FALSE) {
+     saveRDS(obj, file = file.path(testsDir,"test.COTAN.RDS"))
+   }
+ 
+   cells.names.test <- getCells(objCOTAN = obj)[c(1L:10L, 591L:610L, 991L:1000L)]
+   genes.names.test <- getGenes(objCOTAN = 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"))
+ 
+   pcaRaw <- runPCA(x = getRawData(objCOTAN = obj), rank = 10L,
+                    BSPARAM = IrlbaParam(), get.rotation = FALSE)[["x"]]
+ 
+   pca.raw.test <- pcaRaw[genes.names.test, ]
+   saveRDS(pca.raw.test, file.path(testsDir, "pca.raw.test.RDS"))
+ 
+   dispersion.test <- getDispersion(objCOTAN = obj)[genes.names.test]
+   saveRDS(dispersion.test, file.path(testsDir, "dispersion.test.RDS"))
+ 
+   raw.norm.test <-
+     getNuNormData(objCOTAN = obj)[genes.names.test, cells.names.test]
+   saveRDS(raw.norm.test, file.path(testsDir, "raw.norm.test.RDS"))
+ 
+   coex.test <-
+     getGenesCoex(objCOTAN = obj, genes = genes.names.test, zeroDiagonal = FALSE)
+   saveRDS(coex.test, file.path(testsDir, "coex.test.RDS"))
+ 
+   lambda.test <- getLambda(objCOTAN = obj)[genes.names.test]
+   saveRDS(lambda.test, file.path(testsDir, "lambda.test.RDS"))
+ 
+   GDI.test <- calculateGDI(objCOTAN = obj)
+   GDI.test <- GDI.test[genes.names.test, ]
+   saveRDS(GDI.test, file.path(testsDir, "GDI.test.RDS"))
+ 
+   nu.test <- getNu(objCOTAN = obj)[cells.names.test]
+   saveRDS(nu.test, file.path(testsDir, "nu.test.RDS"))
+ 
+   pvalues.test <- calculatePValue(objCOTAN = obj,
+                                   geneSubsetCol = genes.names.test,
+                                   geneSubsetRow = genes.names.test)
+   saveRDS(pvalues.test, file.path(testsDir, "pvalues.test.RDS"))
+ 
+   groupMarkers <- list(G1 = c("g-000010", "g-000020", "g-000030"),
+                        G2 = c("g-000300", "g-000330", "g-000660"),
+                        G3 = c("g-000510", "g-000530", "g-000550",
+                               "g-000570", "g-000590"))
+ 
+   gcsData <- genesCoexSpace(objCOTAN = obj,
+                             primaryMarkers = unlist(groupMarkers),
+                             numGenesPerMarker = 11L)
+ 
+   genes.coex.space.test <- gcsData[["GCS"]]
+   saveRDS(genes.coex.space.test,
+           file.path(testsDir, "genes.coex.space.test.RDS"))
+ 
+   genesClustersData <-
+     establishGenesClusters(objCOTAN = obj,
+                            groupMarkers = groupMarkers,
+                            numGenesPerMarker = 11L,
+                            kCuts = 6L, distance = "cosine",
+                            hclustMethod = "ward.D2")
+ 
+   pca.genes.clusters.test <- genesClustersData[["pca_clusters"]]
+   saveRDS(pca.genes.clusters.test,
+           file.path(testsDir, "pca.genes.clusters.test.RDS"))
+ 
+   # Make it a less strict check as it is only for testing
+   checker <- new("AdvancedGDIUniformityCheck")
+   checker <- shiftCheckerThresholds(checker, 0.1)
+ 
+   initialResolution <- 0.8
+   splitData <- cellsUniformClustering(objCOTAN = obj,
+                                       checker = checker,
+                                       initialResolution = initialResolution,
+                                       useCoexEigen = TRUE,
+                                       dataMethod = "LL",
+                                       numReducedComp = 50L,
+                                       cores = 6L, optimizeForSpeed = TRUE,
+                                       deviceStr = "cuda", saveObj = FALSE)
+ 
+   split.clusters.test <- splitData[["clusters"]]
+   saveRDS(split.clusters.test, file.path(testsDir, "split.clusters.test.RDS"))
+ 
+   obj <- addClusterization(objCOTAN = obj,
+                            clName = "split",
+                            clusters = splitData[["clusters"]],
+                            coexDF = splitData[["coex"]])
+ 
+   coex.clusters.test <- splitData[["coex"]][genes.names.test, ]
+   saveRDS(coex.clusters.test, file.path(testsDir, "coex.clusters.test.RDS"))
+ 
+   pvalDF <- pValueFromDEA(splitData[["coex"]],
+                           getNumCells(objCOTAN = obj),
+                           adjustmentMethod = "none")
+ 
+   pvalues.clusters.test <- pvalDF[genes.names.test, ]
+   saveRDS(pvalues.clusters.test,
+           file.path(testsDir, "pvalues.clusters.test.RDS"))
+ 
+   mergedData <- mergeUniformCellsClusters(objCOTAN = obj,
+                                           clusters = splitData[["clusters"]],
+                                           checkers = checker,
+                                           cores = 6L,
+                                           distance = "cosine",
+                                           hclustMethod = "ward.D2",
+                                           saveObj = FALSE)
+ 
+   merge.clusters.test <- mergedData[["clusters"]]
+   saveRDS(merge.clusters.test, file.path(testsDir, "merge.clusters.test.RDS"))
+ 
+   obj <- addClusterization(objCOTAN = obj,
+                            clName = "merge",
+                            clusters = mergedData[["clusters"]],
+                            coexDF = mergedData[["coex"]])
+ }
> 
> proc.time()
   user  system elapsed 
  0.257   0.094   0.404

COTAN.Rcheck/tests/spelling.Rout


R version 4.5.1 Patched (2025-09-10 r88807) -- "Great Square Root"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin20

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> if (requireNamespace("spelling", quietly = TRUE))
+   spelling::spell_check_test(vignettes = TRUE, error = FALSE,
+                              skip_on_cran = TRUE)
All Done!
> 
> proc.time()
   user  system elapsed 
  0.223   0.098   0.351

COTAN.Rcheck/tests/testthat.Rout


R version 4.5.1 Patched (2025-09-10 r88807) -- "Great Square Root"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin20

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> Sys.setenv(R_TESTS = "")
> library(testthat)
> library(COTAN)
> library(zeallot)
> 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
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.85019588470459
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.903596292963576 | max: 4.67949426963453 | % 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.62690901756287
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.000844955444335938
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00735998153686523
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.877569913864136
Total calculations elapsed time: 3.5126838684082
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 5 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 4 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 5.2992479801178
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.903596292963576 | max: 4.67949426963453 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 7 cells batches
Total calculations elapsed time: 0.858681917190552
Estimate `nu`: DONE
`nu` change (abs) | max: 1.75595238095238 | median:  1.07156808035714 | mean:  1.07156808035714
Estimate `dispersion`/`nu`: START
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.71418499946594
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 1.03734183779544 | max: 4.61112297272667 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.914663076400757
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0295487661210321 | median:  0.016072261901129 | mean:  0.016072261901129
`nu` mean: 1.69786080481514
Marginal errors | max: 1.95799528047999 | median 1.32268477389029 | mean: 1.33563364783916
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.8302218914032
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.0585450475968788 | max: 3.52739274166646 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.904983043670654
Estimate `nu`: DONE
`nu` change (abs) | max: 0.417150556912302 | median:  0.240156504520849 | mean:  0.240156504520849
`nu` mean: 0.823005947608547
Marginal errors | max: 0.837672524335753 | median 0.704691669139061 | mean: 0.646575831699996
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.83638978004456
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.328845855849147 | max: 4.03031691940854 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.926531076431274
Estimate `nu`: DONE
`nu` change (abs) | max: 0.164439840041 | median:  0.0957105992188852 | mean:  0.0957105992188852
`nu` mean: 1.06872924082212
Marginal errors | max: 0.260109847946692 | median 0.21392791155499 | mean: 0.197565106613055
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.76353597640991
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.229469002841045 | max: 3.87206403670079 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.928059101104736
Estimate `nu`: DONE
`nu` change (abs) | max: 0.05539667181205 | median:  0.0321044921874999 | mean:  0.0321044921874999
`nu` mean: 0.97670782037545
Marginal errors | max: 0.095347789718808 | median 0.0796003827042702 | mean: 0.0725509735761708
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.82825613021851
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.263830434255107 | max: 3.92953043775907 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.935965061187744
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0196209299693411 | median:  0.0114037947200925 | mean:  0.0114037947200925
`nu` mean: 1.00821713524925
Marginal errors | max: 0.0329495161656741 | median 0.0273479725376307 | mean: 0.0250296751521313
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.80164790153503
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.251734659898778 | max: 3.9096133246649 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.877485990524292
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0069104540740097 | median:  0.0040109538496062 | mean:  0.0040109538496062
`nu` mean: 0.997100499775596
Marginal errors | max: 0.0116950057438405 | median 0.00972280144874116 | mean: 0.00889105008777609
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.72330403327942
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.256009779031907 | max: 3.91667271750655 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.886433839797974
Estimate `nu`: DONE
`nu` change (abs) | max: 0.00251004753849005 | median:  0.00146484374999994 | mean:  0.00146484374999994
`nu` mean: 1.00104520378849
Marginal errors | max: 0.00430755646112857 | median 0.00358203399504298 | mean: 0.00327684342604062
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.80814099311829
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.254429717857093 | max: 3.91406908073254 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.85191798210144
Estimate `nu`: DONE
`nu` change (abs) | max: 0.00104629156456748 | median:  0.000592910406456731 | mean:  0.000592910406456731
`nu` mean: 0.999546618841889
Marginal errors | max: 0.00162038706091039 | median 0.00134441676098973 | mean: 0.00123030512553246
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.70035696029663
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.255022049769696 | max: 3.91504587939868 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.900820016860962
Estimate `nu`: DONE
`nu` change (abs) | max: 0.000418450901304768 | median:  0.000244140625000028 | mean:  0.000244140625000028
`nu` mean: 1.0001743102763
Marginal errors | max: 0.00071569233034019 | median 0.000597297218559056 | mean: 0.000545906995079548
Total calculations elapsed time: 47.8489689826965
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
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.70960187911987
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.903596292963576 | max: 4.67949426963453 | % 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 available 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.69518709182739
Calculating genes' COEX normalization factor
Fraction of genes' with very low expected contingency tables: 0.181818181818182
Calculate genes' normalization factor elapsed time: 0.00107693672180176
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00317788124084473
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.687229156494141
Total calculations elapsed time: 3.38667106628418
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.34444618225098
Calculating cells' COEX normalization factor
Fraction of cells' with very low expected contingency tables: 0
Calculate cells' normalization factor elapsed time: 0.000842809677124023
Retrieving observed cells' yes/yes contingency table
calculating YY.. done
Observed cells' contingency table elapsed time: 0.00272202491760254
Estimating cells' COEX
Calculate cells' COEX elapsed time: 0.929540157318115
Total calculations elapsed time: 3.27755117416382
Calculate cells' COEX (legacy): DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0120480060577393
Calculate `GDI`: DONE
Total calculations elapsed time: 2.65867304801941
Calculate GDI dataframe: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.00391697883605957
Calculate `GDI`: DONE
Total calculations elapsed time: 2.73789691925049
Calculate GDI dataframe: DONE
Included 1 genes with GDI below 1.3
Selected 10 genes using HGDI selector
Running genes' selection: DONE
Running genes' selection: START
Included 1 genes with GDI below 1.3
Selected 10 genes using HGDI selector
Running genes' selection: DONE
Running genes' selection: START
Included 1 genes with GDI below 1.3
Selected 10 genes using HGDI selector
Running genes' selection: DONE
Running genes' selection: START
Selected 5 genes using HGDI selector
Running genes' selection: DONE
Running genes' selection: 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%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Selected 5 genes using HVG_Seurat selector
Running genes' selection: DONE
Running genes' selection: START
Normalizing layer: counts
Performing log-normalization
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Finding variable features for layer data
Calculating gene means
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Calculating gene variance to mean ratios
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Running genes' selection: START
Given 5 genes as input
Running genes' selection: DONE
Elaborating Reduced dimensionality Data Matrix - START
Running genes' selection: START
Included 1 genes with GDI below 1.3
Selected 10 genes using HGDI selector
Running genes' selection: DONE
Elaborating PCA - START
Elaborating PCA - DONE
Total calculations elapsed time: 1.9897289276123
Elaborating Reduced dimensionality Data Matrix - DONE
Elaborating Reduced dimensionality Data Matrix - START
Running genes' selection: START
Included 1 genes with GDI below 1.3
Selected 10 genes using HGDI selector
Running genes' selection: DONE
Elaborating PCA - START
Elaborating PCA - DONE
Total calculations elapsed time: 1.7896900177002
Elaborating Reduced dimensionality Data Matrix - DONE
Elaborating Reduced dimensionality Data Matrix - START
Running genes' selection: START
Selected 8 genes using HGDI selector
Running genes' selection: DONE
Elaborating PCA - START
Elaborating PCA - DONE
Total calculations elapsed time: 1.80332398414612
Elaborating Reduced dimensionality Data Matrix - DONE
Elaborating Reduced dimensionality Data Matrix - START
Running genes' selection: 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%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Selected 8 genes using HVG_Seurat selector
Running genes' selection: DONE
Elaborating PCA - START
Elaborating PCA - DONE
Total calculations elapsed time: 2.85831809043884
Elaborating Reduced dimensionality Data Matrix - DONE
Elaborating Reduced dimensionality Data Matrix - START
Elaborating COEX Eigen Vectors - START
Elaborating COEX Eigen Vectors - DONE
Total calculations elapsed time: 2.75763392448425
Elaborating Reduced dimensionality Data Matrix - DONE
Elaborating Reduced dimensionality Data Matrix - START
Elaborating COEX Eigen Vectors - START
Elaborating COEX Eigen Vectors - DONE
Total calculations elapsed time: 2.81905817985535
Elaborating Reduced dimensionality Data Matrix - 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
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.92954993247986
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.903596292963576 | max: 4.67949426963453 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.888302087783813
Estimate `nu`: DONE
`nu` change (abs) | max: 1.75719246031746 | median:  1.07229953342014 | mean:  1.07229953342014
`nu` mean: 1.68489292689732
Marginal errors | max: 1.735699617106 | median 1.3799744133366 | mean: 1.32181397451301
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.78014183044434
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.0655494303592343 | max: 3.54400580804152 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.900699138641357
Estimate `nu`: DONE
`nu` change (abs) | max: 0.402649984216273 | median:  0.231868788425666 | mean:  0.231868788425666
`nu` mean: 0.829218804209393
Marginal errors | max: 0.803248939260364 | median 0.677504196374578 | mean: 0.619386446933353
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.82906103134155
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.326019058449519 | max: 4.02612397202046 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.935062885284424
Estimate `nu`: DONE
`nu` change (abs) | max: 0.158004893526231 | median:  0.0919692884670312 | mean:  0.0919692884670312
`nu` mean: 1.0660356050592
Marginal errors | max: 0.250677002875358 | median 0.206241644985115 | mean: 0.190402167377958
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.80829691886902
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.230397316576369 | max: 3.87365659194104 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.912688970565796
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0532774732102337 | median:  0.0308837890624999 | mean:  0.0308837890624999
`nu` mean: 0.977606315852266
Marginal errors | max: 0.0916920322976189 | median 0.0765360853958157 | mean: 0.069767303009408
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.84523487091064
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.263457616228195 | max: 3.928924407347 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.915508985519409
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0187878665732679 | median:  0.0109155550210726 | mean:  0.0109155550210726
`nu` mean: 1.0078723115522
Marginal errors | max: 0.0315207323253404 | median 0.0261635598193228 | mean: 0.0239435763640454
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.79739499092102
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.251890641847532 | max: 3.90987172241566 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.902047157287598
Estimate `nu`: DONE
`nu` change (abs) | max: 0.00649137083358098 | median:  0.00376680890457393 | mean:  0.00376680890457393
`nu` mean: 0.997275438070993
Marginal errors | max: 0.0109748802265024 | median 0.00912386940474885 | mean: 0.00834350038184297
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.77841210365295
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.255902752363649 | max: 3.91649655995187 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.927188158035278
Estimate `nu`: DONE
`nu` change (abs) | max: 0.00230094227000666 | median:  0.00132529141000018 | mean:  0.00132529141000018
`nu` mean: 1.00097565086001
Marginal errors | max: 0.00376864912836616 | median 0.00313421412286008 | mean: 0.00286709489297508
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.74691319465637
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.254520292119112 | max: 3.91421850264948 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.86217999458313
Estimate `nu`: DONE
`nu` change (abs) | max: 0.000837013154366462 | median:  0.000488281250000028 | mean:  0.000488281250000028
`nu` mean: 0.999651268095634
Marginal errors | max: 0.0014405532051498 | median 0.00119489364392944 | mean: 0.00109355439757195
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 1.80145406723022
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.255046800817985 | max: 3.91508667609651 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Total calculations elapsed time: 0.920964002609253
Estimate `nu`: DONE
`nu` change (abs) | max: 0.000418448154490214 | median:  0.000244140624999917 | mean:  0.000244140624999917
`nu` mean: 1.00017430752949
Marginal errors | max: 0.000715683388516908 | median 0.000597284734240588 | mean: 0.000545899916522963
Total calculations elapsed time: 48.7757339477539
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 available 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.79717612266541
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.000794887542724609
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00208210945129395
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.938452959060669
Total calculations elapsed time: 3.73850607872009
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.75935077667236
Calculating cells' COEX normalization factor
Calculate cells' normalization factor elapsed time: 0.000793218612670898
Retrieving observed cells' yes/yes contingency table
calculating YY.. done
Observed cells' contingency table elapsed time: 0.00167679786682129
Estimating cells' COEX
Calculate cells' COEX elapsed time: 0.891630172729492
Total calculations elapsed time: 3.65345096588135
Calculate cells' COEX (legacy): DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.00357794761657715
Calculate `GDI`: DONE
Total calculations elapsed time: 2.8411979675293
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
Total calculations elapsed time: 7.2959189414978
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
Total calculations elapsed time: 7.21833801269531
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
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 3.71356916427612
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.903596292963576 | max: 4.67949426963453 | % 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.70288300514221
Calculating genes' COEX normalization factor
Fraction of genes' with very low expected contingency tables: 0.181818181818182
Calculate genes' normalization factor elapsed time: 0.000760078430175781
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00190377235412598
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.931997060775757
Total calculations elapsed time: 3.63754391670227
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 available 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.70214986801147
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.000482082366943359
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00152492523193359
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.850414037704468
Total calculations elapsed time: 3.55457091331482
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
Total calculations elapsed time: 7.35657691955566
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
Total calculations elapsed time: 7.42557501792908
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
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 3.7520911693573
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.903596292963576 | max: 4.67949426963453 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Total calculations elapsed time: 0.930413961410522
Estimate `nu`: DONE
`nu` change (abs) | max: 1.75719246031746 | median:  1.07229953342014 | mean:  1.07229953342014
`nu` mean: 1.68489292689732
Marginal errors | max: 0.255290533439036 | median 0.0807313765005375 | mean: 0.101968565509113
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 3.59460806846619
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 1.03709166563768 | max: 4.61069655774133 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Total calculations elapsed time: 0.909102916717529
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0273438105507502 | median:  0.0148852611818011 | mean:  0.0148852611818011
`nu` mean: 1.69735147626627
Marginal errors | max: 0.00330148418741771 | median 0.00111713187082163 | mean: 0.00132007186543017
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 3.69048714637756
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 1.03887279473672 | max: 4.60974134092906 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Total calculations elapsed time: 0.944160938262939
Estimate `nu`: DONE
`nu` change (abs) | max: 0 | median:  0 | mean:  0
`nu` mean: 1.69735147626627
Marginal errors | max: 5.14174422185931e-05 | median 1.68182247684001e-07 | mean: 5.57649710817287e-06
Total calculations elapsed time: 22.0896100997925
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.8149847984314
Calculating cells' COEX normalization factor
Fraction of cells' with very low expected contingency tables: 0
Calculate cells' normalization factor elapsed time: 0.00115108489990234
Retrieving observed cells' yes/yes contingency table
calculating YY.. done
Observed cells' contingency table elapsed time: 0.00211310386657715
Estimating cells' COEX
Calculate cells' COEX elapsed time: 0.936187028884888
Total calculations elapsed time: 3.75443601608276
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
Total calculations elapsed time: 7.41419196128845
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
Total calculations elapsed time: 7.30747413635254
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
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 3.59776401519775
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.903596292963576 | max: 4.67949426963453 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Total calculations elapsed time: 0.956881046295166
Estimate `nu`: DONE
`nu` change (abs) | max: 1.75719246031746 | median:  1.07229953342014 | mean:  1.07229953342014
`nu` mean: 1.68489292689732
Marginal errors | max: 0.255290533439036 | median 0.0807313765005375 | mean: 0.101968565509113
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 3.75420308113098
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 1.03709166563768 | max: 4.61069655774133 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Total calculations elapsed time: 1.0006730556488
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0273438105507502 | median:  0.0148852611818011 | mean:  0.0148852611818011
`nu` mean: 1.69735147626627
Marginal errors | max: 0.00330148418741771 | median 0.00111713187082163 | mean: 0.00132007186543017
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Parallel dispersion Newton-Raphson: used up to 5 iterations
Total calculations elapsed time: 3.80162596702576
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 1.03887279473672 | max: 4.60974134092906 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Total calculations elapsed time: 0.882180213928223
Estimate `nu`: DONE
`nu` change (abs) | max: 0 | median:  0 | mean:  0
`nu` mean: 1.69735147626627
Marginal errors | max: 5.14174422185931e-05 | median 1.68182247684001e-07 | mean: 5.57649710817287e-06
Total calculations elapsed time: 22.2828030586243
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.8286030292511
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00083613395690918
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00214195251464844
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.916130065917969
Total calculations elapsed time: 3.74771118164062
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
Total calculations elapsed time: 6.33845210075378
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
Total calculations elapsed time: 0.00247001647949219
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
Total calculations elapsed time: 6.41618394851685
Calculating G: DONE
calculating PValues: START
Get p-values on a set of genes on columns and on a set of genes on rows
Total calculations elapsed time: 6.4174599647522
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0021970272064209
Calculate `GDI`: DONE
Total calculations elapsed time: 2.82330083847046
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
Total calculations elapsed time: 6.51238703727722
Calculating G: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.00329899787902832
Calculate `GDI`: DONE
Total calculations elapsed time: 9.31096291542053
Calculate GDI dataframe: DONE
Calculate `GDI`: START
Calculating S: START
Calculating S: DONE
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.00461101531982422
Calculate `GDI`: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.00415897369384766
Calculate `GDI`: DONE
Initializing `COTAN` meta-data
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Dataset cleaning elapsed time: 0.909078121185303
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.10199785232544
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Model parameter estimation elapsed time: 3.03962588310242
COTAN genes' COEX estimation: START
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 available 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.42653393745422
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.0116019248962402
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0631370544433594
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.941540002822876
Total calculations elapsed time: 4.4428129196167
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.45476889610291
Dataset analysis elapsed time: 9.40347290039062
COTAN dataset analysis: DONE
Initializing `COTAN` meta-data
Condition test
n cells 1200
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Dataset cleaning elapsed time: 0.961167097091675
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.21266412734985
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Model parameter estimation elapsed time: 3.14563703536987
COTAN genes' COEX estimation: START
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 available 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.35911393165588
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.0100259780883789
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0543630123138428
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.898377180099487
Total calculations elapsed time: 4.32188010215759
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.20294690132141
Dataset analysis elapsed time: 9.30975103378296
COTAN dataset analysis: DONE
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
Total calculations elapsed time: 0.0253190994262695
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.03426194190979
Calculate `GDI`: DONE
Total calculations elapsed time: 2.88720798492432
Calculate GDI dataframe: DONE
Initializing `COTAN` meta-data
Condition test
n cells 1200
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Dataset cleaning elapsed time: 0.907362937927246
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.23406386375427
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Model parameter estimation elapsed time: 3.16041612625122
COTAN genes' COEX estimation: START
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 available 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.31560111045837
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00754094123840332
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0513780117034912
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.871222972869873
Total calculations elapsed time: 4.24574303627014
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.13575196266174
Dataset analysis elapsed time: 9.20353102684021
COTAN dataset analysis: DONE
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
Total calculations elapsed time: 0.0151290893554688
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.031757116317749
Calculate `GDI`: DONE
Total calculations elapsed time: 2.72204780578613
Calculate GDI dataframe: DONE
Initializing `COTAN` meta-data
Condition test
n cells 1200
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Dataset cleaning elapsed time: 0.909581899642944
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.24370312690735
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Model parameter estimation elapsed time: 3.21257090568542
COTAN genes' COEX estimation: START
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 available 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.44063019752502
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.0059969425201416
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0478429794311523
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.951380014419556
Total calculations elapsed time: 4.44585013389587
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.355229139328
Dataset analysis elapsed time: 9.47738194465637
COTAN dataset analysis: DONE
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
Total calculations elapsed time: 0.00634002685546875
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0333249568939209
Calculate `GDI`: DONE
Total calculations elapsed time: 2.9143979549408
Calculate GDI dataframe: DONE
Initializing `COTAN` meta-data
Asked to drop 0 genes and 1160 cells
COTAN dataset analysis: START
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
Dataset cleaning elapsed time: 0.953305006027222
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 6 iterations
Total calculations elapsed time: 1.79482388496399
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: -0.212044752677142 | max: 85.4528331044021 | % negative: 62.2144112478032
Model parameter estimation elapsed time: 2.69865894317627
COTAN genes' COEX estimation not requested
Only genes' COEX elapsed time: 0.885313987731934
Dataset analysis elapsed time: 4.53727793693542
COTAN dataset analysis: 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.80587005615234
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00625491142272949
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0091850757598877
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.936659812927246
Total calculations elapsed time: 3.75796985626221
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 available 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.95598697662354
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00706601142883301
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00546503067016602
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.941934108734131
Total calculations elapsed time: 3.91045212745667
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 available 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.79556584358215
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00863409042358398
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00794506072998047
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.901695013046265
Total calculations elapsed time: 3.71384000778198
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
Total calculations elapsed time: 0.00578808784484863
calculating PValues: DONE
Handling genes type: G2
Handling genes type: G3
Initializing `COTAN` meta-data
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Clean plots: START
PCA: START
PCA: DONE
Hierarchical clustering: START
Hierarchical clustering: DONE
Total calculations elapsed time: 2.11943697929382
Clean plots: DONE
Dataset cleaning elapsed time: 4.50619506835938
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.10467791557312
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Model parameter estimation elapsed time: 3.03894996643066
COTAN genes' COEX estimation: START
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 available 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.45270204544067
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00585794448852539
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0427711009979248
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.94084095954895
Total calculations elapsed time: 4.44217205047607
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.38757610321045
Dataset analysis elapsed time: 12.9327211380005
COTAN dataset analysis: DONE
Saving elaborated data locally at: /tmp/RtmpuQrPnJ/test.cotan.RDS
Creating new clusterization: START
Elaborating Reduced dimensionality Data Matrix - START
Elaborating COEX Eigen Vectors - START
Elaborating COEX Eigen Vectors - DONE
Total calculations elapsed time: 3.27352499961853
Elaborating Reduced dimensionality Data Matrix - DONE
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: 41850

Running Louvain algorithm with multilevel refinement...
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Maximum modularity in 10 random starts: 0.5751
Number of communities: 2
Elapsed time: 0 seconds
Number of clusters is too small. Reclustering at resolution higher than: 0.8
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck

Number of nodes: 1200
Number of edges: 41850

Running Louvain algorithm with multilevel refinement...
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Maximum modularity in 10 random starts: 0.3786
Number of communities: 7
Elapsed time: 0 seconds
2 singletons identified. 5 final clusters.
Used resolution for Seurat clusterization is: 1.3
Total calculations elapsed time: 4.61880111694336
Creating new clusterization: DONE
Creating cells' uniform clustering: START
In iteration 1 the number of cells to re-cluster is 1200 cells belonging to 0 clusters
Asked to drop no genes or cells
Creating new clusterization: START
Elaborating Reduced dimensionality Data Matrix - START
Elaborating COEX Eigen Vectors - START
Elaborating COEX Eigen Vectors - DONE
Total calculations elapsed time: 3.1293089389801
Elaborating Reduced dimensionality Data Matrix - DONE
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: 29356

Running Louvain algorithm with multilevel refinement...
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Maximum modularity in 10 random starts: 0.5763
Number of communities: 2
Elapsed time: 0 seconds
Used resolution for Seurat clusterization is: 0.8
Total calculations elapsed time: 4.13958597183228
Creating new clusterization: DONE
Creating PDF UMAP in file:  /tmp/RtmpuQrPnJ/test/reclustering/pdf_umap_1.pdf
UMAP plot: START
Calculating UMAP: START
21:02:11 UMAP embedding parameters a = 0.9922 b = 1.112
21:02:11 Read 1200 rows and found 65 numeric columns
21:02:11 Using Annoy for neighbor search, n_neighbors = 30
21:02:11 Building Annoy index with metric = cosine, n_trees = 50
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
21:02:12 Writing NN index file to temp file /tmp/RtmpuQrPnJ/file38415bb4dab2
21:02:12 Searching Annoy index using 1 thread, search_k = 3000
21:02:12 Annoy recall = 100%
21:02:13 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
21:02:15 Initializing from normalized Laplacian + noise (using RSpectra)
21:02:15 Commencing optimization for 500 epochs, with 48472 positive edges
21:02:15 Using rng type: pcg
Using method 'umap'
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
21:02:18 Optimization finished
Calculating UMAP: DONE
Total calculations elapsed time: 6.44192290306091
UMAP plot: DONE
UMAP plot: START
Calculating UMAP: START
21:02:18 UMAP embedding parameters a = 0.9922 b = 1.112
21:02:18 Read 1200 rows and found 65 numeric columns
21:02:18 Using Annoy for neighbor search, n_neighbors = 30
21:02:18 Building Annoy index with metric = cosine, n_trees = 50
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
21:02:18 Writing NN index file to temp file /tmp/RtmpuQrPnJ/file38414baa4d80
21:02:18 Searching Annoy index using 1 thread, search_k = 3000
21:02:19 Annoy recall = 100%
21:02:20 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
21:02:22 Initializing from normalized Laplacian + noise (using RSpectra)
21:02:22 Commencing optimization for 500 epochs, with 48472 positive edges
21:02:22 Using rng type: pcg
Using method 'umap'
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
21:02:25 Optimization finished
Calculating UMAP: DONE
Total calculations elapsed time: 6.58395504951477
UMAP plot: DONE
* checking uniformity of cluster '0' of 2 clusters
Asked to drop 0 genes and 587 cells
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [613] cells
Dataset cleaning elapsed time: 0.981482028961182
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.10370111465454
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: -0.0334270578203965 | max: 8.68674670630126 | % negative: 4.16666666666667
Model parameter estimation elapsed time: 3.09470582008362
COTAN genes' COEX estimation: START
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 available 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.3027069568634
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00660991668701172
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0427150726318359
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.06948184967041
Total calculations elapsed time: 4.42151379585266
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.48274922370911
Dataset analysis elapsed time: 9.55893707275391
COTAN dataset analysis: DONE
Checking uniformity for the cluster '01_0000' with 613 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0392050743103027
Calculate `GDI`: DONE
Total calculations elapsed time: 3.04662299156189
Calculate GDI dataframe: DONE
Clean plots: START
Total calculations elapsed time: 0.0481059551239014
Clean plots: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0000, with size 613, is uniform
Cluster 01_0000's shift to uniformity: -0.0159698218976119
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.397, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.0333333333333333, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.38103017810239, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.407, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.0233333333333333, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.42514299918547, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.5, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.00333333333333333, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.4503882472624, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 2, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.5015551958585, isUniform = TRUE, clusterSize = 613
cluster 01_0000 is uniform

* checking uniformity of cluster '1' of 2 clusters
Asked to drop 0 genes and 613 cells
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [587] cells
Dataset cleaning elapsed time: 1.0213770866394
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.0963180065155
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: -0.0383377991708709 | max: 19.013469852533 | % negative: 7
Model parameter estimation elapsed time: 3.10998201370239
COTAN genes' COEX estimation: START
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 available 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.27694511413574
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00643706321716309
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0324938297271729
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.0670599937439
Total calculations elapsed time: 4.38293600082397
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.36225700378418
Dataset analysis elapsed time: 9.49361610412598
COTAN dataset analysis: DONE
Checking uniformity for the cluster '01_0001' with 587 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.031013011932373
Calculate `GDI`: DONE
Total calculations elapsed time: 3.04540014266968
Calculate GDI dataframe: DONE
Clean plots: START
Total calculations elapsed time: 0.0472970008850098
Clean plots: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0001, with size 587, is uniform
Cluster 01_0001's shift to uniformity: -0.00930197949808731
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.397, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.0433333333333333, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.38769802050191, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.407, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.0316666666666667, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.42646635133469, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.5, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.00833333333333333, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.46799318315158, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 5, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.52850104542575, isUniform = TRUE, clusterSize = 587
cluster 01_0001 is uniform


Found 2 uniform and  0 non-uniform clusters
Stopping: too few cells left to cluster
Unclustered cell left: 0
The final raw clusterization contains [ 2 ] different clusters: 01_0000, 01_0001
Differential Expression Analysis - START
* DEA on cluster '1' with 613 cells
* DEA on cluster '2' with 587 cells

Total calculations elapsed time: 0.0572969913482666
Differential Expression Analysis - DONE
Applied reordering to clusterization is:
1  ->  1, 2  ->  2
Total calculations elapsed time: 52.4185419082642
Creating cells' uniform clustering: DONE
Applied reordering to clusterization is:
1  ->  1, 2  ->  2
Asked to drop 0 genes and 613 cells
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [587] cells
Dataset cleaning elapsed time: 0.956133127212524
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.04156899452209
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: -0.0383377991708709 | max: 19.013469852533 | % negative: 7
Model parameter estimation elapsed time: 3.01972103118896
COTAN genes' COEX estimation: START
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 available 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.37732291221619
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.0063478946685791
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0339310169219971
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.992995023727417
Total calculations elapsed time: 4.41059684753418
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.43427801132202
Dataset analysis elapsed time: 9.41013216972351
COTAN dataset analysis: DONE
Checking uniformity for the cluster 'Cluster_2' with 587 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0335168838500977
Calculate `GDI`: DONE
Total calculations elapsed time: 3.09170699119568
Calculate GDI dataframe: DONE
Clean plots: START
Total calculations elapsed time: 0.0455141067504883
Clean plots: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster Cluster_2, with size 587, is uniform
Cluster Cluster_2's shift to uniformity: -0.00930197949808731
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.397, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.0433333333333333, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.38769802050191, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.407, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.0316666666666667, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.42646635133469, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.5, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.00833333333333333, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.46799318315158, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 5, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.52850104542575, isUniform = TRUE, clusterSize = 587
Differential Expression Analysis - START
* DEA on cluster '-1' with 50 cells
* DEA on cluster '1' with 611 cells
* DEA on cluster '2' with 539 cells

Total calculations elapsed time: 0.0580549240112305
Differential Expression Analysis - DONE
Applied reordering to clusterization is:
1  ->  2, 2  ->  1, -1  ->  -1
Applied reordering to clusterization is:
1  ->  1, 2  ->  2, -1  ->  -1
Creating cells' uniform clustering: START
In iteration 1 the number of cells to re-cluster is 1200 cells belonging to 0 clusters
Asked to drop no genes or cells
Creating new clusterization: START
Elaborating Reduced dimensionality Data Matrix - START
Running genes' selection: START
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0320789813995361
Calculate `GDI`: DONE
Total calculations elapsed time: 3.12031197547913
Calculate GDI dataframe: DONE
Included 186 genes with GDI below 1.3
Selected 600 genes using HGDI selector
Running genes' selection: DONE
Elaborating PCA - START
Elaborating PCA - DONE
Total calculations elapsed time: 6.1804780960083
Elaborating Reduced dimensionality Data Matrix - DONE
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
Total calculations elapsed time: 7.15892505645752
Creating new clusterization: DONE
Using passed in clusterization
Creating PDF UMAP in file:  /tmp/RtmpuQrPnJ/test/reclustering/pdf_umap_1.pdf
UMAP plot: START
Calculating UMAP: START
21:03:24 UMAP embedding parameters a = 0.9922 b = 1.112
21:03:24 Read 1200 rows and found 40 numeric columns
21:03:24 Using Annoy for neighbor search, n_neighbors = 30
21:03:24 Building Annoy index with metric = cosine, n_trees = 50
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
21:03:24 Writing NN index file to temp file /tmp/RtmpuQrPnJ/file384142b7e1dc
21:03:24 Searching Annoy index using 1 thread, search_k = 3000
21:03:24 Annoy recall = 100%
21:03:25 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
21:03:27 Initializing from normalized Laplacian + noise (using RSpectra)
21:03:27 Commencing optimization for 500 epochs, with 42644 positive edges
21:03:27 Using rng type: pcg
Using method 'umap'
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
21:03:30 Optimization finished
Calculating UMAP: DONE
Total calculations elapsed time: 6.11458015441895
UMAP plot: DONE
UMAP plot: START
Calculating UMAP: START
21:03:30 UMAP embedding parameters a = 0.9922 b = 1.112
21:03:30 Read 1200 rows and found 40 numeric columns
21:03:30 Using Annoy for neighbor search, n_neighbors = 30
21:03:30 Building Annoy index with metric = cosine, n_trees = 50
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
21:03:30 Writing NN index file to temp file /tmp/RtmpuQrPnJ/file3841ef61e7e
21:03:30 Searching Annoy index using 1 thread, search_k = 3000
21:03:31 Annoy recall = 100%
21:03:31 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
21:03:33 Initializing from normalized Laplacian + noise (using RSpectra)
21:03:33 Commencing optimization for 500 epochs, with 42644 positive edges
21:03:33 Using rng type: pcg
Using method 'umap'
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
21:03:36 Optimization finished
Calculating UMAP: DONE
Total calculations elapsed time: 6.0181450843811
UMAP plot: DONE
* checking uniformity of cluster '1' of 2 clusters
Asked to drop 0 genes and 600 cells
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [600] cells
Dataset cleaning elapsed time: 0.970396041870117
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.03861808776855
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: -0.0387079873074399 | max: 19.3963889884382 | % negative: 6.5
Model parameter estimation elapsed time: 3.01619696617126
COTAN genes' COEX estimation: START
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.23048090934753
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.010390043258667
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0446169376373291
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.978039026260376
Total calculations elapsed time: 4.26352691650391
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.22887897491455
Dataset analysis elapsed time: 9.21547198295593
COTAN dataset analysis: DONE
Checking uniformity for the cluster '01_0001' with 600 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.03407883644104
Calculate `GDI`: DONE
Total calculations elapsed time: 2.94425010681152
Calculate GDI dataframe: DONE
Clean plots: START
Total calculations elapsed time: 0.0452220439910889
Clean plots: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0001, with size 600, is uniform
Cluster 01_0001's shift to uniformity: -0.0123041937718968
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.397, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.0383333333333333, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.3846958062281, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.407, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.035, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.42417208332107, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.5, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.00666666666666667, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.45335006989431, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 4, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.51519118633036, isUniform = TRUE, clusterSize = 600
cluster 01_0001 is uniform

* checking uniformity of cluster '2' of 2 clusters
Asked to drop 0 genes and 600 cells
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [600] cells
Dataset cleaning elapsed time: 0.847238063812256
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 1.90716195106506
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: -0.0329292435898063 | max: 10.0824885336881 | % negative: 3.66666666666667
Model parameter estimation elapsed time: 2.71428298950195
COTAN genes' COEX estimation: START
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.13145279884338
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.0132851600646973
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0436749458312988
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.08342003822327
Total calculations elapsed time: 4.27183294296265
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.2344229221344
Dataset analysis elapsed time: 8.79594397544861
COTAN dataset analysis: DONE
Checking uniformity for the cluster '01_0002' with 600 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0361831188201904
Calculate `GDI`: DONE
Total calculations elapsed time: 2.77078819274902
Calculate GDI dataframe: DONE
Clean plots: START
Total calculations elapsed time: 0.0476889610290527
Clean plots: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0002, with size 600, is uniform
Cluster 01_0002's shift to uniformity: -0.0203364702469133
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.397, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.0366666666666667, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.37666352975309, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.407, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.0316666666666667, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.42307834225193, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.5, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.00333333333333333, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.43781655605474, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 2, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.50180130506577, isUniform = TRUE, clusterSize = 600
cluster 01_0002 is uniform


Found 2 uniform and  0 non-uniform clusters
Stopping: too few cells left to cluster
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

Total calculations elapsed time: 0.0470907688140869
Differential Expression Analysis - DONE
Applied reordering to clusterization is:
1  ->  1, 2  ->  2
Total calculations elapsed time: 52.8257308006287
Creating cells' uniform clustering: DONE
findClustersMarkers - START
Log Fold Change Analysis - START
* Analysis of cluster: '1'
* Analysis of cluster: '2'

Total calculations elapsed time: 1.90795683860779
Log Fold Change Analysis - DONE
Total calculations elapsed time: 1.91966485977173
findClustersMarkers - DONE
findClustersMarkers - START
Log Fold Change Analysis - START
* Analysis of cluster: '1'
* Analysis of cluster: '2'

Total calculations elapsed time: 1.94528293609619
Log Fold Change Analysis - DONE
Total calculations elapsed time: 1.95456409454346
findClustersMarkers - DONE
findClustersMarkers - START
Differential Expression Analysis - START
* DEA on cluster '1' with 613 cells
* DEA on cluster '2' with 587 cells

Total calculations elapsed time: 0.0586230754852295
Differential Expression Analysis - DONE
Log Fold Change Analysis - START
* Analysis of cluster: '1'
* Analysis of cluster: '2'

Total calculations elapsed time: 1.93155384063721
Log Fold Change Analysis - DONE
Total calculations elapsed time: 1.99557900428772
findClustersMarkers - DONE
[1] "Tested cluster: 1"
Asked to drop 0 genes and 587 cells
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [613] cells
Dataset cleaning elapsed time: 0.91259503364563
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.08431887626648
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: -0.0334270578203965 | max: 8.68674670630126 | % negative: 4.16666666666667
Model parameter estimation elapsed time: 3.03630185127258
COTAN genes' COEX estimation: START
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 available 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.31905484199524
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00850105285644531
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0376400947570801
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.01555895805359
Total calculations elapsed time: 4.38075494766235
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.34038114547729
Dataset analysis elapsed time: 9.28927803039551
COTAN dataset analysis: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0312027931213379
Calculate `GDI`: DONE
Total calculations elapsed time: 3.08256506919861
Calculate GDI dataframe: DONE
[1] "Tested cluster: 2"
Asked to drop 0 genes and 613 cells
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [587] cells
Dataset cleaning elapsed time: 1.01269888877869
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.13010501861572
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: -0.0383377991708709 | max: 19.013469852533 | % negative: 7
Model parameter estimation elapsed time: 3.20253705978394
COTAN genes' COEX estimation: START
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 available 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.51959300041199
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00623202323913574
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0444750785827637
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.09131598472595
Total calculations elapsed time: 4.66161608695984
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.68692207336426
Dataset analysis elapsed time: 9.90215802192688
COTAN dataset analysis: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0329239368438721
Calculate `GDI`: DONE
Total calculations elapsed time: 3.0064549446106
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
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.98919892311096
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Initializing `COTAN` meta-data
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Dataset cleaning elapsed time: 1.00678420066833
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.12933707237244
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Model parameter estimation elapsed time: 2.87894797325134
COTAN genes' COEX estimation not requested
Only genes' COEX elapsed time: 0.835690021514893
Dataset analysis elapsed time: 4.72142219543457
COTAN dataset analysis: DONE
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-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 
Negative:  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 
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-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 
Negative:  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 
PC_ 5 
Positive:  g-000518, g-000108, g-000186, g-000170, g-000401, g-000337, g-000047, g-000599, g-000432, g-000578 
	   g-000042, g-000065, g-000493, g-000261, g-000533, g-000256, g-000560, g-000596, g-000368, g-000381 
	   g-000535, g-000338, g-000215, g-000159, g-000365, g-000234, g-000173, g-000387, g-000225, g-000272 
Negative:  g-000451, g-000339, g-000295, g-000328, g-000544, g-000061, g-000227, g-000391, g-000556, g-000237 
	   g-000067, g-000165, g-000449, g-000591, g-000087, g-000129, g-000197, g-000203, g-000487, g-000505 
	   g-000333, g-000029, g-000271, g-000064, g-000583, g-000156, g-000448, g-000153, g-000526, g-000393 
Computing nearest neighbor graph
Computing SNN
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck

Number of nodes: 1200
Number of edges: 55940

Running Louvain algorithm with multilevel refinement...
0%   10   20   30   40   50   60   70   80   90   100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Maximum modularity in 10 random starts: 0.5920
Number of communities: 4
Elapsed time: 0 seconds
Clusterizations found: RNA_snn_res.0.8, seurat_clusters
Conditions found: orig.ident
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Dataset cleaning elapsed time: 0.988193988800049
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.11480093002319
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Model parameter estimation elapsed time: 3.01102495193481
COTAN genes' COEX estimation not requested
Only genes' COEX elapsed time: 0.867878913879395
Dataset analysis elapsed time: 4.86709785461426
COTAN dataset analysis: DONE
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Dataset cleaning elapsed time: 1.00654196739197
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.34881901741028
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Model parameter estimation elapsed time: 3.36617994308472
COTAN genes' COEX estimation: START
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 available 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.81864404678345
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.0103309154510498
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0579371452331543
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.08049702644348
Total calculations elapsed time: 4.96740913391113
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 6.01465320587158
Dataset analysis elapsed time: 10.3873751163483
COTAN dataset analysis: DONE
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
Total calculations elapsed time: 0.00832486152648926
calculating PValues: DONE
Number of selected secondary markers: 6
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.00874900817871094
Calculate `GDI`: DONE
Number of columns (V set - secondary markers): 6
Number of rows (U set): 60
Total calculations elapsed time: 1.15077805519104
Calculating gene co-expression space - DONE
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
Total calculations elapsed time: 0.0111119747161865
calculating PValues: DONE
Number of selected secondary markers: 109
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0106050968170166
Calculate `GDI`: DONE
Number of columns (V set - secondary markers): 109
Number of rows (U set): 60
Total calculations elapsed time: 1.0714099407196
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
Total calculations elapsed time: 0.00774788856506348
calculating PValues: DONE
Number of selected secondary markers: 109
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0109260082244873
Calculate `GDI`: DONE
Number of columns (V set - secondary markers): 109
Number of rows (U set): 60
Total calculations elapsed time: 1.09635615348816
Calculating gene co-expression space - DONE
Total calculations elapsed time: 1.19138598442078
Establishing gene clusters - DONE
Initializing `COTAN` meta-data
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Dataset cleaning elapsed time: 0.947901010513306
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.34206199645996
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Model parameter estimation elapsed time: 3.30254602432251
COTAN genes' COEX estimation not requested
Only genes' COEX elapsed time: 1.06635093688965
Dataset analysis elapsed time: 5.31679797172546
COTAN dataset analysis: DONE
Differential Expression Analysis - START
* DEA on cluster '1' with 608 cells
* DEA on cluster '2' with 592 cells

Total calculations elapsed time: 0.0619261264801025
Differential Expression Analysis - DONE
Log Fold Change Analysis - START
* Analysis of cluster: '1'
* Analysis of cluster: '2'

Total calculations elapsed time: 2.08469915390015
Log Fold Change Analysis - DONE
clustersDeltaExpression - START
Handling cluster '1' with mean UDE 1.09087489420394
Handling cluster '2' with mean UDE 0.906669027574328
clustersDeltaExpression - DONE
In group G1 there are 3 detected over 3 genes
In group G2 there are 2 detected over 3 genes
In group G3 there are 5 detected over 5 genes
Merging cells' uniform clustering: START
The merge algorithm will use 2 passes
Start merging nearest clusters - the main threshold is: 1.297
Start merging nearest clusters: iteration 1
Differential Expression Analysis - START
* DEA on cluster '1' with 608 cells
* DEA on cluster '2' with 592 cells

Total calculations elapsed time: 0.0585050582885742
Differential Expression Analysis - DONE
Updating check results for the 0 already tested pairs to align to new checker
1 new clusters pairs to be tested for merging:
c("1", "2")
*1_2-merge
Asked to drop no genes or cells
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Dataset cleaning elapsed time: 0.957567930221558
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.31399917602539
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: 0.219880139527933 | max: 6.08024473714308 | % negative: 0
Model parameter estimation elapsed time: 3.38426685333252
COTAN genes' COEX estimation: START
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 available 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.65379905700684
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00911092758178711
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0537509918212891
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.947988986968994
Total calculations elapsed time: 4.66464996337891
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.65587019920349
Dataset analysis elapsed time: 9.99770498275757
COTAN dataset analysis: DONE
Checking uniformity for the cluster '1_2-merge' with 1200 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0341689586639404
Calculate `GDI`: DONE
Total calculations elapsed time: 3.07646799087524
Calculate GDI dataframe: DONE
Clean plots: START
Total calculations elapsed time: 0.0452840328216553
Clean plots: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 1_2-merge, with size 1200, is not uniform
Cluster 1_2-merge's shift to uniformity: 0.273171158219277
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.297, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.696666666666667, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.57017115821928, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.307, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.665, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.61241976619328, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.4, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.411666666666667, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.64049230721132, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 247, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.66511539543561, isUniform = FALSE, clusterSize = 1200
Clusters 1 and 2 cannot be merged
Finished the first batch - no merges were executed
Loop calculations elapsed time: 18.1778028011322
Start merging nearest clusters: iteration 2
Differential Expression Analysis - START
* DEA on cluster '1' with 608 cells
* DEA on cluster '2' with 592 cells

Total calculations elapsed time: 0.0515658855438232
Differential Expression Analysis - DONE
Updating check results for the 1 already tested pairs to align to new checker
0 new clusters pairs to be tested for merging:

No clusters will be merged
None of the remaining tested cluster pairs could be merged
Executed all merges for threshold 1.297 out of 1 checks
Start merging nearest clusters - the main threshold is: 1.397
Start merging nearest clusters: iteration 3
Differential Expression Analysis - START
* DEA on cluster '1' with 608 cells
* DEA on cluster '2' with 592 cells

Total calculations elapsed time: 0.0499701499938965
Differential Expression Analysis - DONE
Updating check results for the 1 already tested pairs to align to new checker
0 new clusters pairs to be tested for merging:

No clusters will be merged
None of the remaining tested cluster pairs could be merged
Executed all merges for threshold 1.397 out of 1 checks
The final merged clusterization contains [2] different clusters: 1, 2
Differential Expression Analysis - START
* DEA on cluster '1' with 608 cells
* DEA on cluster '2' with 592 cells

Total calculations elapsed time: 0.0471110343933105
Differential Expression Analysis - DONE
Applied reordering to clusterization is:
1  ->  1, 2  ->  2
Total calculations elapsed time: 20.3603990077972
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'

Total calculations elapsed time: 2.09596014022827
Log Fold Change Analysis - DONE
[1] "Tested cluster: 1"
Asked to drop 0 genes and 592 cells
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [608] cells
Dataset cleaning elapsed time: 1.02784895896912
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.14790201187134
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: -0.0337582853418583 | max: 10.9273420652533 | % negative: 4.33333333333333
Model parameter estimation elapsed time: 3.12254500389099
COTAN genes' COEX estimation: START
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 available 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.23764896392822
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00642800331115723
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0407869815826416
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.995471000671387
Total calculations elapsed time: 4.28033494949341
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 5.30439710617065
Dataset analysis elapsed time: 9.45479106903076
COTAN dataset analysis: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0339498519897461
Calculate `GDI`: DONE
Total calculations elapsed time: 2.81459498405457
Calculate GDI dataframe: DONE
[1] "Tested cluster: 2"
Asked to drop 0 genes and 608 cells
COTAN dataset analysis: START
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [592] cells
Dataset cleaning elapsed time: 1.00355386734009
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
Parallel dispersion Newton-Raphson: used up to 4 iterations
Total calculations elapsed time: 2.04840087890625
Estimate `dispersion`: DONE
Estimate `dispersion`: DONE
`dispersion` | min: -0.0390500444479481 | max: 19.2061656068405 | % negative: 7.33333333333333
Model parameter estimation elapsed time: 2.85454106330872
COTAN genes' COEX estimation: START
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 available 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: 5.00582313537598
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.260423898696899
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.185270071029663
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.919557094573975
Total calculations elapsed time: 6.37107419967651
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 7.53525996208191
Dataset analysis elapsed time: 11.3933548927307
COTAN dataset analysis: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Total calculations elapsed time: 0.0534458160400391
Calculate `GDI`: DONE
Total calculations elapsed time: 3.38691186904907
Calculate GDI dataframe: DONE
Dispersion bisection: straddling used 2 iterations
Dispersion bisection: used 8 iterations
Dispersion bisection: straddling used 1 iterations
Dispersion bisection: used 9 iterations
Dispersion Newton-Raphson: used up to 4 iterations
Dispersion Newton-Raphson: used up to 4 iterations
parallel dispersion bisection: straddling used up to 2 iterations
Parallel dispersion bisection: used up to 9 iterations
Parallel dispersion Newton-Raphson: used up to 4 iterations
[ FAIL 0 | WARN 3 | SKIP 0 | PASS 572 ]

[ FAIL 0 | WARN 3 | SKIP 0 | PASS 572 ]
> 
> proc.time()
   user  system elapsed 
653.796  10.881 671.557

Example timings

COTAN.Rcheck/COTAN-Ex.timings

name	user	system	elapsed
COTAN_Legacy	0.002	0.001	0.004
COTAN_ObjectCreation	10.690	0.304	11.070
CalculatingCOEX	55.689	1.145	57.105
ClustersList	0.009	0.003	0.012
Conversions	5.021	0.128	5.175
GenesStatistics	11.139	0.309	11.360
HandleMetaData	0.093	0.009	0.105
HandlingClusterizations	43.884	1.174	44.691
HandlingConditions	0.114	0.007	0.122
HeatmapPlots	34.730	0.941	35.883
Installing_torch	0.000	0.001	0.001
LoggingFunctions	0.002	0.001	0.004
ParametersEstimations	23.568	0.737	24.414
RawDataCleaning	5.099	0.332	5.268
RawDataGetters	0.091	0.010	0.103
UniformClusters	117.537	2.626	119.597
getColorsVector	0.001	0.015	0.001