From bcc64bccd23293e5bdbf9d290827fd87358ea9c5 Mon Sep 17 00:00:00 2001 From: "Verweij, Stefan" <stefan.verweij@wur.nl> Date: Tue, 25 Apr 2017 10:52:14 +0200 Subject: [PATCH] Add new file --- Command_Lines_R | 22 ++++++++++++++++++++++ 1 file changed, 22 insertions(+) create mode 100644 Command_Lines_R diff --git a/Command_Lines_R b/Command_Lines_R new file mode 100644 index 0000000..b35f9dd --- /dev/null +++ b/Command_Lines_R @@ -0,0 +1,22 @@ +#Command Lines used to plot the bimodal distribution: + #The chosen input file contained the GC and GC3 values against their percentage in the CDS +>GC3 <- read.table(choose.files(), h=T) +>ggplot(GC3, aes(GC3, CDS)) ++geom_point() ++geom_lline(color=”red”) ++ggtitle(“B”) + +#Command Lines used for randomization: + #The chosen input file contained all the loci of the CDS and their GC3 values +>my.data <- as.numeric(as.vector(read.table(choose.files())[,1])) +>replicate(1000, mean(sample(my.data, 590))) + +#Command Lines used to plot normal distribution + observed mean: + #The chosen input files were 1) ND: data of randomization and 2) NDobs: data containing observed mean +>ND <- read.table(choose.files(), h=T) +>NDobs <- read.table(choose.files(), h=T) +>ggplot(ND,aes(GC3,Appearance)) ++geom_point() ++geom_smooth(color="red",size=3) ++stat_function(fun=dnorm,color="grey") ++geom_point(data=NDobs,shape=21,color="black",fill="yellow",size=3,stroke=3) -- GitLab