VEC = scRNAsub@reductions$umap@cell.embeddings
rownames(VEC) = colnames(scRNAsub)
PCA = scRNAsub@reductions$pca@cell.embeddings
source('https://raw.githubusercontent.com/jumphone/Vector/master/Vector.R')
Remove quantile-based colinearity among PCs (new feature in VECTOR 0.0.3):
PCA=vector.rankPCA(PCA)
rownames(PCA) = colnames(scRNAsub)
source('https://raw.githubusercontent.com/jumphone/Vector/master/Vector.R')
Define pixel
OUT=vector.buildGrid(VEC, N=30,SHOW=TRUE)
Build network
OUT=vector.buildNet(OUT, CUT=1, SHOW=TRUE)
Calculate Quantile Polarization (QP) score
OUT=vector.getValue(OUT, PCA, SHOW=TRUE)
Get pixel's QP score
OUT=vector.gridValue(OUT,SHOW=TRUE)
Find starting point
OUT=vector.autoCenter(OUT,UP=0.9,SHOW=TRUE)
Infer vector
OUT=vector.drawArrow(OUT,P=0.9,SHOW=TRUE, COL=OUT$COL, SHOW.SUMMIT=TRUE)
错误于(function (cond) :
在为函数“t”选择方法时计算参数“x”时出错：不是所有的is.atomic(y)都是TRUE
Called from: (function (cond)
.Internal(C_tryCatchHelper(addr, 1L, cond)))(list(message = "不是所有的is.atomic(y)都是TRUE",
call = var(x)))
