NOTE: A preprint is coming soon! If you choose to try this, any constructive feedback is greatly appreciated!
This package simulates RNA-Seq experiments with an explicit definition of the absolute RNA counts.
This package extends the polyester
package by simulating an RNA-Seq using absolute RNA counts before
converting the result into relative proportions (in TPMs and expected
fragments) that are then used as input for polyester.
The basic workflow is the following:
- Take an estimate of TPM values from a sleuth object.
- Conceptually shift the unit from TPM to copy numbers per cell.
- Simulate differential expression on the copy numbers
- Convert the new copy numbers back to TPMs
- Convert TPM values to expected reads per transcript
- Simulate reads using the polyester (bioconductor page) package.
See the vignette for full details!