library(RaceID)
library(Seurat)
## Attaching SeuratObject
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
library(ggplot2)
pbmc <- readRDS("/mnt/nectar_volume/home/eraz0001/pbmc_human/pbmc_tutorial_final.rds")
ndata <- pbmc@assays$RNA@counts[,]
sc <- SCseq(ndata)
The first step is the application of filtering for the purpose of quality control. Cells with a relatively low total number of transcripts are discarded.
sc <- filterdata(sc,mintotal=2000)
In this example, we filter out cells with <2,000 transcipts. The function allows further filtering of genes by choosing the input parameters minexpr and minnumber, i.e. only genes with at least minexpr transcripts in at least minnumber cells are retained. The filtered and normalized expression matrix (normalized to the minimum total transcript count across all cells retained after filtering) can be retrieved by the getfdata function:
fdata <- getfdata(sc)
sc <- compdist(sc,metric="pearson")
sc <- clustexp(sc)
## Clustering k = 1,2,..., K.max (= 30): ..
## k = 1
k = 2
k = 3
k = 4
k = 5
k = 6
k = 7
k = 8
k = 9
k = 10
k = 11
k = 12
k = 13
k = 14
k = 15
k = 16
k = 17
k = 18
k = 19
k = 20
k = 21
k = 22
k = 23
k = 24
k = 25
k = 26
k = 27
k = 28
k = 29
k = 30
## done.
## boot 1
## boot 2
## boot 3
## boot 4
## boot 5
## boot 6
## boot 7
## boot 8
## boot 9
## boot 10
## boot 11
## boot 12
## boot 13
## boot 14
## boot 15
## boot 16
## boot 17
## boot 18
## boot 19
## boot 20
## boot 21
## boot 22
## boot 23
## boot 24
## boot 25
## boot 26
## boot 27
## boot 28
## boot 29
## boot 30
## boot 31
## boot 32
## boot 33
## boot 34
## boot 35
## boot 36
## boot 37
## boot 38
## boot 39
## boot 40
## boot 41
## boot 42
## boot 43
## boot 44
## boot 45
## boot 46
## boot 47
## boot 48
## boot 49
## boot 50
plotsaturation(sc,disp=FALSE)
plotsaturation(sc,disp=TRUE)
plotjaccard(sc)
Here, the automated criterion overestimated the cluster number leading
to instability as indicated by low Jaccard’s similarity. Based on visual
inspection of the average within-cluster dispersion as a function of the
cluster number, we manually set the cluster number to 10 without
recomputing the saturation behaviour.
sc <- clustexp(sc,cln=10,sat=FALSE)
## boot 1
## boot 2
## boot 3
## boot 4
## boot 5
## boot 6
## boot 7
## boot 8
## boot 9
## boot 10
## boot 11
## boot 12
## boot 13
## boot 14
## boot 15
## boot 16
## boot 17
## boot 18
## boot 19
## boot 20
## boot 21
## boot 22
## boot 23
## boot 24
## boot 25
## boot 26
## boot 27
## boot 28
## boot 29
## boot 30
## boot 31
## boot 32
## boot 33
## boot 34
## boot 35
## boot 36
## boot 37
## boot 38
## boot 39
## boot 40
## boot 41
## boot 42
## boot 43
## boot 44
## boot 45
## boot 46
## boot 47
## boot 48
## boot 49
## boot 50
This function perform k-medoids clustering by default. K-means clustering or hierarchical clustering can be chosen with the FUNcluster argument. For very large datasets, hierarchical clustering leads to significantly smaller run time.
Subsequently, outliers in the initial k-medoids clusters are identified based on an internally computed background model for the expected gene expression variability and the assumption that transcript counts follow a negative binomial distribution defined by the mean and the variance of the expression of each gene in a cluster. Outlier genes will be in the tail of this distribution at a p-value defined by the probthr parameter (1e-3 by default), and outlier cells require the presence of a number of outlier genes defined by the outlg parameter (2 by default).
sc <- findoutliers(sc)
## find outliers in cluster 1
find outliers in cluster 2
find outliers in cluster 3
find outliers in cluster 4
find outliers in cluster 5
find outliers in cluster 6
find outliers in cluster 7
find outliers in cluster 8
find outliers in cluster 9
find outliers in cluster 10
## merging outliers 1
merging outliers 2
merging outliers 3
merging outliers 4
merging outliers 5
merging outliers 6
merging outliers 7
merging outliers 8
merging outliers 9
merging outliers 10
merging outliers 11
merging outliers 12
merging outliers 13
merging outliers 14
merging outliers 15
merging outliers 16
merging outliers 17
merging outliers 18
merging outliers 19
merging outliers 20
merging outliers 21
merging outliers 22
merging outliers 23
merging outliers 24
merging outliers 25
merging outliers 26
merging outliers 27
merging outliers 28
merging outliers 29
merging outliers 30
merging outliers 31
merging outliers 32
merging outliers 33
merging outliers 34
merging outliers 35
merging outliers 36
merging outliers 37
merging outliers 38
merging outliers 39
merging outliers 40
merging outliers 41
merging outliers 42
merging outliers 43
## determine final clustering partition 1
determine final clustering partition 2
determine final clustering partition 3
determine final clustering partition 4
determine final clustering partition 5
determine final clustering partition 6
determine final clustering partition 7
determine final clustering partition 8
determine final clustering partition 9
determine final clustering partition 10
determine final clustering partition 11
determine final clustering partition 12
determine final clustering partition 13
determine final clustering partition 14
determine final clustering partition 15
determine final clustering partition 16
determine final clustering partition 17
determine final clustering partition 18
determine final clustering partition 19
determine final clustering partition 20
determine final clustering partition 21
determine final clustering partition 22
determine final clustering partition 23
determine final clustering partition 24
determine final clustering partition 25
determine final clustering partition 26
determine final clustering partition 27
determine final clustering partition 28
determine final clustering partition 29
determine final clustering partition 30
determine final clustering partition 31
determine final clustering partition 32
determine final clustering partition 33
determine final clustering partition 34
determine final clustering partition 35
determine final clustering partition 36
determine final clustering partition 37
determine final clustering partition 38
determine final clustering partition 39
determine final clustering partition 40
determine final clustering partition 41
determine final clustering partition 42
determine final clustering partition 43
determine final clustering partition 44
determine final clustering partition 45
determine final clustering partition 46
determine final clustering partition 47
determine final clustering partition 48
determine final clustering partition 49
determine final clustering partition 50
determine final clustering partition 51
determine final clustering partition 52
determine final clustering partition 53
The background noise model can be inspected:
plotbackground(sc)
plotsensitivity(sc)
Another way of checking the presence of outliers is the inspection of a
barplot of p-values across all cells in each cluster:
plotoutlierprobs(sc)
A heatmap of cell-to-cell distances grouped be the final clusters
inferred based on the original clusters and the outliers allows visual
inspection of the clusters:
clustheatmap(sc)
## [1] 52 21 53 50 49 27 48 36 47 41 9 51 40 46 35 15 24 4 7 12 6 13 14 45 34
## [26] 33 31 42 26 44 37 30 18 29 22 32 28 19 38 20 23 43 16 39 2 17 1 3 5 11
## [51] 10 8 25
sc <- compumap(sc)
Parameters for umap can be given as additional argument
plotmap(sc)
sc <- comptsne(sc)
sc <- compfr(sc,knn=10)
The t-SNE map can be plotted by
plotmap(sc)
The Fruchterman-Rheingold layout can be plotted by
plotmap(sc,fr=TRUE)
umap presentation
plotmap(sc,um=TRUE)
The dimensional reduction maps can be inspected, e.g., for the
localization of (a subset of) samples included in the analysis:
types <- sub("(\\_\\d+)$","", colnames(sc@ndata))
plotsymbolsmap(sc,types,fr=TRUE)
plotexpmap(sc,"CCR7",logsc=TRUE,fr=TRUE)
g <- c("S100A4", "IL7R")
plotexpmap(sc,g,n="Memory CD4+",logsc=TRUE,fr=TRUE)
For instance, to obtain differentially expressed genes within cluster 4
compared to all other cells:
d <- clustdiffgenes(sc,4,pvalue=.01)
dg <- d$dg
head(dg,25)
## mean.ncl mean.cl fc pv padj
## FTH1 9.6290373 31.6791749 3.28996284 0.000000e+00 0.000000e+00
## MALAT1 34.1778602 11.6575998 0.34108630 0.000000e+00 0.000000e+00
## LYZ 4.3617687 21.0657130 4.82962628 0.000000e+00 0.000000e+00
## CST3 1.8530394 8.1435953 4.39472320 0.000000e+00 0.000000e+00
## TYROBP 1.4462505 6.7012763 4.63355152 0.000000e+00 0.000000e+00
## FTL 12.0404304 44.5472409 3.69980469 0.000000e+00 0.000000e+00
## LGALS2 0.2856779 2.8015730 9.80675286 3.000210e-113 2.880202e-111
## S100A9 2.7875662 9.3448063 3.35231726 1.342336e-96 1.127563e-94
## S100A8 1.2979021 5.0767657 3.91151676 2.534042e-76 1.892085e-74
## FCN1 0.6709750 3.2285731 4.81176373 6.592657e-69 4.430266e-67
## LGALS1 1.3402498 4.8105813 3.58931696 3.380739e-65 2.065324e-63
## SAT1 1.2152254 4.3952305 3.61680268 1.351832e-61 7.570261e-60
## S100A6 1.9930820 5.7346055 2.87725523 1.492331e-52 7.714201e-51
## LTB 2.4241102 0.1524501 0.06288908 1.647321e-51 7.907139e-50
## TYMP 0.5048875 2.2602897 4.47681872 4.300544e-48 1.926644e-46
## GSTP1 0.5807220 2.4348709 4.19283390 1.940096e-47 8.148402e-46
## HLA-DRA 2.4502621 6.3635437 2.59708697 4.146877e-47 1.639236e-45
## S100A4 3.6110289 8.4502015 2.34010909 1.380971e-46 5.155626e-45
## AIF1 1.0270824 3.3575036 3.26897201 3.613774e-43 1.278135e-41
## RPS27A 7.0126593 2.3495410 0.33504280 5.909077e-40 1.985450e-38
## HLA-DRB1 1.0089787 3.1903187 3.16192880 1.331112e-39 4.259557e-38
## CTSS 0.9612186 3.0628646 3.18643914 4.248990e-39 1.297873e-37
## COTL1 1.4458032 3.9831465 2.75497146 2.510319e-37 7.334498e-36
## CCL3 0.1648204 1.0547455 6.39936369 4.113219e-37 1.151701e-35
## MS4A6A 0.1126701 0.8670669 7.69562594 1.880493e-36 5.054764e-35
types <- sub("(\\_\\d+)$","", colnames(sc@ndata))
genes <- head(rownames(dg)[dg$fc>1],10)
plotmarkergenes(sc,genes,samples=types)
plotmarkergenes(sc,genes,cl=c(2,6,7,8,10),samples=types,order.cells=TRUE)
fractDotPlot(sc, genes, cluster=c(2,6,7,8,10), zsc=TRUE)
A <- names(sc@cpart)[sc@cpart %in% c(2,4)]
B <- names(sc@cpart)[sc@cpart %in% c(3)]
x <- diffexpnb(getfdata(sc,n=c(A,B)), A=A, B=B )
plotdiffgenesnb(x,pthr=.05,lthr=.5,mthr=-1,Aname="Cl.2",Bname="Cl.3,5",show_names=TRUE,padj=TRUE)
ltr <- Ltree(sc)
ltr <- compentropy(ltr)
ltr <- projcells(ltr,cthr=5,nmode=FALSE,fr=TRUE)
ltr <- projback(ltr,pdishuf=100)
## pdishuffle: 1
pdishuffle: 2
pdishuffle: 3
pdishuffle: 4
pdishuffle: 5
pdishuffle: 6
pdishuffle: 7
pdishuffle: 8
pdishuffle: 9
pdishuffle: 10
pdishuffle: 11
pdishuffle: 12
pdishuffle: 13
pdishuffle: 14
pdishuffle: 15
pdishuffle: 16
pdishuffle: 17
pdishuffle: 18
pdishuffle: 19
pdishuffle: 20
pdishuffle: 21
pdishuffle: 22
pdishuffle: 23
pdishuffle: 24
pdishuffle: 25
pdishuffle: 26
pdishuffle: 27
pdishuffle: 28
pdishuffle: 29
pdishuffle: 30
pdishuffle: 31
pdishuffle: 32
pdishuffle: 33
pdishuffle: 34
pdishuffle: 35
pdishuffle: 36
pdishuffle: 37
pdishuffle: 38
pdishuffle: 39
pdishuffle: 40
pdishuffle: 41
pdishuffle: 42
pdishuffle: 43
pdishuffle: 44
pdishuffle: 45
pdishuffle: 46
pdishuffle: 47
pdishuffle: 48
pdishuffle: 49
pdishuffle: 50
pdishuffle: 51
pdishuffle: 52
pdishuffle: 53
pdishuffle: 54
pdishuffle: 55
pdishuffle: 56
pdishuffle: 57
pdishuffle: 58
pdishuffle: 59
pdishuffle: 60
pdishuffle: 61
pdishuffle: 62
pdishuffle: 63
pdishuffle: 64
pdishuffle: 65
pdishuffle: 66
pdishuffle: 67
pdishuffle: 68
pdishuffle: 69
pdishuffle: 70
pdishuffle: 71
pdishuffle: 72
pdishuffle: 73
pdishuffle: 74
pdishuffle: 75
pdishuffle: 76
pdishuffle: 77
pdishuffle: 78
pdishuffle: 79
pdishuffle: 80
pdishuffle: 81
pdishuffle: 82
pdishuffle: 83
pdishuffle: 84
pdishuffle: 85
pdishuffle: 86
pdishuffle: 87
pdishuffle: 88
pdishuffle: 89
pdishuffle: 90
pdishuffle: 91
pdishuffle: 92
pdishuffle: 93
pdishuffle: 94
pdishuffle: 95
pdishuffle: 96
pdishuffle: 97
pdishuffle: 98
pdishuffle: 99
pdishuffle: 100
pdishuffle:done.
ltr <- lineagegraph(ltr)
## Building tree: 1
Building tree: 2
Building tree: 3
Building tree: 4
Building tree: 5
Building tree: 6
Building tree: 7
Building tree: 8
Building tree: 9
Building tree: 10
Building tree: 11
Building tree: 12
Building tree: 13
Building tree: 14
Building tree: 15
Building tree: 16
Building tree: 17
Building tree: 18
Building tree: 19
Building tree: 20
Building tree: 21
Building tree: 22
Building tree: 23
Building tree: 24
Building tree: 25
Building tree: 26
Building tree: 27
Building tree: 28
Building tree: 29
Building tree: 30
Building tree: 31
Building tree: 32
Building tree: 33
Building tree: 34
Building tree: 35
Building tree: 36
Building tree: 37
Building tree: 38
Building tree: 39
Building tree: 40
Building tree: 41
Building tree: 42
Building tree: 43
Building tree: 44
Building tree: 45
Building tree: 46
Building tree: 47
Building tree: 48
Building tree: 49
Building tree: 50
Building tree: 51
Building tree: 52
Building tree: 53
Building tree: 54
Building tree: 55
Building tree: 56
Building tree: 57
Building tree: 58
Building tree: 59
Building tree: 60
Building tree: 61
Building tree: 62
Building tree: 63
Building tree: 64
Building tree: 65
Building tree: 66
Building tree: 67
Building tree: 68
Building tree: 69
Building tree: 70
Building tree: 71
Building tree: 72
Building tree: 73
Building tree: 74
Building tree: 75
Building tree: 76
Building tree: 77
Building tree: 78
Building tree: 79
Building tree: 80
Building tree: 81
Building tree: 82
Building tree: 83
Building tree: 84
Building tree: 85
Building tree: 86
Building tree: 87
Building tree: 88
Building tree: 89
Building tree: 90
Building tree: 91
Building tree: 92
Building tree: 93
Building tree: 94
Building tree: 95
Building tree: 96
Building tree: 97
Building tree: 98
Building tree: 99
Building tree: 100
Building tree: 101
Building tree: 102
Building tree: 103
Building tree: 104
Building tree: 105
Building tree: 106
Building tree: 107
Building tree: 108
Building tree: 109
Building tree: 110
Building tree: 111
Building tree: 112
Building tree: 113
Building tree: 114
Building tree: 115
Building tree: 116
Building tree: 117
Building tree: 118
Building tree: 119
Building tree: 120
Building tree: 121
Building tree: 122
Building tree: 123
Building tree: 124
Building tree: 125
Building tree: 126
Building tree: 127
Building tree: 128
Building tree: 129
Building tree: 130
Building tree: 131
Building tree: 132
Building tree: 133
Building tree: 134
Building tree: 135
Building tree: 136
Building tree: 137
Building tree: 138
Building tree: 139
Building tree: 140
Building tree: 141
Building tree: 142
Building tree: 143
Building tree: 144
Building tree: 145
Building tree: 146
Building tree: 147
Building tree: 148
Building tree: 149
Building tree: 150
Building tree: 151
Building tree: 152
Building tree: 153
Building tree: 154
Building tree: 155
Building tree: 156
Building tree: 157
Building tree: 158
Building tree: 159
Building tree: 160
Building tree: 161
Building tree: 162
Building tree: 163
Building tree: 164
Building tree: 165
Building tree: 166
Building tree: 167
Building tree: 168
Building tree: 169
Building tree: 170
Building tree: 171
Building tree: 172
Building tree: 173
Building tree: 174
Building tree: 175
Building tree: 176
Building tree: 177
Building tree: 178
Building tree: 179
Building tree: 180
Building tree: 181
Building tree: 182
Building tree: 183
Building tree: 184
Building tree: 185
Building tree: 186
Building tree: 187
Building tree: 188
Building tree: 189
Building tree: 190
Building tree: 191
Building tree: 192
Building tree: 193
Building tree: 194
Building tree: 195
Building tree: 196
Building tree: 197
Building tree: 198
Building tree: 199
Building tree: 200
Building tree: 201
Building tree: 202
Building tree: 203
Building tree: 204
Building tree: 205
Building tree: 206
Building tree: 207
Building tree: 208
Building tree: 209
Building tree: 210
Building tree: 211
Building tree: 212
Building tree: 213
Building tree: 214
Building tree: 215
Building tree: 216
Building tree: 217
Building tree: 218
Building tree: 219
Building tree: 220
Building tree: 221
Building tree: 222
Building tree: 223
Building tree: 224
Building tree: 225
Building tree: 226
Building tree: 227
Building tree: 228
Building tree: 229
Building tree: 230
Building tree: 231
Building tree: 232
Building tree: 233
Building tree: 234
Building tree: 235
Building tree: 236
Building tree: 237
Building tree: 238
Building tree: 239
Building tree: 240
Building tree: 241
Building tree: 242
Building tree: 243
Building tree: 244
Building tree: 245
Building tree: 246
Building tree: 247
Building tree: 248
Building tree: 249
Building tree: 250
Building tree: 251
Building tree: 252
Building tree: 253
Building tree: 254
Building tree: 255
Building tree: 256
Building tree: 257
Building tree: 258
Building tree: 259
Building tree: 260
Building tree: 261
Building tree: 262
Building tree: 263
Building tree: 264
Building tree: 265
Building tree: 266
Building tree: 267
Building tree: 268
Building tree: 269
Building tree: 270
Building tree: 271
Building tree: 272
Building tree: 273
Building tree: 274
Building tree: 275
Building tree: 276
Building tree: 277
Building tree: 278
Building tree: 279
Building tree: 280
Building tree: 281
Building tree: 282
Building tree: 283
Building tree: 284
Building tree: 285
Building tree: 286
Building tree: 287
Building tree: 288
Building tree: 289
Building tree: 290
Building tree: 291
Building tree: 292
Building tree: 293
Building tree: 294
Building tree: 295
Building tree: 296
Building tree: 297
Building tree: 298
Building tree: 299
Building tree: 300
Building tree: 301
Building tree: 302
Building tree: 303
Building tree: 304
Building tree: 305
Building tree: 306
Building tree: 307
Building tree: 308
Building tree: 309
Building tree: 310
Building tree: 311
Building tree: 312
Building tree: 313
Building tree: 314
Building tree: 315
Building tree: 316
Building tree: 317
Building tree: 318
Building tree: 319
Building tree: 320
Building tree: 321
Building tree: 322
Building tree: 323
Building tree: 324
Building tree: 325
Building tree: 326
Building tree: 327
Building tree: 328
Building tree: 329
Building tree: 330
Building tree: 331
Building tree: 332
Building tree: 333
Building tree: 334
Building tree: 335
Building tree: 336
Building tree: 337
Building tree: 338
Building tree: 339
Building tree: 340
Building tree: 341
Building tree: 342
Building tree: 343
Building tree: 344
Building tree: 345
Building tree: 346
Building tree: 347
Building tree: 348
Building tree: 349
Building tree: 350
Building tree: 351
Building tree: 352
Building tree: 353
Building tree: 354
Building tree: 355
Building tree: 356
Building tree: 357
Building tree: 358
Building tree: 359
Building tree: 360
Building tree: 361
Building tree: 362
Building tree: 363
Building tree: 364
Building tree: 365
Building tree: 366
Building tree: 367
Building tree: 368
Building tree: 369
Building tree: 370
Building tree: 371
Building tree: 372
Building tree: 373
Building tree: 374
Building tree: 375
Building tree: 376
Building tree: 377
Building tree: 378
Building tree: 379
Building tree: 380
Building tree: 381
Building tree: 382
Building tree: 383
Building tree: 384
Building tree: 385
Building tree: 386
Building tree: 387
Building tree: 388
Building tree: 389
Building tree: 390
Building tree: 391
Building tree: 392
Building tree: 393
Building tree: 394
Building tree: 395
Building tree: 396
Building tree: 397
Building tree: 398
Building tree: 399
Building tree: 400
Building tree: 401
Building tree: 402
Building tree: 403
Building tree: 404
Building tree: 405
Building tree: 406
Building tree: 407
Building tree: 408
Building tree: 409
Building tree: 410
Building tree: 411
Building tree: 412
Building tree: 413
Building tree: 414
Building tree: 415
Building tree: 416
Building tree: 417
Building tree: 418
Building tree: 419
Building tree: 420
Building tree: 421
Building tree: 422
Building tree: 423
Building tree: 424
Building tree: 425
Building tree: 426
Building tree: 427
Building tree: 428
Building tree: 429
Building tree: 430
Building tree: 431
Building tree: 432
Building tree: 433
Building tree: 434
Building tree: 435
Building tree: 436
Building tree: 437
Building tree: 438
Building tree: 439
Building tree: 440
Building tree: 441
Building tree: 442
Building tree: 443
Building tree: 444
Building tree: 445
Building tree: 446
Building tree: 447
Building tree: 448
Building tree: 449
Building tree: 450
Building tree: 451
Building tree: 452
Building tree: 453
Building tree: 454
Building tree: 455
Building tree: 456
Building tree: 457
Building tree: 458
Building tree: 459
Building tree: 460
Building tree: 461
Building tree: 462
Building tree: 463
Building tree: 464
Building tree: 465
Building tree: 466
Building tree: 467
Building tree: 468
Building tree: 469
Building tree: 470
Building tree: 471
Building tree: 472
Building tree: 473
Building tree: 474
Building tree: 475
Building tree: 476
Building tree: 477
Building tree: 478
Building tree: 479
Building tree: 480
Building tree: 481
Building tree: 482
Building tree: 483
Building tree: 484
Building tree: 485
Building tree: 486
Building tree: 487
Building tree: 488
Building tree: 489
Building tree: 490
Building tree: 491
Building tree: 492
Building tree: 493
Building tree: 494
Building tree: 495
Building tree: 496
Building tree: 497
Building tree: 498
Building tree: 499
Building tree: 500
Building tree: 501
Building tree: 502
Building tree: 503
Building tree: 504
Building tree: 505
Building tree: 506
Building tree: 507
Building tree: 508
Building tree: 509
Building tree: 510
Building tree: 511
Building tree: 512
Building tree: 513
Building tree: 514
Building tree: 515
Building tree: 516
Building tree: 517
Building tree: 518
Building tree: 519
Building tree: 520
Building tree: 521
Building tree: 522
Building tree: 523
Building tree: 524
Building tree: 525
Building tree: 526
Building tree: 527
Building tree: 528
Building tree: 529
Building tree: 530
Building tree: 531
Building tree: 532
Building tree: 533
Building tree: 534
Building tree: 535
Building tree: 536
Building tree: 537
Building tree: 538
Building tree: 539
Building tree: 540
Building tree: 541
Building tree: 542
Building tree: 543
Building tree: 544
Building tree: 545
Building tree: 546
Building tree: 547
Building tree: 548
Building tree: 549
Building tree: 550
Building tree: 551
Building tree: 552
Building tree: 553
Building tree: 554
Building tree: 555
Building tree: 556
Building tree: 557
Building tree: 558
Building tree: 559
Building tree: 560
Building tree: 561
Building tree: 562
Building tree: 563
Building tree: 564
Building tree: 565
Building tree: 566
Building tree: 567
Building tree: 568
Building tree: 569
Building tree: 570
Building tree: 571
Building tree: 572
Building tree: 573
Building tree: 574
Building tree: 575
Building tree: 576
Building tree: 577
Building tree: 578
Building tree: 579
Building tree: 580
Building tree: 581
Building tree: 582
Building tree: 583
Building tree: 584
Building tree: 585
Building tree: 586
Building tree: 587
Building tree: 588
Building tree: 589
Building tree: 590
Building tree: 591
Building tree: 592
Building tree: 593
Building tree: 594
Building tree: 595
Building tree: 596
Building tree: 597
Building tree: 598
Building tree: 599
Building tree: 600
Building tree: 601
Building tree: 602
Building tree: 603
Building tree: 604
Building tree: 605
Building tree: 606
Building tree: 607
Building tree: 608
Building tree: 609
Building tree: 610
Building tree: 611
Building tree: 612
Building tree: 613
Building tree: 614
Building tree: 615
Building tree: 616
Building tree: 617
Building tree: 618
Building tree: 619
Building tree: 620
Building tree: 621
Building tree: 622
Building tree: 623
Building tree: 624
Building tree: 625
Building tree: 626
Building tree: 627
Building tree: 628
Building tree: 629
Building tree: 630
Building tree: 631
Building tree: 632
Building tree: 633
Building tree: 634
Building tree: 635
Building tree: 636
Building tree: 637
Building tree: 638
Building tree: 639
Building tree: 640
Building tree: 641
Building tree: 642
Building tree: 643
Building tree: 644
Building tree: 645
Building tree: 646
Building tree: 647
Building tree: 648
Building tree: 649
Building tree: 650
Building tree: 651
Building tree: 652
Building tree: 653
Building tree: 654
Building tree: 655
Building tree: 656
Building tree: 657
Building tree: 658
Building tree: 659
Building tree: 660
Building tree: 661
Building tree: 662
Building tree: 663
Building tree: 664
Building tree: 665
Building tree: 666
Building tree: 667
Building tree: 668
Building tree: 669
Building tree: 670
Building tree: 671
Building tree: 672
Building tree: 673
Building tree: 674
Building tree: 675
Building tree: 676
Building tree: 677
Building tree: 678
Building tree: 679
Building tree: 680
Building tree: 681
Building tree: 682
Building tree: 683
Building tree: 684
Building tree: 685
Building tree: 686
Building tree: 687
Building tree: 688
Building tree: 689
Building tree: 690
Building tree: 691
Building tree: 692
Building tree: 693
Building tree: 694
Building tree: 695
Building tree: 696
Building tree: 697
Building tree: 698
Building tree: 699
Building tree: 700
Building tree: 701
Building tree: 702
Building tree: 703
Building tree: 704
Building tree: 705
Building tree: 706
Building tree: 707
Building tree: 708
Building tree: 709
Building tree: 710
Building tree: 711
Building tree: 712
Building tree: 713
Building tree: 714
Building tree: 715
Building tree: 716
Building tree: 717
Building tree: 718
Building tree: 719
Building tree: 720
Building tree: 721
Building tree: 722
Building tree: 723
Building tree: 724
Building tree: 725
Building tree: 726
Building tree: 727
Building tree: 728
Building tree: 729
Building tree: 730
Building tree: 731
Building tree: 732
Building tree: 733
Building tree: 734
Building tree: 735
Building tree: 736
Building tree: 737
Building tree: 738
Building tree: 739
Building tree: 740
Building tree: 741
Building tree: 742
Building tree: 743
Building tree: 744
Building tree: 745
Building tree: 746
Building tree: 747
Building tree: 748
Building tree: 749
Building tree: 750
Building tree: 751
Building tree: 752
Building tree: 753
Building tree: 754
Building tree: 755
Building tree: 756
Building tree: 757
Building tree: 758
Building tree: 759
Building tree: 760
Building tree: 761
Building tree: 762
Building tree: 763
Building tree: 764
Building tree: 765
Building tree: 766
Building tree: 767
Building tree: 768
Building tree: 769
Building tree: 770
Building tree: 771
Building tree: 772
Building tree: 773
Building tree: 774
Building tree: 775
Building tree: 776
Building tree: 777
Building tree: 778
Building tree: 779
Building tree: 780
Building tree: 781
Building tree: 782
Building tree: 783
Building tree: 784
Building tree: 785
Building tree: 786
Building tree: 787
Building tree: 788
Building tree: 789
Building tree: 790
Building tree: 791
Building tree: 792
Building tree: 793
Building tree: 794
Building tree: 795
Building tree: 796
Building tree: 797
Building tree: 798
Building tree: 799
Building tree: 800
Building tree: 801
Building tree: 802
Building tree: 803
Building tree: 804
Building tree: 805
Building tree: 806
Building tree: 807
Building tree: 808
Building tree: 809
Building tree: 810
Building tree: 811
Building tree: 812
Building tree: 813
Building tree: 814
Building tree: 815
Building tree: 816
Building tree: 817
Building tree: 818
Building tree: 819
Building tree: 820
Building tree: 821
Building tree: 822
Building tree: 823
Building tree: 824
Building tree: 825
Building tree: 826
Building tree: 827
Building tree: 828
Building tree: 829
Building tree: 830
Building tree: 831
Building tree: 832
Building tree: 833
Building tree: 834
Building tree: 835
Building tree: 836
Building tree: 837
Building tree: 838
Building tree: 839
Building tree: 840
Building tree: 841
Building tree: 842
Building tree: 843
Building tree: 844
Building tree: 845
Building tree: 846
Building tree: 847
Building tree: 848
Building tree: 849
Building tree: 850
Building tree: 851
Building tree: 852
Building tree: 853
Building tree: 854
Building tree: 855
Building tree: 856
Building tree: 857
Building tree: 858
Building tree: 859
Building tree: 860
Building tree: 861
Building tree: 862
Building tree: 863
Building tree: 864
Building tree: 865
Building tree: 866
Building tree: 867
Building tree: 868
Building tree: 869
Building tree: 870
Building tree: 871
Building tree: 872
Building tree: 873
Building tree: 874
Building tree: 875
Building tree: 876
Building tree: 877
Building tree: 878
Building tree: 879
Building tree: 880
Building tree: 881
Building tree: 882
Building tree: 883
Building tree: 884
Building tree: 885
Building tree: 886
Building tree: 887
Building tree: 888
Building tree: 889
Building tree: 890
Building tree: 891
Building tree: 892
Building tree: 893
Building tree: 894
Building tree: 895
Building tree: 896
Building tree: 897
Building tree: 898
Building tree: 899
Building tree: 900
Building tree: 901
Building tree: 902
Building tree: 903
Building tree: 904
Building tree: 905
Building tree: 906
Building tree: 907
Building tree: 908
Building tree: 909
Building tree: 910
Building tree: 911
Building tree: 912
Building tree: 913
Building tree: 914
Building tree: 915
Building tree: 916
Building tree: 917
Building tree: 918
Building tree: 919
Building tree: 920
Building tree: 921
Building tree: 922
Building tree: 923
Building tree: 924
Building tree: 925
Building tree: 926
Building tree: 927
Building tree: 928
Building tree: 929
Building tree: 930
Building tree: 931
Building tree: 932
Building tree: 933
Building tree: 934
Building tree: 935
Building tree: 936
Building tree: 937
Building tree: 938
Building tree: 939
Building tree: 940
Building tree: 941
Building tree: 942
Building tree: 943
Building tree: 944
Building tree: 945
Building tree: 946
Building tree: 947
Building tree: 948
Building tree: 949
Building tree: 950
Building tree: 951
Building tree: 952
Building tree: 953
Building tree: 954
Building tree: 955
Building tree: 956
Building tree: 957
Building tree: 958
Building tree: 959
Building tree: 960
Building tree: 961
Building tree: 962
Building tree: 963
Building tree: 964
Building tree: 965
Building tree: 966
Building tree: 967
Building tree: 968
Building tree: 969
Building tree: 970
Building tree: 971
Building tree: 972
Building tree: 973
Building tree: 974
Building tree: 975
Building tree: 976
Building tree: 977
Building tree: 978
Building tree: 979
Building tree: 980
Building tree: 981
Building tree: 982
Building tree: 983
Building tree: 984
Building tree: 985
Building tree: 986
Building tree: 987
Building tree: 988
Building tree: 989
Building tree: 990
Building tree: 991
Building tree: 992
Building tree: 993
Building tree: 994
Building tree: 995
Building tree: 996
Building tree: 997
Building tree: 998
Building tree: 999
Building tree: 1000
Building tree: 1001
Building tree: 1002
Building tree: 1003
Building tree: 1004
Building tree: 1005
Building tree: 1006
Building tree: 1007
Building tree: 1008
Building tree: 1009
Building tree: 1010
Building tree: 1011
Building tree: 1012
Building tree: 1013
Building tree: 1014
Building tree: 1015
Building tree: 1016
Building tree: 1017
Building tree: 1018
Building tree: 1019
Building tree: 1020
Building tree: 1021
Building tree: 1022
Building tree: 1023
Building tree: 1024
Building tree: 1025
Building tree: 1026
Building tree: 1027
Building tree: 1028
Building tree: 1029
Building tree: 1030
Building tree: 1031
Building tree: 1032
Building tree: 1033
Building tree: 1034
Building tree: 1035
Building tree: 1036
Building tree: 1037
Building tree: 1038
Building tree: 1039
Building tree: 1040
Building tree: 1041
Building tree: 1042
Building tree: 1043
Building tree: 1044
Building tree: 1045
Building tree: 1046
Building tree: 1047
Building tree: 1048
Building tree: 1049
Building tree: 1050
Building tree: 1051
Building tree: 1052
Building tree: 1053
Building tree: 1054
Building tree: 1055
Building tree: 1056
Building tree: 1057
Building tree: 1058
Building tree: 1059
Building tree: 1060
Building tree: 1061
Building tree: 1062
Building tree: 1063
Building tree: 1064
Building tree: 1065
Building tree: 1066
Building tree: 1067
Building tree: 1068
Building tree: 1069
Building tree: 1070
Building tree: 1071
Building tree: 1072
Building tree: 1073
Building tree: 1074
Building tree: 1075
Building tree: 1076
Building tree: 1077
Building tree: 1078
Building tree: 1079
Building tree: 1080
Building tree: 1081
Building tree: 1082
Building tree: 1083
Building tree: 1084
Building tree: 1085
Building tree: 1086
Building tree: 1087
Building tree: 1088
Building tree: 1089
Building tree: 1090
Building tree: 1091
Building tree: 1092
Building tree: 1093
Building tree: 1094
Building tree: 1095
Building tree: 1096
Building tree: 1097
Building tree: 1098
Building tree: 1099
Building tree: 1100
Building tree: 1101
Building tree: 1102
Building tree: 1103
Building tree: 1104
Building tree: 1105
Building tree: 1106
Building tree: 1107
Building tree: 1108
Building tree: 1109
Building tree: 1110
Building tree: 1111
Building tree: 1112
Building tree: 1113
Building tree: 1114
Building tree: 1115
Building tree: 1116
Building tree: 1117
Building tree: 1118
Building tree: 1119
Building tree: 1120
Building tree: 1121
Building tree: 1122
Building tree: 1123
Building tree: 1124
Building tree: 1125
Building tree: 1126
Building tree: 1127
Building tree: 1128
Building tree: 1129
Building tree: 1130
Building tree: 1131
Building tree: 1132
Building tree: 1133
Building tree: 1134
Building tree: 1135
Building tree: 1136
Building tree: 1137
Building tree: 1138
Building tree: 1139
Building tree: 1140
Building tree: 1141
Building tree: 1142
Building tree: 1143
Building tree: 1144
Building tree: 1145
Building tree: 1146
Building tree: 1147
Building tree: 1148
Building tree: 1149
Building tree: 1150
Building tree: 1151
Building tree: 1152
Building tree: 1153
Building tree: 1154
Building tree: 1155
Building tree: 1156
Building tree: 1157
Building tree: 1158
Building tree: 1159
Building tree: 1160
Building tree: 1161
Building tree: 1162
Building tree: 1163
Building tree: 1164
Building tree: 1165
Building tree: 1166
Building tree: 1167
Building tree: 1168
Building tree: 1169
Building tree: 1170
Building tree: 1171
Building tree: 1172
Building tree: 1173
Building tree: 1174
Building tree: 1175
Building tree: 1176
Building tree: 1177
Building tree: 1178
Building tree: 1179
Building tree: 1180
Building tree: 1181
Building tree: 1182
Building tree: 1183
Building tree: 1184
Building tree: 1185
Building tree: 1186
Building tree: 1187
Building tree: 1188
Building tree: 1189
Building tree: 1190
Building tree: 1191
Building tree: 1192
Building tree: 1193
Building tree: 1194
Building tree: 1195
Building tree: 1196
Building tree: 1197
Building tree: 1198
Building tree: 1199
Building tree: 1200
Building tree: 1201
Building tree: 1202
Building tree: 1203
Building tree: 1204
Building tree: 1205
Building tree: 1206
Building tree: 1207
Building tree: 1208
Building tree: 1209
Building tree: 1210
Building tree: 1211
Building tree: 1212
Building tree: 1213
Building tree: 1214
Building tree: 1215
Building tree: 1216
Building tree: 1217
Building tree: 1218
Building tree: 1219
Building tree: 1220
Building tree: 1221
Building tree: 1222
Building tree: 1223
Building tree: 1224
Building tree: 1225
Building tree: 1226
Building tree: 1227
Building tree: 1228
Building tree: 1229
Building tree: 1230
Building tree: 1231
Building tree: 1232
Building tree: 1233
Building tree: 1234
Building tree: 1235
Building tree: 1236
Building tree: 1237
Building tree: 1238
Building tree: 1239
Building tree: 1240
Building tree: 1241
Building tree: 1242
Building tree: 1243
Building tree: 1244
Building tree: 1245
Building tree: 1246
Building tree: 1247
Building tree: 1248
Building tree: 1249
Building tree: 1250
Building tree: 1251
Building tree: 1252
Building tree: 1253
Building tree: 1254
Building tree: 1255
Building tree: 1256
Building tree: 1257
Building tree: 1258
Building tree: 1259
Building tree: 1260
Building tree: 1261
Building tree: 1262
Building tree: 1263
Building tree: 1264
Building tree: 1265
Building tree: 1266
Building tree: 1267
Building tree: 1268
Building tree: 1269
Building tree: 1270
Building tree: 1271
Building tree: 1272
Building tree: 1273
Building tree: 1274
Building tree: 1275
Building tree: 1276
Building tree: 1277
Building tree: 1278
Building tree: 1279
Building tree: 1280
Building tree: 1281
Building tree: 1282
Building tree: 1283
Building tree: 1284
Building tree: 1285
Building tree: 1286
Building tree: 1287
Building tree: 1288
Building tree: 1289
Building tree: 1290
Building tree: 1291
Building tree: 1292
Building tree: 1293
Building tree: 1294
Building tree: 1295
Building tree: 1296
Building tree: 1297
Building tree: 1298
Building tree: 1299
Building tree: 1300
Building tree: 1301
Building tree: 1302
Building tree: 1303
Building tree: 1304
Building tree: 1305
Building tree: 1306
Building tree: 1307
Building tree: 1308
Building tree: 1309
Building tree: 1310
Building tree: 1311
Building tree: 1312
Building tree: 1313
Building tree: 1314
Building tree: 1315
Building tree: 1316
Building tree: 1317
Building tree: 1318
Building tree: 1319
Building tree: 1320
Building tree: 1321
Building tree: 1322
Building tree: 1323
Building tree: 1324
Building tree: 1325
Building tree: 1326
Building tree: 1327
Building tree: 1328
Building tree: 1329
Building tree: 1330
Building tree: 1331
Building tree: 1332
Building tree: 1333
Building tree: 1334
Building tree: 1335
Building tree: 1336
Building tree: 1337
Building tree: 1338
Building tree: 1339
Building tree: 1340
Building tree: 1341
Building tree: 1342
Building tree: 1343
Building tree: 1344
Building tree: 1345
Building tree: 1346
Building tree: 1347
Building tree: 1348
Building tree: 1349
Building tree: 1350
Building tree: 1351
Building tree: 1352
Building tree: 1353
Building tree: 1354
Building tree: 1355
Building tree: 1356
Building tree: 1357
Building tree: 1358
Building tree: 1359
Building tree: 1360
Building tree: 1361
Building tree: 1362
Building tree: 1363
Building tree: 1364
Building tree: 1365
Building tree: 1366
Building tree: 1367
Building tree: 1368
Building tree: 1369
Building tree: 1370
Building tree: 1371
Building tree: 1372
Building tree: 1373
Building tree: 1374
Building tree: 1375
Building tree: 1376
Building tree: 1377
Building tree: 1378
Building tree: 1379
Building tree: 1380
Building tree: 1381
Building tree: 1382
Building tree: 1383
Building tree: 1384
Building tree: 1385
Building tree: 1386
Building tree: 1387
Building tree: 1388
Building tree: 1389
Building tree: 1390
Building tree: 1391
Building tree: 1392
Building tree: 1393
Building tree: 1394
Building tree: 1395
Building tree: 1396
Building tree: 1397
Building tree: 1398
Building tree: 1399
Building tree: 1400
Building tree: 1401
Building tree: 1402
Building tree: 1403
Building tree: 1404
Building tree: 1405
Building tree: 1406
Building tree: 1407
Building tree: 1408
Building tree: 1409
Building tree: 1410
Building tree: 1411
Building tree: 1412
Building tree: 1413
Building tree: 1414
Building tree: 1415
Building tree: 1416
Building tree: 1417
Building tree: 1418
Building tree: 1419
Building tree: 1420
Building tree: 1421
Building tree: 1422
Building tree: 1423
Building tree: 1424
Building tree: 1425
Building tree: 1426
Building tree: 1427
Building tree: 1428
Building tree: 1429
Building tree: 1430
Building tree: 1431
Building tree: 1432
Building tree: 1433
Building tree: 1434
Building tree: 1435
Building tree: 1436
Building tree: 1437
Building tree: 1438
Building tree: 1439
Building tree: 1440
Building tree: 1441
Building tree: 1442
Building tree: 1443
Building tree: 1444
Building tree: 1445
Building tree: 1446
Building tree: 1447
Building tree: 1448
Building tree: 1449
Building tree: 1450
Building tree: 1451
Building tree: 1452
Building tree: 1453
Building tree: 1454
Building tree: 1455
Building tree: 1456
Building tree: 1457
Building tree: 1458
Building tree: 1459
Building tree: 1460
Building tree: 1461
Building tree: 1462
Building tree: 1463
Building tree: 1464
Building tree: 1465
Building tree: 1466
Building tree: 1467
Building tree: 1468
Building tree: 1469
Building tree: 1470
Building tree: 1471
Building tree: 1472
Building tree: 1473
Building tree: 1474
Building tree: 1475
Building tree: 1476
Building tree: 1477
Building tree: 1478
Building tree: 1479
Building tree: 1480
Building tree: 1481
Building tree: 1482
Building tree: 1483
Building tree: 1484
Building tree: 1485
Building tree: 1486
Building tree: 1487
Building tree: 1488
Building tree: 1489
Building tree: 1490
Building tree: 1491
Building tree: 1492
Building tree: 1493
Building tree: 1494
Building tree: 1495
Building tree: 1496
Building tree: 1497
Building tree: 1498
Building tree: 1499
Building tree: 1500
Building tree: 1501
Building tree: 1502
Building tree: 1503
Building tree: 1504
Building tree: 1505
Building tree: 1506
Building tree: 1507
Building tree: 1508
Building tree: 1509
Building tree: 1510
Building tree: 1511
Building tree: 1512
Building tree: 1513
Building tree: 1514
Building tree: 1515
Building tree: 1516
Building tree: 1517
Building tree: 1518
Building tree: 1519
Building tree: 1520
Building tree: 1521
Building tree: 1522
Building tree: 1523
Building tree: 1524
Building tree: 1525
Building tree: 1526
Building tree: 1527
Building tree: 1528
Building tree: 1529
Building tree: 1530
Building tree: 1531
Building tree: 1532
Building tree: 1533
Building tree: 1534
Building tree: 1535
Building tree: 1536
Building tree: 1537
Building tree: 1538
Building tree: 1539
Building tree: 1540
Building tree: 1541
Building tree: 1542
Building tree: 1543
Building tree: 1544
Building tree: 1545
Building tree: 1546
Building tree: 1547
Building tree: 1548
Building tree: 1549
Building tree: 1550
Building tree: 1551
Building tree: 1552
Building tree: 1553
Building tree: 1554
Building tree: 1555
Building tree: 1556
Building tree: 1557
Building tree: 1558
Building tree: 1559
Building tree: 1560
Building tree: 1561
Building tree: 1562
Building tree: 1563
Building tree: 1564
Building tree: 1565
Building tree: 1566
Building tree: 1567
Building tree: 1568
Building tree: 1569
Building tree: 1570
Building tree: 1571
Building tree: 1572
Building tree: 1573
Building tree: 1574
Building tree: 1575
Building tree: 1576
Building tree: 1577
Building tree: 1578
Building tree: 1579
Building tree: 1580
Building tree: 1581
Building tree: 1582
Building tree: 1583
Building tree: 1584
Building tree: 1585
Building tree: 1586
Building tree: 1587
Building tree: 1588
Building tree: 1589
Building tree: 1590
Building tree: 1591
Building tree: 1592
Building tree: 1593
Building tree: 1594
Building tree: 1595
Building tree: done.
ltr <- comppvalue(ltr,pthr=0.1)
plotgraph(ltr,scthr=0.2,showCells=FALSE,showMap=TRUE)
To predict the stem cell, the StemID score can be computed and
visualized:
x <- compscore(ltr,scthr=0.2)
plotdistanceratio(ltr)
plotspantree(ltr)
The cell projections onto the links can be directly compared with this
minimal spanning tree:
plotspantree(ltr,projections=TRUE)
plotlinkscore(ltr)
projenrichment(ltr)
x <- getproj(ltr,i=3)
sessionInfo()
## R version 4.0.2 (2020-06-22)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 20.04.4 LTS
##
## Matrix products: default
## BLAS: /mnt/nectar_volume/software/apps/R/4.0.2/lib/R/lib/libRblas.so
## LAPACK: /mnt/nectar_volume/software/apps/R/4.0.2/lib/R/lib/libRlapack.so
##
## locale:
## [1] LC_CTYPE=en_AU.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_AU.UTF-8 LC_COLLATE=en_AU.UTF-8
## [5] LC_MONETARY=en_AU.UTF-8 LC_MESSAGES=en_AU.UTF-8
## [7] LC_PAPER=en_AU.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_AU.UTF-8 LC_IDENTIFICATION=C
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] ggplot2_3.3.6 dplyr_1.0.8 SeuratObject_4.0.4 Seurat_4.1.0
## [5] RaceID_0.2.3
##
## loaded via a namespace (and not attached):
## [1] plyr_1.8.6 igraph_1.2.11 lazyeval_0.2.2
## [4] splines_4.0.2 listenv_0.8.0 scattermore_0.7
## [7] lle_1.1 digest_0.6.29 htmltools_0.5.2
## [10] fansi_1.0.2 magrittr_2.0.2 tensor_1.5
## [13] cluster_2.1.3 ROCR_1.0-11 globals_0.14.0
## [16] matrixStats_0.61.0 askpass_1.1 spatstat.sparse_2.1-0
## [19] princurve_2.1.6 colorspace_2.0-2 ggrepel_0.9.1
## [22] xfun_0.30 crayon_1.4.2 jsonlite_1.7.3
## [25] spatstat.data_2.1-2 survival_3.1-12 zoo_1.8-9
## [28] glue_1.6.1 polyclip_1.10-0 gtable_0.3.0
## [31] leiden_0.3.9 coop_0.6-3 kernlab_0.9-29
## [34] future.apply_1.8.1 DEoptimR_1.0-11 prabclus_2.3-2
## [37] abind_1.4-5 scales_1.1.1 pheatmap_1.0.12
## [40] DBI_1.1.2 som_0.3-5.1 miniUI_0.1.1.1
## [43] Rcpp_1.0.8 viridisLite_0.4.0 xtable_1.8-4
## [46] reticulate_1.24 spatstat.core_2.3-2 mclust_5.4.9
## [49] stats4_4.0.2 umap_0.2.7.0 propr_4.2.6
## [52] htmlwidgets_1.5.4 httr_1.4.2 FNN_1.1.3
## [55] RColorBrewer_1.1-2 fpc_2.2-9 modeltools_0.2-23
## [58] ellipsis_0.3.2 ica_1.0-2 farver_2.1.0
## [61] flexmix_2.3-17 pkgconfig_2.0.3 nnet_7.3-14
## [64] sass_0.4.0 uwot_0.1.11 deldir_1.0-6
## [67] locfit_1.5-9.4 utf8_1.2.2 labeling_0.4.2
## [70] tidyselect_1.1.2 rlang_1.0.2 reshape2_1.4.4
## [73] later_1.3.0 munsell_0.5.0 tools_4.0.2
## [76] cli_3.2.0 generics_0.1.2 ggridges_0.5.3
## [79] evaluate_0.15 stringr_1.4.0 fastmap_1.1.0
## [82] yaml_2.2.2 goftest_1.2-3 knitr_1.39
## [85] fitdistrplus_1.1-6 robustbase_0.95-0 purrr_0.3.4
## [88] randomForest_4.6-14 RANN_2.6.1 pbapply_1.5-0
## [91] future_1.23.0 nlme_3.1-148 mime_0.12
## [94] compiler_4.0.2 rstudioapi_0.13 plotly_4.10.0
## [97] png_0.1-7 spatstat.utils_2.3-0 tibble_3.1.6
## [100] bslib_0.3.1 stringi_1.7.6 highr_0.9
## [103] RSpectra_0.16-0 lattice_0.20-41 Matrix_1.4-0
## [106] vegan_2.6-2 permute_0.9-7 FateID_0.2.1
## [109] vctrs_0.3.8 pillar_1.7.0 lifecycle_1.0.1
## [112] spatstat.geom_2.3-1 lmtest_0.9-39 jquerylib_0.1.4
## [115] RcppAnnoy_0.0.19 data.table_1.14.2 cowplot_1.1.1
## [118] irlba_2.3.5 httpuv_1.6.5 patchwork_1.1.1
## [121] R6_2.5.1 promises_1.2.0.1 KernSmooth_2.23-17
## [124] gridExtra_2.3 parallelly_1.30.0 codetools_0.2-16
## [127] MASS_7.3-51.6 assertthat_0.2.1 openssl_1.4.6
## [130] withr_2.4.3 sctransform_0.3.3 diptest_0.76-0
## [133] mgcv_1.8-31 parallel_4.0.2 quadprog_1.5-8
## [136] grid_4.0.2 rpart_4.1-15 class_7.3-17
## [139] tidyr_1.2.0 rmarkdown_2.14 Rtsne_0.15
## [142] scatterplot3d_0.3-41 shiny_1.7.1 snowfall_1.84-6.1