Mercurial > repos > laurenmarazzi > netisce_test
diff tools/myTools/bin/pertanalysis.R @ 1:7e5c71b2e71f draft default tip
Uploaded
author | laurenmarazzi |
---|---|
date | Wed, 22 Dec 2021 16:00:34 +0000 |
parents | |
children |
line wrap: on
line diff
--- /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