Mercurial > repos > iuc > raceid_inspectclusters
diff scripts/pseudotemporal.R @ 0:9fec5dd8fbb9 draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/raceid3 commit f880060c478d42202df5b78a81329f8af56b1138
author | iuc |
---|---|
date | Thu, 22 Nov 2018 04:45:41 -0500 |
parents | |
children | 41f34e925bd5 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/scripts/pseudotemporal.R Thu Nov 22 04:45:41 2018 -0500 @@ -0,0 +1,58 @@ +#!/usr/bin/env R +VERSION = "0.1" + +args = commandArgs(trailingOnly = T) + +if (length(args) != 1){ + message(paste("VERSION:", VERSION)) + stop("Please provide the config file") +} + +suppressWarnings(suppressPackageStartupMessages(require(RaceID))) +source(args[1]) + +test <- list() +test$side = 3 +test$line = 3 +second <- test +second$cex = 0.5 +second$line = 2.5 + + +do.pseudotemp <- function(sc){ + pdf(out.pdf) + ltr <- Ltree(sc) + ltr <- compentropy(ltr) + ltr <- do.call(projcells, c(ltr, pstc.projc)) + ltr <- do.call(projback, c(ltr, pstc.projb)) + ltr <- lineagegraph(ltr) + ltr <- do.call(comppvalue, c(ltr, pstc.comppval)) + x <- do.call(compscore, c(ltr, pstc.compscore)) + print(do.call(mtext, c("Compute Score", test))) + print(do.call(mtext, c("No. of inter-cluster links / Delta median entropy of each cluster / StemID2 score (combination of both)", second))) + plotdistanceratio(ltr) + print(do.call(mtext, c("Cell-to-Cell Distance Ratio", test))) + print(do.call(mtext, c("Original vs High-dimensional Embedded Space", second))) + do.call(plotgraph, c(ltr, pstc.plotgraph)) + print(do.call(mtext, c("Lineage Trajectories ", test))) + print(do.call(mtext, c("Colour = Level of Significance, Width = Link Score ", second))) + plotspantree(ltr) + print(do.call(mtext, c("Minimum Spanning Tree", test))) + plotprojections(ltr) + print(do.call(mtext, c("Minimum Spanning Tree", test))) + print(do.call(mtext, c("Cells Projected onto Links", second))) + test$side = 4 + test$line = 0 + plotlinkscore(ltr) + print(do.call(mtext, c("Link Score", test))) + projenrichment(ltr) + print(do.call(mtext, c("Enrichment Ratios", test))) + dev.off() + + return(ltr) +} + +ltr <- do.pseudotemp(in.rdat) + + +saveRDS(ltr, out.rdat)