Description Usage Arguments Details Value Author(s) References Examples
Perform model-based clustering on expression values
1 | exprmclust(data, clusternum = 2:9, modelNames = "VVV", reduce = T)
|
data |
The raw single_cell data, which is a numeric matrix or data.frame. Rows represent genes/features and columns represent single cells. |
clusternum |
An integer vector specifying all possible cluster numbers. The best cluster number will be picked using BIC. The minimum value should be two other |
modelNames |
model to be used in model-based clustering. By default "ellipsoidal, varying volume, shape, and orientation" is used. |
reduce |
Whether to perform the PCA on the expression data. |
By default, this function first uses principal component analysis (PCA) to reduce dimensionality of original data. It then performs model-based clustering on the transformed expression values. A minimum-spanning-tree is constructed to link the cluster centers. The clustering results will be used for TSCAN ordering.
if more than one cluster detected, a list containing
pcareduceres Numeric matrix containing the transformed expression values after PCA.
MSTtree igraph object which is the result of constructing MST.
clusterid A named vector specifying which cluster the cells belong to.
clucenter Numeric matrix of the cluster centers.
if only one cluster detected, a list containing
pcareduceres Numeric matrix containing the transformed expression values after PCA.
Zhicheng Ji, Hongkai Ji <zji4@zji4.edu>
Fraley, C., & Raftery, A. E. (2002). Model-based clustering, discriminant analysis, and density estimation. Journal of the American Statistical Association, 97(458), 611-631.
1 2 3 | data(lpsdata)
procdata <- preprocess(lpsdata)
exprmclust(procdata)
|
$pcareduceres
PC1 PC2 PC3 PC4 PC5
Unstimulated_S15 -4.92638297 -8.8529185 3.22733252 -2.530967511 -2.32018858
Unstimulated_S16 -6.53424105 2.7790245 -0.07747889 0.828806525 1.38991796
Unstimulated_S17 -4.37822709 -7.1209402 3.75669380 -1.880575370 -0.20894592
Unstimulated_S18 17.81566572 -4.9989393 5.89433717 0.498551168 0.37654921
Unstimulated_S23 4.63958402 -7.0623736 -0.16541202 1.884503386 -4.20403502
Unstimulated_S24 -5.11519374 -7.8997135 3.39796650 -3.705876957 -0.72953511
Unstimulated_S26 -5.61797356 -7.0595907 0.59760040 -0.830618710 -4.50418081
Unstimulated_S27 -16.00062653 -3.1312633 -0.71870458 1.113940481 -2.63587791
Unstimulated_S28 -7.66518422 -8.1902394 1.26091879 -2.130027451 -5.02037825
Unstimulated_S2 -7.01346590 -5.4934888 -1.37784833 1.466572277 4.56611985
Unstimulated_S31 -2.96507054 -3.5346277 -0.64206453 5.027971580 4.34014341
Unstimulated_S32 -2.61022408 -4.4616316 -1.10533938 0.434571738 4.00092481
Unstimulated_S34 -7.21965614 -4.0424538 -4.26155077 2.719748564 1.25508809
Unstimulated_S37 16.95835092 -5.4442503 5.10034345 1.603287505 -0.74659474
Unstimulated_S38 -11.22916670 -6.1224529 -1.21729258 -0.267079854 -1.73711700
Unstimulated_S39 -8.59346505 -6.6919448 0.07142350 0.549971696 1.19148416
Unstimulated_S3 -9.23322744 -3.5455576 -2.70894287 3.510278152 0.33302224
Unstimulated_S41 -2.38737204 -4.3388880 0.33955341 -2.645652053 10.18993650
Unstimulated_S46 -9.94907597 -4.8907295 -2.99509140 2.355224934 -1.41274142
Unstimulated_S47 -0.31980414 -5.5103084 -2.42844969 -0.440886113 1.38752524
Unstimulated_S48 -4.85597695 -6.6271405 -0.78445250 0.464599295 -6.21643850
Unstimulated_S49 -8.52832342 -7.0757567 -0.23212922 -1.232910859 -0.04610125
Unstimulated_S4 -1.67313748 -4.0589471 -1.76452598 0.711142629 -0.32216581
Unstimulated_S51 -6.05022337 -5.2222314 -9.85843013 0.032263153 1.44143540
Unstimulated_S52 -6.46344442 -5.8367555 0.48176676 3.486264194 0.60760780
Unstimulated_S53 -5.63582158 -3.3383028 -0.05760304 1.476669275 2.06340185
Unstimulated_S54 -0.84628345 -6.1867255 -5.47061553 0.209137854 -2.77552085
Unstimulated_S55 11.32263613 -7.0996464 3.91329851 0.307710659 0.50135403
Unstimulated_S5 0.35007434 -2.5303971 -0.69994892 -1.999123649 4.25006301
Unstimulated_S60 -8.38574750 -8.0671125 1.35441912 -1.448207726 -0.12567796
Unstimulated_S61 -8.97361103 -5.9023228 0.67988538 4.960491262 -1.11922858
Unstimulated_S62 -0.54632004 -7.5029215 2.73790070 -2.444717075 -2.65743490
Unstimulated_S63 -9.00060740 -4.8634030 -0.71524850 2.479860400 0.85320474
Unstimulated_S66 -6.96045778 -7.5509841 -0.45159996 0.734743680 -3.27054592
Unstimulated_S6 -3.71662075 -8.8188517 3.01350227 -1.390360390 -1.58972241
Unstimulated_S70 -6.31449548 -4.7744833 -2.53975012 0.903786856 0.39503315
Unstimulated_S71 -5.20887527 -6.7001860 -0.54456262 1.400067578 2.55685691
Unstimulated_S72 0.06053991 -3.8589266 -1.38365411 6.023712402 3.16191871
Unstimulated_S74 -1.21926738 -7.3115827 1.94698254 -1.931832975 -0.81712112
Unstimulated_S75 -2.42607251 -7.3594098 -0.78034205 -0.386519407 2.42879469
Unstimulated_S76 -1.14472728 -5.7108653 -3.92242548 1.964556081 4.86660657
Unstimulated_S77 -1.91827615 -4.9925797 2.10799225 1.867632369 5.41470502
Unstimulated_S78 -9.11918319 -6.6336481 0.50797940 -3.363948549 0.57338160
Unstimulated_S79 4.71193963 -7.0725255 0.84396847 -3.113547137 2.84111610
Unstimulated_S7 15.02717140 -3.1812831 -20.37863291 -0.195088580 -3.81874672
Unstimulated_S80 -6.16676261 -4.3178415 -4.46734816 4.330139633 3.38227484
Unstimulated_S81 14.68838230 -4.8070524 5.95065699 3.154664077 2.03882286
Unstimulated_S82 -11.18519764 -5.1538116 1.38012271 -0.546715656 -2.40312946
Unstimulated_S83 7.86295942 -8.4072488 1.72028689 -2.583946521 0.95359758
Unstimulated_S84 12.45244947 -7.3521724 6.14909952 -0.644095292 3.28545509
Unstimulated_S85 17.77990759 -5.6904395 6.05606450 0.429767132 -0.08124407
Unstimulated_S8 0.48902448 -3.1708219 -8.27172906 -2.165200456 5.66737518
Unstimulated_S90 -7.04643419 -0.6512254 -0.51337113 0.004834302 3.04673465
Unstimulated_S92 -1.76995672 -5.1857664 0.50378468 6.754768765 1.37016863
Unstimulated_S93 -6.32462213 -6.8513650 -2.26678468 4.944123477 -3.39052517
Unstimulated_S95 15.90686277 -5.9030986 4.81258389 1.951157858 -1.48268050
Unstimulated_S96 -0.52027265 -5.8988891 2.93362818 -0.119240743 -4.71252497
LPS_6h_S10 3.41059912 4.8026164 0.85047833 1.780947630 2.24298684
LPS_6h_S11 -7.00637482 3.0573865 3.72876703 -4.809485380 -1.01143877
LPS_6h_S12 -4.95769809 3.8583591 0.58008114 -4.185700849 1.35934163
LPS_6h_S14 -0.07711724 4.3329526 3.15257157 -3.697925011 -1.57703436
LPS_6h_S15 -5.99004365 8.5170997 -1.98223087 -1.112486090 1.72521617
LPS_6h_S16 3.25273365 6.9835752 -2.50910046 -3.692120681 6.61555100
LPS_6h_S17 0.16567101 6.5309662 -1.76889325 8.012914146 3.22717809
LPS_6h_S18 3.91965487 4.4234775 0.25898982 -1.530347415 2.43127779
LPS_6h_S19 12.97789991 -0.5000554 -25.02491519 -1.781497016 -2.97902446
LPS_6h_S1 -7.03719726 7.3093473 -0.77765713 5.999730375 -0.53922083
LPS_6h_S20 -2.46467845 6.5408224 0.65057878 -3.523567289 -1.27990616
LPS_6h_S21 1.70142697 3.9383284 3.90496265 -2.314455582 0.58854748
LPS_6h_S22 -1.07460769 5.2207599 1.27816169 1.794528181 1.68747620
LPS_6h_S23 3.10603550 2.8995015 0.87847805 7.187948276 -2.40253086
LPS_6h_S24 2.75229554 6.4929678 4.28272143 1.379262311 2.22417765
LPS_6h_S25 -5.69985751 5.4299413 1.84199357 -5.745291220 0.76557472
LPS_6h_S26 2.71074604 7.2928612 -3.27872848 0.033622977 -1.24818749
LPS_6h_S27 2.90714330 6.2549945 1.06585578 4.290505200 2.31777167
LPS_6h_S28 -7.89316669 9.1771654 -1.93875068 6.521291494 -0.78088003
LPS_6h_S29 -4.46761743 6.4024421 -1.46099382 -4.947388517 1.68933704
LPS_6h_S2 -3.70134238 4.2744424 2.18699427 -3.738905322 -2.27810844
LPS_6h_S30 -7.67553517 8.8862112 -0.65398568 -1.856313382 -0.16399463
LPS_6h_S34 -3.55916101 6.2490548 -0.35526252 5.723292570 -2.82033325
LPS_6h_S35 10.26778553 0.4457759 -14.57019357 -3.176624180 -1.57986357
LPS_6h_S37 6.02404404 2.3526877 0.41900957 6.941185763 3.36490615
LPS_6h_S38 -2.48745627 3.0196229 1.72433835 -5.331153776 -0.78868036
LPS_6h_S39 -3.33697250 5.0115097 -0.47750058 -5.776854123 -3.53991962
LPS_6h_S3 -0.83089207 5.5145241 -1.01037500 -5.480371858 3.06046678
LPS_6h_S40 -1.81799956 8.1188865 0.03516849 3.021234920 -3.55901507
LPS_6h_S42 -1.65975486 6.9848698 0.04203313 -2.974717382 -0.76239455
LPS_6h_S43 1.33163332 0.5941597 4.86591317 -4.589921089 -3.17288891
LPS_6h_S44 3.85282419 3.4886721 1.01623599 -3.487203830 3.33621142
LPS_6h_S46 9.05482707 -1.0029173 5.80845865 -0.571412817 -0.78218627
LPS_6h_S48 -6.31474913 8.4861760 -0.86959683 0.815151097 -1.84176737
LPS_6h_S49 2.16789966 7.4756361 1.63731632 4.872575491 0.91741317
LPS_6h_S4 4.12396129 4.3565106 5.66052678 3.763670326 -4.61965416
LPS_6h_S50 5.05947785 3.1909071 -1.11006714 0.842027130 1.73581780
LPS_6h_S51 -3.60413100 5.9866181 0.66599037 1.456320190 -1.55938739
LPS_6h_S53 -1.87468034 3.9185788 0.87665262 -3.886568930 3.06134877
LPS_6h_S55 3.05305193 3.6942515 4.64567130 -1.659525273 0.17408672
LPS_6h_S56 16.02356893 -2.4976457 6.04171666 -0.778524362 -2.12019207
LPS_6h_S57 1.40436177 3.6939251 1.82268348 2.165641163 -3.37870649
LPS_6h_S58 -1.06997022 4.6632060 -3.01245754 -6.021806669 2.54405641
LPS_6h_S59 4.83412988 0.8066678 3.08974965 -1.002525091 -0.25579360
LPS_6h_S5 6.72565183 3.6389907 5.81248089 2.443312906 -2.09759812
LPS_6h_S60 8.58558030 0.3053119 5.96519081 0.754855349 0.98161175
LPS_6h_S61 8.62886105 0.9265787 5.58806183 1.165773221 -2.88207983
LPS_6h_S63 -3.71412623 6.4200202 0.15583612 1.784154397 1.78659200
LPS_6h_S64 3.70134435 5.2299460 1.10664240 3.591449839 -3.97923918
LPS_6h_S65 5.39827669 3.2897982 -2.30315801 -3.283495666 0.48639437
LPS_6h_S66 2.72151929 1.0535344 1.93643449 -8.168103819 3.23690978
LPS_6h_S67 13.72681695 0.1066333 -14.98690823 -0.121052195 -3.67368564
LPS_6h_S68 -4.29003096 6.2523365 3.01941856 -1.582325433 2.52132017
LPS_6h_S6 4.65411117 4.1693852 -1.97612448 -0.511589227 -3.11966973
LPS_6h_S70 1.65116754 3.5729276 6.06804155 -0.470158640 0.96624273
LPS_6h_S74 -0.03292703 1.9839559 3.27426863 -5.816737301 -1.17881452
LPS_6h_S76 -2.99857718 6.1842686 1.21836077 -4.259185554 1.04919203
LPS_6h_S77 4.71109312 6.8728790 -1.72215350 3.925529129 -4.45211994
LPS_6h_S78 12.97375958 -0.1385897 -14.76922619 -4.623364748 -2.15527185
LPS_6h_S79 14.81506369 -2.1037154 0.94989966 -2.204786844 -0.75277072
LPS_6h_S7 4.12414975 0.7361231 4.93784192 4.797098269 -2.78452929
LPS_6h_S80 -11.32182828 5.7084215 -2.80129888 4.149089105 -0.03603544
LPS_6h_S81 -1.27659653 1.1660312 1.41473961 -6.802943862 -5.72519432
LPS_6h_S83 16.36402421 -1.6575890 6.21153834 1.427235737 -0.73814919
LPS_6h_S84 -1.83424737 6.3887328 1.35240135 2.220242728 -4.49141816
LPS_6h_S86 6.56470875 4.0209304 -5.36589766 1.891070782 2.41680209
LPS_6h_S87 6.93158472 4.5905606 -4.89029618 5.173998039 5.14783421
LPS_6h_S88 -9.68164837 6.6423318 -1.42270777 -1.702127235 -3.04629747
LPS_6h_S89 2.65038377 8.6424984 -0.13397975 3.829378276 1.11653735
LPS_6h_S8 -3.68096328 4.3332544 1.53220666 -6.333910685 2.89330302
LPS_6h_S90 -8.58848657 9.5353756 0.65084652 -1.374325026 -1.21007512
LPS_6h_S91 9.01267053 0.5864734 1.08695356 -6.216677798 1.57159708
LPS_6h_S93 -6.55617226 4.9763087 4.31312777 -1.008142282 -5.76542630
LPS_6h_S9 4.22759613 7.1684539 3.62730613 5.808200530 1.50418649
$MSTtree
IGRAPH 37b136c UNW- 3 2 --
+ attr: name (v/c), weight (e/n)
+ edges from 37b136c (vertex names):
[1] 1--2 2--3
$clusterid
Unstimulated_S15 Unstimulated_S16 Unstimulated_S17 Unstimulated_S18
1 3 1 3
Unstimulated_S23 Unstimulated_S24 Unstimulated_S26 Unstimulated_S27
1 1 1 1
Unstimulated_S28 Unstimulated_S2 Unstimulated_S31 Unstimulated_S32
1 1 1 1
Unstimulated_S34 Unstimulated_S37 Unstimulated_S38 Unstimulated_S39
1 3 1 1
Unstimulated_S3 Unstimulated_S41 Unstimulated_S46 Unstimulated_S47
1 1 1 1
Unstimulated_S48 Unstimulated_S49 Unstimulated_S4 Unstimulated_S51
1 1 1 3
Unstimulated_S52 Unstimulated_S53 Unstimulated_S54 Unstimulated_S55
1 1 1 3
Unstimulated_S5 Unstimulated_S60 Unstimulated_S61 Unstimulated_S62
3 1 1 1
Unstimulated_S63 Unstimulated_S66 Unstimulated_S6 Unstimulated_S70
1 1 1 1
Unstimulated_S71 Unstimulated_S72 Unstimulated_S74 Unstimulated_S75
1 1 1 1
Unstimulated_S76 Unstimulated_S77 Unstimulated_S78 Unstimulated_S79
1 1 1 1
Unstimulated_S7 Unstimulated_S80 Unstimulated_S81 Unstimulated_S82
3 1 3 1
Unstimulated_S83 Unstimulated_S84 Unstimulated_S85 Unstimulated_S8
1 3 3 3
Unstimulated_S90 Unstimulated_S92 Unstimulated_S93 Unstimulated_S95
3 1 1 3
Unstimulated_S96 LPS_6h_S10 LPS_6h_S11 LPS_6h_S12
1 2 2 2
LPS_6h_S14 LPS_6h_S15 LPS_6h_S16 LPS_6h_S17
2 2 2 3
LPS_6h_S18 LPS_6h_S19 LPS_6h_S1 LPS_6h_S20
2 3 3 2
LPS_6h_S21 LPS_6h_S22 LPS_6h_S23 LPS_6h_S24
2 2 3 2
LPS_6h_S25 LPS_6h_S26 LPS_6h_S27 LPS_6h_S28
2 2 2 3
LPS_6h_S29 LPS_6h_S2 LPS_6h_S30 LPS_6h_S34
2 2 2 3
LPS_6h_S35 LPS_6h_S37 LPS_6h_S38 LPS_6h_S39
3 3 2 2
LPS_6h_S3 LPS_6h_S40 LPS_6h_S42 LPS_6h_S43
2 2 2 2
LPS_6h_S44 LPS_6h_S46 LPS_6h_S48 LPS_6h_S49
2 2 2 2
LPS_6h_S4 LPS_6h_S50 LPS_6h_S51 LPS_6h_S53
2 2 2 2
LPS_6h_S55 LPS_6h_S56 LPS_6h_S57 LPS_6h_S58
2 3 2 2
LPS_6h_S59 LPS_6h_S5 LPS_6h_S60 LPS_6h_S61
2 2 2 2
LPS_6h_S63 LPS_6h_S64 LPS_6h_S65 LPS_6h_S66
2 2 2 2
LPS_6h_S67 LPS_6h_S68 LPS_6h_S6 LPS_6h_S70
3 2 2 2
LPS_6h_S74 LPS_6h_S76 LPS_6h_S77 LPS_6h_S78
2 2 2 3
LPS_6h_S79 LPS_6h_S7 LPS_6h_S80 LPS_6h_S81
3 3 3 2
LPS_6h_S83 LPS_6h_S84 LPS_6h_S86 LPS_6h_S87
3 2 3 3
LPS_6h_S88 LPS_6h_S89 LPS_6h_S8 LPS_6h_S90
2 2 2 2
LPS_6h_S91 LPS_6h_S93 LPS_6h_S9
2 2 2
$clucenter
[,1] [,2] [,3] [,4] [,5]
[1,] -4.701208 -6.0001699 -0.2304098 0.6541102 0.1480404
[2,] 0.165610 4.8079024 1.4012641 -1.3636384 -0.2311658
[3,] 6.580446 -0.3347654 -2.3244674 1.6315513 0.2220892
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.