# NOT RUN {
set.seed(1)
Counts <- matrix(rpois(50*10, 2), ncol = 10)
rownames(Counts) <- c(paste0('Gene', 1:40), paste0('Spike', 1:10))
Tech <- c(rep(FALSE,40),rep(TRUE,10))
set.seed(2)
SpikeInput <- rgamma(10,1,1)
SpikeInfo <- data.frame('SpikeID' = paste0('Spike', 1:10),
'SpikeInput' = SpikeInput)
Filter = BASiCS_Filter(Counts, Tech, SpikeInput,
MinTotalCountsPerCell = 2,
MinTotalCountsPerGene = 2,
MinCellsWithExpression = 2,
MinAvCountsPerCellsWithExpression = 2)
SpikeInfoFilter <- SpikeInfo[SpikeInfo$SpikeID %in% rownames(Filter$Counts),]
FilterData <- newBASiCS_Data(Filter$Counts, Filter$Tech, SpikeInfoFilter)
# }
Run the code above in your browser using DataLab