R/dimsumms__library_bottleneck.R
In lehner-lab/dimsumms: Analysis scripts to reproduce the figures and results from DiMSum manuscript
Defines functions
dimsumms__library_bottleneck
Documented in
dimsumms__library_bottleneck
#' dimsumms__library_bottleneck
#'
#' Simulate a library bottleneck.
#'
#' @param input_file path to input file (required)
#' @param outpath output path for plots and saved objects (required)
#' @param bottleneck.alpha bottleneck factor
#' @param Sequencing.Error sequencing error (default:0.02)
#'
#' @return Nothing
#' @export
#' @import data.table
dimsumms__library_bottleneck <- function(
input_file,
outpath,
bottleneck.alpha,
Sequencing.Error = 0.02
){
original.data <- as.data.frame(fread(input_file))
# subsample input 1 and use that as the "real library
original.library <- (original.data$input1_e1_s0_bNA_count + original.data$input2_e2_s0_bNA_count + original.data$input3_e3_s0_bNA_count + original.data$input4_e4_s0_bNA_count) / 4
bottleneck.factor <- rpois(n = length(original.library),
lambda = original.library * bottleneck.alpha) / original.library
bottleneck.factor[is.na(bottleneck.factor)] <- 0
set.seed(123)
### Input
New.Input.1 <- sample(x = 1:length(bottleneck.factor),
size = sum(original.data$input1_e1_s0_bNA_count),
prob = original.data$input1_e1_s0_bNA_count * bottleneck.factor,
replace = TRUE)
New.Input.1 <- factor(New.Input.1, levels = 1:length(bottleneck.factor))
New.Input.1 <- as.numeric(table(New.Input.1))
New.Input.2 <- sample(x = 1:length(bottleneck.factor),
size = sum(original.data$input2_e2_s0_bNA_count),
prob = original.data$input2_e2_s0_bNA_count * bottleneck.factor,
replace = TRUE)
New.Input.2 <- factor(New.Input.2, levels = 1:length(bottleneck.factor))
New.Input.2 <- as.numeric(table(New.Input.2))
New.Input.3 <- sample(x = 1:length(bottleneck.factor),
size = sum(original.data$input3_e3_s0_bNA_count),
prob = original.data$input3_e3_s0_bNA_count * bottleneck.factor,
replace = TRUE)
New.Input.3 <- factor(New.Input.3, levels = 1:length(bottleneck.factor))
New.Input.3 <- as.numeric(table(New.Input.3))
New.Input.4 <- sample(x = 1:length(bottleneck.factor),
size = sum(original.data$input4_e4_s0_bNA_count),
prob = original.data$input4_e4_s0_bNA_count * bottleneck.factor,
replace = TRUE)
New.Input.4 <- factor(New.Input.4, levels = 1:length(bottleneck.factor))
New.Input.4 <- as.numeric(table(New.Input.4))
### Output
New.Output.1A <- sample(x = 1:nrow(original.data),
size = sum(original.data$output1A_e1_s1_b1_count),
prob = original.data$output1A_e1_s1_b1_count * bottleneck.factor,
replace = TRUE)
New.Output.1A <- factor(New.Output.1A, levels = 1:nrow(original.data))
New.Output.1A <- as.numeric(table(New.Output.1A))
New.Output.1B <- sample(x = 1:nrow(original.data),
size = sum(original.data$output1B_e1_s1_b2_count),
prob = original.data$output1B_e1_s1_b2_count * bottleneck.factor,
replace = TRUE)
New.Output.1B <- factor(New.Output.1B, levels = 1:nrow(original.data))
New.Output.1B <- as.numeric(table(New.Output.1B))
New.Output.2A <- sample(x = 1:nrow(original.data),
size = sum(original.data$output2A_e2_s1_b1_count),
prob = original.data$output2A_e2_s1_b1_count * bottleneck.factor,
replace = TRUE)
New.Output.2A <- factor(New.Output.2A, levels = 1:nrow(original.data))
New.Output.2A <- as.numeric(table(New.Output.2A))
New.Output.2B <- sample(x = 1:nrow(original.data),
size = sum(original.data$output2B_e2_s1_b2_count),
prob = original.data$output2B_e2_s1_b2_count * bottleneck.factor,
replace = TRUE)
New.Output.2B <- factor(New.Output.2B, levels = 1:nrow(original.data))
New.Output.2B <- as.numeric(table(New.Output.2B))
New.Output.3A <- sample(x = 1:nrow(original.data),
size = sum(original.data$output3A_e3_s1_b1_count),
prob = original.data$output3A_e3_s1_b1_count * bottleneck.factor,
replace = TRUE)
New.Output.3A <- factor(New.Output.3A, levels = 1:nrow(original.data))
New.Output.3A <- as.numeric(table(New.Output.3A))
New.Output.3B <- sample(x = 1:nrow(original.data),
size = sum(original.data$output3B_e3_s1_b2_count),
prob = original.data$output3B_e3_s1_b2_count * bottleneck.factor,
replace = TRUE)
New.Output.3B <- factor(New.Output.3B, levels = 1:nrow(original.data))
New.Output.3B <- as.numeric(table(New.Output.3B))
New.Output.4A <- sample(x = 1:nrow(original.data),
size = sum(original.data$output4A_e4_s1_b1_count),
prob = original.data$output4A_e4_s1_b1_count * bottleneck.factor,
replace = TRUE)
New.Output.4A <- factor(New.Output.4A, levels = 1:nrow(original.data))
New.Output.4A <- as.numeric(table(New.Output.4A))
New.Output.4B <- sample(x = 1:nrow(original.data),
size = sum(original.data$output4B_e4_s1_b2_count),
prob = original.data$output4B_e4_s1_b2_count * bottleneck.factor,
replace = TRUE)
New.Output.4B <- factor(New.Output.4B, levels = 1:nrow(original.data))
New.Output.4B <- as.numeric(table(New.Output.4B))
Bottlenecked.Data <- original.data
Bottlenecked.Data$input1_e1_s0_bNA_count <- New.Input.1
Bottlenecked.Data$input2_e2_s0_bNA_count <- New.Input.2
Bottlenecked.Data$input3_e3_s0_bNA_count <- New.Input.3
Bottlenecked.Data$input4_e4_s0_bNA_count <- New.Input.4
Bottlenecked.Data$output1A_e1_s1_b1_count <- New.Output.1A
Bottlenecked.Data$output2A_e2_s1_b1_count <- New.Output.2A
Bottlenecked.Data$output3A_e3_s1_b1_count <- New.Output.3A
Bottlenecked.Data$output4A_e4_s1_b1_count <- New.Output.4A
Bottlenecked.Data$output1B_e1_s1_b2_count <- New.Output.1B
Bottlenecked.Data$output2B_e2_s1_b2_count <- New.Output.2B
Bottlenecked.Data$output3B_e3_s1_b2_count <- New.Output.3B
Bottlenecked.Data$output4B_e4_s1_b2_count <- New.Output.4B
Idx.WT <- which(Bottlenecked.Data$Nham_nt == 0)
Idx.Singles <- which(Bottlenecked.Data$Nham_nt == 1)
Idx.Doubles <- which(Bottlenecked.Data$Nham_nt == 2)
Idx.Triples <- which(Bottlenecked.Data$Nham_nt == 3)
Idx.Quadruples <- which(Bottlenecked.Data$Nham_nt == 4)
# sequencing errors input 1
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Quadruples] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.WT] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors input 2
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Quadruples] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.WT] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors input 3
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Quadruples] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.WT] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors input 4
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Quadruples] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.WT] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 1
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Quadruples] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.WT] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 2
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Quadruples] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.WT] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 3
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Quadruples] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.WT] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 4
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Quadruples] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.WT] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 1
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Quadruples] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.WT] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 2
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Quadruples] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.WT] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 3
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Quadruples] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.WT] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 4
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Quadruples] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.WT] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.WT] - Counts.To.Subtract.From.WT
#Finalise
variant_data_merge <- as.data.table(Bottlenecked.Data)
#Replace negative counts with 0
for(j in grep("_count$", names(variant_data_merge))){
set(variant_data_merge, i=which(variant_data_merge[[j]]<0), j=j, value=0)
}
#Remove variants with zero counts in all samples
variant_data_merge <- variant_data_merge[apply(variant_data_merge[,.SD,.SDcols = names(variant_data_merge)[grep("_count$", names(variant_data_merge))]], 1, sum)!=0,]
#Reformat
variant_data_merge <- variant_data_merge[,.SD,,.SDcols = grep("nt_seq|_count$", names(variant_data_merge))]
names(variant_data_merge) <- sapply(strsplit(names(variant_data_merge), "_e"), '[', 1)
#Save
OutFileName <- file.path(outpath, paste0("BB_TARDBP_290_2017-06-13_LibraryBottleneck_", substr(as.character(bottleneck.alpha), 3, 4), "_t0_variant_data_merge.tsv"))
write.table(variant_data_merge, OutFileName, sep = "\t", quote = F, row.names = F)
}
lehner-lab/dimsumms documentation built on June 23, 2020, 5:48 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Defines functions dimsumms__library_bottleneck
Documented in dimsumms__library_bottleneck
#' dimsumms__library_bottleneck
#'
#' Simulate a library bottleneck.
#'
#' @param input_file path to input file (required)
#' @param outpath output path for plots and saved objects (required)
#' @param bottleneck.alpha bottleneck factor
#' @param Sequencing.Error sequencing error (default:0.02)
#'
#' @return Nothing
#' @export
#' @import data.table
dimsumms__library_bottleneck <- function(
input_file,
outpath,
bottleneck.alpha,
Sequencing.Error = 0.02
){
original.data <- as.data.frame(fread(input_file))
# subsample input 1 and use that as the "real library
original.library <- (original.data$input1_e1_s0_bNA_count + original.data$input2_e2_s0_bNA_count + original.data$input3_e3_s0_bNA_count + original.data$input4_e4_s0_bNA_count) / 4
bottleneck.factor <- rpois(n = length(original.library),
lambda = original.library * bottleneck.alpha) / original.library
bottleneck.factor[is.na(bottleneck.factor)] <- 0
set.seed(123)
### Input
New.Input.1 <- sample(x = 1:length(bottleneck.factor),
size = sum(original.data$input1_e1_s0_bNA_count),
prob = original.data$input1_e1_s0_bNA_count * bottleneck.factor,
replace = TRUE)
New.Input.1 <- factor(New.Input.1, levels = 1:length(bottleneck.factor))
New.Input.1 <- as.numeric(table(New.Input.1))
New.Input.2 <- sample(x = 1:length(bottleneck.factor),
size = sum(original.data$input2_e2_s0_bNA_count),
prob = original.data$input2_e2_s0_bNA_count * bottleneck.factor,
replace = TRUE)
New.Input.2 <- factor(New.Input.2, levels = 1:length(bottleneck.factor))
New.Input.2 <- as.numeric(table(New.Input.2))
New.Input.3 <- sample(x = 1:length(bottleneck.factor),
size = sum(original.data$input3_e3_s0_bNA_count),
prob = original.data$input3_e3_s0_bNA_count * bottleneck.factor,
replace = TRUE)
New.Input.3 <- factor(New.Input.3, levels = 1:length(bottleneck.factor))
New.Input.3 <- as.numeric(table(New.Input.3))
New.Input.4 <- sample(x = 1:length(bottleneck.factor),
size = sum(original.data$input4_e4_s0_bNA_count),
prob = original.data$input4_e4_s0_bNA_count * bottleneck.factor,
replace = TRUE)
New.Input.4 <- factor(New.Input.4, levels = 1:length(bottleneck.factor))
New.Input.4 <- as.numeric(table(New.Input.4))
### Output
New.Output.1A <- sample(x = 1:nrow(original.data),
size = sum(original.data$output1A_e1_s1_b1_count),
prob = original.data$output1A_e1_s1_b1_count * bottleneck.factor,
replace = TRUE)
New.Output.1A <- factor(New.Output.1A, levels = 1:nrow(original.data))
New.Output.1A <- as.numeric(table(New.Output.1A))
New.Output.1B <- sample(x = 1:nrow(original.data),
size = sum(original.data$output1B_e1_s1_b2_count),
prob = original.data$output1B_e1_s1_b2_count * bottleneck.factor,
replace = TRUE)
New.Output.1B <- factor(New.Output.1B, levels = 1:nrow(original.data))
New.Output.1B <- as.numeric(table(New.Output.1B))
New.Output.2A <- sample(x = 1:nrow(original.data),
size = sum(original.data$output2A_e2_s1_b1_count),
prob = original.data$output2A_e2_s1_b1_count * bottleneck.factor,
replace = TRUE)
New.Output.2A <- factor(New.Output.2A, levels = 1:nrow(original.data))
New.Output.2A <- as.numeric(table(New.Output.2A))
New.Output.2B <- sample(x = 1:nrow(original.data),
size = sum(original.data$output2B_e2_s1_b2_count),
prob = original.data$output2B_e2_s1_b2_count * bottleneck.factor,
replace = TRUE)
New.Output.2B <- factor(New.Output.2B, levels = 1:nrow(original.data))
New.Output.2B <- as.numeric(table(New.Output.2B))
New.Output.3A <- sample(x = 1:nrow(original.data),
size = sum(original.data$output3A_e3_s1_b1_count),
prob = original.data$output3A_e3_s1_b1_count * bottleneck.factor,
replace = TRUE)
New.Output.3A <- factor(New.Output.3A, levels = 1:nrow(original.data))
New.Output.3A <- as.numeric(table(New.Output.3A))
New.Output.3B <- sample(x = 1:nrow(original.data),
size = sum(original.data$output3B_e3_s1_b2_count),
prob = original.data$output3B_e3_s1_b2_count * bottleneck.factor,
replace = TRUE)
New.Output.3B <- factor(New.Output.3B, levels = 1:nrow(original.data))
New.Output.3B <- as.numeric(table(New.Output.3B))
New.Output.4A <- sample(x = 1:nrow(original.data),
size = sum(original.data$output4A_e4_s1_b1_count),
prob = original.data$output4A_e4_s1_b1_count * bottleneck.factor,
replace = TRUE)
New.Output.4A <- factor(New.Output.4A, levels = 1:nrow(original.data))
New.Output.4A <- as.numeric(table(New.Output.4A))
New.Output.4B <- sample(x = 1:nrow(original.data),
size = sum(original.data$output4B_e4_s1_b2_count),
prob = original.data$output4B_e4_s1_b2_count * bottleneck.factor,
replace = TRUE)
New.Output.4B <- factor(New.Output.4B, levels = 1:nrow(original.data))
New.Output.4B <- as.numeric(table(New.Output.4B))
Bottlenecked.Data <- original.data
Bottlenecked.Data$input1_e1_s0_bNA_count <- New.Input.1
Bottlenecked.Data$input2_e2_s0_bNA_count <- New.Input.2
Bottlenecked.Data$input3_e3_s0_bNA_count <- New.Input.3
Bottlenecked.Data$input4_e4_s0_bNA_count <- New.Input.4
Bottlenecked.Data$output1A_e1_s1_b1_count <- New.Output.1A
Bottlenecked.Data$output2A_e2_s1_b1_count <- New.Output.2A
Bottlenecked.Data$output3A_e3_s1_b1_count <- New.Output.3A
Bottlenecked.Data$output4A_e4_s1_b1_count <- New.Output.4A
Bottlenecked.Data$output1B_e1_s1_b2_count <- New.Output.1B
Bottlenecked.Data$output2B_e2_s1_b2_count <- New.Output.2B
Bottlenecked.Data$output3B_e3_s1_b2_count <- New.Output.3B
Bottlenecked.Data$output4B_e4_s1_b2_count <- New.Output.4B
Idx.WT <- which(Bottlenecked.Data$Nham_nt == 0)
Idx.Singles <- which(Bottlenecked.Data$Nham_nt == 1)
Idx.Doubles <- which(Bottlenecked.Data$Nham_nt == 2)
Idx.Triples <- which(Bottlenecked.Data$Nham_nt == 3)
Idx.Quadruples <- which(Bottlenecked.Data$Nham_nt == 4)
# sequencing errors input 1
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Quadruples] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.WT] <- Bottlenecked.Data$input1_e1_s0_bNA_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors input 2
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Quadruples] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.WT] <- Bottlenecked.Data$input2_e2_s0_bNA_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors input 3
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Quadruples] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.WT] <- Bottlenecked.Data$input3_e3_s0_bNA_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors input 4
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Quadruples] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.WT] <- Bottlenecked.Data$input4_e4_s0_bNA_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 1
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Quadruples] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))))
Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.WT] <- Bottlenecked.Data$output1A_e1_s1_b1_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 2
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Quadruples] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.WT] <- Bottlenecked.Data$output2A_e2_s1_b1_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 3
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Quadruples] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.WT] <- Bottlenecked.Data$output3A_e3_s1_b1_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 4
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Quadruples] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.WT] <- Bottlenecked.Data$output4A_e4_s1_b1_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 1
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Quadruples] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))))
Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.WT] <- Bottlenecked.Data$output1B_e1_s1_b2_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 2
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Quadruples] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.WT] <- Bottlenecked.Data$output2B_e2_s1_b2_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 3
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Quadruples] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.WT] <- Bottlenecked.Data$output3B_e3_s1_b2_count[Idx.WT] - Counts.To.Subtract.From.WT
# sequencing errors output 4
Counts.To.Add.To.Quadruples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Quadruples))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Quadruples] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Quadruples] + Counts.To.Add.To.Quadruples
Counts.To.Subtract.From.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples] - Counts.To.Subtract.From.Triples
Counts.To.Add.To.Triples <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Triples))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Triples] + Counts.To.Add.To.Triples
Counts.To.Subtract.From.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles] - Counts.To.Subtract.From.Doubles
Counts.To.Add.To.Doubles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Doubles))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Doubles] + Counts.To.Add.To.Doubles
Counts.To.Subtract.From.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles] - Counts.To.Subtract.From.Singles
Counts.To.Add.To.Singles <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.Singles))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.Singles] + Counts.To.Add.To.Singles
Counts.To.Subtract.From.WT <- as.vector(rmultinom(n = 1,
size = sum(Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.WT]) * Sequencing.Error,
prob = rep(1, times = length(Idx.WT))
))
Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.WT] <- Bottlenecked.Data$output4B_e4_s1_b2_count[Idx.WT] - Counts.To.Subtract.From.WT
#Finalise
variant_data_merge <- as.data.table(Bottlenecked.Data)
#Replace negative counts with 0
for(j in grep("_count$", names(variant_data_merge))){
set(variant_data_merge, i=which(variant_data_merge[[j]]<0), j=j, value=0)
}
#Remove variants with zero counts in all samples
variant_data_merge <- variant_data_merge[apply(variant_data_merge[,.SD,.SDcols = names(variant_data_merge)[grep("_count$", names(variant_data_merge))]], 1, sum)!=0,]
#Reformat
variant_data_merge <- variant_data_merge[,.SD,,.SDcols = grep("nt_seq|_count$", names(variant_data_merge))]
names(variant_data_merge) <- sapply(strsplit(names(variant_data_merge), "_e"), '[', 1)
#Save
OutFileName <- file.path(outpath, paste0("BB_TARDBP_290_2017-06-13_LibraryBottleneck_", substr(as.character(bottleneck.alpha), 3, 4), "_t0_variant_data_merge.tsv"))
write.table(variant_data_merge, OutFileName, sep = "\t", quote = F, row.names = F)
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.