# NOT RUN {
data(demo.sp);
demo.sp = scaleCountMatrix(
obj=demo.sp,
mat="gmat",
cov=SnapATAC::rowSums(demo.sp, mat="bmat"),
method="RPM"
)
gene.names = c(
"Prdm14", "E330040D14Rik", "Gm17971",
"Gm17970", "Defb44-ps", "Gm7357",
"Gm37265", "Kctd18", "Gm37143"
)
plotGene(
obj=demo.sp,
gene.names=gene.names,
viz.method="tsne"
)
# }
Run the code above in your browser using DataLab