diff --git a/src/tutorial_htrsim.R b/src/tutorial_htrsim.R index d31e7b88d0dd3dc6469bc33cbafe16dd04de8e11..def0def94177b97942e3d65ed1a7c0ba6b8f0463 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")