From aab1671b7e76bb245eff81cb54e05fce902284d3 Mon Sep 17 00:00:00 2001 From: aduvermy <arnaud.duvermy@ens-lyon.fr> Date: Tue, 15 Feb 2022 12:11:38 +0000 Subject: [PATCH] modif tuto htrsim --- src/tutorial_htrsim.R | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) diff --git a/src/tutorial_htrsim.R b/src/tutorial_htrsim.R index d31e7b8..def0def 100644 --- a/src/tutorial_htrsim.R +++ b/src/tutorial_htrsim.R @@ -1,5 +1,7 @@ ################################### Getting started ######################################################## - +### Required +library(tidyverse) +library(data.table) ############################################################################### @@ -24,14 +26,14 @@ set_gene_name = paste0('gene', 1:N_genes) set_gene_name ### Fill INPUT Dataframe -genes_NB_params <- samples$name %>% - map(~(list(name=., #sample_name - n_replicates = sample(1:max_N_replicates, 1), #random int between 1 & max_N_replicates +#samples$name %>% +genes_NB_params <- map2(.x= samples$name, .y = samples$n_rep, ~(list(name=.x, #sample_name + n_replicates = .y,#sample(1:max_N_replicates, 1), #random int between 1 & max_N_replicates name_gene = set_gene_name, # gene_name mu = runif(100,1000, n = N_genes), #mu(ij) alpha = set_alpha_per_gene))) %>% # alpha(i) rbindlist(.) %>% data.frame() ## convert to dtf - +genes_NB_params ## Use filter to understand our dtf ## Notice that alpha is equal for equivalent gene between sample genes_NB_params %>% filter(name_gene == "gene2") -- GitLab