#A short script plotting values of stat1 gene
#read the input data stat1_interactions_new <- read.csv("string_interactions.tsv",sep ="\t") #extract the first two and the last column statsmallnew <- stat1_interactions_new[c(1,2,15)] #extract the interactions with more than 0.5 score statsmallnewfiltered = subset(statsmallnew, combined_score >= 0.5) #write the filtered interactions in to a file write.table(statsmallnewfiltered,"statsmallfiltered.txt",sep="\t",quote = FALSE, row.names = FALSE) #print the score distribution graph png('statsmallnew.png') plot(statsmallnewfiltered$combined_score) dev.off()
message("script ran successfully, all outputs are in the working directory")