diff tools/myTools/bin/pertanalysis.R @ 1:7e5c71b2e71f draft default tip

Uploaded

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tools/myTools/bin/pertanalysis.R	Wed Dec 22 16:00:34 2021 +0000
@@ -0,0 +1,24 @@
+#!/usr/bin/env Rscript 
+.libPaths('~/R_libs_36')
+args = commandArgs(trailingOnly=TRUE)
+
+library("dplyr", warn.conflicts = FALSE)
+library(ggplot2)
+
+
+perts_lab<-read.delim(args[1],sep=" ",row.names = 1) 
+
+perts_lab<-replace(perts_lab,perts_lab==-1,"down")
+perts_lab<-replace(perts_lab,perts_lab==1,"up")
+perts_lab<-replace(perts_lab,perts_lab==0,"nochange")
+
+df1<-data.frame(perts_lab)
+df1$up <- apply(df1, 1, function(x) length(which(x=="up")))
+df1$down<-apply(df1, 1, function(x) length(which(x=="down")))
+df1$total<-df1$up+df1$down
+df1<-df1[order(df1$total),]
+write.table(df1,"successful_controlnode_perturbations.txt",quote=F)
+
+# oh<-t(perts_lab) %>% reshape2::melt()
+# pdf(file="pertrubation_trends.pdf")
+# ggplot(oh, aes(x=value)) + facet_wrap(~Var1,scales = "free_x",shrink=FALSE) + geom_bar(aes(y = (..count..)/nrow(perts_lab))) + scale_y_continuous(labels=scales::percent)
\ No newline at end of file