OUT=vector.drawArrow(OUT,P=0.9,SHOW=TRUE, COL=OUT$COL, SHOW.SUMMIT=TRUE)

[WANGROUOYI]

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)))