diff --git a/src/bin/kmerclust_clust.R b/src/bin/kmerclust_clust.R index 346e6a87094fd4401a81fb5013635a1439216aa7..c1153022a698572a5e2174ef8b83a33377f16956 100644 --- a/src/bin/kmerclust_clust.R +++ b/src/bin/kmerclust_clust.R @@ -5,16 +5,6 @@ print(args) load(file = paste0(args[1], ".Rdata")) -params = list( - A = list(kappa = 0.4, l1 = mean(count$count_m), l2 = mean(count$count_m), l3 = mean(count$count_m)), -) -if (args[2] %in% c("XY", "XO")) { - params$X <- list(kappa = 0.3, l1 = 1, l2 = mean(count$count_m), l3 = mean(count$count_m)) -} -if (args[2] == "XY") { - params$Y <- list(kappa = 0.3, l1 = mean(count$count_m), l2 = 1, l3 = mean(count$count_m)) -} - res <- count %>% dplyr::select(count_m, count_f) %>% mutate(