simulate_dropseq_experiment_countmat: Simulator: simulate_dropseq_experiment_countmat Simulate...

Description Usage Arguments Value

View source: R/simulate_dropseq_experiment_countmat.R

Description

create FASTA files containing RNA-seq reads simulated from provided transcripts, given a countmat matrix

Usage

1
2
simulate_dropseq_experiment_countmat(fasta = NULL, readmat,
  polyAnum = 15, bias = "empirical", model = NULL, outdir = ".")

Arguments

fasta

path to FASTA file containing transcripts from which to simulate reads.

readmat

a gene-cell matrix, each entry represents number of reads to simulate.

polyAnum

minimum number of 'A's in a polyA region, integer. A region should contain at least n continous 'A's to be considered as a polyA region.

bias

polyA sampling bias, one of 'empirical', 'custom', 'naive' (default 'empirical'). If 'empirical', the built-in model is adopted. This model is trained on data sets from the Drop-seq experiments described in the Cell paper (Macosko et al, 2015). If 'custom', polyA sampling bias is captured from user input model. The path of user input model is specified in model. If 'naive', sampling probability is calculated as (1/distance to polyA region).

model

path to polyA sampling bias model if bias set to 'custom'. Model format should be a data frame with 2 columns. Column 1 is distance to polyA region, while Column 2 is probability of getting a fragment at that position. Column 2 sums to 1.

outdir

path to folder where simulated reads should be written. By default, reads written to the working directory.

Value

No return, but simulated reads are written to outdir.


Helena-L/DropseqSimulator documentation built on May 14, 2019, 9:37 a.m.