Galfrè Silvia Giulia
| Snapshot Date: 2025-01-23 13:40 -0500 (Thu, 23 Jan 2025) |
| git_url: https://git.bioconductor.org/packages/COTAN |
| git_branch: devel |
| git_last_commit: ff57026 |
| git_last_commit_date: 2024-12-09 09:29:06 -0500 (Mon, 09 Dec 2024) |
| Back to Multiple platform build/check report for BioC 3.21: simplified long |
|
This page was generated on 2025-01-24 11:41 -0500 (Fri, 24 Jan 2025).
| Hostname | OS | Arch (*) | R version | Installed pkgs |
|---|---|---|---|---|
| nebbiolo1 | Linux (Ubuntu 24.04.1 LTS) | x86_64 | R Under development (unstable) (2025-01-20 r87609) -- "Unsuffered Consequences" | 4609 |
| palomino7 | Windows Server 2022 Datacenter | x64 | R Under development (unstable) (2025-01-21 r87610 ucrt) -- "Unsuffered Consequences" | 4393 |
| lconway | macOS 12.7.1 Monterey | x86_64 | R Under development (unstable) (2025-01-20 r87609) -- "Unsuffered Consequences" | 3839 |
| kjohnson3 | macOS 13.7.1 Ventura | arm64 | R Under development (unstable) (2025-01-20 r87609) -- "Unsuffered Consequences" | 3835 |
| kunpeng2 | Linux (openEuler 22.03 LTS-SP1) | aarch64 | R Under development (unstable) (2024-11-24 r87369) -- "Unsuffered Consequences" | 4408 |
| 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 440/2286 | Hostname | OS / Arch | INSTALL | BUILD | CHECK | BUILD BIN | ||||||||
| COTAN 2.7.2 (landing page) Galfrè Silvia Giulia
| nebbiolo1 | Linux (Ubuntu 24.04.1 LTS) / x86_64 | OK | OK | OK | | ||||||||
| palomino7 | Windows Server 2022 Datacenter / x64 | OK | OK | OK | OK | | ||||||||
| lconway | macOS 12.7.1 Monterey / x86_64 | OK | OK | OK | OK | | ||||||||
| kjohnson3 | macOS 13.7.1 Ventura / arm64 | OK | OK | OK | OK | | ||||||||
| kunpeng2 | Linux (openEuler 22.03 LTS-SP1) / 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.7.2 |
| Command: /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:COTAN.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings COTAN_2.7.2.tar.gz |
| StartedAt: 2025-01-23 18:23:14 -0500 (Thu, 23 Jan 2025) |
| EndedAt: 2025-01-23 18:48:56 -0500 (Thu, 23 Jan 2025) |
| EllapsedTime: 1542.0 seconds |
| RetCode: 0 |
| Status: OK |
| CheckDir: COTAN.Rcheck |
| Warnings: 0 |
##############################################################################
##############################################################################
###
### Running command:
###
### /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:COTAN.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings COTAN_2.7.2.tar.gz
###
##############################################################################
##############################################################################
* using log directory ‘/Users/biocbuild/bbs-3.21-bioc/meat/COTAN.Rcheck’
* using R Under development (unstable) (2025-01-20 r87609)
* using platform: x86_64-apple-darwin20
* R was compiled by
Apple clang version 14.0.0 (clang-1400.0.29.202)
GNU Fortran (GCC) 14.2.0
* running under: macOS Monterey 12.7.6
* using session charset: UTF-8
* using option ‘--no-vignettes’
* checking for file ‘COTAN/DESCRIPTION’ ... OK
* checking extension type ... Package
* this is package ‘COTAN’ version ‘2.7.2’
* package encoding: UTF-8
* checking package namespace information ... OK
* checking package dependencies ... INFO
Imports includes 28 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 dependencies in R code ... NOTE
'library' or 'require' call to ‘torch’ in package code.
Please use :: or requireNamespace() instead.
See section 'Suggested packages' in the 'Writing R Extensions' manual.
Unexported object imported by a ':::' call: ‘ggplot2:::ggname’
See the note in ?`:::` about the use of this operator.
* checking S3 generic/method consistency ... OK
* checking replacement functions ... OK
* checking foreign function calls ... OK
* checking R code for possible problems ... NOTE
mergeUniformCellsClusters : fromMergedName: warning in
vapply(currentClNames, function(clName, mergedName) {: partial
argument match of 'FUN.VAL' to 'FUN.VALUE'
mergeUniformCellsClusters : fromMergedName: warning in
return(str_detect(mergedName, clName)): partial argument match of
'FUN.VAL' to 'FUN.VALUE'
mergeUniformCellsClusters : fromMergedName: warning in }, FUN.VAL =
logical(1L), mergedClName): partial argument match of 'FUN.VAL' to
'FUN.VALUE'
ECDPlot: no visible binding for global variable ‘.’
GDIPlot: no visible binding for global variable ‘sum.raw.norm’
GDIPlot: no visible binding for global variable ‘GDI’
UMAPPlot: no visible binding for global variable ‘x’
UMAPPlot: no visible binding for global variable ‘y’
UMAPPlot: no visible binding for global variable ‘types’
calculateCoex_Legacy: no visible binding for global variable
‘expectedNN’
calculateCoex_Legacy: no visible binding for global variable
‘expectedNY’
calculateCoex_Legacy: no visible binding for global variable
‘expectedYN’
calculateCoex_Legacy: no visible binding for global variable
‘expectedYY’
calculateCoex_Legacy: no visible binding for global variable
‘observedYY’
calculateCoex_Legacy: no visible binding for global variable ‘.’
calculateG: no visible binding for global variable ‘observedNN’
calculateG: no visible binding for global variable ‘observedNY’
calculateG: no visible binding for global variable ‘observedYN’
calculateG: no visible binding for global variable ‘observedYY’
calculateG: no visible binding for global variable ‘expectedNN’
calculateG: no visible binding for global variable ‘expectedNY’
calculateG: no visible binding for global variable ‘expectedYN’
calculateG: no visible binding for global variable ‘expectedYY’
calculatePartialCoex: no visible binding for global variable
‘expectedNN’
calculatePartialCoex: no visible binding for global variable
‘expectedNY’
calculatePartialCoex: no visible binding for global variable
‘expectedYN’
calculatePartialCoex: no visible binding for global variable
‘expectedYY’
calculatePartialCoex: no visible binding for global variable
‘observedYY’
calculatePartialCoex: no visible binding for global variable ‘.’
cellSizePlot: no visible binding for global variable ‘.’
cellsUniformClustering: no visible binding for global variable
‘objSeurat’
cellsUniformClustering: no visible binding for global variable
‘usedMaxResolution’
cellsUniformClustering: no visible binding for global variable
‘permMap’
checkClusterUniformity: no visible binding for global variable
‘pcaCellsPlot’
checkClusterUniformity: no visible binding for global variable ‘.’
checkClusterUniformity: no visible binding for global variable
‘bGroupGenesPlot’
checkClusterUniformity: no visible binding for global variable
‘UDEPLot’
checkClusterUniformity: no visible binding for global variable ‘nuPlot’
checkClusterUniformity: no visible binding for global variable
‘zoomedNuPlot’
cleanPlots: no visible binding for global variable ‘PC1’
cleanPlots: no visible binding for global variable ‘PC2’
cleanPlots: no visible binding for global variable ‘n’
cleanPlots: no visible binding for global variable ‘means’
cleanPlots: no visible binding for global variable ‘nu’
clustersMarkersHeatmapPlot: no visible binding for global variable
‘condName’
clustersMarkersHeatmapPlot: no visible binding for global variable
‘conditions’
clustersSummaryPlot: no visible binding for global variable ‘keys’
clustersSummaryPlot: no visible binding for global variable ‘values’
clustersSummaryPlot: no visible binding for global variable
‘CellNumber’
clustersSummaryPlot: no visible binding for global variable ‘ExpGenes’
clustersSummaryPlot: no visible binding for global variable ‘Cluster’
clustersSummaryPlot: no visible binding for global variable ‘Condition’
establishGenesClusters: no visible binding for global variable
‘secondaryMarkers’
establishGenesClusters: no visible binding for global variable ‘GCS’
establishGenesClusters: no visible binding for global variable
‘rankGenes’
expectedContingencyTables: no visible binding for global variable
‘expectedN’
expectedPartialContingencyTables: no visible binding for global
variable ‘expectedNN’
expectedPartialContingencyTables: no visible binding for global
variable ‘expectedN’
genesSizePlot: no visible binding for global variable ‘.’
geom_flat_violin : <anonymous>: no visible binding for global variable
‘group’
geom_flat_violin : <anonymous>: no visible binding for global variable
‘y’
geom_flat_violin : <anonymous>: no visible binding for global variable
‘x’
geom_flat_violin : <anonymous>: no visible binding for global variable
‘width’
geom_flat_violin : <anonymous>: no visible binding for global variable
‘violinwidth’
geom_flat_violin : <anonymous>: no visible binding for global variable
‘xmax’
geom_flat_violin : <anonymous>: no visible binding for global variable
‘xminv’
geom_flat_violin : <anonymous>: no visible binding for global variable
‘xmaxv’
heatmapPlot: no visible binding for global variable ‘obj’
heatmapPlot: no visible binding for global variable ‘type’
heatmapPlot: no visible binding for global variable ‘g2’
heatmapPlot: no visible binding for global variable ‘coex’
mergeUniformCellsClusters : fromMergedName: no visible binding for
global variable ‘clName1’
mergeUniformCellsClusters : fromMergedName: no visible binding for
global variable ‘clName2’
mergeUniformCellsClusters : testPairListMerge: no visible binding for
global variable ‘cl1’
mergeUniformCellsClusters : testPairListMerge: no visible binding for
global variable ‘cl2’
mergeUniformCellsClusters : mergeAllClusters : <anonymous>: no visible
binding for global variable ‘clName1’
mergeUniformCellsClusters : mergeAllClusters : <anonymous>: no visible
binding for global variable ‘clName2’
mergeUniformCellsClusters : mergeAllClusters: no visible binding for
global variable ‘clName1’
mergeUniformCellsClusters : mergeAllClusters: no visible binding for
global variable ‘clName2’
mergeUniformCellsClusters: no visible binding for global variable
‘permMap’
mitochondrialPercentagePlot: no visible binding for global variable ‘.’
mitochondrialPercentagePlot: no visible binding for global variable
‘mit.percentage’
observedContingencyTables: no visible binding for global variable
‘observedY’
observedPartialContingencyTables: no visible binding for global
variable ‘observedYY’
observedPartialContingencyTables: no visible binding for global
variable ‘observedY’
scatterPlot: no visible binding for global variable ‘.’
scatterPlot: no visible binding for global variable ‘.x’
calculateCoex,COTAN: no visible binding for global variable ‘useTorch’
coerce,COTAN-scCOTAN: no visible binding for global variable ‘rawNorm’
coerce,COTAN-scCOTAN: no visible binding for global variable ‘nu’
coerce,COTAN-scCOTAN: no visible binding for global variable ‘lambda’
coerce,COTAN-scCOTAN: no visible binding for global variable ‘a’
coerce,COTAN-scCOTAN: no visible binding for global variable ‘hk’
coerce,COTAN-scCOTAN: no visible binding for global variable ‘clusters’
coerce,COTAN-scCOTAN: no visible binding for global variable
‘clusterData’
Undefined global functions or variables:
. .x CellNumber Cluster Condition ExpGenes GCS GDI PC1 PC2 UDEPLot a
bGroupGenesPlot cl1 cl2 clName1 clName2 clusterData clusters coex
condName conditions expectedN expectedNN expectedNY expectedYN
expectedYY g2 group hk keys lambda means mit.percentage n nu nuPlot
obj objSeurat observedNN observedNY observedY observedYN observedYY
pcaCellsPlot permMap rankGenes rawNorm secondaryMarkers sum.raw.norm
type types useTorch usedMaxResolution values violinwidth width x xmax
xmaxv xminv y zoomedNuPlot
* checking Rd files ... OK
* checking Rd metadata ... OK
* checking Rd cross-references ... NOTE
Found the following Rd file(s) with Rd \link{} targets missing package
anchors:
Conversions.Rd: SingleCellExperiment-class, Assays-class,
SingleCellExperiment
HandlingClusterizations.Rd: umap.defaults
getColorsVector.Rd: brewer.pal, brewer.pal.info
Please provide package anchors for all Rd \link{} targets not in the
package itself and the base packages.
* 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 222.538 3.946 226.495
CalculatingCOEX 58.574 1.312 60.279
HeatmapPlots 39.935 1.011 41.239
HandlingClusterizations 31.896 1.480 32.975
ParametersEstimations 27.417 0.857 28.466
GenesStatistics 11.943 0.423 12.199
COTAN_ObjectCreation 10.732 0.333 11.137
Conversions 6.167 0.148 6.347
RawDataCleaning 4.988 0.329 5.113
* checking for unstated dependencies in ‘tests’ ... OK
* checking tests ...
Running ‘outputTestDatasetCreation.R’
Running ‘spelling.R’
Running ‘testthat.R’
OK
* checking for unstated dependencies in vignettes ... OK
* checking package vignettes ... OK
* checking running R code from vignettes ... SKIPPED
* checking re-building of vignette outputs ... SKIPPED
* checking PDF version of manual ... OK
* DONE
Status: 3 NOTEs
See
‘/Users/biocbuild/bbs-3.21-bioc/meat/COTAN.Rcheck/00check.log’
for details.
COTAN.Rcheck/00install.out
############################################################################## ############################################################################## ### ### Running command: ### ### /Library/Frameworks/R.framework/Resources/bin/R CMD INSTALL COTAN ### ############################################################################## ############################################################################## * installing to library ‘/Library/Frameworks/R.framework/Versions/4.5-x86_64/Resources/library’ * installing *source* package ‘COTAN’ ... ** this is package ‘COTAN’ version ‘2.7.2’ ** 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 Under development (unstable) (2025-01-20 r87609) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin20
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
>
> # Creates the files to be reloaded by the tests for comparisons
> 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.223 0.088 0.302
COTAN.Rcheck/tests/spelling.Rout
R Under development (unstable) (2025-01-20 r87609) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin20
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> if (requireNamespace("spelling", quietly = TRUE))
+ spelling::spell_check_test(vignettes = TRUE, error = FALSE,
+ skip_on_cran = TRUE)
NULL
>
> proc.time()
user system elapsed
0.213 0.085 0.301
COTAN.Rcheck/tests/testthat.Rout
R Under development (unstable) (2025-01-20 r87609) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin20
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> Sys.setenv(R_TESTS = "")
> library(testthat)
> library(COTAN)
> test_check("COTAN")
Setting new log level to 4
Initializing `COTAN` meta-data
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.9033203125 | max: 4.6796875 | % negative: 10
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.49200296401978
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.000841140747070312
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00788998603820801
Estimating genes' coex
Calculate genes' coex elapsed time: 0.808491945266724
Total calculations elapsed time: 3.30922603607178
Calculate genes' coex (legacy): DONE
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 4 genes batches from [1:2] to [7:8]
Executing 1 genes batches from [9:10] to [9:10]
Estimate dispersion: DONE
dispersion | min: 0.9033203125 | max: 4.6796875 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 4 cells batches from [1:3] to [10:11]
Executing 3 cells batches from [12:14] to [18:20]
Estimate nu: DONE
nu change (abs) | max: 1.75595238095238 | median: 1.07174634176587 | mean: 1.07174634176587
Estimate 'dispersion'/'nu': START
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 1.0362548828125 | max: 4.60986328125 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.0265938895089288 | median: 0.0144680038331048 | mean: 0.0144680038331048
Nu mean: 1.69633192486233
Marginal errors | max: 1.95570586131367 | median 1.32068160171502 | mean: 1.33375826507259
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.058837890625 | max: 3.528076171875 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.416683423613994 | median: 0.239880630367975 | mean: 0.239880630367975
Nu mean: 0.823197206753982
Marginal errors | max: 0.836359531101206 | median 0.703684202571891 | mean: 0.645537958989614
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.32879638671875 | max: 4.0302734375 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.164237872898673 | median: 0.0955985184389135 | mean: 0.0955985184389135
Nu mean: 1.06863935445976
Marginal errors | max: 0.259872988828242 | median 0.213703042752633 | mean: 0.197386407582083
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.2294921875 | max: 3.8720703125 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.055185575120883 | median: 0.0319991762044448 | mean: 0.0319991762044448
Nu mean: 0.976813601083562
Marginal errors | max: 0.0951586919577032 | median 0.079429703709466 | mean: 0.0724140148396648
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.2637939453125 | max: 3.929443359375 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.0196211148938294 | median: 0.01138609597457 | mean: 0.01138609597457
Nu mean: 1.00823501891926
Marginal errors | max: 0.0327747321002292 | median 0.0272104747849538 | mean: 0.0248963830312038
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.25177001953125 | max: 3.90966796875 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.00670099066960717 | median: 0.00388888266671264 | mean: 0.00388888266671264
Nu mean: 0.997187891997105
Marginal errors | max: 0.0114324509186883 | median 0.0094232649770607 | mean: 0.00863113610779571
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.25592041015625 | max: 3.91650390625 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.00230093811689414 | median: 0.00132529122122246 | mean: 0.00132529122122246
Nu mean: 1.00097564689567
Marginal errors | max: 0.00387133150664631 | median 0.0031091017608853 | mean: 0.00286071175800213
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.2545166015625 | max: 3.914306640625 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.000837011646904529 | median: 0.000470837393363011 | mean: 0.000470837393363011
Nu mean: 0.999633825746458
Marginal errors | max: 0.00122501723202006 | median 0.00102126435760308 | mean: 0.000943992659051318
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.2550048828125 | max: 3.9150390625 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.000209227351122054 | median: 0.000122070312500028 | mean: 0.000122070312500028
Nu mean: 1.00008715703862
Marginal errors | max: 0.000364602956581805 | median 0.000313956936819793 | mean: 0.000282899574318485
Estimate dispersion/nu: DONE
Estimate 'dispersion'/'nu': START
Initializing `COTAN` meta-data
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.9033203125 | max: 4.6796875 | % negative: 10
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.7362699508667
Calculating genes' coex normalization factor
Fraction of genes' with very low expected contingency tables: 0.181818181818182
Calculate genes' normalization factor elapsed time: 0.00118517875671387
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00497078895568848
Estimating genes' coex
Calculate genes' coex elapsed time: 0.949783086776733
Total calculations elapsed time: 3.69220900535583
Calculate genes' coex (legacy): DONE
Calculate cells' coex (legacy): START
Retrieving expected cells' contingency table
calculating NN.. done
calculating YN..NY..YY..t().. done
Expected cells' contingency table elapsed time: 2.73819899559021
Calculating cells' coex normalization factor
Fraction of cells' with very low expected contingency tables: 0
Calculate cells' normalization factor elapsed time: 0.00111293792724609
Retrieving observed cells' yes/yes contingency table
calculating YY.. done
Observed cells' contingency table elapsed time: 0.00314998626708984
Estimating cells' coex
Calculate cells' coex elapsed time: 0.884376049041748
Total calculations elapsed time: 3.62683796882629
Calculate cells' coex (legacy): DONE
Initializing `COTAN` meta-data
Initializing `COTAN` meta-data
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
Estimate 'dispersion'/'nu': START
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.903564453125 | max: 4.679443359375 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 1.75719246031746 | median: 1.07229953342014 | mean: 1.07229953342014
Nu mean: 1.68489292689732
Marginal errors | max: 1.73564890252257 | median 1.37996360874076 | mean: 1.32180348113228
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.0655517578125 | max: 3.5439453125 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.402649984216273 | median: 0.231868788425666 | mean: 0.231868788425666
Nu mean: 0.829218804209393
Marginal errors | max: 0.80321315986594 | median 0.677497553540581 | mean: 0.619375430892821
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.3260498046875 | max: 4.026123046875 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.158004893526231 | median: 0.0919692884670312 | mean: 0.0919692884670312
Nu mean: 1.0660356050592
Marginal errors | max: 0.250724014302326 | median 0.206232152124436 | mean: 0.190425623677198
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.23040771484375 | max: 3.8736572265625 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.0532774732102337 | median: 0.0308837890624999 | mean: 0.0308837890624999
Nu mean: 0.977606315852266
Marginal errors | max: 0.0916983669060123 | median 0.0765266929824948 | mean: 0.0697593208689684
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.26348876953125 | max: 3.928955078125 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.0189966206463044 | median: 0.0110199320575908 | mean: 0.0110199320575908
Nu mean: 1.00797668858871
Marginal errors | max: 0.0317151207459254 | median 0.0262702142278251 | mean: 0.0240886952086962
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.2518310546875 | max: 3.9097900390625 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.00670088501353994 | median: 0.00388888101583662 | mean: 0.00388888101583662
Nu mean: 0.997187996002297
Marginal errors | max: 0.0113693316356223 | median 0.00939669372836249 | mean: 0.00860715734056932
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.2559814453125 | max: 3.9166259765625 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.00251007446998996 | median: 0.00144735987374958 | mean: 0.00144735987374958
Nu mean: 1.00106271459624
Marginal errors | max: 0.00406746973787442 | median 0.00343393462175534 | mean: 0.00313496757119527
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.25445556640625 | max: 3.9140625 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.000837027590858019 | median: 0.000488281249999889 | mean: 0.000488281249999889
Nu mean: 0.999651253659142
Marginal errors | max: 0.00143433714371533 | median 0.00116636244706747 | mean: 0.00109289166947839
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:10] to [1:10]
Estimate dispersion: DONE
dispersion | min: 0.2550048828125 | max: 3.9150390625 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 1 cells batches from [1:20] to [1:20]
Estimate nu: DONE
nu change (abs) | max: 0.000209227688885871 | median: 0.0001220703125 | mean: 0.0001220703125
Nu mean: 1.00008715737639
Marginal errors | max: 0.000379524846207957 | median 0.000325685250686547 | mean: 0.000295532844332236
Estimate dispersion/nu: DONE
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.54849195480347
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.000970125198364258
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00201797485351562
Estimating genes' coex
Calculate genes' coex elapsed time: 0.827430963516235
Total calculations elapsed time: 3.37891101837158
Calculate genes' coex (legacy): DONE
Calculate cells' coex (legacy): START
Retrieving expected cells' contingency table
calculating NN.. done
calculating YN..NY..YY..t().. done
Expected cells' contingency table elapsed time: 2.49949908256531
Calculating cells' coex normalization factor
Calculate cells' normalization factor elapsed time: 0.000559091567993164
Retrieving observed cells' yes/yes contingency table
calculating YY.. done
Observed cells' contingency table elapsed time: 0.00143098831176758
Estimating cells' coex
Calculate cells' coex elapsed time: 0.862774848937988
Total calculations elapsed time: 3.36426401138306
Calculate cells' coex (legacy): DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
Calculate genes' partial coex: START
Retrieving expected genes' partial contingency table
calculating partial NN.. done
calculating partial NY..YN..YY.. done
Calculating genes' partial coex normalization factor
Fraction of genes with very low expected contingency tables: 0.325
Retrieving observed genes' yes/yes partial contingency table
calculating partial YY.. done
Estimating genes' partial coex
Calculate genes' partial coex: DONE
Calculate cells' partial coex: START
Retrieving expected cells' partial contingency table
calculating partial NN.. done
calculating partial YN..NY..YY.. done
Calculating cells' partial coex normalization factor
Fraction of genes with very low expected contingency tables: 0
Retrieving observed cells' yes/yes partial contingency table
calculating partial YY.. done
Estimating cells' partial coex
Calculate cells' partial coex: DONE
Asked to drop 2 genes and 0 cells
Asked to drop 0 genes and 4 cells
Asked to drop 2 genes and 2 cells
Attaching package: 'rlang'
The following objects are masked from 'package:testthat':
is_false, is_null, is_true
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
calculating YY.. done
calculating YY.. done
calculating YN..NY..NN..t().. done
Estimate dispersion: START
Effective number of cores used: 1
Executing 3 genes batches from [1:3] to [8:10]
Estimate dispersion: DONE
dispersion | min: 0.9033203125 | max: 4.6796875 | % negative: 10
calculating NN.. done
calculating NN.. done
calculating NY..YN..YY..t().. done
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.73847389221191
Calculating genes' coex normalization factor
Fraction of genes' with very low expected contingency tables: 0.181818181818182
Calculate genes' normalization factor elapsed time: 0.000669956207275391
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00147104263305664
Estimating genes' coex
Calculate genes' coex elapsed time: 0.908044099807739
Total calculations elapsed time: 3.64865899085999
Calculate genes' coex (legacy): DONE
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.8790819644928
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.000809907913208008
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00224804878234863
Estimating genes' coex
Calculate genes' coex elapsed time: 0.964878082275391
Total calculations elapsed time: 3.84701800346375
Calculate genes' coex (legacy): DONE
Calculate genes' partial coex: START
Retrieving expected genes' partial contingency table
calculating partial NN.. done
calculating partial NY..YN..YY.. done
Calculating genes' partial coex normalization factor
Fraction of genes with very low expected contingency tables: 0.1
Retrieving observed genes' yes/yes partial contingency table
calculating partial YY.. done
Estimating genes' partial coex
Calculate genes' partial coex: DONE
Calculate genes' partial coex: START
Retrieving expected genes' partial contingency table
calculating partial NN.. done
calculating partial NY..YN..YY.. done
Calculating genes' partial coex normalization factor
Fraction of genes with very low expected contingency tables: 0.4
Retrieving observed genes' yes/yes partial contingency table
calculating partial YY.. done
Estimating genes' partial coex
Calculate genes' partial coex: DONE
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
calculating YY.. done
calculating YY.. done
calculating NY..YN..NN..t().. done
Estimate 'dispersion'/'nu': START
Estimate dispersion: START
Effective number of cores used: 1
Executing 3 genes batches from [1:3] to [8:10]
Estimate dispersion: DONE
dispersion | min: 0.903564453125 | max: 4.679443359375 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 4 cells batches from [1:4] to [13:16]
Executing 1 cells batches from [17:20] to [17:20]
Estimate nu: DONE
nu change (abs) | max: 1.75719246031746 | median: 1.07229953342014 | mean: 1.07229953342014
Nu mean: 1.68489292689732
Marginal errors | max: 0.25535328937316 | median 0.0807577993228135 | mean: 0.101980750205762
Estimate dispersion: START
Effective number of cores used: 1
Executing 3 genes batches from [1:3] to [8:10]
Estimate dispersion: DONE
dispersion | min: 1.037109375 | max: 4.6107177734375 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 4 cells batches from [1:4] to [13:16]
Executing 1 cells batches from [17:20] to [17:20]
Estimate nu: DONE
nu change (abs) | max: 0.0273438105507502 | median: 0.0148852611818011 | mean: 0.0148852611818011
Nu mean: 1.69735147626627
Marginal errors | max: 0.00326864580272002 | median 0.00111524657842743 | mean: 0.00131556083122533
Estimate dispersion: START
Effective number of cores used: 1
Executing 3 genes batches from [1:3] to [8:10]
Estimate dispersion: DONE
dispersion | min: 1.03887939453125 | max: 4.6097412109375 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 4 cells batches from [1:4] to [13:16]
Executing 1 cells batches from [17:20] to [17:20]
Estimate nu: DONE
nu change (abs) | max: 0 | median: 0 | mean: 0
Nu mean: 1.69735147626627
Marginal errors | max: 7.56328383637594e-05 | median 1.72948087246994e-05 | mean: 2.99252342145451e-05
Estimate dispersion/nu: DONE
calculating NN.. done
calculating NN.. done
calculating YN..NY..YY..t().. done
Calculate cells' coex (legacy): START
Retrieving expected cells' contingency table
calculating NN.. done
calculating YN..NY..YY..t().. done
Expected cells' contingency table elapsed time: 3.30222702026367
Calculating cells' coex normalization factor
Fraction of cells' with very low expected contingency tables: 0
Calculate cells' normalization factor elapsed time: 0.000974893569946289
Retrieving observed cells' yes/yes contingency table
calculating YY.. done
Observed cells' contingency table elapsed time: 0.00207996368408203
Estimating cells' coex
Calculate cells' coex elapsed time: 0.939470052719116
Total calculations elapsed time: 4.24475193023682
Calculate cells' coex (legacy): DONE
Calculate cells' partial coex: START
Retrieving expected cells' partial contingency table
calculating partial NN.. done
calculating partial YN..NY..YY.. done
Calculating cells' partial coex normalization factor
Fraction of genes with very low expected contingency tables: 0
Retrieving observed cells' yes/yes partial contingency table
calculating partial YY.. done
Estimating cells' partial coex
Calculate cells' partial coex: DONE
Calculate cells' partial coex: START
Retrieving expected cells' partial contingency table
calculating partial NN.. done
calculating partial YN..NY..YY.. done
Calculating cells' partial coex normalization factor
Fraction of genes with very low expected contingency tables: 0
Retrieving observed cells' yes/yes partial contingency table
calculating partial YY.. done
Estimating cells' partial coex
Calculate cells' partial coex: DONE
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
Estimate 'dispersion'/'nu': START
Estimate dispersion: START
Effective number of cores used: 1
Executing 3 genes batches from [1:3] to [8:10]
Estimate dispersion: DONE
dispersion | min: 0.903564453125 | max: 4.679443359375 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 4 cells batches from [1:4] to [13:16]
Executing 1 cells batches from [17:20] to [17:20]
Estimate nu: DONE
nu change (abs) | max: 1.75719246031746 | median: 1.07229953342014 | mean: 1.07229953342014
Nu mean: 1.68489292689732
Marginal errors | max: 0.25535328937316 | median 0.0807577993228135 | mean: 0.101980750205762
Estimate dispersion: START
Effective number of cores used: 1
Executing 3 genes batches from [1:3] to [8:10]
Estimate dispersion: DONE
dispersion | min: 1.037109375 | max: 4.6107177734375 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 4 cells batches from [1:4] to [13:16]
Executing 1 cells batches from [17:20] to [17:20]
Estimate nu: DONE
nu change (abs) | max: 0.0273438105507502 | median: 0.0148852611818011 | mean: 0.0148852611818011
Nu mean: 1.69735147626627
Marginal errors | max: 0.00326864580272002 | median 0.00111524657842743 | mean: 0.00131556083122533
Estimate dispersion: START
Effective number of cores used: 1
Executing 3 genes batches from [1:3] to [8:10]
Estimate dispersion: DONE
dispersion | min: 1.03887939453125 | max: 4.6097412109375 | % negative: 10
Estimate nu: START
Effective number of cores used: 1
Executing 4 cells batches from [1:4] to [13:16]
Executing 1 cells batches from [17:20] to [17:20]
Estimate nu: DONE
nu change (abs) | max: 0 | median: 0 | mean: 0
Nu mean: 1.69735147626627
Marginal errors | max: 7.56328383637594e-05 | median 1.72948087246994e-05 | mean: 2.99252342145451e-05
Estimate dispersion/nu: DONE
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.66298794746399
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.000825881958007812
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00202608108520508
Estimating genes' coex
Calculate genes' coex elapsed time: 0.865792989730835
Total calculations elapsed time: 3.53163290023804
Calculate genes' coex (legacy): DONE
Calculating S: START
Calculating S: DONE
Calculating G: START
calculating YY.. done
calculating YN..NY..NN..t().. done
calculating NN.. done
calculating NY..YN..YY..t().. done
Estimating G
Calculating G: DONE
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values genome wide on columns and genome wide on rows
calculating PValues: DONE
Using G
Calculating G: START
calculating YY.. done
calculating YN..NY..NN..t().. done
calculating NN.. done
calculating NY..YN..YY..t().. done
Estimating G
Calculating G: DONE
calculating PValues: START
Get p-values on a set of genes on columns and on a set of genes on rows
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
Calculate GDI dataframe: START
Using G
Calculating G: START
calculating YY.. done
calculating YN..NY..NN..t().. done
calculating NN.. done
calculating NY..YN..YY..t().. done
Estimating G
Calculating G: DONE
S matrix sorted
Calculate GDI dataframe: DONE
Calculating S: START
Calculating S: DONE
S matrix sorted
S matrix sorted
Initializing `COTAN` meta-data
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 3.3619019985199
Cotan genes' coex estimation started
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.23350214958191
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00951790809631348
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0453200340270996
Estimating genes' coex
Calculate genes' coex elapsed time: 0.84571099281311
Total calculations elapsed time: 4.13405108451843
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.03733801841736
Total elapsed time: 10.0579619407654
Initializing `COTAN` meta-data
Condition test
n cells 1200
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 3.39166712760925
Cotan genes' coex estimation started
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.22714710235596
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00900101661682129
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0505838394165039
Estimating genes' coex
Calculate genes' coex elapsed time: 0.928417205810547
Total calculations elapsed time: 4.21514916419983
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.07456994056702
Total elapsed time: 10.28191614151
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and genome wide on rows
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
Initializing `COTAN` meta-data
Condition test
n cells 1200
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 3.53620004653931
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.34893679618835
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00876712799072266
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.047882080078125
Estimating genes' coex
Calculate genes' coex elapsed time: 0.913444995880127
Total calculations elapsed time: 4.31903100013733
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.1916868686676
Total elapsed time: 10.5891377925873
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and genome wide on rows
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
Initializing `COTAN` meta-data
Condition test
n cells 1200
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 3.31586694717407
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.43107080459595
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00836706161499023
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0437519550323486
Estimating genes' coex
Calculate genes' coex elapsed time: 0.937448978424072
Total calculations elapsed time: 4.42063879966736
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.32046699523926
Total elapsed time: 10.4479939937592
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and genome wide on rows
calculating PValues: DONE
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
Initializing `COTAN` meta-data
Asked to drop 0 genes and 1160 cells
Cotan analysis functions started
Asked to drop 31 genes and 0 cells
Genes/cells selection done: dropped [31] genes and [0] cells
Working on [569] genes and [40] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:569] to [1:569]
Estimate dispersion: DONE
dispersion | min: -0.2119140625 | max: 85.5 | % negative: 62.2144112478032
Cotan genes' coex estimation not requested
Total elapsed time: 4.60278415679932
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.52506613731384
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.0172770023345947
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0258188247680664
Estimating genes' coex
Calculate genes' coex elapsed time: 0.844638109207153
Total calculations elapsed time: 3.41280007362366
Calculate genes' coex (legacy): DONE
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.4552628993988
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00545215606689453
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00685691833496094
Estimating genes' coex
Calculate genes' coex elapsed time: 0.877334117889404
Total calculations elapsed time: 3.34490609169006
Calculate genes' coex (legacy): DONE
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.698646068573
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00577902793884277
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.00680804252624512
Estimating genes' coex
Calculate genes' coex elapsed time: 0.955452919006348
Total calculations elapsed time: 3.66668605804443
Calculate genes' coex (legacy): DONE
Hangling COTAN object with condition: test
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and on a set of genes on rows
calculating PValues: DONE
Handling genes type: G2
Handling genes type: G3
Initializing `COTAN` meta-data
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
PCA: START
PCA: DONE
Hierarchical clustering: START
Hierarchical clustering: DONE
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Cotan genes' coex estimation not requested
Total elapsed time: 9.50845813751221
Saving elaborated data locally at: /tmp/RtmpJTHXF4/test.cotan.RDS
Creating cells' uniform clustering: START
In iteration 1 the number of cells to re-cluster is 1200 cells belonging to 0 clusters
Creating Seurat object: START
Normalizing layer: counts
Performing log-normalization
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Finding variable features for layer counts
Calculating gene variances
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Calculating feature variances of standardized and clipped values
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Centering and scaling data matrix
|
| | 0%
|
|======================================================================| 100%
PC_ 1
Positive: g-000558, g-000570, g-000499, g-000504, g-000546, g-000506, g-000503, g-000517, g-000596, g-000528
g-000527, g-000580, g-000592, g-000578, g-000509, g-000488, g-000555, g-000577, g-000534, g-000583
g-000598, g-000535, g-000512, g-000554, g-000519, g-000525, g-000548, g-000544, g-000502, g-000541
Negative: g-000133, g-000007, g-000074, g-000141, g-000057, g-000235, g-000170, g-000019, g-000195, g-000140
g-000183, g-000031, g-000046, g-000178, g-000177, g-000161, g-000157, g-000139, g-000011, g-000135
g-000125, g-000208, g-000061, g-000085, g-000204, g-000104, g-000237, g-000004, g-000038, g-000128
PC_ 2
Positive: g-000039, g-000050, g-000175, g-000078, g-000116, g-000189, g-000135, g-000047, g-000072, g-000087
g-000063, g-000235, g-000066, g-000109, g-000018, g-000074, g-000231, g-000136, g-000034, g-000207
g-000128, g-000167, g-000171, g-000049, g-000182, g-000013, g-000054, g-000062, g-000240, g-000158
Negative: g-000584, g-000583, g-000544, g-000519, g-000575, g-000516, g-000585, g-000486, g-000489, g-000539
g-000484, g-000502, g-000523, g-000595, g-000305, g-000574, g-000599, g-000589, g-000509, g-000538
g-000526, g-000551, g-000579, g-000590, g-000445, g-000556, g-000543, g-000501, g-000504, g-000570
PC_ 3
Positive: g-000015, g-000575, g-000483, g-000316, g-000025, g-000364, g-000050, g-000278, g-000443, g-000360
g-000332, g-000124, g-000212, g-000387, g-000536, g-000252, g-000251, g-000321, g-000501, g-000470
g-000582, g-000106, g-000455, g-000368, g-000081, g-000104, g-000437, g-000288, g-000386, g-000317
Negative: g-000211, g-000337, g-000129, g-000185, g-000397, g-000403, g-000253, g-000098, g-000390, g-000303
g-000052, g-000088, g-000463, g-000468, g-000236, g-000209, g-000005, g-000375, g-000342, g-000262
g-000388, g-000091, g-000413, g-000285, g-000003, g-000095, g-000142, g-000205, g-000432, g-000241
PC_ 4
Positive: g-000379, g-000193, g-000212, g-000434, g-000593, g-000513, g-000177, g-000223, g-000069, g-000131
g-000162, g-000345, g-000462, g-000484, g-000448, g-000229, g-000365, g-000302, g-000010, g-000366
g-000051, g-000535, g-000269, g-000270, g-000155, g-000529, g-000373, g-000008, g-000393, g-000306
Negative: g-000334, g-000398, g-000292, g-000095, g-000097, g-000202, g-000382, g-000195, g-000007, g-000079
g-000086, g-000240, g-000263, g-000317, g-000576, g-000557, g-000160, g-000154, g-000214, g-000228
g-000313, g-000053, g-000524, g-000374, g-000568, g-000188, g-000358, g-000528, g-000362, g-000150
PC_ 5
Positive: g-000451, g-000339, g-000295, g-000328, g-000544, g-000061, g-000227, g-000391, g-000556, g-000237
g-000067, g-000165, g-000449, g-000591, g-000087, g-000129, g-000197, g-000203, g-000487, g-000505
g-000333, g-000029, g-000271, g-000064, g-000583, g-000156, g-000448, g-000153, g-000526, g-000393
Negative: g-000518, g-000108, g-000186, g-000170, g-000401, g-000337, g-000047, g-000599, g-000432, g-000578
g-000042, g-000065, g-000493, g-000261, g-000533, g-000256, g-000560, g-000596, g-000368, g-000381
g-000535, g-000338, g-000215, g-000159, g-000365, g-000234, g-000173, g-000387, g-000225, g-000272
Computing nearest neighbor graph
Computing SNN
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 1200
Number of edges: 55489
Running Louvain algorithm with multilevel refinement...
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Maximum modularity in 10 random starts: 0.5973
Number of communities: 4
Elapsed time: 0 seconds
Used resolution for Seurat clusterization is: 0.8
18:41:40 UMAP embedding parameters a = 0.9922 b = 1.112
18:41:40 Read 1200 rows and found 50 numeric columns
18:41:40 Using Annoy for neighbor search, n_neighbors = 30
18:41:40 Building Annoy index with metric = cosine, n_trees = 50
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
18:41:40 Writing NN index file to temp file /tmp/RtmpJTHXF4/file104557f5c2608
18:41:40 Searching Annoy index using 1 thread, search_k = 3000
18:41:41 Annoy recall = 100%
18:41:42 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
18:41:44 Initializing from normalized Laplacian + noise (using RSpectra)
18:41:44 Commencing optimization for 500 epochs, with 42270 positive edges
Using method 'umap'
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
18:41:46 Optimization finished
Creating PDF UMAP in file: /tmp/RtmpJTHXF4/test/reclustering/pdf_umap_1.pdf
Creating Seurat object: DONE
* checking uniformity of cluster '0' of 4 clusters
Asked to drop 0 genes and 847 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [353] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.038818359375 | max: 11.109375 | % negative: 5
Only analysis elapsed time: 3.0869300365448
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.16611790657043
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00688409805297852
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.040895938873291
Estimating genes' coex
Calculate genes' coex elapsed time: 1.07380104064941
Total calculations elapsed time: 4.28769898414612
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.17698693275452
Total elapsed time: 10.1495878696442
Checking uniformity for the cluster '01_0000' with 353 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0000, with size 353, is uniform
check.isCheckAbove = 0, 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 = 1, 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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.0501708984375 | max: 14.515625 | % negative: 8.5
Only analysis elapsed time: 2.97205305099487
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.0971839427948
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00718998908996582
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0351228713989258
Estimating genes' coex
Calculate genes' coex elapsed time: 0.931315183639526
Total calculations elapsed time: 4.07081198692322
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.02519392967224
Total elapsed time: 9.95584893226624
Checking uniformity for the cluster '01_0001' with 321 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0001, with size 321, is uniform
check.isCheckAbove = 0, 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 = 1, 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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.054443359375 | max: 104 | % negative: 36.5
Only analysis elapsed time: 2.92716503143311
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.9569981098175
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.0103728771209717
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.016334056854248
Estimating genes' coex
Calculate genes' coex elapsed time: 1.0110650062561
Total calculations elapsed time: 3.99477005004883
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 4.95594382286072
Total elapsed time: 9.76996397972107
Checking uniformity for the cluster '01_0002' with 295 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0002, with size 295, is uniform
check.isCheckAbove = 0, 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 = 1, 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 from [1:599] to [1:599]
Estimate dispersion: DONE
dispersion | min: -0.0662841796875 | max: 82.5 | % negative: 33.889816360601
Only analysis elapsed time: 2.88409399986267
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.01539993286133
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.0213508605957031
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0202090740203857
Estimating genes' coex
Calculate genes' coex elapsed time: 0.967530965805054
Total calculations elapsed time: 4.02449083328247
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 4.87022399902344
Total elapsed time: 9.58199286460876
Checking uniformity for the cluster '01_0003' with 231 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0003, with size 231, is uniform
check.isCheckAbove = 0, 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 = 1, clusterSize = 231
cluster 01_0003 is uniform
Found 4 uniform and 0 non-uniform clusters
NO new possible uniform clusters! Unclustered cell left: 0
The final raw clusterization contains [ 4 ] different clusters: 01_0000, 01_0001, 01_0002, 01_0003
Differential Expression Analysis - START
* DEA on cluster '1' with 353 cells
* DEA on cluster '2' with 321 cells
* DEA on cluster '3' with 295 cells
* DEA on cluster '4' with 231 cells
Differential Expression Analysis - DONE
Applied reordering to clusterization is:
1 -> 3, 2 -> 1, 3 -> 2, 4 -> 4
Creating cells' uniform clustering: DONE
Applied reordering to clusterization is:
1 -> 1, 2 -> 2, 3 -> 3, 4 -> 4
Asked to drop 0 genes and 905 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [295] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.054443359375 | max: 104 | % negative: 36.5
Only analysis elapsed time: 2.88822102546692
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.97659587860107
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00971603393554688
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0174720287322998
Estimating genes' coex
Calculate genes' coex elapsed time: 0.981085062026978
Total calculations elapsed time: 3.9848690032959
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 4.92576003074646
Total elapsed time: 9.7654390335083
Checking uniformity for the cluster 'Cluster_2' with 295 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster Cluster_2, with size 295, is uniform
firstCheck.isCheckAbove = 0, 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 = 1, 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 = 0, 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 = 0, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 5, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.4279762218014, isUniform = 1, 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
Creating Seurat object: START
Normalizing layer: counts
Performing log-normalization
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 3.54986000061035
Cotan genes' coex estimation started
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.51986980438232
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.0063471794128418
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0655770301818848
Estimating genes' coex
Calculate genes' coex elapsed time: 0.871279001235962
Total calculations elapsed time: 4.46307301521301
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.46463799476624
Total elapsed time: 10.8578388690948
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
Centering and scaling data matrix
|
| | 0%
|
|======================================================================| 100%
PC_ 1
Positive: g-000558, g-000570, g-000499, g-000504, g-000546, g-000506, g-000503, g-000517, g-000596, g-000528
g-000527, g-000580, g-000592, g-000578, g-000509, g-000488, g-000555, g-000577, g-000534, g-000583
g-000598, g-000535, g-000512, g-000554, g-000519, g-000525, g-000548, g-000544, g-000502, g-000541
Negative: g-000133, g-000007, g-000074, g-000141, g-000057, g-000235, g-000170, g-000019, g-000195, g-000140
g-000183, g-000031, g-000046, g-000178, g-000177, g-000161, g-000157, g-000139, g-000011, g-000135
g-000125, g-000208, g-000061, g-000085, g-000204, g-000104, g-000237, g-000004, g-000038, g-000128
PC_ 2
Positive: g-000039, g-000050, g-000175, g-000078, g-000116, g-000189, g-000135, g-000047, g-000072, g-000087
g-000063, g-000235, g-000066, g-000109, g-000018, g-000074, g-000231, g-000136, g-000034, g-000207
g-000128, g-000167, g-000171, g-000049, g-000182, g-000013, g-000054, g-000062, g-000240, g-000158
Negative: g-000584, g-000583, g-000544, g-000519, g-000575, g-000516, g-000585, g-000486, g-000489, g-000539
g-000484, g-000502, g-000523, g-000595, g-000305, g-000574, g-000599, g-000589, g-000509, g-000538
g-000526, g-000551, g-000579, g-000590, g-000445, g-000556, g-000543, g-000501, g-000504, g-000570
PC_ 3
Positive: g-000015, g-000575, g-000483, g-000316, g-000025, g-000364, g-000050, g-000278, g-000443, g-000360
g-000332, g-000124, g-000212, g-000387, g-000536, g-000252, g-000251, g-000321, g-000501, g-000470
g-000582, g-000106, g-000455, g-000368, g-000081, g-000104, g-000437, g-000288, g-000386, g-000317
Negative: g-000211, g-000337, g-000129, g-000185, g-000397, g-000403, g-000253, g-000098, g-000390, g-000303
g-000052, g-000088, g-000463, g-000468, g-000236, g-000209, g-000005, g-000375, g-000342, g-000262
g-000388, g-000091, g-000413, g-000285, g-000003, g-000095, g-000142, g-000205, g-000432, g-000241
PC_ 4
Positive: g-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: 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
18:43:29 UMAP embedding parameters a = 0.9922 b = 1.112
18:43:29 Read 1200 rows and found 50 numeric columns
18:43:29 Using Annoy for neighbor search, n_neighbors = 30
18:43:29 Building Annoy index with metric = cosine, n_trees = 50
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
18:43:29 Writing NN index file to temp file /tmp/RtmpJTHXF4/file1045573512c66
18:43:29 Searching Annoy index using 1 thread, search_k = 3000
18:43:29 Annoy recall = 100%
18:43:30 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
18:43:32 Initializing from normalized Laplacian + noise (using RSpectra)
18:43:32 Commencing optimization for 500 epochs, with 42246 positive edges
Using method 'umap'
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
18:43:35 Optimization finished
Creating PDF UMAP in file: /tmp/RtmpJTHXF4/test/reclustering/pdf_umap_1.pdf
Creating Seurat object: DONE
Using passed in clusterization
* checking uniformity of cluster '1' of 2 clusters
Asked to drop 0 genes and 600 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [600] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.0386962890625 | max: 19.40625 | % negative: 6.5
Only analysis elapsed time: 3.16951012611389
Cotan genes' coex estimation started
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.1760458946228
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00792503356933594
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0327839851379395
Estimating genes' coex
Calculate genes' coex elapsed time: 0.940655946731567
Total calculations elapsed time: 4.15741086006165
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.07613086700439
Total elapsed time: 10.1218140125275
Checking uniformity for the cluster '01_0001' with 600 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0001, with size 600, is uniform
firstCheck.isCheckAbove = 0, 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 = 1, 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 = 0, 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 = 0, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 4, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.51521594648007, isUniform = 1, 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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.032958984375 | max: 10.0859375 | % negative: 3.66666666666667
Only analysis elapsed time: 3.08177399635315
Cotan genes' coex estimation started
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.16224193572998
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00791406631469727
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0398180484771729
Estimating genes' coex
Calculate genes' coex elapsed time: 0.92964506149292
Total calculations elapsed time: 4.13961911201477
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.06401109695435
Total elapsed time: 9.91395998001099
Checking uniformity for the cluster '01_0002' with 600 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 01_0002, with size 600, is uniform
firstCheck.isCheckAbove = 0, 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 = 1, 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 = 0, 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 = 0, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 2, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.50182972335853, isUniform = 1, clusterSize = 600
cluster 01_0002 is uniform
Found 2 uniform and 0 non-uniform clusters
NO new possible uniform clusters! Unclustered cell left: 0
The final raw clusterization contains [ 2 ] different clusters: 01_0001, 01_0002
Differential Expression Analysis - START
* DEA on cluster '1' with 600 cells
* DEA on cluster '2' with 600 cells
Differential Expression Analysis - DONE
Applied reordering to clusterization is:
1 -> 1, 2 -> 2
Creating cells' uniform clustering: DONE
findClustersMarkers - START
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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.0501708984375 | max: 14.515625 | % negative: 8.5
Only analysis elapsed time: 2.98437404632568
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.08042192459106
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00674700736999512
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0371239185333252
Estimating genes' coex
Calculate genes' coex elapsed time: 0.928807020187378
Total calculations elapsed time: 4.05309987068176
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.02966594696045
Total elapsed time: 9.93330502510071
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
[1] "Tested cluster: 3"
Asked to drop 0 genes and 847 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [353] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.038818359375 | max: 11.109375 | % negative: 5
Only analysis elapsed time: 2.90614008903503
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 2.9804060459137
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00711178779602051
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0302600860595703
Estimating genes' coex
Calculate genes' coex elapsed time: 0.882358074188232
Total calculations elapsed time: 3.90013599395752
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 4.78937005996704
Total elapsed time: 9.51611113548279
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [10] genes and [20] cells
Initializing `COTAN` meta-data
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Initializing `COTAN` meta-data
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Cotan genes' coex estimation not requested
Total elapsed time: 6.19049286842346
Creating Seurat object: START
Normalizing layer: counts
Performing log-normalization
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Finding variable features for layer counts
Calculating gene variances
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Calculating feature variances of standardized and clipped values
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Centering and scaling data matrix
|
| | 0%
|
|======================================================================| 100%
PC_ 1
Positive: g-000558, g-000570, g-000499, g-000504, g-000546, g-000506, g-000503, g-000517, g-000596, g-000528
g-000527, g-000580, g-000592, g-000578, g-000509, g-000488, g-000555, g-000577, g-000534, g-000583
g-000598, g-000535, g-000512, g-000554, g-000519, g-000525, g-000548, g-000544, g-000502, g-000541
Negative: g-000133, g-000007, g-000074, g-000141, g-000057, g-000235, g-000170, g-000019, g-000195, g-000140
g-000183, g-000031, g-000046, g-000178, g-000177, g-000161, g-000157, g-000139, g-000011, g-000135
g-000125, g-000208, g-000061, g-000085, g-000204, g-000104, g-000237, g-000004, g-000038, g-000128
PC_ 2
Positive: g-000039, g-000050, g-000175, g-000078, g-000116, g-000189, g-000135, g-000047, g-000072, g-000087
g-000063, g-000235, g-000066, g-000109, g-000018, g-000074, g-000231, g-000136, g-000034, g-000207
g-000128, g-000167, g-000171, g-000049, g-000182, g-000013, g-000054, g-000062, g-000240, g-000158
Negative: g-000584, g-000583, g-000544, g-000519, g-000575, g-000516, g-000585, g-000486, g-000489, g-000539
g-000484, g-000502, g-000523, g-000595, g-000305, g-000574, g-000599, g-000589, g-000509, g-000538
g-000526, g-000551, g-000579, g-000590, g-000445, g-000556, g-000543, g-000501, g-000504, g-000570
PC_ 3
Positive: g-000015, g-000575, g-000483, g-000316, g-000025, g-000364, g-000050, g-000278, g-000443, g-000360
g-000332, g-000124, g-000212, g-000387, g-000536, g-000252, g-000251, g-000321, g-000501, g-000470
g-000582, g-000106, g-000455, g-000368, g-000081, g-000104, g-000437, g-000288, g-000386, g-000317
Negative: g-000211, g-000337, g-000129, g-000185, g-000397, g-000403, g-000253, g-000098, g-000390, g-000303
g-000052, g-000088, g-000463, g-000468, g-000236, g-000209, g-000005, g-000375, g-000342, g-000262
g-000388, g-000091, g-000413, g-000285, g-000003, g-000095, g-000142, g-000205, g-000432, g-000241
PC_ 4
Positive: g-000379, g-000193, g-000212, g-000434, g-000593, g-000513, g-000177, g-000223, g-000069, g-000131
g-000162, g-000345, g-000462, g-000484, g-000448, g-000229, g-000365, g-000302, g-000010, g-000366
g-000051, g-000535, g-000269, g-000270, g-000155, g-000529, g-000373, g-000008, g-000393, g-000306
Negative: g-000334, g-000398, g-000292, g-000095, g-000097, g-000202, g-000382, g-000195, g-000007, g-000079
g-000086, g-000240, g-000263, g-000317, g-000576, g-000557, g-000160, g-000154, g-000214, g-000228
g-000313, g-000053, g-000524, g-000374, g-000568, g-000188, g-000358, g-000528, g-000362, g-000150
PC_ 5
Positive: g-000451, g-000339, g-000295, g-000328, g-000544, g-000061, g-000227, g-000391, g-000556, g-000237
g-000067, g-000165, g-000449, g-000591, g-000087, g-000129, g-000197, g-000203, g-000487, g-000505
g-000333, g-000029, g-000271, g-000064, g-000583, g-000156, g-000448, g-000153, g-000526, g-000393
Negative: g-000518, g-000108, g-000186, g-000170, g-000401, g-000337, g-000047, g-000599, g-000432, g-000578
g-000042, g-000065, g-000493, g-000261, g-000533, g-000256, g-000560, g-000596, g-000368, g-000381
g-000535, g-000338, g-000215, g-000159, g-000365, g-000234, g-000173, g-000387, g-000225, g-000272
Computing nearest neighbor graph
Computing SNN
Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
Number of nodes: 1200
Number of edges: 55489
Running Louvain algorithm with multilevel refinement...
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Maximum modularity in 10 random starts: 0.5973
Number of communities: 4
Elapsed time: 0 seconds
Used resolution for Seurat clusterization is: 0.8
18:45:04 UMAP embedding parameters a = 0.9922 b = 1.112
18:45:04 Read 1200 rows and found 50 numeric columns
18:45:04 Using Annoy for neighbor search, n_neighbors = 30
18:45:04 Building Annoy index with metric = cosine, n_trees = 50
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
18:45:04 Writing NN index file to temp file /tmp/RtmpJTHXF4/file104555857db9e
18:45:04 Searching Annoy index using 1 thread, search_k = 3000
18:45:04 Annoy recall = 100%
18:45:05 Commencing smooth kNN distance calibration using 1 thread with target n_neighbors = 30
18:45:07 Initializing from normalized Laplacian + noise (using RSpectra)
18:45:07 Commencing optimization for 500 epochs, with 42270 positive edges
Using method 'umap'
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
18:45:10 Optimization finished
Creating Seurat object: DONE
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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Cotan genes' coex estimation not requested
Total elapsed time: 6.31002116203308
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 3.46503090858459
Cotan genes' coex estimation started
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.38273000717163
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00861501693725586
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0420129299163818
Estimating genes' coex
Calculate genes' coex elapsed time: 0.988624095916748
Total calculations elapsed time: 4.42198204994202
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.32884407043457
Total elapsed time: 10.6703410148621
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
S matrix sorted
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
S matrix sorted
Number of columns (V set - secondary markers): 109
Number of rows (U set): 60
Calculating gene co-expression space - DONE
Establishing gene clusters - START
Calculating gene co-expression space - START
Using S
Calculating S: START
Calculating S: DONE
calculating PValues: START
Get p-values on a set of genes on columns and genome wide on rows
calculating PValues: DONE
Number of selected secondary markers: 109
Calculating S: START
Calculating S: DONE
S matrix sorted
Number of columns (V set - secondary markers): 109
Number of rows (U set): 60
Calculating gene co-expression space - DONE
Establishing gene clusters - DONE
Initializing `COTAN` meta-data
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [1200] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Cotan genes' coex estimation not requested
Total elapsed time: 6.40424299240112
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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.0340576171875 | max: 10.4375 | % negative: 4.33333333333333
Only analysis elapsed time: 3.13513207435608
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.13448286056519
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00781607627868652
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0384979248046875
Estimating genes' coex
Calculate genes' coex elapsed time: 0.949808120727539
Total calculations elapsed time: 4.1306049823761
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.04301691055298
Total elapsed time: 10.0877709388733
Checking uniformity for the cluster '1_2-merge' with 574 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 1_2-merge, with size 574, is not uniform
firstCheck.isCheckAbove = 0, 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 = 1, 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 = 0, 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 = 0, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 20, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.5067755752749, isUniform = 0, 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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.03839111328125 | max: 17.40625 | % negative: 6.66666666666667
Only analysis elapsed time: 3.15704393386841
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.18448805809021
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00696992874145508
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0337111949920654
Estimating genes' coex
Calculate genes' coex elapsed time: 0.937991857528687
Total calculations elapsed time: 4.16316103935242
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.09931993484497
Total elapsed time: 10.1389260292053
Checking uniformity for the cluster '3_4-merge' with 626 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 3_4-merge, with size 626, is not uniform
firstCheck.isCheckAbove = 0, 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 = 1, 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 = 0, 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 = 0, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 22, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.52592182249635, isUniform = 0, 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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.0390625 | max: 12.59375 | % negative: 9.16666666666667
Only analysis elapsed time: 3.11973285675049
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.14291501045227
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00905394554138184
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0278501510620117
Estimating genes' coex
Calculate genes' coex elapsed time: 0.965594053268433
Total calculations elapsed time: 4.1454131603241
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.1109881401062
Total elapsed time: 10.18323802948
Checking uniformity for the cluster '2_3-merge' with 536 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 2_3-merge, with size 536, is not uniform
firstCheck.isCheckAbove = 0, 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 = 1, 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 = 0, 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 = 0, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 30, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.54477108741147, isUniform = 0, 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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.04754638671875 | max: 10.234375 | % negative: 5
Only analysis elapsed time: 3.11333394050598
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.09881687164307
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00654911994934082
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0313849449157715
Estimating genes' coex
Calculate genes' coex elapsed time: 0.957501888275146
Total calculations elapsed time: 4.09425282478333
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.04372215270996
Total elapsed time: 10.035590171814
Checking uniformity for the cluster '2_4-merge' with 532 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 2_4-merge, with size 532, is not uniform
firstCheck.isCheckAbove = 0, 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 = 1, 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 = 0, 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 = 0, fourthCheck.GDIThreshold = 1.4, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 109, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.68475283375668, isUniform = 0, 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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.046630859375 | max: 9.03515625 | % negative: 1.33333333333333
Only analysis elapsed time: 3.35406589508057
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.26157593727112
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00758218765258789
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0391058921813965
Estimating genes' coex
Calculate genes' coex elapsed time: 0.940958976745605
Total calculations elapsed time: 4.24922299385071
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.23888397216797
Total elapsed time: 10.4745888710022
Checking uniformity for the cluster '1_3-merge' with 668 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 1_3-merge, with size 668, is not uniform
firstCheck.isCheckAbove = 0, 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 = 1, 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 = 0, 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 = 0, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 77, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.68881338289977, isUniform = 0, 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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.0164794921875 | max: 5.5595703125 | % negative: 0.166666666666667
Only analysis elapsed time: 3.27414083480835
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.21557211875916
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00650596618652344
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0410559177398682
Estimating genes' coex
Calculate genes' coex elapsed time: 0.994947195053101
Total calculations elapsed time: 4.25808119773865
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.22615098953247
Total elapsed time: 10.4477868080139
Checking uniformity for the cluster '1_4-merge' with 664 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 1_4-merge, with size 664, is not uniform
firstCheck.isCheckAbove = 0, 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 = 1, 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 = 0, 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 = 0, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 59, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.64605433203111, isUniform = 0, 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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: 0.2197265625 | max: 6.08056640625 | % negative: 0
Only analysis elapsed time: 3.59679794311523
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.57318115234375
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00766777992248535
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.0499651432037354
Estimating genes' coex
Calculate genes' coex elapsed time: 0.947123050689697
Total calculations elapsed time: 4.57793712615967
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.50250506401062
Total elapsed time: 11.0805320739746
Checking uniformity for the cluster '1_2-merge_3_4-merge-merge' with 1200 cells
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
GDI plot
Removed 0 low GDI genes (such as the fully-expressed) in GDI plot
Cluster 1_2-merge_3_4-merge-merge, with size 1200, is not uniform
firstCheck.isCheckAbove = 0, 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 = 1, 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 = 0, 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 = 0, fourthCheck.GDIThreshold = 1.5, fourthCheck.maxRatioBeyond = NaN, fourthCheck.maxRankBeyond = 2, fourthCheck.fractionBeyond = NaN, fourthCheck.thresholdRank = 90, fourthCheck.quantileAtRatio = NaN, fourthCheck.quantileAtRank = 1.66512085964758, isUniform = 0, 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 from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.0340576171875 | max: 10.4375 | % negative: 4.33333333333333
Only analysis elapsed time: 3.09929704666138
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.08357095718384
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.00797605514526367
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.039078950881958
Estimating genes' coex
Calculate genes' coex elapsed time: 0.936607122421265
Total calculations elapsed time: 4.06723308563232
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.02627086639404
Total elapsed time: 9.96630883216858
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
Asked to drop 0 genes and 574 cells
Cotan analysis functions started
Genes/cells selection done: dropped [0] genes and [0] cells
Working on [600] genes and [626] cells
Estimate dispersion: START
Effective number of cores used: 1
Executing 1 genes batches from [1:600] to [1:600]
Estimate dispersion: DONE
dispersion | min: -0.03839111328125 | max: 17.40625 | % negative: 6.66666666666667
Only analysis elapsed time: 3.19249296188354
Cotan genes' coex estimation started
While trying to load the `torch` library Error in doTryCatch(return(expr), name, parentenv, handler): The `torch` library is installed but the required additional libraries are not avalable yet
Calculate genes' coex (legacy): START
Retrieving expected genes' contingency table
calculating NN.. done
calculating NY..YN..YY..t().. done
Expected genes' contingency table elapsed time: 3.65121793746948
Calculating genes' coex normalization factor
Calculate genes' normalization factor elapsed time: 0.0865349769592285
Retrieving observed genes' yes/yes contingency table
calculating YY.. done
Observed genes' contingency table elapsed time: 0.082050085067749
Estimating genes' coex
Calculate genes' coex elapsed time: 0.987390995025635
Total calculations elapsed time: 4.80719399452209
Calculate genes' coex (legacy): DONE
Only genes' coex elapsed time: 5.91460108757019
Total elapsed time: 11.0437519550323
Calculate GDI dataframe: START
Using S
Calculating S: START
Calculating S: DONE
S matrix sorted
Calculate GDI dataframe: DONE
[ FAIL 0 | WARN 0 | SKIP 0 | PASS 498 ]
>
> proc.time()
user system elapsed
748.167 12.965 775.534
COTAN.Rcheck/COTAN-Ex.timings
| name | user | system | elapsed | |
| COTAN_Legacy | 0.003 | 0.002 | 0.004 | |
| COTAN_ObjectCreation | 10.732 | 0.333 | 11.137 | |
| CalculatingCOEX | 58.574 | 1.312 | 60.279 | |
| ClustersList | 0.007 | 0.003 | 0.010 | |
| Conversions | 6.167 | 0.148 | 6.347 | |
| GenesStatistics | 11.943 | 0.423 | 12.199 | |
| HandleMetaData | 0.089 | 0.007 | 0.096 | |
| HandlingClusterizations | 31.896 | 1.480 | 32.975 | |
| HandlingConditions | 0.127 | 0.007 | 0.136 | |
| HeatmapPlots | 39.935 | 1.011 | 41.239 | |
| LoggingFunctions | 0.002 | 0.001 | 0.004 | |
| ParametersEstimations | 27.417 | 0.857 | 28.466 | |
| RawDataCleaning | 4.988 | 0.329 | 5.113 | |
| RawDataGetters | 0.080 | 0.005 | 0.086 | |
| UniformClusters | 222.538 | 3.946 | 226.495 | |
| getColorsVector | 0.001 | 0.016 | 0.001 | |