Galfrè Silvia Giulia
| Snapshot Date: 2025-09-21 13:45 -0400 (Sun, 21 Sep 2025) |
| git_url: https://git.bioconductor.org/packages/COTAN |
| git_branch: devel |
| git_last_commit: bb9a730 |
| git_last_commit_date: 2025-07-15 13:06:23 -0400 (Tue, 15 Jul 2025) |
| Back to Multiple platform build/check report for BioC 3.22: simplified long |
|
This page was generated on 2025-09-22 12:03 -0400 (Mon, 22 Sep 2025).
| 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" | 4814 |
| lconway | macOS 12.7.1 Monterey | x86_64 | 4.5.1 Patched (2025-09-10 r88807) -- "Great Square Root" | 4603 |
| kjohnson3 | macOS 13.7.7 Ventura | arm64 | 4.5.1 Patched (2025-09-10 r88807) -- "Great Square Root" | 4547 |
| taishan | Linux (openEuler 24.03 LTS) | aarch64 | 4.5.0 (2025-04-11) -- "How About a Twenty-Six" | 4553 |
| Click on any hostname to see more info about the system (e.g. compilers) (*) as reported by 'uname -p', except on Windows and Mac OS X | ||||
| Package 453/2333 | Hostname | OS / Arch | INSTALL | BUILD | CHECK | BUILD BIN | ||||||||
| COTAN 2.9.5 (landing page) Galfrè Silvia Giulia
| nebbiolo2 | Linux (Ubuntu 24.04.3 LTS) / x86_64 | OK | OK | TIMEOUT | |||||||||
| lconway | macOS 12.7.1 Monterey / x86_64 | OK | OK | OK | OK | | ||||||||
| kjohnson3 | macOS 13.7.7 Ventura / arm64 | OK | OK | OK | OK | | ||||||||
| taishan | Linux (openEuler 24.03 LTS) / aarch64 | OK | OK | OK | ||||||||||
|
To the developers/maintainers of the COTAN package: - Allow up to 24 hours (and sometimes 48 hours) for your latest push to git@git.bioconductor.org:packages/COTAN.git to reflect on this report. See Troubleshooting Build Report for more information. - Use the following Renviron settings to reproduce errors and warnings. - If 'R CMD check' started to fail recently on the Linux builder(s) over a missing dependency, add the missing dependency to 'Suggests:' in your DESCRIPTION file. See Renviron.bioc for more information. |
| Package: COTAN |
| Version: 2.9.5 |
| Command: /home/biocbuild/bbs-3.22-bioc/R/bin/R CMD check --install=check:COTAN.install-out.txt --library=/home/biocbuild/bbs-3.22-bioc/R/site-library --timings COTAN_2.9.5.tar.gz |
| StartedAt: 2025-09-21 22:28:27 -0400 (Sun, 21 Sep 2025) |
| EndedAt: 2025-09-21 23:08:27 -0400 (Sun, 21 Sep 2025) |
| EllapsedTime: 2400.1 seconds |
| RetCode: None |
| Status: TIMEOUT |
| CheckDir: COTAN.Rcheck |
| Warnings: NA |
##############################################################################
##############################################################################
###
### Running command:
###
### /home/biocbuild/bbs-3.22-bioc/R/bin/R CMD check --install=check:COTAN.install-out.txt --library=/home/biocbuild/bbs-3.22-bioc/R/site-library --timings COTAN_2.9.5.tar.gz
###
##############################################################################
##############################################################################
* using log directory ‘/home/biocbuild/bbs-3.22-bioc/meat/COTAN.Rcheck’
* using R version 4.5.1 Patched (2025-08-23 r88802)
* using platform: x86_64-pc-linux-gnu
* R was compiled by
gcc (Ubuntu 13.3.0-6ubuntu2~24.04) 13.3.0
GNU Fortran (Ubuntu 13.3.0-6ubuntu2~24.04) 13.3.0
* running under: Ubuntu 24.04.3 LTS
* using session charset: UTF-8
* checking for file ‘COTAN/DESCRIPTION’ ... OK
* checking extension type ... Package
* this is package ‘COTAN’ version ‘2.9.5’
* 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 loading without being on the library search path ... 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 4041.270 17.554 410.295
CalculatingCOEX 1490.455 2.851 116.696
HeatmapPlots 1078.592 2.341 76.217
HandlingClusterizations 387.745 2.204 61.365
GenesStatistics 364.379 0.983 26.849
COTAN_ObjectCreation 354.796 0.870 23.444
ParametersEstimations 31.564 1.013 32.105
RawDataCleaning 7.509 0.021 7.690
Conversions 6.521 0.137 6.899
* checking for unstated dependencies in ‘tests’ ... OK
* checking tests ...
Running ‘outputTestDatasetCreation.R’
Running ‘spelling.R’
Running ‘testthat.R’
COTAN.Rcheck/00install.out
############################################################################## ############################################################################## ### ### Running command: ### ### /home/biocbuild/bbs-3.22-bioc/R/bin/R CMD INSTALL COTAN ### ############################################################################## ############################################################################## * installing to library ‘/home/biocbuild/bbs-3.22-bioc/R/site-library’ * installing *source* package ‘COTAN’ ... ** this is package ‘COTAN’ version ‘2.9.5’ ** 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)
COTAN.Rcheck/tests/outputTestDatasetCreation.Rout
R version 4.5.1 Patched (2025-08-23 r88802) -- "Great Square Root"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
>
> # Creates the files to be reloaded by the tests for comparisons
> library(zeallot)
>
> outputTestDatasetCreation <- function(testsDir = file.path("tests",
+ "testthat")) {
+ utils::data("test.dataset", package = "COTAN")
+ options(parallelly.fork.enable = TRUE)
+
+ obj <- COTAN(raw = test.dataset)
+ obj <- initializeMetaDataset(obj, GEO = " ",
+ sequencingMethod = "artificial",
+ sampleCondition = "test")
+
+ obj <- proceedToCoex(obj, cores = 6L, saveObj = FALSE)
+ #saveRDS(obj, file = file.path(testsDir,"temp.RDS"))
+
+ cells.names.test <- getCells(obj)[c(1L:10L, 591L:610L, 991L:1000L)]
+ genes.names.test <- getGenes(obj)[c(1L:10L, 291L:310L, 591L: 600L)]
+ saveRDS(cells.names.test, file.path(testsDir, "cells.names.test.RDS"))
+ saveRDS(genes.names.test, file.path(testsDir, "genes.names.test.RDS"))
+
+ dispersion.test <- getDispersion(obj)[genes.names.test]
+ saveRDS(dispersion.test, file.path(testsDir, "dispersion.test.RDS"))
+
+ raw.norm.test <- getNuNormData(obj)[genes.names.test, cells.names.test]
+ saveRDS(raw.norm.test, file.path(testsDir, "raw.norm.test.RDS"))
+
+ coex.test <- getGenesCoex(obj, genes = genes.names.test, zeroDiagonal = FALSE)
+ saveRDS(coex.test, file.path(testsDir, "coex.test.RDS"))
+
+ lambda.test <- getLambda(obj)[genes.names.test]
+ saveRDS(lambda.test, file.path(testsDir, "lambda.test.RDS"))
+
+ GDI.test <- calculateGDI(obj)
+ GDI.test <- GDI.test[genes.names.test, ]
+ saveRDS(GDI.test, file.path(testsDir, "GDI.test.RDS"))
+
+ nu.test <- getNu(obj)[cells.names.test]
+ saveRDS(nu.test, file.path(testsDir, "nu.test.RDS"))
+
+ pval.test <- calculatePValue(obj, geneSubsetCol = genes.names.test)
+ saveRDS(pval.test, file.path(testsDir, "pval.test.RDS"))
+
+ GDIThreshold <- 1.46
+ initialResolution <- 0.8
+ clusters <- cellsUniformClustering(obj, GDIThreshold = GDIThreshold,
+ initialResolution = initialResolution,
+ cores = 6L, saveObj = FALSE)[["clusters"]]
+ saveRDS(clusters, file.path(testsDir, "clusters1.RDS"))
+
+ coexDF <- DEAOnClusters(obj, clusters = clusters)
+ obj <- addClusterization(obj, clName = "clusters",
+ clusters = clusters, coexDF = coexDF)
+
+ saveRDS(coexDF[genes.names.test, ],
+ file.path(testsDir, "coex.test.cluster1.RDS"))
+
+ pvalDF <- pValueFromDEA(coexDF, getNumCells(obj), adjustmentMethod = "none")
+
+ saveRDS(pvalDF[genes.names.test, ],
+ file.path(testsDir, "pval.test.cluster1.RDS"))
+
+ c(mergedClusters, mCoexDF) %<-%
+ mergeUniformCellsClusters(objCOTAN = obj,
+ clusters = NULL,
+ GDIThreshold = GDIThreshold,
+ cores = 6L,
+ distance = "cosine",
+ hclustMethod = "ward.D2",
+ saveObj = FALSE)
+
+ saveRDS(mergedClusters[genes.names.test],
+ file.path(testsDir, "cluster_data_merged.RDS"))
+ }
>
> proc.time()
user system elapsed
0.184 0.045 0.213
COTAN.Rcheck/tests/spelling.Rout
R version 4.5.1 Patched (2025-08-23 r88802) -- "Great Square Root"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu
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.132 0.046 0.164
COTAN.Rcheck/tests/testthat.Rout
R version 4.5.1 Patched (2025-08-23 r88802) -- "Great Square Root"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> Sys.setenv(R_TESTS = "")
> library(testthat)
> library(COTAN)
> test_check("COTAN")
Setting new log level to 4
Initializing `COTAN` meta-data
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.9033203125 | max: 4.6796875 | % negative: 10
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.53793787956238
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00101447105407715
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00616192817687988
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.10036563873291
Total calculations elapsed time: 4.64547991752625
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
`dispersion` | min: 0.9033203125 | max: 4.6796875 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 7 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 1.75595238095238 | median: 1.07174634176587 | mean: 1.07174634176587
Estimate `dispersion`/`nu`: START
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 1.0362548828125 | max: 4.60986328125 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0265938895089288 | median: 0.0144680038331048 | mean: 0.0144680038331048
`nu` mean: 1.69633192486233
Marginal errors | max: 1.95570586131367 | median 1.32068160171502 | mean: 1.33375826507259
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.058837890625 | max: 3.528076171875 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.416683423613994 | median: 0.239880630367975 | mean: 0.239880630367975
`nu` mean: 0.823197206753982
Marginal errors | max: 0.836359531101206 | median 0.703684202571891 | mean: 0.645537958989614
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.32879638671875 | max: 4.0302734375 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.164237872898673 | median: 0.0955985184389135 | mean: 0.0955985184389135
`nu` mean: 1.06863935445976
Marginal errors | max: 0.259872988828242 | median 0.213703042752633 | mean: 0.197386407582083
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2294921875 | max: 3.8720703125 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.055185575120883 | median: 0.0319991762044448 | mean: 0.0319991762044448
`nu` mean: 0.976813601083562
Marginal errors | max: 0.0951586919577032 | median 0.079429703709466 | mean: 0.0724140148396648
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2637939453125 | max: 3.929443359375 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0196211148938294 | median: 0.01138609597457 | mean: 0.01138609597457
`nu` mean: 1.00823501891926
Marginal errors | max: 0.0327747321002292 | median 0.0272104747849538 | mean: 0.0248963830312038
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.25177001953125 | max: 3.90966796875 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.00670099066960717 | median: 0.00388888266671264 | mean: 0.00388888266671264
`nu` mean: 0.997187891997105
Marginal errors | max: 0.0114324509186883 | median 0.0094232649770607 | mean: 0.00863113610779571
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.25592041015625 | max: 3.91650390625 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.00230093811689414 | median: 0.00132529122122246 | mean: 0.00132529122122246
`nu` mean: 1.00097564689567
Marginal errors | max: 0.00387133150664631 | median 0.0031091017608853 | mean: 0.00286071175800213
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2545166015625 | max: 3.914306640625 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.000837011646904529 | median: 0.000470837393363011 | mean: 0.000470837393363011
`nu` mean: 0.999633825746458
Marginal errors | max: 0.00122501723202006 | median 0.00102126435760308 | mean: 0.000943992659051318
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2550048828125 | max: 3.9150390625 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.000209227351122054 | median: 0.000122070312500028 | mean: 0.000122070312500028
`nu` mean: 1.00008715703862
Marginal errors | max: 0.000364602956581805 | median 0.000313956936819793 | mean: 0.000282899574318485
Estimate `dispersion`/`nu`: DONE
Estimate `dispersion`/`nu`: START
Initializing `COTAN` meta-data
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.9033203125 | max: 4.6796875 | % negative: 10
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.87749671936035
Calculating genes' COEX normalization factor
Fraction of genes' with very low expected contingency tables: 0.181818181818182
Calculate genes' normalization factor elapsed time: 0.000800371170043945
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0036618709564209
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.26865077018738
Total calculations elapsed time: 5.15060973167419
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: 3.88468337059021
Calculating cells' COEX normalization factor
Fraction of cells' with very low expected contingency tables: 0
Calculate cells' normalization factor elapsed time: 0.000895023345947266
Retrieving observed cells' yes/yes contingency table
calculating YY.. done
Observed cells' contingency table elapsed time: 0.00264501571655273
Estimating cells' COEX
Calculate cells' COEX elapsed time: 1.09895944595337
Total calculations elapsed time: 4.98718285560608
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
Calculate `GDI`: DONE
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
Calculate `GDI`: DONE
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
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
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
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
Elaborating Reduced dimensionality Data Matrix - DONE
Elaborating Reduced dimensionality Data Matrix - START
Elaborating COEX Eigen Vectors - START
Elaborating COEX Eigen Vectors - DONE
Elaborating Reduced dimensionality Data Matrix - DONE
Elaborating Reduced dimensionality Data Matrix - START
Elaborating COEX Eigen Vectors - START
Elaborating COEX Eigen Vectors - DONE
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
`dispersion` | min: 0.903564453125 | max: 4.679443359375 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 1.75719246031746 | median: 1.07229953342014 | mean: 1.07229953342014
`nu` mean: 1.68489292689732
Marginal errors | max: 1.73564890252257 | median 1.37996360874076 | mean: 1.32180348113228
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.0655517578125 | max: 3.5439453125 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.402649984216273 | median: 0.231868788425666 | mean: 0.231868788425666
`nu` mean: 0.829218804209393
Marginal errors | max: 0.80321315986594 | median 0.677497553540581 | mean: 0.619375430892821
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.3260498046875 | max: 4.026123046875 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.158004893526231 | median: 0.0919692884670312 | mean: 0.0919692884670312
`nu` mean: 1.0660356050592
Marginal errors | max: 0.250724014302326 | median 0.206232152124436 | mean: 0.190425623677198
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.23040771484375 | max: 3.8736572265625 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0532774732102337 | median: 0.0308837890624999 | mean: 0.0308837890624999
`nu` mean: 0.977606315852266
Marginal errors | max: 0.0916983669060123 | median 0.0765266929824948 | mean: 0.0697593208689684
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.26348876953125 | max: 3.928955078125 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0189966206463044 | median: 0.0110199320575908 | mean: 0.0110199320575908
`nu` mean: 1.00797668858871
Marginal errors | max: 0.0317151207459254 | median 0.0262702142278251 | mean: 0.0240886952086962
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2518310546875 | max: 3.9097900390625 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.00670088501353994 | median: 0.00388888101583662 | mean: 0.00388888101583662
`nu` mean: 0.997187996002297
Marginal errors | max: 0.0113693316356223 | median 0.00939669372836249 | mean: 0.00860715734056932
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2559814453125 | max: 3.9166259765625 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.00251007446998996 | median: 0.00144735987374958 | mean: 0.00144735987374958
`nu` mean: 1.00106271459624
Marginal errors | max: 0.00406746973787442 | median 0.00343393462175534 | mean: 0.00313496757119527
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.25445556640625 | max: 3.9140625 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.000837027590858019 | median: 0.000488281249999889 | mean: 0.000488281249999889
`nu` mean: 0.999651253659142
Marginal errors | max: 0.00143433714371533 | median 0.00116636244706747 | mean: 0.00109289166947839
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2550048828125 | max: 3.9150390625 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 1 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.000209227688885871 | median: 0.0001220703125 | mean: 0.0001220703125
`nu` mean: 1.00008715737639
Marginal errors | max: 0.000379524846207957 | median 0.000325685250686547 | mean: 0.000295532844332236
Estimate `dispersion`/`nu`: DONE
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.64940142631531
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.000613212585449219
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00160312652587891
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.11861443519592
Total calculations elapsed time: 4.77023220062256
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: 3.67435765266418
Calculating cells' COEX normalization factor
Calculate cells' normalization factor elapsed time: 0.000790834426879883
Retrieving observed cells' yes/yes contingency table
calculating YY.. done
Observed cells' contingency table elapsed time: 0.00211930274963379
Estimating cells' COEX
Calculate cells' COEX elapsed time: 1.06359767913818
Total calculations elapsed time: 4.74086546897888
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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
Calculate genes' partial COEX: START
Retrieving expected genes' partial contingency table
calculating partial NN.. done
calculating partial NY..YN..YY.. done
Calculating genes' partial COEX normalization factor
Fraction of genes with very low expected contingency tables: 0.325
Retrieving observed genes' yes/yes partial contingency table
calculating partial YY.. done
Estimating genes' partial COEX
Calculate genes' partial COEX: DONE
Calculate cells' partial COEX: START
Retrieving expected cells' partial contingency table
calculating partial NN.. done
calculating partial YN..NY..YY.. done
Calculating cells' partial COEX normalization factor
Fraction of genes with very low expected contingency tables: 0
Retrieving observed cells' yes/yes partial contingency table
calculating partial YY.. done
Estimating cells' partial COEX
Calculate cells' partial COEX: DONE
Asked to drop 2 genes and 0 cells
Asked to drop 0 genes and 4 cells
Asked to drop 2 genes and 2 cells
Attaching package: 'rlang'
The following objects are masked from 'package:testthat':
is_false, is_null, is_true
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
calculating YY.. done
calculating YY.. done
calculating YN..NY..NN..t().. done
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
`dispersion` | min: 0.9033203125 | max: 4.6796875 | % negative: 10
calculating NN.. done
calculating NN.. done
calculating NY..YN..YY..t().. done
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.6039776802063
Calculating genes' COEX normalization factor
Fraction of genes' with very low expected contingency tables: 0.181818181818182
Calculate genes' normalization factor elapsed time: 0.00118780136108398
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00240230560302734
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.07955050468445
Total calculations elapsed time: 4.68711829185486
Calculate genes' COEX (legacy): DONE
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.51546931266785
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.000799417495727539
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00204181671142578
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.06448531150818
Total calculations elapsed time: 4.58279585838318
Calculate genes' COEX (legacy): DONE
Calculate genes' partial COEX: START
Retrieving expected genes' partial contingency table
calculating partial NN.. done
calculating partial NY..YN..YY.. done
Calculating genes' partial COEX normalization factor
Fraction of genes with very low expected contingency tables: 0.1
Retrieving observed genes' yes/yes partial contingency table
calculating partial YY.. done
Estimating genes' partial COEX
Calculate genes' partial COEX: DONE
Calculate genes' partial COEX: START
Retrieving expected genes' partial contingency table
calculating partial NN.. done
calculating partial NY..YN..YY.. done
Calculating genes' partial COEX normalization factor
Fraction of genes with very low expected contingency tables: 0.4
Retrieving observed genes' yes/yes partial contingency table
calculating partial YY.. done
Estimating genes' partial COEX
Calculate genes' partial COEX: DONE
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
calculating YY.. done
calculating YY.. done
calculating NY..YN..NN..t().. done
Estimate `dispersion`/`nu`: START
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
`dispersion` | min: 0.903564453125 | max: 4.679443359375 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 1.75719246031746 | median: 1.07229953342014 | mean: 1.07229953342014
`nu` mean: 1.68489292689732
Marginal errors | max: 0.25535328937316 | median 0.0807577993228135 | mean: 0.101980750205762
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
`dispersion` | min: 1.037109375 | max: 4.6107177734375 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0273438105507502 | median: 0.0148852611818011 | mean: 0.0148852611818011
`nu` mean: 1.69735147626627
Marginal errors | max: 0.00326864580272002 | median 0.00111524657842743 | mean: 0.00131556083122533
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
`dispersion` | min: 1.03887939453125 | max: 4.6097412109375 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0 | median: 0 | mean: 0
`nu` mean: 1.69735147626627
Marginal errors | max: 7.56328383637594e-05 | median 1.72948087246994e-05 | mean: 2.99252342145451e-05
Estimate `dispersion`/`nu`: DONE
calculating NN.. done
calculating NN.. done
calculating YN..NY..YY..t().. done
Calculate cells' COEX (legacy): START
Retrieving expected cells' contingency table
calculating NN.. done
calculating YN..NY..YY..t().. done
Expected cells' contingency table elapsed time: 3.4207022190094
Calculating cells' COEX normalization factor
Fraction of cells' with very low expected contingency tables: 0
Calculate cells' normalization factor elapsed time: 0.00117874145507812
Retrieving observed cells' yes/yes contingency table
calculating YY.. done
Observed cells' contingency table elapsed time: 0.00237369537353516
Estimating cells' COEX
Calculate cells' COEX elapsed time: 0.922086238861084
Total calculations elapsed time: 4.3463408946991
Calculate cells' COEX (legacy): DONE
Calculate cells' partial COEX: START
Retrieving expected cells' partial contingency table
calculating partial NN.. done
calculating partial YN..NY..YY.. done
Calculating cells' partial COEX normalization factor
Fraction of genes with very low expected contingency tables: 0
Retrieving observed cells' yes/yes partial contingency table
calculating partial YY.. done
Estimating cells' partial COEX
Calculate cells' partial COEX: DONE
Calculate cells' partial COEX: START
Retrieving expected cells' partial contingency table
calculating partial NN.. done
calculating partial YN..NY..YY.. done
Calculating cells' partial COEX normalization factor
Fraction of genes with very low expected contingency tables: 0
Retrieving observed cells' yes/yes partial contingency table
calculating partial YY.. done
Estimating cells' partial COEX
Calculate cells' partial COEX: DONE
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
Estimate `dispersion`/`nu`: START
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
`dispersion` | min: 0.903564453125 | max: 4.679443359375 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 1.75719246031746 | median: 1.07229953342014 | mean: 1.07229953342014
`nu` mean: 1.68489292689732
Marginal errors | max: 0.25535328937316 | median 0.0807577993228135 | mean: 0.101980750205762
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
`dispersion` | min: 1.037109375 | max: 4.6107177734375 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0.0273438105507502 | median: 0.0148852611818011 | mean: 0.0148852611818011
`nu` mean: 1.69735147626627
Marginal errors | max: 0.00326864580272002 | median 0.00111524657842743 | mean: 0.00131556083122533
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 3 genes batches
`dispersion` | min: 1.03887939453125 | max: 4.6097412109375 | % negative: 10
Estimate `nu`: START
Effective number of cores used: 1
Executing 5 cells batches
Estimate `nu`: DONE
`nu` change (abs) | max: 0 | median: 0 | mean: 0
`nu` mean: 1.69735147626627
Marginal errors | max: 7.56328383637594e-05 | median 1.72948087246994e-05 | mean: 2.99252342145451e-05
Estimate `dispersion`/`nu`: DONE
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.08433938026428
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00078272819519043
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0020909309387207
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.854230880737305
Total calculations elapsed time: 3.9414439201355
Calculate genes' COEX (legacy): DONE
Calculating S: START
Calculating S: DONE
Calculating G: START
calculating YY.. done
calculating YN..NY..NN..t().. done
calculating NN.. done
calculating NY..YN..YY..t().. done
Estimating G
Calculating G: DONE
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values genome wide on columns and genome wide on rows
calculating PValues: DONE
Using G
Calculating G: START
calculating YY.. done
calculating YN..NY..NN..t().. done
calculating NN.. done
calculating NY..YN..YY..t().. done
Estimating G
Calculating G: DONE
calculating PValues: START
Get p-values on a set of genes on columns and on a set of genes on rows
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
Calculate GDI dataframe: START
Using G
Calculating G: START
calculating YY.. done
calculating YN..NY..NN..t().. done
calculating NN.. done
calculating NY..YN..YY..t().. done
Estimating G
Calculating G: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Calculate `GDI`: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Calculate `GDI`: DONE
Initializing `COTAN` meta-data
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 3.6418514251709
Cotan genes' COEX estimation started
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 16.1512460708618
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.0060126781463623
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0692641735076904
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.02027463912964
Total calculations elapsed time: 17.2467975616455
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 18.2421979904175
Total elapsed time: 23.6339087486267
Initializing `COTAN` meta-data
Condition test
n cells 1200
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 3.66778469085693
Cotan genes' COEX estimation started
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 13.0719764232635
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00855326652526855
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0702660083770752
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.00134682655334
Total calculations elapsed time: 14.1521425247192
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 15.211886882782
Total elapsed time: 20.8286173343658
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and genome wide on rows
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
Initializing `COTAN` meta-data
Condition test
n cells 1200
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 4.73668909072876
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 11.9553077220917
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00842714309692383
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0797066688537598
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.971833229064941
Total calculations elapsed time: 13.0152747631073
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 13.9685442447662
Total elapsed time: 21.4547843933105
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and genome wide on rows
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
Initializing `COTAN` meta-data
Condition test
n cells 1200
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 4.61577248573303
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 15.1869695186615
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00617766380310059
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0517704486846924
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.891450643539429
Total calculations elapsed time: 16.1363682746887
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 17.0348889827728
Total elapsed time: 24.5100457668304
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and genome wide on rows
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
Initializing `COTAN` meta-data
Asked to drop 0 genes and 1160 cells
Cotan analysis functions started
Asked to drop 31 genes and 0 cells
Genes/cells selection done: dropped [31] genes and [0] cells
Working on [569] genes and [40] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.2119140625 | max: 85.5 | % negative: 62.2144112478032
Cotan genes' COEX estimation not requested
Total elapsed time: 5.58332300186157
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: 15.0861020088196
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00670504570007324
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0093696117401123
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.06074070930481
Total calculations elapsed time: 16.1629173755646
Calculate genes' COEX (legacy): DONE
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 13.9471039772034
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00652813911437988
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00851869583129883
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.25555658340454
Total calculations elapsed time: 15.2177073955536
Calculate genes' COEX (legacy): DONE
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 14.9142985343933
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00764656066894531
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0107033252716064
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.997379541397095
Total calculations elapsed time: 15.930027961731
Calculate genes' COEX (legacy): DONE
Hangling COTAN object with condition: test
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and on a set of genes on rows
calculating PValues: DONE
Handling genes type: G2
Handling genes type: G3
Initializing `COTAN` meta-data
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
PCA: START
Elaborating PCA - START
Elaborating PCA - DONE
PCA: DONE
Hierarchical clustering: START
Hierarchical clustering: DONE
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 3.90823745727539
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 11.3671288490295
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00678825378417969
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0665769577026367
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.13510346412659
Total calculations elapsed time: 12.5755975246429
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 13.6686375141144
Total elapsed time: 23.5497903823853
Saving elaborated data locally at: /tmp/RtmpmhC9id/test.cotan.RDS
Creating new clusterization: START
Elaborating Reduced dimensionality Data Matrix - START
Elaborating COEX Eigen Vectors - START
Elaborating COEX Eigen Vectors - DONE
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: 37723
Running Louvain algorithm with multilevel refinement...
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Maximum modularity in 10 random starts: 0.5827
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: 37723
Running Louvain algorithm with multilevel refinement...
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Maximum modularity in 10 random starts: 0.3954
Number of communities: 5
Elapsed time: 0 seconds
Used resolution for Seurat clusterization is: 1.3
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
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 600 genes using HVG_Seurat selector
Running genes' selection: DONE
Elaborating PCA - START
Elaborating PCA - DONE
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
Creating new clusterization: DONE
Creating PDF UMAP in file: /tmp/RtmpmhC9id/test/reclustering/pdf_umap_1.pdf
UMAP plot
Calculating UMAP: START
22:57:22 UMAP embedding parameters a = 0.9922 b = 1.112
22:57:22 Read 1200 rows and found 40 numeric columns
22:57:22 Using Annoy for neighbor search, n_neighbors = 30
22:57:22 Building Annoy index with metric = cosine, n_trees = 50
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
22:57:23 Writing NN index file to temp file /tmp/RtmpmhC9id/file31be685bc06080
22:57:23 Searching Annoy index using 1 thread, search_k = 3000
22:57:23 Annoy recall = 100%
22:57:24 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
22:57:27 Initializing from normalized Laplacian + noise (using RSpectra)
22:57:27 Commencing optimization for 500 epochs, with 42674 positive edges
22:57:27 Using rng type: pcg
Using method 'umap'
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
22:57:29 Optimization finished
Calculating UMAP: DONE
UMAP plot
Calculating UMAP: START
22:57:30 UMAP embedding parameters a = 0.9922 b = 1.112
22:57:30 Read 1200 rows and found 40 numeric columns
22:57:30 Using Annoy for neighbor search, n_neighbors = 30
22:57:30 Building Annoy index with metric = cosine, n_trees = 50
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
22:57:30 Writing NN index file to temp file /tmp/RtmpmhC9id/file31be681d303ea1
22:57:30 Searching Annoy index using 1 thread, search_k = 3000
22:57:30 Annoy recall = 100%
22:57:31 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
22:57:33 Initializing from normalized Laplacian + noise (using RSpectra)
22:57:33 Commencing optimization for 500 epochs, with 42674 positive edges
22:57:33 Using rng type: pcg
Using method 'umap'
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
22:57:36 Optimization finished
Calculating UMAP: DONE
* checking uniformity of cluster '0' of 4 clusters
Asked to drop 0 genes and 847 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [353] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.038818359375 | max: 11.109375 | % negative: 5
Only analysis elapsed time: 3.49287295341492
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 17.090616941452
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.0108528137207031
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0529975891113281
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.29897928237915
Total calculations elapsed time: 18.4534466266632
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 19.7869777679443
Total elapsed time: 25.403742313385
Checking uniformity for the cluster '01_0000' with 353 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0000, with size 353, is uniform
Cluster 01_0000's shift to uniformity: -0.0160990564073067
class = SimpleGDIUniformityCheck, check.isCheckAbove = FALSE, check.GDIThreshold = 1.46, check.maxRatioBeyond = 0.01, check.maxRankBeyond = 0, check.fractionBeyond = 0.00833333333333333, check.thresholdRank = 0, check.quantileAtRatio = 1.44390094359269, check.quantileAtRank = NaN, isUniform = TRUE, clusterSize = 353
cluster 01_0000 is uniform
* checking uniformity of cluster '1' of 4 clusters
Asked to drop 0 genes and 879 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [321] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.0501708984375 | max: 14.515625 | % negative: 8.5
Only analysis elapsed time: 4.21849989891052
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 16.1392428874969
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.0081791877746582
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0490689277648926
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.17887210845947
Total calculations elapsed time: 17.375363111496
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 18.7105095386505
Total elapsed time: 25.8504505157471
Checking uniformity for the cluster '01_0001' with 321 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0001, with size 321, is uniform
Cluster 01_0001's shift to uniformity: -0.029501740796509
class = SimpleGDIUniformityCheck, check.isCheckAbove = FALSE, check.GDIThreshold = 1.46, check.maxRatioBeyond = 0.01, check.maxRankBeyond = 0, check.fractionBeyond = 0.005, check.thresholdRank = 0, check.quantileAtRatio = 1.43049825920349, check.quantileAtRank = NaN, isUniform = TRUE, clusterSize = 321
cluster 01_0001 is uniform
* checking uniformity of cluster '2' of 4 clusters
Asked to drop 0 genes and 905 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [295] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.054443359375 | max: 104 | % negative: 36.5
Only analysis elapsed time: 4.18443918228149
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 15.5011703968048
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00893115997314453
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0241382122039795
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.42386603355408
Total calculations elapsed time: 16.958105802536
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 18.7126989364624
Total elapsed time: 25.3948168754578
Checking uniformity for the cluster '01_0002' with 295 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0002, with size 295, is uniform
Cluster 01_0002's shift to uniformity: -0.0675248136510636
class = SimpleGDIUniformityCheck, check.isCheckAbove = FALSE, check.GDIThreshold = 1.46, check.maxRatioBeyond = 0.01, check.maxRankBeyond = 0, check.fractionBeyond = 0, check.thresholdRank = 0, check.quantileAtRatio = 1.39247518634894, check.quantileAtRank = NaN, isUniform = TRUE, clusterSize = 295
cluster 01_0002 is uniform
* checking uniformity of cluster '3' of 4 clusters
Asked to drop 0 genes and 969 cells
Cotan analysis functions started
Asked to drop 1 genes and 0 cells
Genes/cells selection done: dropped [1] genes and [0] cells
Working on [599] genes and [231] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.0662841796875 | max: 82.5 | % negative: 33.889816360601
Only analysis elapsed time: 3.58606839179993
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 16.3066072463989
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00890374183654785
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0246157646179199
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.33383011817932
Total calculations elapsed time: 17.6739568710327
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 19.0272307395935
Total elapsed time: 25.1301443576813
Checking uniformity for the cluster '01_0003' with 231 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0003, with size 231, is uniform
Cluster 01_0003's shift to uniformity: -0.0869783653733056
class = SimpleGDIUniformityCheck, check.isCheckAbove = FALSE, check.GDIThreshold = 1.46, check.maxRatioBeyond = 0.01, check.maxRankBeyond = 0, check.fractionBeyond = 0, check.thresholdRank = 0, check.quantileAtRatio = 1.37302163462669, check.quantileAtRank = NaN, isUniform = TRUE, clusterSize = 231
cluster 01_0003 is uniform
Found 4 uniform and 0 non-uniform clusters
Stopping: too few cells left to cluster
Unclustered cell left: 0
The final raw clusterization contains [ 4 ] different clusters: 01_0000, 01_0001, 01_0002, 01_0003
Differential Expression Analysis - START
* DEA on cluster '1' with 353 cells
* DEA on cluster '2' with 321 cells
* DEA on cluster '3' with 295 cells
* DEA on cluster '4' with 231 cells
Differential Expression Analysis - DONE
Applied reordering to clusterization is:
1 -> 3, 2 -> 1, 3 -> 2, 4 -> 4
Creating cells' uniform clustering: DONE
Applied reordering to clusterization is:
1 -> 1, 2 -> 2, 3 -> 3, 4 -> 4
Asked to drop 0 genes and 905 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [295] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.054443359375 | max: 104 | % negative: 36.5
Only analysis elapsed time: 3.92374062538147
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 15.1475224494934
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.0065457820892334
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0185856819152832
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.16374492645264
Total calculations elapsed time: 16.3363988399506
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 17.6223003864288
Total elapsed time: 24.499748468399
Checking uniformity for the cluster 'Cluster_2' with 295 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster Cluster_2, with size 295, is uniform
Cluster Cluster_2's shift to uniformity: -0.00560790972802017
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.297, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.0416666666666667, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.29139209027198, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.307, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.035, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.34700400911787, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.4, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.00833333333333333, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.39247518634894, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 5, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.4279762218014, isUniform = TRUE, clusterSize = 295
Differential Expression Analysis - START
* DEA on cluster '-1' with 50 cells
* DEA on cluster '1' with 321 cells
* DEA on cluster '2' with 245 cells
* DEA on cluster '3' with 353 cells
* DEA on cluster '4' with 231 cells
Differential Expression Analysis - DONE
Applied reordering to clusterization is:
1 -> 2, 2 -> 1, 3 -> 4, 4 -> 3, -1 -> -1
Applied reordering to clusterization is:
1 -> 2, 2 -> 3, 3 -> 1, 4 -> 4, -1 -> -1
Creating cells' uniform clustering: START
In iteration 1 the number of cells to re-cluster is 1200 cells belonging to 0 clusters
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
Calculate `GDI`: DONE
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
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
Creating new clusterization: DONE
Using passed in clusterization
Creating PDF UMAP in file: /tmp/RtmpmhC9id/test/reclustering/pdf_umap_1.pdf
UMAP plot
Calculating UMAP: START
23:00:41 UMAP embedding parameters a = 0.9922 b = 1.112
23:00:41 Read 1200 rows and found 40 numeric columns
23:00:41 Using Annoy for neighbor search, n_neighbors = 30
23:00:41 Building Annoy index with metric = cosine, n_trees = 50
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
23:00:41 Writing NN index file to temp file /tmp/RtmpmhC9id/file31be683d345a47
23:00:41 Searching Annoy index using 1 thread, search_k = 3000
23:00:42 Annoy recall = 100%
23:00:43 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
23:00:45 Initializing from normalized Laplacian + noise (using RSpectra)
23:00:45 Commencing optimization for 500 epochs, with 42644 positive edges
23:00:45 Using rng type: pcg
Using method 'umap'
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
23:00:48 Optimization finished
Calculating UMAP: DONE
UMAP plot
Calculating UMAP: START
23:00:48 UMAP embedding parameters a = 0.9922 b = 1.112
23:00:48 Read 1200 rows and found 40 numeric columns
23:00:48 Using Annoy for neighbor search, n_neighbors = 30
23:00:48 Building Annoy index with metric = cosine, n_trees = 50
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
23:00:49 Writing NN index file to temp file /tmp/RtmpmhC9id/file31be681e806e2b
23:00:49 Searching Annoy index using 1 thread, search_k = 3000
23:00:49 Annoy recall = 100%
23:00:50 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
23:00:53 Initializing from normalized Laplacian + noise (using RSpectra)
23:00:53 Commencing optimization for 500 epochs, with 42644 positive edges
23:00:53 Using rng type: pcg
Using method 'umap'
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
23:00:56 Optimization finished
Calculating UMAP: DONE
* checking uniformity of cluster '1' of 2 clusters
Asked to drop 0 genes and 600 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [600] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.0386962890625 | max: 19.40625 | % negative: 6.5
Only analysis elapsed time: 3.64788103103638
Cotan genes' COEX estimation started
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 15.3565781116486
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00681972503662109
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.047415018081665
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.05965304374695
Total calculations elapsed time: 16.4704658985138
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 17.7009620666504
Total elapsed time: 23.5562527179718
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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0001, with size 600, is uniform
Cluster 01_0001's shift to uniformity: -0.012256083873327
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.38474391612667, 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.42418110276821, 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.45338035249414, 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.51521594648007, 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 analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [600] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.032958984375 | max: 10.0859375 | % negative: 3.66666666666667
Only analysis elapsed time: 3.94116163253784
Cotan genes' COEX estimation started
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 15.2431845664978
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00767803192138672
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0571169853210449
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.1251335144043
Total calculations elapsed time: 16.4331130981445
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 17.6050734519958
Total elapsed time: 24.0618622303009
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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0002, with size 600, is uniform
Cluster 01_0002's shift to uniformity: -0.0203445505533575
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.37665544944664, 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.42305229982216, 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.43783135828069, 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.50182972335853, 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
Differential Expression Analysis - DONE
Applied reordering to clusterization is:
1 -> 1, 2 -> 2
Creating cells' uniform clustering: DONE
findClustersMarkers - START
Log Fold Change Analysis - START
* Analysis of cluster: '1'
* Analysis of cluster: '2'
* Analysis of cluster: '3'
* Analysis of cluster: '4'
Log Fold Change Analysis - DONE
findClustersMarkers - DONE
findClustersMarkers - START
Log Fold Change Analysis - START
* Analysis of cluster: '1'
* Analysis of cluster: '2'
* Analysis of cluster: '3'
* Analysis of cluster: '4'
Log Fold Change Analysis - DONE
findClustersMarkers - DONE
findClustersMarkers - START
Differential Expression Analysis - START
* DEA on cluster '1' with 321 cells
* DEA on cluster '2' with 295 cells
* DEA on cluster '3' with 353 cells
* DEA on cluster '4' with 231 cells
Differential Expression Analysis - DONE
Log Fold Change Analysis - START
* Analysis of cluster: '1'
* Analysis of cluster: '2'
* Analysis of cluster: '3'
* Analysis of cluster: '4'
Log Fold Change Analysis - DONE
findClustersMarkers - DONE
[1] "Tested cluster: 1"
Asked to drop 0 genes and 879 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [321] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.0501708984375 | max: 14.515625 | % negative: 8.5
Only analysis elapsed time: 4.20996904373169
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 14.0603575706482
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00826644897460938
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0449700355529785
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.16757702827454
Total calculations elapsed time: 15.2811710834503
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 16.4935922622681
Total elapsed time: 23.4543120861053
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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
[1] "Tested cluster: 4"
Asked to drop 0 genes and 969 cells
Cotan analysis functions started
Asked to drop 1 genes and 0 cells
Genes/cells selection done: dropped [1] genes and [0] cells
Working on [599] genes and [231] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.0662841796875 | max: 82.5 | % negative: 33.889816360601
Only analysis elapsed time: 3.64061498641968
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 15.3296434879303
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00747776031494141
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0219783782958984
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.00737977027893
Total calculations elapsed time: 16.3664793968201
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 17.4122223854065
Total elapsed time: 23.3660519123077
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
Calculate `GDI`: DONE
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
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Initializing `COTAN` meta-data
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Cotan genes' COEX estimation not requested
Total elapsed time: 9.46793818473816
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-000451, g-000339, g-000295, g-000328, g-000544, g-000061, g-000227, g-000391, g-000556, g-000237
g-000067, g-000165, g-000449, g-000591, g-000087, g-000129, g-000197, g-000203, g-000487, g-000505
g-000333, g-000029, g-000271, g-000064, g-000583, g-000156, g-000448, g-000153, g-000526, g-000393
Negative: g-000518, g-000108, g-000186, g-000170, g-000401, g-000337, g-000047, g-000599, g-000432, g-000578
g-000042, g-000065, g-000493, g-000261, g-000533, g-000256, g-000560, g-000596, g-000368, g-000381
g-000535, g-000338, g-000215, g-000159, g-000365, g-000234, g-000173, g-000387, g-000225, g-000272
Computing nearest neighbor graph
Computing SNN
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 1200
Number of edges: 55884
Running Louvain algorithm with multilevel refinement...
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Maximum modularity in 10 random starts: 0.5927
Number of communities: 4
Elapsed time: 0 seconds
Clusterizations found: RNA_snn_res.0.8, seurat_clusters
Conditions found: orig.ident
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Cotan genes' COEX estimation not requested
Total elapsed time: 7.44602608680725
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 4.18468546867371
Cotan genes' COEX estimation started
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 15.5476019382477
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00792288780212402
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0700626373291016
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.998740196228027
Total calculations elapsed time: 16.6243276596069
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 17.8419427871704
Total elapsed time: 24.2710494995117
Calculating gene co-expression space - START
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and genome wide on rows
calculating PValues: DONE
Number of selected secondary markers: 6
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Calculate `GDI`: DONE
Number of columns (V set - secondary markers): 6
Number of rows (U set): 60
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
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
Calculate `GDI`: DONE
Number of columns (V set - secondary markers): 109
Number of rows (U set): 60
Calculating gene co-expression space - DONE
Establishing gene clusters - START
Calculating gene co-expression space - START
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and genome wide on rows
calculating PValues: DONE
Number of selected secondary markers: 109
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Calculate `GDI`: DONE
Number of columns (V set - secondary markers): 109
Number of rows (U set): 60
Calculating gene co-expression space - DONE
Establishing gene clusters - DONE
Initializing `COTAN` meta-data
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Cotan genes' COEX estimation not requested
Total elapsed time: 7.10424709320068
Differential Expression Analysis - START
* DEA on cluster '1' with 353 cells
* DEA on cluster '2' with 221 cells
* DEA on cluster '3' with 315 cells
* DEA on cluster '4' with 311 cells
Differential Expression Analysis - DONE
Log Fold Change Analysis - START
* Analysis of cluster: '1'
* Analysis of cluster: '2'
* Analysis of cluster: '3'
* Analysis of cluster: '4'
Log Fold Change Analysis - DONE
clustersDeltaExpression - START
Handling cluster '1' with mean UDE 1.43530796540674
Handling cluster '2' with mean UDE 0.640931107489519
Handling cluster '3' with mean UDE 0.546546914955575
Handling cluster '4' with mean UDE 1.22034802329657
clustersDeltaExpression - DONE
In group G1 there are 3 detected over 3 genes
In group G2 there are 2 detected over 2 genes
In group G3 there are 5 detected over 5 genes
Merging cells' uniform clustering: START
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 353 cells
* DEA on cluster '2' with 221 cells
* DEA on cluster '3' with 315 cells
* DEA on cluster '4' with 311 cells
Differential Expression Analysis - DONE
Updating check results for the 0 already tested pairs to align to new checker
2 new clusters pairs to be tested for merging:
c("1", "2") c("3", "4")
*1_2-merge
Asked to drop 0 genes and 626 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [574] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.0340576171875 | max: 10.4375 | % negative: 4.33333333333333
Only analysis elapsed time: 4.072509765625
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 15.0558400154114
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00873231887817383
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0580601692199707
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.22391963005066
Total calculations elapsed time: 16.3465521335602
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 17.5184328556061
Total elapsed time: 23.9833772182465
Checking uniformity for the cluster '1_2-merge' with 574 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 1_2-merge, with size 574, is not uniform
Cluster 1_2-merge's shift to uniformity: 0.0820221369127518
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.297, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.326666666666667, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.37902213691275, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.307, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.28, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.41724738266283, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.4, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.0333333333333333, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.42971854140312, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 20, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.5067755752749, isUniform = FALSE, clusterSize = 574
Clusters 1 and 2 cannot be merged
*3_4-merge
Asked to drop 0 genes and 574 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [626] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.03839111328125 | max: 17.40625 | % negative: 6.66666666666667
Only analysis elapsed time: 3.82468748092651
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 17.1670649051666
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00606393814086914
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0430517196655273
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.05011606216431
Total calculations elapsed time: 18.2662966251373
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 19.3683197498322
Total elapsed time: 25.5337083339691
Checking uniformity for the cluster '3_4-merge' with 626 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 3_4-merge, with size 626, is not uniform
Cluster 3_4-merge's shift to uniformity: 0.0914611744346074
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.297, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.248333333333333, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.38846117443461, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.307, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.218333333333333, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.42066795244751, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.4, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.0366666666666667, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.45442688129939, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 22, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.52592182249635, isUniform = FALSE, clusterSize = 626
Clusters 3 and 4 cannot be merged
Finished the first batch - no merges were executed
Start merging nearest clusters: iteration 2
Differential Expression Analysis - START
* DEA on cluster '1' with 353 cells
* DEA on cluster '2' with 221 cells
* DEA on cluster '3' with 315 cells
* DEA on cluster '4' with 311 cells
Differential Expression Analysis - DONE
Updating check results for the 2 already tested pairs to align to new checker
2 new clusters pairs to be tested for merging:
c("2", "3") c("2", "4")
*2_3-merge
Asked to drop 0 genes and 664 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [536] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.0390625 | max: 12.59375 | % negative: 9.16666666666667
Only analysis elapsed time: 3.58826160430908
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 13.7845740318298
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00853729248046875
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0317196846008301
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.11960101127625
Total calculations elapsed time: 14.9444320201874
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 16.0143435001373
Total elapsed time: 21.7427451610565
Checking uniformity for the cluster '2_3-merge' with 536 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 2_3-merge, with size 536, is not uniform
Cluster 2_3-merge's shift to uniformity: 0.102213698003777
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.297, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.29, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.39921369800378, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.307, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.251666666666667, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.45276697088286, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.4, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.05, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.49016449535422, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 30, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.54477108741147, isUniform = FALSE, clusterSize = 536
Clusters 2 and 3 cannot be merged
*2_4-merge
Asked to drop 0 genes and 668 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [532] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.04754638671875 | max: 10.234375 | % negative: 5
Only analysis elapsed time: 3.27253389358521
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 10.0250902175903
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00606799125671387
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0397226810455322
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.27024388313293
Total calculations elapsed time: 11.3411247730255
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 12.6845469474792
Total elapsed time: 18.0383830070496
Checking uniformity for the cluster '2_4-merge' with 532 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 2_4-merge, with size 532, is not uniform
Cluster 2_4-merge's shift to uniformity: 0.213129217682233
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.297, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.39, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.51012921768223, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.307, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.361666666666667, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.55599610798341, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.4, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.181666666666667, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.59088431608469, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 109, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.68475283375668, isUniform = FALSE, clusterSize = 532
Clusters 2 and 4 cannot be merged
No clusters will be merged
None of the remaining tested cluster pairs could be merged
Executed all merges for threshold 1.297 out of 4 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 353 cells
* DEA on cluster '2' with 221 cells
* DEA on cluster '3' with 315 cells
* DEA on cluster '4' with 311 cells
Differential Expression Analysis - DONE
Updating check results for the 4 already tested pairs to align to new checker
2 new clusters pairs to be tested for merging:
c("1", "3") c("1", "4")
*1_3-merge
Asked to drop 0 genes and 532 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [668] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.046630859375 | max: 9.03515625 | % negative: 1.33333333333333
Only analysis elapsed time: 3.54457831382751
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 9.81931734085083
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00573992729187012
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0436253547668457
Estimating genes' COEX
Calculate genes' COEX elapsed time: 0.936333656311035
Total calculations elapsed time: 10.8050162792206
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 11.8075966835022
Total elapsed time: 17.7869513034821
Checking uniformity for the cluster '1_3-merge' with 668 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 1_3-merge, with size 668, is not uniform
Cluster 1_3-merge's shift to uniformity: 0.177923707405303
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.397, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.33, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.5749237074053, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.407, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.31, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.6368755707341, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.5, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.128333333333333, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.65704996173771, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 77, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.68881338289977, isUniform = FALSE, clusterSize = 668
Clusters 1 and 3 cannot be merged
*1_4-merge
Asked to drop 0 genes and 536 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [664] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.0164794921875 | max: 5.5595703125 | % negative: 0.166666666666667
Only analysis elapsed time: 3.32131290435791
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 9.67338252067566
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00630593299865723
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0465245246887207
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.08942317962646
Total calculations elapsed time: 10.8156361579895
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 11.9177119731903
Total elapsed time: 17.2770166397095
Checking uniformity for the cluster '1_4-merge' with 664 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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 1_4-merge, with size 664, is not uniform
Cluster 1_4-merge's shift to uniformity: 0.138277355686157
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.397, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.348333333333333, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.53527735568616, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.407, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.316666666666667, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.55951216335001, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.5, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.0983333333333333, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.59105696847609, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 59, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.64605433203111, isUniform = FALSE, clusterSize = 664
Clusters 1 and 4 cannot be merged
Clusters 2 and 1 will be merged
Clusters 3 and 4 will be merged
Executed 2 merges
Start merging nearest clusters: iteration 4
Differential Expression Analysis - START
* DEA on cluster '1_2-merge' with 574 cells
* DEA on cluster '3_4-merge' with 626 cells
Differential Expression Analysis - DONE
Updating check results for the 6 already tested pairs to align to new checker
1 new clusters pairs to be tested for merging:
c("1_2-merge", "3_4-merge")
*1_2-merge_3_4-merge-merge
Asked to drop no genes or cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 4.30208849906921
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 10.3464066982269
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00925612449645996
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0530517101287842
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.13837623596191
Total calculations elapsed time: 11.5470907688141
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 12.6159687042236
Total elapsed time: 19.441073179245
Checking uniformity for the cluster '1_2-merge_3_4-merge-merge' with 1200 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
Calculate `GDI`: START
Effective number of cores used: 1
Executing 1 genes batches
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 1_2-merge_3_4-merge-merge, with size 1200, is not uniform
Cluster 1_2-merge_3_4-merge-merge's shift to uniformity: 0.173174985379327
class = AdvancedGDIUniformityCheck, firstCheck.isCheckAbove = FALSE, firstCheck.GDIThreshold = 1.397, firstCheck.maxRatioBeyond = 0.05, firstCheck.maxRankBeyond = 0, firstCheck.fractionBeyond = 0.415, firstCheck.thresholdRank = 0, firstCheck.quantileAtRatio = 1.57017498537933, firstCheck.quantileAtRank = NaN, secondCheck.isCheckAbove = TRUE, secondCheck.GDIThreshold = 1.407, secondCheck.maxRatioBeyond = 0.02, secondCheck.maxRankBeyond = 0, secondCheck.fractionBeyond = 0.388333333333333, secondCheck.thresholdRank = 0, secondCheck.quantileAtRatio = 1.61240081401504, secondCheck.quantileAtRank = NaN, thirdCheck.isCheckAbove = FALSE, thirdCheck.GDIThreshold = 1.5, thirdCheck.maxRatioBeyond = 0.01, thirdCheck.maxRankBeyond = 0, thirdCheck.fractionBeyond = 0.15, thirdCheck.thresholdRank = 0, thirdCheck.quantileAtRatio = 1.64046660675367, thirdCheck.quantileAtRank = NaN, fourthCheck.isCheckAbove = FALSE, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 90, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.66512085964758, isUniform = FALSE, clusterSize = 1200
Clusters 1_2-merge and 3_4-merge cannot be merged
No clusters will be merged
None of the remaining tested cluster pairs could be merged
Executed all merges for threshold 1.397 out of 7 checks
The final merged clusterization contains [2] different clusters: 1_2-merge, 3_4-merge
Differential Expression Analysis - START
* DEA on cluster '1' with 574 cells
* DEA on cluster '2' with 626 cells
Differential Expression Analysis - DONE
Applied reordering to clusterization is:
1 -> 1, 2 -> 2
Merging cells' uniform clustering: DONE
Applied reordering to clusterization is:
1 -> 1, 2 -> 2
Log Fold Change Analysis - START
* Analysis of cluster: '1'
* Analysis of cluster: '2'
Log Fold Change Analysis - DONE
Asked to drop 0 genes and 626 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [574] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.0340576171875 | max: 10.4375 | % negative: 4.33333333333333
Only analysis elapsed time: 3.84807872772217
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 10.6958250999451
Calculating genes' COEX normalization factor
Calculate genes' normalization factor elapsed time: 0.00773119926452637
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.043071985244751
Estimating genes' COEX
Calculate genes' COEX elapsed time: 1.23074555397034
Total calculations elapsed time: 11.9773738384247
Calculate genes' COEX (legacy): DONE
Only genes' COEX elapsed time: 13.4372475147247
Total elapsed time: 19.6171412467957
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
Calculate `GDI`: DONE
Calculate GDI dataframe: DONE
Asked to drop 0 genes and 574 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [626] cells
Estimate `dispersion`: START
Effective number of cores used: 1
Executing 1 genes batches
`dispersion` | min: -0.03839111328125 | max: 17.40625 | % negative: 6.66666666666667
Only analysis elapsed time: 4.52861595153809
Cotan genes' COEX estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' COEX (legacy): START
Retrieving expected genes' contingency table
calculating NN..
COTAN.Rcheck/COTAN-Ex.timings
| name | user | system | elapsed | |
| COTAN_Legacy | 0.002 | 0.000 | 0.009 | |
| COTAN_ObjectCreation | 354.796 | 0.870 | 23.444 | |
| CalculatingCOEX | 1490.455 | 2.851 | 116.696 | |
| ClustersList | 0.009 | 0.000 | 0.009 | |
| Conversions | 6.521 | 0.137 | 6.899 | |
| GenesStatistics | 364.379 | 0.983 | 26.849 | |
| HandleMetaData | 0.080 | 0.020 | 0.111 | |
| HandlingClusterizations | 387.745 | 2.204 | 61.365 | |
| HandlingConditions | 0.146 | 0.001 | 0.147 | |
| HeatmapPlots | 1078.592 | 2.341 | 76.217 | |
| LoggingFunctions | 0.001 | 0.001 | 0.002 | |
| ParametersEstimations | 31.564 | 1.013 | 32.105 | |
| RawDataCleaning | 7.509 | 0.021 | 7.690 | |
| RawDataGetters | 0.103 | 0.005 | 0.109 | |
| UniformClusters | 4041.270 | 17.554 | 410.295 | |
| getColorsVector | 0.001 | 0.000 | 0.001 | |