data(sampleReadCounts)
data(referenceReadCounts)
## Gene names should be same size as row columns
geneNames <- row.names(referenceReadCounts)
ampliconNames <- NULL
normalizedReadCounts <- CombinedNormalizedCounts(sampleReadCounts,
referenceReadCounts,
ampliconNames = ampliconNames)
# After normalization data sets need to be splitted again to perform bootstrap
samplesNormalizedReadCounts = normalizedReadCounts["samples"][[1]]
referenceNormalizedReadCounts = normalizedReadCounts["reference"][[1]]
# Should be used values above 10000
replicates <- 10
# Perform the bootstrap based analysis
bootList <- BootList(geneNames,
samplesNormalizedReadCounts,
referenceNormalizedReadCounts,
replicates = replicates)
backgroundNoise <- Background(geneNames,
samplesNormalizedReadCounts,
referenceNormalizedReadCounts,
bootList,
replicates = replicates)
reportTables <- ReportTables(geneNames,
samplesNormalizedReadCounts,
referenceNormalizedReadCounts,
bootList,
backgroundNoise)
PlotBootstrapDistributions(bootList, reportTables, save = FALSE)
Run the code above in your browser using DataLab