Skip to content
Snippets Groups Projects
Commit a1cb82de authored by Arnaud Duvermy's avatar Arnaud Duvermy
Browse files

update setup_cntsGene

parent da30f316
No related branches found
No related tags found
No related merge requests found
...@@ -66,7 +66,6 @@ setup_countGener <- function(sample_names = "my_first_lib", n_rep = NULL , gene_ ...@@ -66,7 +66,6 @@ setup_countGener <- function(sample_names = "my_first_lib", n_rep = NULL , gene_
######### HANDLE EXCEPTION ####### ######### HANDLE EXCEPTION #######
setup = handle_except(sample_names, n_rep , gene_names , gene_dispersion, n_genes) setup = handle_except(sample_names, n_rep , gene_names , gene_dispersion, n_genes)
######## HANDLE TYPE MU ########## ######## HANDLE TYPE MU ##########
if(is.null(mu)) mu = .mu_generator # default function to generate mu if(is.null(mu)) mu = .mu_generator # default function to generate mu
...@@ -88,7 +87,7 @@ setup_countGener <- function(sample_names = "my_first_lib", n_rep = NULL , gene_ ...@@ -88,7 +87,7 @@ setup_countGener <- function(sample_names = "my_first_lib", n_rep = NULL , gene_
######## BUILD AN INPUT DTF FOR count_generator ############ ######## BUILD AN INPUT DTF FOR count_generator ############
nBinom_params <- purrr::map2(.x= setup$samples, .y = setup$rep, nBinom_params <- purrr::map2(.x= setup$samples, .y = setup$rep,
~(list(name=.x, #sample_name ~(list(name=.x, #sample_name
n_replicates = .y, # random int between 1 & max_N_replicates n_replicates = .y, # number replicates
name_gene = setup$genes, # gene_name name_gene = setup$genes, # gene_name
mu = mu.dtf %>% dplyr::select(all_of(.x)) %>% unlist() , #mu(ij) mu = mu.dtf %>% dplyr::select(all_of(.x)) %>% unlist() , #mu(ij)
alpha = setup$alpha))) %>% # alpha(i) alpha = setup$alpha))) %>% # alpha(i)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment