From df60a11fd49b44697213f957680c943eda3e66ef Mon Sep 17 00:00:00 2001 From: aduvermy <arnaud.duvermy@ens-lyon.fr> Date: Mon, 7 Feb 2022 09:50:56 +0000 Subject: [PATCH] add src repo --- src/counts_matrix_generator.R | 34 ++++++++++++++++++++++++++++++++++ 1 file changed, 34 insertions(+) create mode 100644 src/counts_matrix_generator.R diff --git a/src/counts_matrix_generator.R b/src/counts_matrix_generator.R new file mode 100644 index 0000000..4d42d01 --- /dev/null +++ b/src/counts_matrix_generator.R @@ -0,0 +1,34 @@ +#seed +library(tidyverse) +library(data.table) + + +# fix seed +set.seed(123) + +#params +N_gene = 6000 + + + +## function generator +rnbinom_generator <- function(n_gene, mu_simul, alpha){ + + my_counts <- rnbinom(n = n_gene, mu = mu_simul, size = alpha) + return(my_counts) + + } + + + +## build my lib counts +env1 <- rnbinom_generator(N_gene, 50, 5) +env2 <- rnbinom_generator(N_gene, 100, 5) + + + +## export my lib count +names(env1) <- sprintf("gene%d", 1:N_gene) +names(env2) <- sprintf("gene%d", 1:N_gene) +env1 %>% data.frame(.) %>% fwrite(., "~/mydatalocal/counts_simulation/env1.csv", row.names = T, col.names = F, sep = "\t") +env2 %>% data.frame(.) %>% fwrite(., "~/mydatalocal/counts_simulation/env2.csv", row.names = T, col.names = F, sep = "\t") -- GitLab